Changeset 6116 for LMDZ6/trunk/libf/phylmd/phyetat0_mod.f90

Timestamp:

Mar 21, 2026, 7:18:46 PM (2 weeks ago)

Author:

Laurent Fairhead

Message:

1. SZ98 interpolation (Sheng and Zwiers, Climate Dynamics (1998)

Mean-preserving linear interpolation of monthly boundary conditions to daily values

2. Concentration-Driven mode

Prescribed atmospheric CO2 concentrations without atmospheric CO2 transport

3. CO2 Fossil fuel emissions

Reading of monthly or daily emissions with SZ98 interpolation when needed

4. Flux corrections

Implementation of scalar (1D) ans spatially varying (2D) carbon flux corrections for piControl simulations.

PCadule

File:

• LMDZ6/trunk/libf/phylmd/phyetat0_mod.f90 (modified) (4 diffs)

LMDZ6/trunk/libf/phylmd/phyetat0_mod.f90

@@ -49 +49 @@
   USE traclmdz_mod,     ONLY: traclmdz_from_restart
   USE carbon_cycle_mod, ONLY: carbon_cycle_init, carbon_cycle_cpl, carbon_cycle_tr, carbon_cycle_rad, co2_send, RCO2_glo
+! << PC
+  USE carbon_cycle_mod, ONLY: carbon_cycle_conc_driven
+  USE mod_phys_lmdz_mpi_data, ONLY: is_mpi_root
+! >> PC  
   USE indice_sol_mod,   ONLY: nbsrf, is_ter, epsfra, is_lic, is_oce, is_sic
   !GG USE ocean_slab_mod,   ONLY: nslay, tslab, seaice, tice, ocean_slab_init
@@ -145 +149 @@
   ! 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
+! << PC  
+!  ! co2_ppm0 : initial value of atmospheric CO2 (from create_etat0_limit.e .def)
+!  co2_ppm0 = 284.32
+  ! co2_ppm0 is read from config.def by conf_phys_m
+  ! Only fall back to a default if it is not set to a value
+  IF (co2_ppm0 <= 0.0) co2_ppm0 = 284.32
+! >> PC  
+
   ! when no initial value is available e.g., from a restart 
   ! this variable must be set  in a .def file which will then be 
@@ -153 +163 @@
   ! 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
+! << PC  
+  IF (carbon_cycle_tr .OR. carbon_cycle_cpl .OR. carbon_cycle_conc_driven) THEN
+
+     IF (.NOT. carbon_cycle_conc_driven) THEN
+        ! Emission-driven: restart controls the initial CO2 value
+        co2_ppm = tab_cntrl(3)
+     ELSE
+        ! Concentration-driven: co2_ppm must be prescribed by the forcing
+        IF (co2_ppm <= 0.0) co2_ppm = tab_cntrl(3)
+     END IF
+
+     RCO2 = co2_ppm * 1.0e-06 * RMCO2 / RMD
+
+     IF (carbon_cycle_conc_driven) THEN
+        ! Concentration-driven: radiative CO2 must follow prescribed co2_ppm
+        RCO2_glo = co2_ppm * 1.0e-06 * RMCO2 / RMD
+     ELSE
+        ! Emission-driven: may restore from restart
+        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
+        ELSE
+           RCO2_glo = co2_ppm0 * 1.0e-06 * RMCO2 / RMD
+        END IF
+     END IF
+
+  END IF
+
+! >> PC
 
   solaire_etat0      = tab_cntrl(4)
@@ -602 +630 @@
 !===========================================
 
-!--OB now this is for co2i - ThL: and therefore also for inco
+!--PC+OB now this is for co2i - ThL: and therefore also for inco
   IF (ANY(type_trac == ['co2i','inco'])) THEN
-     IF (carbon_cycle_cpl) THEN
-        ALLOCATE(co2_send(klon), stat=ierr)
-        IF (ierr /= 0) CALL abort_physic('phyetat0', 'pb allocation co2_send', 1)
-        found=phyetat0_get(co2_send,"co2_send","co2 send",co2_ppm0)
-     ENDIF
+! << PC
+    ! Always allocate co2_send
+    IF (.NOT. ALLOCATED(co2_send)) ALLOCATE(co2_send(klon)) 
+
+    IF (carbon_cycle_conc_driven) THEN
+       ! In concentration driven mode, we do NOT read co2_send from restart
+       ! We rely on tracco2i to set it to co2_ppm 
+          PRINT*, 'phyetat0: Carbon cycle is Concentration Driven'
+          PRINT*, 'phyetat0: co2_send will be initialized to co2_ppm in tracco2i'
+          co2_send(:) = co2_ppm
+          IF (is_mpi_root) PRINT*, 'phyetat0: Conc-Driven. co2_send init to ', co2_ppm
+    ELSE
+       ! Emission driven: Standard behavior (alloc and read from restart)
+       IF (.NOT. ALLOCATED(co2_send)) ALLOCATE(co2_send(klon))
+       found = phyetat0_get(co2_send, "co2_send", "co2 send", co2_ppm0)
+       IF (.NOT.found) co2_send(:) = co2_ppm0
+    END IF
+! >> PC             
+    
   ELSE IF (type_trac == 'lmdz') THEN
      it = 0