Changeset 3865 for LMDZ6/trunk/libf/phylmd/physiq_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/physiq_mod.F90 (modified) (5 diffs)

LMDZ6/trunk/libf/phylmd/physiq_mod.F90

@@ -39 +39 @@
     USE ioipsl_getin_p_mod, ONLY : getin_p
     USE indice_sol_mod
-    USE infotrac_phy, ONLY: nqtot, nbtr, nqo, type_trac
+    USE infotrac_phy, ONLY: nqtot, nbtr, nqo, type_trac, nqINCA
     USE iophy
     USE limit_read_mod, ONLY : init_limit_read
@@ -1412 +1412 @@
        tau_overturning_th(:)=0.
 
-       IF (type_trac == 'inca') THEN
+       IF (type_trac == 'inca' .OR. type_trac == 'inco') THEN
           ! jg : initialisation jusqu'au ces variables sont dans restart
           ccm(:,:,:) = 0.
@@ -1968 +1968 @@
           cg_aero(:,:,:,:) = init_cginca
 !          
+
+
+          CALL VTe(VTinca)
+          CALL VTb(VTphysiq)
+#endif
+       ELSEIF (type_trac == 'inco') THEN
+#ifdef INCA
+          CALL VTe(VTphysiq)
+          CALL VTb(VTinca)
+          calday = REAL(days_elapsed) + jH_cur
+          WRITE(lunout,*) 'initial time chemini', days_elapsed, calday
+
+          CALL chemini(  &
+               rg, &
+               ra, &
+               cell_area, &
+               latitude_deg, &
+               longitude_deg, &
+               presnivs, &
+               calday, &
+               klon, &
+               nqtot, &
+               nqo+1, &   ! Note ThL: diff is here with case 'inca'
+               pdtphys, &
+               annee_ref, &
+               year_cur, &
+               day_ref,  &
+               day_ini, &
+               start_time, &
+               itau_phy, &
+               date0, &
+               io_lon, &
+               io_lat, &
+               chemistry_couple, &
+               init_source, &
+               init_tauinca, &
+               init_pizinca, &
+               init_cginca, &
+               init_ccminca)
+
+
+          ! initialisation des variables depuis le restart de inca
+          ccm(:,:,:) = init_ccminca
+          tau_aero(:,:,:,:) = init_tauinca
+          piz_aero(:,:,:,:) = init_pizinca
+          cg_aero(:,:,:,:) = init_cginca
+!
 
 
@@ -3792 +3839 @@
        CALL VTe(VTinca)
        CALL VTb(VTphysiq)
-#endif 
-    ENDIF !type_trac = inca
+#endif
+    ELSEIF (type_trac == 'inco') THEN
+#ifdef INCA
+       CALL VTe(VTphysiq)
+       CALL VTb(VTinca)
+       calday = REAL(days_elapsed + 1) + jH_cur
+
+       CALL chemtime(itap+itau_phy-1, date0, phys_tstep, itap)
+       CALL AEROSOL_METEO_CALC( &
+            calday,pdtphys,pplay,paprs,t,pmflxr,pmflxs, &
+            prfl,psfl,pctsrf,cell_area, &
+            latitude_deg,longitude_deg,u10m,v10m)
+
+       zxsnow_dummy(:) = 0.0
+
+       CALL chemhook_begin (calday, &
+            days_elapsed+1, &
+            jH_cur, &
+            pctsrf(1,1), &
+            latitude_deg, &
+            longitude_deg, &
+            cell_area, &
+            paprs, &
+            pplay, &
+            coefh(1:klon,1:klev,is_ave), &
+            pphi, &
+            t_seri, &
+            u, &
+            v, &
+            rot, &
+            wo(:, :, 1), &
+            q_seri, &
+            zxtsol, &
+            zt2m, &
+            zxsnow_dummy, &
+            solsw, &
+            albsol1, &
+            rain_fall, &
+            snow_fall, &
+            itop_con, &
+            ibas_con, &
+            cldfra, &
+            nbp_lon, &
+            nbp_lat-1, &
+            tr_seri(:,:,2:nbtr), &  ! Note ThL: diff is here with case 'inca'
+            ftsol, &
+            paprs, &
+            cdragh, &
+            cdragm, &
+            pctsrf, &
+            pdtphys, &
+            itap)
+
+       CALL VTe(VTinca)
+       CALL VTb(VTphysiq)
+#endif
+    ENDIF !type_trac = inca or inco
     IF (type_trac == 'repr') THEN
 #ifdef REPROBUS
@@ -4945 +5047 @@
        CALL VTb(VTphysiq)
 #endif
+    ELSEIF (type_trac == 'inco') THEN
+#ifdef INCA
+       CALL VTe(VTphysiq)
+       CALL VTb(VTinca)
+
+       CALL chemhook_end ( &
+            phys_tstep, &
+            pplay, &
+            t_seri, &
+            tr_seri(:,:,2:nbtr), & ! Note ThL: diff is here with case 'inca'
+            nbtr, &
+            paprs, &
+            q_seri, &
+            cell_area, &
+            pphi, &
+            pphis, &
+            zx_rh, &
+            aps, bps, ap, bp)
+
+       CALL VTe(VTinca)
+       CALL VTb(VTphysiq)
+#endif
     ENDIF