Changeset 4013 for LMDZ6/branches/Ocean_skin/libf/phylmd/phytrac_mod.F90

Timestamp:

Nov 19, 2021, 4:58:59 PM (3 years ago)

Author:

lguez

Message:

Sync latest trunk changes to Ocean_skin

Location:

LMDZ6/branches/Ocean_skin

Files:

• . (modified) (1 prop)
• libf/phylmd/phytrac_mod.F90 (modified) (13 diffs)

LMDZ6/branches/Ocean_skin/libf/phylmd/phytrac_mod.F90

@@ -56 +56 @@
   SUBROUTINE phytrac_init()
     USE dimphy
-    USE infotrac_phy, ONLY: nbtr, type_trac
+    USE infotrac_phy, ONLY: nbtr, nqCO2, type_trac
     USE tracco2i_mod, ONLY: tracco2i_init
     IMPLICIT NONE
@@ -81 +81 @@
     CASE('co2i')
        !   -- CO2 interactif --
+       CALL tracco2i_init()
+    CASE('inco')
        CALL tracco2i_init()
     END SELECT
@@ -122 +124 @@
     USE phys_cal_mod, only : hour
     USE dimphy
-    USE infotrac_phy, ONLY: nbtr, type_trac, conv_flg, solsym, pbl_flg
+    USE infotrac_phy, ONLY: nbtr, nqCO2, type_trac, conv_flg, solsym, pbl_flg
     USE mod_grid_phy_lmdz
     USE mod_phys_lmdz_para
@@ -176 +178 @@
     REAL,DIMENSION(klon,klev),INTENT(IN)   :: sh      ! humidite specifique
     REAL,DIMENSION(klon,klev),INTENT(IN)   :: rh      ! humidite relative
-    REAL,DIMENSION(klon,klev),INTENT(IN)   :: ch      ! eau liquide
+    REAL,DIMENSION(klon,klev),INTENT(IN)   :: ch      ! eau liquide (+ glace si le traceur existe) 
     REAL,DIMENSION(klon,klev+1),INTENT(IN) :: paprs   ! pression pour chaque inter-couche (en Pa)
     REAL,DIMENSION(klon,klev),INTENT(IN)   :: pplay   ! pression pour le mileu de chaque couche (en Pa)
@@ -305 +307 @@
 
     !
-    !Entrees/Sorties: (cf ini_histrac.h et write_histrac.h)  
+    !Entrees/Sorties:
     !---------------
     INTEGER                   :: iiq, ierr
@@ -332 +334 @@
     !---------- 
     REAL,DIMENSION(klon,klev,nbtr) :: flestottr ! flux de lessivage dans chaque couche 
-    REAL,DIMENSION(klon,klev)      :: zmasse    ! densité atmosphérique Kg/m2
+    REAL,DIMENSION(klon,klev)      :: zmasse    ! densite atmospherique Kg/m2
     REAL,DIMENSION(klon,klev)      :: ztra_th
     !PhH
@@ -505 +507 @@
           iflag_vdf_trac= 1
           iflag_con_trac= 1
+       CASE('inco')
+          source(:,1:nqCO2) = 0.                          ! from CO2i   ModThL
+          source(:,nqCO2+1:nbtr)=init_source(:,:)         ! from INCA   ModThL
+          aerosol(1:nqCO2) = .FALSE.                      ! from CO2i   ModThL
+          CALL tracinca_init(aerosol(nqCO2+1:nbtr),lessivage)     ! from INCA   ModThL
+          pbl_flg(1:nqCO2) = 1              ! From CO2i         ModThL
+          iflag_the_trac= 1           ! From CO2i
+          iflag_vdf_trac= 1           ! From CO2i
+          iflag_con_trac= 1           ! From CO2i
 #ifdef CPP_StratAer
        CASE('coag')
@@ -571 +582 @@
                 !--co2 tracers are not scavenged
                 flag_cvltr(it)=.FALSE.
-
+             CASE('inco')     ! Add ThL
+                flag_cvltr(it)=.FALSE.
 #ifdef CPP_StratAer
              CASE('coag')
@@ -590 +602 @@
           flag_cvltr(:) = .FALSE. 
        ENDIF
-       !
-       ! Initialize diagnostic output
-       ! ----------------------------
-#ifdef CPP_IOIPSL
-       !     INCLUDE "ini_histrac.h"
-#endif
        !
        ! print out all tracer flags 
@@ -614 +620 @@
        write(lunout,*)  'flag_cvltr    = ', flag_cvltr
 
-       IF (lessivage .AND. type_trac .EQ. 'inca') THEN
+       IF (lessivage .AND. (type_trac .EQ. 'inca' .OR. type_trac .EQ. 'inco')) THEN     ! Mod ThL
           CALL abort_physic('phytrac', 'lessivage=T config_inca=inca impossible',1)
 !          STOP
@@ -666 +672 @@
        !   -- sign convention : positive into the atmosphere
 
+       CALL tracco2i(pdtphys, debutphy, &
+            xlat, xlon, pphis, pphi, &
+            t_seri, pplay, paprs, tr_seri, source)
+    CASE('inco')      ! Add ThL
        CALL tracco2i(pdtphys, debutphy, &
             xlat, xlon, pphis, pphi, &
@@ -1092 +1102 @@
 
     !    -- CHIMIE INCA  config_inca = aero or chem --
-    IF (type_trac == 'inca') THEN
+    IF (type_trac == 'inca' .OR. type_trac == 'inco') THEN  ! ModThL
 
        CALL tracinca(&
@@ -1104 +1114 @@
             tau_aero, piz_aero, cg_aero,        ccm,       &
             rfname,                                        &
-            tr_seri,  source)      
-       
-       
+            tr_seri(:,:,1+nqCO2:nbtr),  source(:,1+nqCO2:nbtr))  ! ModThL  
     ENDIF
-    !=============================================================
-    !   Ecriture des sorties
-    !=============================================================
-#ifdef CPP_IOIPSL
-    ! INCLUDE "write_histrac.h"
-#endif
 
   END SUBROUTINE phytrac