Changeset 7709 for branches/ORCHIDEE_2_2/ORCHIDEE/src_sechiba/diffuco.f90

Timestamp:

2022-07-20T11:30:43+02:00 (2 years ago)

Author:

josefine.ghattas

Message:

Integrated new irrigation scheme developed by Pedro Arboleda. See ticket #857
This corresponds to revsion 7708 of version pedro.arboleda/ORCHIDEE. Following differences were made but were not made on the pedro.arboleda/ORCHIDEE :

• argumet place in call to routing_wrapper_intialize changed order
• lines with only change in space were not taken
• some indentation changed
• set irrigation output as enalbled false if not do_irrigation

File:

• branches/ORCHIDEE_2_2/ORCHIDEE/src_sechiba/diffuco.f90 (modified) (3 diffs)

branches/ORCHIDEE_2_2/ORCHIDEE/src_sechiba/diffuco.f90

@@ -332 +332 @@
     REAL(r_std),DIMENSION (kjpindex,nvm)     :: cim                !! Intercellular CO2 over nlai 
 
+    !! Based on ZunYin for correct transport
+    REAL(r_std),DIMENSION (kjpindex,nvm)     :: humrel_tmp         !! Soil moisture stress (within range 0 to 1)
+    REAL(r_std),DIMENSION (kjpindex,nvm)     :: vbeta3_tmp         !! Beta for transpiration (-)
+    REAL(r_std),DIMENSION (kjpindex,nvm)     :: rveget_tmp         !! Stomatal resistance for the whole canopy (s m^{-1})
+    REAL(r_std),DIMENSION (kjpindex,nvm)     :: rstruct_tmp        !! Structural resistance for the vegetation
+    REAL(r_std),DIMENSION (kjpindex,nvm)     :: cimean_tmp         !! Mean leaf Ci (ppm)
+    REAL(r_std),DIMENSION (kjpindex,nvm)     :: gsmean_tmp         !! Mean stomatal conductance to CO2 (umol m-2 s-1)
+    REAL(r_std),DIMENSION (kjpindex,nvm)     :: gpp_tmp            !! Assimilation ((gC m^{-2} s^{-1}), total area)
+    REAL(r_std),DIMENSION (kjpindex,nvm)     :: cim_tmp            !! Intercellular CO2 over nlai
+
 !_ ================================================================================================================================
     
@@ -368 +378 @@
 
   !! 5. Calculate the beta coefficient for transpiration
+    !! Dummy variables to compute vbeta3pot (humrel_temp = 1, no water stress)
+    humrel_tmp(:,:) = un
+    vbeta3_tmp = vbeta3
+    rveget_tmp = rveget
+    rstruct_tmp = rstruct
+    cimean_tmp = cimean
+    gsmean_tmp = gsmean
+    gpp_tmp = gpp
+    cim_tmp = cim
 
     CALL diffuco_trans_co2 (kjpindex, swdown, pb, qsurf, qair, temp_air, temp_growth, rau, u, v, q_cdrag, humrel, &
@@ -374 +393 @@
          rveget, rstruct, cimean, gsmean, gpp, &
          vbeta23, hist_id, indexveg, indexlai, index, kjit, cim)
+
+   !! New resistance calc. with temporary (dummy) variables, just for vbeta3pot
+    CALL diffuco_trans_co2 (kjpindex, swdown, pb, qsurf, qair, temp_air, temp_growth, rau, u, v, q_cdrag, humrel_tmp, &
+         assim_param, ccanopy, &
+         veget, veget_max, lai, qsintveg, qsintmax, vbeta3_tmp, vbeta3pot, &
+         rveget_tmp, rstruct_tmp, cimean_tmp, gsmean_tmp, gpp_tmp, &
+         vbeta23, hist_id, indexveg, indexlai, index, kjit, cim_tmp)
 
     !