Changeset 8418 for branches/ORCHIDEE_2_2/ORCHIDEE/src_sechiba/sechiba.f90

Timestamp:

2024-02-12T19:49:09+01:00 (5 months ago)

Author:

bertrand.guenet

Message:

The Moyano function describing the soil moisture effect on OM decomposition is added

File:

• branches/ORCHIDEE_2_2/ORCHIDEE/src_sechiba/sechiba.f90 (modified) (7 diffs)

branches/ORCHIDEE_2_2/ORCHIDEE/src_sechiba/sechiba.f90

@@ -140 +140 @@
                                                                      !! by thermosoil.f90 (unitless, 0-1)
 !$OMP THREADPRIVATE(shumdiag)
+  REAL(r_std), ALLOCATABLE, SAVE, DIMENSION (:,:)  :: shumdiagSAT    !! Mean relative soil moisture (0-1, unitless) 
+                                                                     !! with respect to(mcs-mcw)
+!$OMP THREADPRIVATE(shumdiagSAT) 
   REAL(r_std), ALLOCATABLE, SAVE, DIMENSION (:,:)  :: shumdiag_perma !! Saturation degree of the soil 
 !$OMP THREADPRIVATE(shumdiag_perma)
@@ -146 +149 @@
   REAL(r_std), ALLOCATABLE, SAVE, DIMENSION (:)    :: litterhumdiag  !! Litter dryness factor (unitless, 0-1)
 !$OMP THREADPRIVATE(litterhumdiag)
+  REAL(r_std), ALLOCATABLE, SAVE, DIMENSION (:)    :: litterhumdiagSAT!! Litter dryness factor (unitless, 0-1)
+!$OMP THREADPRIVATE(litterhumdiagSAT)
   REAL(r_std), ALLOCATABLE, SAVE, DIMENSION (:,:)  :: stempdiag      !! Temperature which controls canopy evolution (K)
 !$OMP THREADPRIVATE(stempdiag)
@@ -765 +770 @@
          & tot_melt, transpir, precip_rain, precip_snow, returnflow, reinfiltration, irrigation, &
          & humrel, vegstress, drysoil_frac, evapot, evapot_corr, evap_bare_lim, evap_bare_lim_ns, flood_frac, flood_res, &
-         & shumdiag,shumdiag_perma, k_litt, litterhumdiag, soilcap, soiltile, fraclut, reinf_slope_soil,&
+         & shumdiagSAT,litterhumdiagSAT,shumdiag,shumdiag_perma, k_litt, litterhumdiag, soilcap, soiltile, fraclut, reinf_slope,&
          & rest_id, hist_id, hist2_id,&
          & contfrac, stempdiag, &
@@ -820 +825 @@
     CALL slowproc_main (kjit, kjpij, kjpindex, &
          index, indexveg, lalo, neighbours, resolution, contfrac, soiltile, fraclut, nwdFraclut, &
-         temp_air, temp_sol, stempdiag, &
+         temp_air, temp_sol, stempdiag,shumdiagSAT,litterhumdiagSAT, &
          vegstress, shumdiag, litterhumdiag, precip_rain, precip_snow, gpp, &
          deadleaf_cover, &
@@ -1765 +1770 @@
     ALLOCATE (shumdiag(kjpindex,nslm),stat=ier)
     IF (ier /= 0) CALL ipslerr_p(3,'sechiba_init','Pb in alloc for shumdiag','','')
+
+    ALLOCATE (shumdiagSAT(kjpindex,nslm),stat=ier)
+    IF (ier /= 0) CALL ipslerr_p(3,'sechiba_init','Pb in alloc for shumdiagSAT','','')
     
     ALLOCATE (shumdiag_perma(kjpindex,nslm),stat=ier)
@@ -1771 +1779 @@
     ALLOCATE (litterhumdiag(kjpindex),stat=ier)
     IF (ier /= 0) CALL ipslerr_p(3,'sechiba_init','Pb in alloc for litterhumdiag','','')
+
+    ALLOCATE (litterhumdiagSAT(kjpindex),stat=ier)
+    IF (ier /= 0) CALL ipslerr_p(3,'sechiba_init','Pb in alloc for litterhumdiagSAT','','')
 
     ALLOCATE (ptnlev1(kjpindex),stat=ier)
@@ -2021 +2032 @@
     IF ( ALLOCATED (co2_flux)) DEALLOCATE (co2_flux)
     IF ( ALLOCATED (shumdiag)) DEALLOCATE (shumdiag)
+    IF ( ALLOCATED (shumdiagSAT)) DEALLOCATE (shumdiagSAT)
     IF ( ALLOCATED (shumdiag_perma)) DEALLOCATE (shumdiag_perma)
     IF ( ALLOCATED (litterhumdiag)) DEALLOCATE (litterhumdiag)
+    IF ( ALLOCATED (litterhumdiagSAT)) DEALLOCATE (litterhumdiagSAT)
     IF ( ALLOCATED (ptnlev1)) DEALLOCATE (ptnlev1)
     IF ( ALLOCATED (k_litt)) DEALLOCATE (k_litt)