Changeset 3865 for LMDZ6/trunk/libf/phylmd/phytrac_mod.F90

Timestamp:

Mar 23, 2021, 4:14:07 PM (3 years ago)

Author:

lmdz-users

Message:

Modifications from Thibaut to create an ESM with interactive CO2 + INCA aerosols

File:

• LMDZ6/trunk/libf/phylmd/phytrac_mod.F90 (modified) (8 diffs)

LMDZ6/trunk/libf/phylmd/phytrac_mod.F90

@@ -56 +56 @@
   SUBROUTINE phytrac_init()
     USE dimphy
-    USE infotrac_phy, ONLY: nbtr, type_trac
+    USE infotrac_phy, ONLY: nbtr, nqINCA, 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, nqINCA, type_trac, conv_flg, solsym, pbl_flg
     USE mod_grid_phy_lmdz
     USE mod_phys_lmdz_para
@@ -505 +507 @@
           iflag_vdf_trac= 1
           iflag_con_trac= 1
+       CASE('inco')
+          source(:,1) = 0.                          ! from CO2i
+          source(:,2:nbtr)=init_source(:,:)         ! from INCA
+          aerosol(1) = .FALSE.                      ! from CO2i
+          CALL tracinca_init(aerosol(2:nbtr),lessivage)     ! from INCA
+          pbl_flg(1) = 1              ! From CO2i
+          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')
@@ -614 +626 @@
        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 +678 @@
        !   -- 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, &
@@ -1104 +1120 @@
             tau_aero, piz_aero, cg_aero,        ccm,       &
             rfname,                                        &
-            tr_seri,  source)      
-       
-       
+            tr_seri,  source)
+    ELSEIF (type_trac == 'inco') THEN       ! Add ThL
+       CALL tracinca(&
+            nstep,    julien,   gmtime,         lafin,     &
+            pdtphys,  t_seri,   paprs,          pplay,     &
+            pmfu,     upwd,     ftsol,  pctsrf, pphis,     &
+            pphi,     albsol,   sh,    ch,     rh,        &
+            cldfra,   rneb,     diafra,         cldliq,    &
+            itop_con, ibas_con, pmflxr,         pmflxs,    &
+            prfl,     psfl,     aerosol_couple, flxmass_w, &
+            tau_aero, piz_aero, cg_aero,        ccm,       &
+            rfname,                                        &
+            tr_seri(:,:,2:nbtr),  source(:,2:nbtr))     ! Difference with case 'inca'  
     ENDIF
     !=============================================================