Changeset 4298 for LMDZ6/trunk/libf/phylmd/phyetat0.F90

Timestamp:

Oct 17, 2022, 10:15:06 AM (2 years ago)

Author:

pcadule

Message:

modifications to ensure restartability of the model when CO2 tracer is passed to radiative code, and to add diagnostics variables

File:

• LMDZ6/trunk/libf/phylmd/phyetat0.F90 (modified) (4 diffs)

LMDZ6/trunk/libf/phylmd/phyetat0.F90

@@ -25 +25 @@
   USE infotrac_phy,     ONLY: nqtot, nbtr, type_trac, types_trac, tracers
   USE traclmdz_mod,     ONLY: traclmdz_from_restart
-  USE carbon_cycle_mod, ONLY: carbon_cycle_tr, carbon_cycle_cpl, co2_send
+  USE carbon_cycle_mod, ONLY: carbon_cycle_init, carbon_cycle_cpl, carbon_cycle_tr, carbon_cycle_rad, co2_send, RCO2_glo
   USE indice_sol_mod,   ONLY: nbsrf, is_ter, epsfra, is_lic, is_oce, is_sic
   USE ocean_slab_mod,   ONLY: nslay, tslab, seaice, tice, ocean_slab_init
@@ -113 +113 @@
 
   ! co2_ppm : value from the previous time step
+
+  ! co2_ppm0 : initial value of atmospheric CO2 (from create_etat0_limit.e .def)
+  co2_ppm0 = 284.32
+  ! when no initial value is available e.g., from a restart 
+  ! this variable must be set  in a .def file which will then be 
+  ! used by the conf_phys_m.F90 routine.
+  ! co2_ppm0 = 284.32 (illustrative example on how to set the variable in .def
+  ! file, for a pre-industrial CO2 concentration value)
+
   IF (carbon_cycle_tr .OR. carbon_cycle_cpl) THEN
      co2_ppm = tab_cntrl(3)
      RCO2    = co2_ppm * 1.0e-06 * RMCO2 / RMD
+     IF (tab_cntrl(17) > 0. .AND. carbon_cycle_rad) THEN
+           RCO2_glo = tab_cntrl(17)
+       ELSE
+           RCO2_glo    = co2_ppm0 * 1.0e-06 * RMCO2 / RMD
+     ENDIF
      ! ELSE : keep value from .def
   ENDIF
-
-! co2_ppm0 : initial value of atmospheric CO2 (from create_etat0_limit.e .def)
-! co2_ppm0   = tab_cntrl(16)
-! initial value for interactive CO2 run when there is no tracer field for CO2 in restart
-  co2_ppm0=284.32
 
   solaire_etat0      = tab_cntrl(4)
@@ -145 +154 @@
   clesphy0(7)=tab_cntrl( 11 )
   clesphy0(8)=tab_cntrl( 12 )
+  clesphy0(9)=tab_cntrl( 17 )
 
   ! set time iteration
@@ -454 +464 @@
         ALLOCATE(co2_send(klon), stat=ierr)
         IF (ierr /= 0) CALL abort_physic('phyetat0', 'pb allocation co2_send', 1)
-        found=phyetat0_get(1,co2_send,"co2_send","co2 send",co2_ppm)
+        !found=phyetat0_get(1,co2_send,"co2_send","co2 send",co2_ppm)
+        found=phyetat0_get(1,co2_send,"co2_send","co2 send",co2_ppm0)
      ENDIF
   ELSE IF (type_trac == 'lmdz') THEN