Changeset 8220

Timestamp:

2023-10-10T15:50:31+02:00 (7 months ago)

Author:

anne.cozic

Message:

Add initialization of ptnlev1 used in chemistry (specific case of VOC coupling between atm and surf) before it was calculated by thermosoil_main

Location:

branches/ORCHIDEE_2_2/ORCHIDEE/src_sechiba

Files:

• sechiba.f90 (modified) (1 diff)
• thermosoil.f90 (modified) (3 diffs)

branches/ORCHIDEE_2_2/ORCHIDEE/src_sechiba/sechiba.f90

@@ -517 +517 @@
          mc_layh,  mcl_layh,   soilmoist,       njsc ,     &
          frac_snow_veg,frac_snow_nobio,totfrac_nobio,     &
-         snowdz, snowrho, snowtemp, lambda_snow, cgrnd_snow, dgrnd_snow, pb)
+         snowdz, snowrho, snowtemp, lambda_snow, cgrnd_snow, dgrnd_snow, pb, &
+         ptnlev1)
 
 

branches/ORCHIDEE_2_2/ORCHIDEE/src_sechiba/thermosoil.f90

@@ -235 +235 @@
                                     mc_layh,       mcl_layh,   tmc_layh,        njsc,     &
                                     frac_snow_veg,frac_snow_nobio,totfrac_nobio, &
-                                    snowdz, snowrho, snowtemp, lambda_snow, cgrnd_snow, dgrnd_snow, pb)
+                                    snowdz, snowrho, snowtemp, lambda_snow, cgrnd_snow, dgrnd_snow, pb, &
+                                    ptnlev1)
 
     !! 0. Variable and parameter declaration
@@ -265 +266 @@
     REAL(r_std),DIMENSION (kjpindex,ngrnd), INTENT (out)  :: ftempdiag        !! temperature profile on full depth for stream temperature (K)
     REAL(r_std),DIMENSION (kjpindex),INTENT(out)          :: gtemp            !! First soil layer temperature
+    REAL(r_std),DIMENSION (kjpindex), INTENT (out)        :: ptnlev1          !! 1st level soil temperature   
 
     !! 0.3 Modified variables
@@ -521 +523 @@
     stempdiag(:,:) = ptn(:,1:nslm)
     ftempdiag(:,:) = ptn(:,1:ngrnd)
-    
+
+    !! Initialize output arguments to be used in sechiba
+    ptnlev1(:) = ptn(:,1)
 
     !! 2.3. Computes cgrnd, dgrnd, soilflx and soilcap coefficients only if they were not found in restart file.