Changeset 6116 for LMDZ6/trunk/libf/phylmd/pbl_surface_subsrf_mod.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/pbl_surface_subsrf_mod.F90 (modified) (3 diffs)

LMDZ6/trunk/libf/phylmd/pbl_surface_subsrf_mod.F90

@@ -175 +175 @@
 ! fluxlat_tersrf--output-R- latent heat flux of continental sub-surfaces
     USE carbon_cycle_mod,   ONLY : carbon_cycle_cpl  
-    USE carbon_cycle_mod,   ONLY : co2_send, nbcf_out, fields_out, yfields_out
+    USE carbon_cycle_mod,   ONLY : co2_send, nbcf_out, fields_out, yfields_out, yfields_in
     use hbtm_mod, only: hbtm
     USE indice_sol_mod
@@ -800 +800 @@
 !ORCHIDEE grid and as such should be defined in yfields_in(knon,nbcf_in) but the
 !knon variable is not known at that level of pbl_surface_mod
-  yfields_out(:,:) = 0.
+  IF (ALLOCATED(yfields_out)) yfields_out(:,:) = 0.
+  IF (ALLOCATED(yfields_in))  yfields_in(:,:) = 0.
+! >> PC
 
   ypphi = 0.0  
@@ -891 +893 @@
 #endif
        ENDDO
+! >> PC
 !--compressing fields_out onto ORCHIDEE grid
-!--these fields are shared and used directly surf_land_orchidee_mod
-       DO n = 1, nbcf_out
-         DO j = 1, knon
-           i = ni(j)
-           yfields_out(j,n) = fields_out(i,n)
-         ENDDO
-       ENDDO
+!--these fields are shared and used directly in surf_land_orchidee_mod
+        IF (nbcf_out > 0 .AND. ALLOCATED(yfields_out) .AND. ALLOCATED(fields_out)) THEN
+           DO n = 1, nbcf_out
+              DO j = 1, knon
+                 i = ni(j)
+                 yfields_out(j,n) = fields_out(i,n)
+              ENDDO
+           ENDDO
+        END IF
+! << PC
+
 
        DO k = 1, klev