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

Timestamp:

2022-07-20T13:09:05+02:00 (2 years ago)

Author:

josefine.ghattas

Message:

Integration of temperature of water in highres routing scheme. Done by Jan Polcer and Lucia Rinchiuso

File:

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

branches/ORCHIDEE_2_2/ORCHIDEE/src_sechiba/sechiba.f90

@@ -13 +13 @@
 !! processes as well.
 !!
-!!\n DESCRIPTION  : :: shumdiag, :: litterhumdiag and :: stempdiag are not 
+!!\n DESCRIPTION  : :: shumdiag, :: litterhumdiag and :: stempdiag :: ftempdiag are not 
 !! saved in the restart file because at the first time step because they 
 !! are recalculated. However, they must be saved as they are in slowproc 
@@ -148 +148 @@
   REAL(r_std), ALLOCATABLE, SAVE, DIMENSION (:,:)  :: stempdiag      !! Temperature which controls canopy evolution (K)
 !$OMP THREADPRIVATE(stempdiag)
+  REAL(r_std), ALLOCATABLE, SAVE, DIMENSION (:,:)  :: ftempdiag      !! Temperature over the full soil column for river temperature (K)
+!$OMP THREADPRIVATE(ftempdiag)
   REAL(r_std), ALLOCATABLE, SAVE, DIMENSION (:,:)  :: qsintveg       !! Water on vegetation due to interception 
                                                                      !! @tex $(kg m^{-2})$ @endtex
@@ -511 +513 @@
     CALL thermosoil_initialize (kjit, kjpindex, rest_id, mcs,  &
          temp_sol_new, snow,       shumdiag_perma,        &
-         soilcap,      soilflx,    stempdiag,             &
+         soilcap,      soilflx,    stempdiag, ftempdiag,             &
          gtemp,               &
          mc_layh,  mcl_layh,   soilmoist,       njsc ,     &
@@ -525 +527 @@
             kjit,        kjpindex,       index,                 &
             rest_id,     hist_id,        hist2_id,   lalo,      &
-            neighbours,  resolution,     contfrac,   stempdiag, &
+            neighbours,  resolution,     contfrac,   stempdiag, ftempdiag, &
             soiltile,    irrig_frac,     veget_max,  irrigated_next, &
             returnflow,  reinfiltration, irrigation, riverflow, &
@@ -786 +788 @@
          index, indexgrnd, mcs, &
          temp_sol_new, snow, soilcap, soilflx, &
-         shumdiag_perma, stempdiag, ptnlev1, rest_id, hist_id, hist2_id, &
+         shumdiag_perma, stempdiag, ftempdiag, ptnlev1, rest_id, hist_id, hist2_id, &
          snowdz,snowrho,snowtemp,gtemp,pb,&
          mc_layh, mcl_layh, soilmoist, njsc,frac_snow_veg,frac_snow_nobio,totfrac_nobio,temp_sol_add, &
@@ -797 +799 @@
        CALL routing_wrapper_main (kjit, kjpindex, index, &
             & lalo, neighbours, resolution, contfrac, totfrac_nobio, veget_max, floodout, runoff, &
-            & drainage, transpot, precip_rain, humrel, k_litt, flood_frac, flood_res, &
-            & stempdiag, reinf_slope, returnflow, reinfiltration, irrigation, riverflow, coastalflow, rest_id, hist_id, hist2_id, &
+            & drainage, transpot, precip_rain, humrel, k_litt, flood_frac, flood_res, stempdiag, &
+            & ftempdiag, reinf_slope, returnflow, reinfiltration, irrigation, riverflow, coastalflow, rest_id, hist_id, hist2_id, &
             & soiltile, root_deficit, irrigated_next, irrig_frac, fraction_aeirrig_sw)
     ELSE
@@ -1753 +1755 @@
     IF (ier /= 0) CALL ipslerr_p(3,'sechiba_init','Pb in alloc for stempdiag','','')
 
+    ALLOCATE (ftempdiag(kjpindex, ngrnd),stat=ier)
+    IF (ier /= 0) CALL ipslerr_p(3,'sechiba_init','Pb in alloc for ftempdiag','','')
+
     ALLOCATE (co2_flux(kjpindex,nvm),stat=ier)
     IF (ier /= 0) CALL ipslerr_p(3,'sechiba_init','Pb in alloc for co2_flux','','')
@@ -2012 +2017 @@
     IF ( ALLOCATED (deadleaf_cover)) DEALLOCATE (deadleaf_cover)
     IF ( ALLOCATED (stempdiag)) DEALLOCATE (stempdiag)
+    IF ( ALLOCATED (ftempdiag)) DEALLOCATE (ftempdiag)
     IF ( ALLOCATED (co2_flux)) DEALLOCATE (co2_flux)
     IF ( ALLOCATED (shumdiag)) DEALLOCATE (shumdiag)