Changeset 4298 for LMDZ6/trunk/libf/phylmdiso

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

Location:

LMDZ6/trunk/libf/phylmdiso

Files:

• phyetat0.F90 (modified) (4 diffs)
• phyredem.F90 (modified) (2 diffs)
• physiq_mod.F90 (modified) (1 diff)

LMDZ6/trunk/libf/phylmdiso/phyetat0.F90

@@ -33 +33 @@
   USE infotrac_phy,     ONLY: nqtot, nbtr, 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
@@ -134 +134 @@
 
   ! 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)
@@ -166 +175 @@
   clesphy0(7)=tab_cntrl( 11 )
   clesphy0(8)=tab_cntrl( 12 )
+  clesphy0(9)=tab_cntrl( 17 )
 
   ! set time iteration
@@ -473 +483 @@
         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_ppm0)
      ENDIF
   ELSE IF (ANY(types_trac == 'lmdz')) THEN

LMDZ6/trunk/libf/phylmdiso/phyredem.F90

@@ -45 +45 @@
 #endif
 #endif
-  USE carbon_cycle_mod, ONLY : carbon_cycle_cpl, co2_send
+  USE carbon_cycle_mod, ONLY : carbon_cycle_cpl, co2_send, carbon_cycle_rad, RCO2_glo
   USE indice_sol_mod, ONLY: nbsrf, is_oce, is_sic, is_ter, is_lic, epsfra
   USE surface_data, ONLY: type_ocean, version_ocean
@@ -130 +130 @@
 
   ! co2_ppm0 : initial value of atmospheric CO2
-  tab_cntrl(16) = co2_ppm0
+  ! tab_cntrl(16) = co2_ppm0
+
+  !  PC -- initial value of RCO2 for the radiation scheme
+  !  tab_cntrl(17) = co2_ppm * 1.0e-06 * RMCO2 / RMD
+  IF (carbon_cycle_rad) tab_cntrl(17) = RCO2_glo
+  !PRINT*, "PC : phyredem RCO2_glo =",RCO2_glo
 
   DO pass=1,2   ! pass=1 netcdf definition ; pass=2 netcdf write

LMDZ6/trunk/libf/phylmdiso/physiq_mod.F90

@@ -5497 +5497 @@
           ! 
           !--interactive CO2 in ppm from carbon cycle
-          IF (carbon_cycle_rad.AND..NOT.debut) THEN
-            RCO2=RCO2_glo
-          ENDIF
+          IF (carbon_cycle_rad) RCO2=RCO2_glo
           !
           IF (prt_level .GE.10) THEN