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

Timestamp:

Mar 21, 2026, 7:18:46 PM (3 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/conf_phys_m.f90 (modified) (6 diffs)

LMDZ6/trunk/libf/phylmd/conf_phys_m.f90

@@ -70 +70 @@
     USE phys_cal_mod
     USE carbon_cycle_mod,  ONLY: carbon_cycle_tr, carbon_cycle_cpl, carbon_cycle_rad, level_coupling_esm
+! << PC
+    USE carbon_cycle_mod,  ONLY: read_daily_co2ff
+    USE carbon_cycle_mod,  ONLY: carbon_cycle_conc_driven
+! >> PC        
     USE carbon_cycle_mod,  ONLY: read_fco2_ocean_cor,var_fco2_ocean_cor
     USE carbon_cycle_mod,  ONLY: read_fco2_land_cor,var_fco2_land_cor
+    USE carbon_cycle_mod,  ONLY: read_fco2_ocean_cor_2d, read_fco2_land_cor_2d
     USE chemistry_cycle_mod, ONLY: dms_cycle_cpl, n2o_cycle_cpl, ndp_cycle_cpl, nh3_cycle_cpl
     USE mod_grid_phy_lmdz, ONLY: klon_glo
@@ -1190 +1195 @@
     CALL getin_p('carbon_cycle_rad',carbon_cycle_rad)
 
-    read_fco2_ocean_cor=.FALSE.
-    CALL getin_p('read_fco2_ocean_cor',read_fco2_ocean_cor)
+    !Patricia
+    read_daily_co2ff = .FALSE.
+    CALL getin_p('read_daily_co2ff', read_daily_co2ff)
 
     var_fco2_ocean_cor=0. ! default value
@@ -1201 +1207 @@
     var_fco2_land_cor=0. ! default value
     CALL getin_p('var_fco2_land_cor',var_fco2_land_cor)
+
+    read_fco2_ocean_cor_2d=.FALSE.
+    CALL getin_p('read_fco2_ocean_cor_2d',read_fco2_ocean_cor_2d)
+    
+    read_fco2_land_cor_2d=.FALSE.
+    CALL getin_p('read_fco2_land_cor_2d',read_fco2_land_cor_2d)
 
     ! level_coupling_esm : level of coupling of the biogeochemical fields between LMDZ, ORCHIDEE and NEMO
@@ -1215 +1227 @@
     CALL getin_p('level_coupling_esm',level_coupling_esm)
 
-
+    ! Patricia
+    carbon_cycle_conc_driven=.FALSE.
+    CALL getin_p('carbon_cycle_conc_driven', carbon_cycle_conc_driven)
+
+    ! Check: If concentration-driven, ensure co2_ppm is valid
+    IF (carbon_cycle_conc_driven) THEN
+       IF (co2_ppm <= 0.0) THEN
+          co2_ppm = co2_ppm0
+          WRITE(lunout,*) 'conf_phys: WARNING: co2_ppm was <= 0 in concentration-driven mode'
+          WRITE(lunout,*) '           Resetting to initial value co2_ppm0 = ', co2_ppm
+       END IF
+       WRITE(lunout,*) 'conf_phys: Carbon Cycle is CONCENTRATION-DRIVEN'
+       WRITE(lunout,*) '           Atmosphere forced by co2_ppm'
+    ELSE
+       WRITE(lunout,*) 'conf_phys: Carbon Cycle is EMISSION-DRIVEN (Interactive)'
+    END IF
+
+    
     dms_cycle_cpl=.FALSE.
     CALL getin_p('transm_dms_oa',dms_cycle_cpl)
+
     !
 
@@ -1382 +1412 @@
     ENDIF
 
+! << PC
+    ! Consistency checks for CO2 driving philosophy
+    ! If CO2 is concentration-driven for the surface components, the radiation
+    ! should not use the drifting tracer CO2 (avoid inconsistent forcing)
+    IF (carbon_cycle_conc_driven .AND. carbon_cycle_rad) THEN
+       !CALL abort_physic('conf_phys', &
+       !     'carbon_cycle_rad must be FALSE when carbon_cycle_conc_driven is TRUE', 1)
+       WRITE(lunout,*) 'CO2 concentration-driven: radiation uses prescribed co2_ppm: atmospheric tracer is diagnostic only.'
+    ENDIF
+! >> PC
 
     ! Test on chemistry cycle 
@@ -1565 +1605 @@
     WRITE(lunout,*) 'conf_phys :  carbon_cycle_rad = ', carbon_cycle_rad
     WRITE(lunout,*) 'conf_phys :  level_coupling_esm = ', level_coupling_esm
+    WRITE(lunout,*) 'conf_phys :  carbon_cycle_conc_driven = ', carbon_cycle_conc_driven
+    WRITE(lunout,*) 'conf_phys :  read_daily_co2ff = ', read_daily_co2ff
+    WRITE(lunout,*) 'conf_phys :  read_fco2_ocean_cor_2d = ', read_fco2_ocean_cor_2d    
+    WRITE(lunout,*) 'conf_phys :  read_fco2_land_cor_2d = ', read_fco2_land_cor_2d
     WRITE(lunout,*) 'conf_phys :  read_fco2_ocean_cor = ', read_fco2_ocean_cor
     WRITE(lunout,*) 'conf_phys :  var_fco2_ocean_cor = ', var_fco2_ocean_cor