TEB open source v2.4406.

Manually constructed from the official project release. Archive downloaded and extracted as-is from https://opensource.umr-cnrm.fr/attachments/download/1555/TEB_open_source_v2.4406.tar.gz
TEB-model · Apr 14, 2019 · 4cc9760 · 4cc9760
1 parent 76afbc5
commit 4cc9760
Show file tree

Hide file tree

Showing 113 changed files with 194,389 additions and 193,834 deletions.
diff --git a/output_ref/HVAC_COOL.txt b/output_ref/HVAC_COOL.txt
diff --git a/output_ref/HVAC_HEAT.txt b/output_ref/HVAC_HEAT.txt
diff --git a/output_ref/H_TOWN.txt b/output_ref/H_TOWN.txt
diff --git a/output_ref/LE_TOWN.txt b/output_ref/LE_TOWN.txt
diff --git a/output_ref/Q_CANYON.txt b/output_ref/Q_CANYON.txt
diff --git a/output_ref/RN_TOWN.txt b/output_ref/RN_TOWN.txt
diff --git a/output_ref/TI_BLD.txt b/output_ref/TI_BLD.txt
diff --git a/output_ref/T_CANYON.txt b/output_ref/T_CANYON.txt
diff --git a/output_ref/T_ROAD1.txt b/output_ref/T_ROAD1.txt
diff --git a/output_ref/T_ROOF1.txt b/output_ref/T_ROOF1.txt
diff --git a/output_ref/T_WALLA1.txt b/output_ref/T_WALLA1.txt
diff --git a/output_ref/T_WALLB1.txt b/output_ref/T_WALLB1.txt
diff --git a/output_ref/U_CANYON.txt b/output_ref/U_CANYON.txt
@@ -1673,11 +1673,11 @@
   0.24256437568892866     
   0.32635934183601312     
   0.10143601165173380     
-  7.49744433947597683E-002
+   7.4974443394759768E-002
   0.25579515981741563     
   0.33517986458833782     
   0.18964123917498057     
-  9.26154888994091124E-002
+   9.2615488899409112E-002
   0.25579515981741563     
   0.40133378523077290     
   0.37928247834996115     
@@ -11226,7 +11226,7 @@
   0.34896545959794606     
   0.29290398108872084     
   0.21398279478843338     
-  9.97415892310122287E-002
+   9.9741589231012229E-002
   0.32508257116761408     
   0.45990690759372249     
   0.29668842637560583     
@@ -11985,13 +11985,13 @@
   0.67476999055283782     
   0.41897483073542224     
   0.21610280743195459     
-  7.49744433947597683E-002
-  7.49744433947597683E-002
+   7.4974443394759768E-002
+   7.4974443394759768E-002
   0.10143601165173383     
   0.10584627302789613     
   0.21610280743195459     
-  6.17436592662727601E-002
-  9.26154888994091402E-002
+   6.1743659266272760E-002
+   9.2615488899409140E-002
   0.21610280743195459     
   0.21610280743195459     
   0.18964123917498057     
@@ -12610,7 +12610,7 @@
   0.17200019367033123     
   0.22492333018427932     
   0.19405150055114292     
-  9.70257502755714624E-002
+   9.7025750275571462E-002
   0.36605169422147421     
   0.40133378523077290     
   0.46307744449704563     
@@ -12627,7 +12627,7 @@
   0.39251326247844825     
   0.20287202330346765     
   0.13671810266103251     
-  6.61539206424350962E-002
+   6.6153920642435096E-002
   0.31753881908368842     
   0.59097502440575345     
   0.87764201385630547     
@@ -13239,8 +13239,8 @@
   0.26853861006383678     
   0.22276009697927171     
   0.10171077093546872     
-  5.73333978901104171E-002
-  5.81294500685077126E-002
+   5.7333397890110417E-002
+   5.8129450068507713E-002
   0.12011831986729508     
   0.19983555834597985     
   0.31185972550960411     
@@ -15397,7 +15397,7 @@
   0.37801849743955301     
   0.33949839252787595     
   0.16699145982542368     
-  9.22207705062425970E-002
+   9.2220770506242597E-002
   0.15341226504821984     
   0.16735883459805798     
   0.45205884747484160     

diff --git a/src_driver/driver.F90 b/src_driver/driver.F90
@@ -48,6 +48,7 @@ PROGRAM DRIVER
 USE MODD_SURF_PAR, ONLY: XUNDEF
 USE MODD_TYPE_DATE_SURF
 USE MODE_THERMOS
+USE MODD_REPROD_OPER, ONLY : CQSAT
 !
 USE MODI_SUNPOS
 USE MODI_OL_READ_ATM
@@ -574,10 +575,10 @@ PROGRAM DRIVER
 LGARDEN = .FALSE.
 ! Green roofs
 LGREENROOF = .FALSE. ! greenroof activation
-ZFRAC_GR   = 0.      ! greenroof fraction on roofs
+ZFRAC_GR   = 0.0      ! greenroof fraction on roofs
 ! Solar panels
 LSOLAR_PANEL = .FALSE. ! solar panels activation
-ZFRAC_PANEL  = 0.      ! solar panels fraction on roofs
+ZFRAC_PANEL  = 0.0     ! solar panels fraction on roofs
 ! Road watering
 LPAR_RD_IRRIG= .TRUE.
 ! Natural Ventilation
@@ -922,6 +923,10 @@ PROGRAM DRIVER
   print*, 'Garden     option   is not activated but a non-zero garden    fraction is given'
   STOP
 END IF
+IF ( ZBLD(1)+ZGARDEN(1)>=1.) THEN
+  print*, 'The sum of garden and building fraction is larger than one, so road fraction is negative. Please check their values.'
+  STOP
+END IF
 IF ( (.NOT. LSOLAR_PANEL) .AND. ZFRAC_PANEL(1)>0.) THEN
   print*, 'Solar panels option is not activated but a non-zero solar panels fraction is given'
   STOP
@@ -939,6 +944,8 @@ PROGRAM DRIVER
 !
 CALL INI_CSTS
 !
+CQSAT='OLD' ! saturation is computed relative to water above 0°C, and relative to ice below 0°C
+!
 !* various thresholds
 !
 XCISMIN = 0.5     ! minimum wind shear

diff --git a/src_driver/modd_reprod_oper.F90 b/src_driver/modd_reprod_oper.F90
@@ -0,0 +1,102 @@
+!SFX_LIC Copyright 1994-2014 CNRS, Meteo-France and Universite Paul Sabatier
+!SFX_LIC This is part of the SURFEX software governed by the CeCILL-C licence
+!SFX_LIC version 1. See LICENSE, CeCILL-C_V1-en.txt and CeCILL-C_V1-fr.txt  
+!SFX_LIC for details. version 1.
+!     ######################
+      MODULE MODD_REPROD_OPER
+!     ######################
+!
+!!****  *MODD_REPROD_OPER* - declaration of ISBA parameters
+!!
+!!    PURPOSE
+!!    -------
+!       The purpose of this declarative module is to specify temporary
+!       old parameters related to the surface parameterization ISBA 
+!       to ensure reproductibility with previous oper cycle
+!
+!!
+!!      
+!!
+!!    AUTHOR
+!!    ------
+!!      B. Decharme   *Meteo France*
+!!
+!!    MODIFICATIONS
+!!    -------------
+!!      Original       04/2013                      
+!-------------------------------------------------------------------------------
+!
+!*       0.   DECLARATIONS
+!             ------------
+!
+IMPLICIT NONE
+!
+!------------------------------------------------------------------------------------
+! Old global ISBA param temporary activated in NAM_SURF_REPROD_OPER (for reproductibility)
+!------------------------------------------------------------------------------------
+!
+! * Tropical evergreen forest parameter
+!
+!XEVERG_RSMIN : old = 250. (Manzi 1993) but observations range 
+!               from 140 to 180. According to Delire et al. (1997) and 
+!               new tests over 6 local sites, 175. is recommended
+!               Should be the default after check with AROME/ALADIN
+!
+REAL             :: XEVERG_RSMIN
+!
+!XEVERG_VEG : old = 0.99 (Manzi 1993) but according to Delire et al. (1997) and 
+!             new tests over 6 local sites, 1.0 is recommended because 0.99
+!             induces unrealistic bare soil evaporation for Tropical forest
+!             Should be the default after check with AROME/ALADIN
+!
+REAL             :: XEVERG_VEG
+!
+! * Soil depth average
+!
+!CDGAVG : old         = 'ARI' Arithmetic average for all depths 
+!         recommended = 'INV' Harmonic average for all depths
+!
+ CHARACTER(LEN=3) :: CDGAVG
+!
+! * Soil depth with ISBA-DF
+!
+!CDGDIF : old         = 'SOIL' Total soil depth (d3) in Ecoclimap
+!         recommended = 'ROOT' Root depth (d2) in Ecoclimap
+!
+ CHARACTER(LEN=4) :: CDGDIF
+!
+! * wind implicitation
+!
+ CHARACTER(LEN=3) :: CIMPLICIT_WIND ! wind implicitation option
+!                                  ! 'OLD' = direct
+!                                  ! 'NEW' = Taylor serie, order 1 (recommended)
+!
+! * qsat computation
+!
+ CHARACTER(LEN=3) :: CQSAT ! qsat computation option
+!                         ! 'OLD' = do not depend on temperature
+!                         ! 'NEW' = qsat and qsati merged (recommended)
+!
+! * Charnock parameter
+!
+ CHARACTER(LEN=3) :: CCHARNOCK ! Charnock parameter option
+!                             ! 'OLD' = constant equal to XVCHRNK
+!                             ! 'NEW' = vary between 0.011 et 0.018 according
+!                             !         to Chris Fairall's data as in coare3.0
+!                             !         (recommended)
+!
+!--------------------------------------------------------------------------------
+!
+END MODULE MODD_REPROD_OPER
+
+
+
+
+
+
+
+
+
+
+
+
diff --git a/src_proxi_SVAT/PROXI_OLD/garden.F90 b/src_proxi_SVAT/PROXI_OLD/garden.F90
@@ -0,0 +1,169 @@
+!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
+! Copyright 1998-2013 Meteo-France
+! This is part of the TEB software governed by the CeCILL-C licence version 1.
+! See LICENCE, CeCILL-C_V1-en.txt and CeCILL-C_V1-fr.txt for details.
+! http://www.cecill.info/licences/Licence_CeCILL-C_V1-en.txt
+! http://www.cecill.info/licences/Licence_CeCILL-C_V1-fr.txt
+! The CeCILL-C licence is compatible with L-GPL
+!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
+!     #########
+    SUBROUTINE GARDEN(HIMPLICIT_WIND, TPTIME, PTSUN, PPEW_A_COEF, PPEW_B_COEF,       &
+                PPET_A_COEF, PPEQ_A_COEF, PPET_B_COEF, PPEQ_B_COEF,                  &
+                PTSTEP, PZ_LOWCAN,                                                   &
+                PT_LOWCAN, PQ_LOWCAN, PEXNS, PRHOA, PCO2, PPS, PRR, PSR, PZENITH,    &
+                PSW, PLW, PU_LOWCAN,                                                 &
+                PRN_GARDEN,PH_GARDEN,PLE_GARDEN,PGFLUX_GARDEN,PSFCO2,                &
+                PEVAP_GARDEN, PUW_GARDEN, PRUNOFF_GARDEN,                            &
+                PAC_GARDEN,PQSAT_GARDEN,PTS_GARDEN,                                  &
+                PAC_AGG_GARDEN, PHU_AGG_GARDEN, PDRAIN_GARDEN, PIRRIG_GARDEN         )  
+!   ##########################################################################
+!
+!!****  *GARDEN*  
+!!
+!!    PURPOSE
+!!    -------
+!
+!!call  a proxi of vegetation scheme inside TEB
+!
+!!========================================================================
+!!========================================================================
+!!========================================================================
+!!
+!! ==> YOU ARE (MORE THAN) WELCOME TO USE YOUR OWN VEGETATION SCHEME HERE
+!!
+!!========================================================================
+!!========================================================================
+!!========================================================================
+!     
+!!**  METHOD
+!     ------
+!
+!
+!!    EXTERNAL
+!!    --------
+!!
+!!
+!!    IMPLICIT ARGUMENTS
+!!    ------------------
+!!
+!!      
+!!    REFERENCE
+!!    ---------
+!!
+!!      
+!!    AUTHOR
+!!    ------
+!!
+!!	A. Lemonsu          * Meteo-France *
+!!
+!!    MODIFICATIONS
+!!    -------------
+!!    Original    05/2009
+!-------------------------------------------------------------------------------
+!
+!*       0.     DECLARATIONS
+!               ------------
+!
+USE MODD_CSTS, ONLY : XLVTT , &   ! Latent heat constant for evaporation
+                      XKARMAN     ! Von Karman constant
+USE MODE_THERMOS                  ! Function to compute humidity at saturation
+USE MODD_TYPE_DATE_SURF,    ONLY: DATE_TIME
+!
+IMPLICIT NONE
+!
+!*      0.1    Declarations of arguments
+!
+ CHARACTER(LEN=*),     INTENT(IN)  :: HIMPLICIT_WIND   ! wind implicitation option
+!                                                     ! 'OLD' = direct
+!                                                     ! 'NEW' = Taylor serie, order 1
+TYPE(DATE_TIME)     , INTENT(IN)    :: TPTIME             ! current date and time from teb
+REAL, DIMENSION(:)  , INTENT(IN)    :: PTSUN              ! solar time      (s from midnight)
+REAL, DIMENSION(:)  , INTENT(IN)    :: PPEW_A_COEF        ! implicit coefficients
+REAL, DIMENSION(:)  , INTENT(IN)    :: PPEW_B_COEF        ! for wind coupling
+REAL, DIMENSION(:)  , INTENT(IN)    :: PPEQ_A_COEF        ! implicit coefficients
+REAL, DIMENSION(:)  , INTENT(IN)    :: PPEQ_B_COEF        ! for humidity
+REAL, DIMENSION(:)  , INTENT(IN)    :: PPET_A_COEF        ! implicit coefficients
+REAL, DIMENSION(:)  , INTENT(IN)    :: PPET_B_COEF        ! for temperature
+REAL                , INTENT(IN)    :: PTSTEP             ! time step
+REAL, DIMENSION(:)  , INTENT(IN)    :: PZ_LOWCAN          ! height of atm. var. near the road
+REAL, DIMENSION(:)  , INTENT(IN)    :: PT_LOWCAN          ! temp. near the road
+REAL, DIMENSION(:)  , INTENT(IN)    :: PQ_LOWCAN          ! hum. near the road
+REAL, DIMENSION(:)  , INTENT(IN)    :: PPS                ! pressure at the surface
+REAL, DIMENSION(:)  , INTENT(IN)    :: PEXNS              ! surface exner function
+REAL, DIMENSION(:)  , INTENT(IN)    :: PRHOA              ! air density at the lowest level
+REAL, DIMENSION(:)  , INTENT(IN)    :: PCO2               ! CO2 concentration in the air    (kg/m3)
+REAL, DIMENSION(:)  , INTENT(IN)    :: PRR                ! rain rate
+REAL, DIMENSION(:)  , INTENT(IN)    :: PSR                ! snow rate
+REAL, DIMENSION(:)  , INTENT(IN)    :: PZENITH            ! solar zenithal angle
+REAL, DIMENSION(:),   INTENT(IN)    :: PSW                ! incoming total solar rad on an horizontal surface
+REAL, DIMENSION(:)  , INTENT(IN)    :: PLW                ! atmospheric infrared radiation
+REAL, DIMENSION(:)  , INTENT(IN)    :: PU_LOWCAN          ! wind near the road
+
+REAL, DIMENSION(:)  , INTENT(OUT)   :: PRN_GARDEN         ! net radiation over green areas
+REAL, DIMENSION(:)  , INTENT(OUT)   :: PH_GARDEN          ! sensible heat flux over green areas
+REAL, DIMENSION(:)  , INTENT(OUT)   :: PLE_GARDEN         ! latent heat flux over green areas
+REAL, DIMENSION(:)  , INTENT(OUT)   :: PGFLUX_GARDEN      ! flux through the green areas
+REAL, DIMENSION(:)  , INTENT(OUT)   :: PSFCO2             ! flux of CO2 positive toward the atmosphere (kg/m2/s)
+REAL, DIMENSION(:)  , INTENT(OUT)   :: PEVAP_GARDEN       ! total evaporation over gardens (kg/m2/s)
+REAL, DIMENSION(:)  , INTENT(OUT)   :: PUW_GARDEN         ! friction flux (m2/s2)
+REAL, DIMENSION(:)  , INTENT(OUT)   :: PRUNOFF_GARDEN     ! runoff over garden (kg/m2/s)
+REAL, DIMENSION(:)  , INTENT(OUT)   :: PAC_GARDEN         ! aerodynamical conductance
+REAL, DIMENSION(:)  , INTENT(OUT)   :: PQSAT_GARDEN       ! saturation humidity
+REAL, DIMENSION(:)  , INTENT(INOUT) :: PTS_GARDEN         ! radiative surface temp. (snow free)
+REAL, DIMENSION(:)  , INTENT(OUT)   :: PAC_AGG_GARDEN     ! aggreg. aeodynamic resistance for green areas for latent heat flux
+REAL, DIMENSION(:)  , INTENT(OUT)   :: PHU_AGG_GARDEN     ! aggreg. relative humidity for green areas for latent heat flux
+REAL, DIMENSION(:)  , INTENT(OUT)   :: PDRAIN_GARDEN      ! garden total (vertical) drainage
+REAL, DIMENSION(:)  , INTENT(OUT)   :: PIRRIG_GARDEN      ! garden summer irrigation rate
+!
+!
+!*      0.2    Declarations of local variables
+!
+!
+!-------------------------------------------------------------------------------
+!
+!*      1.     Proxi model based on a fixed Bowen ratio
+!              ----------------------------------------
+!
+!* albedo fixed to 0.15
+PRN_GARDEN(:) = (1.-0.15) * PSW(:)
+!
+!* Bowen ratio fixed to 0.25
+PH_GARDEN (:) = 0.2 * PRN_GARDEN(:)
+PLE_GARDEN(:) = 0.8 * PRN_GARDEN(:)
+!
+!* Conduction heat flux is neglected
+PGFLUX_GARDEN(:) = 0.
+!
+!* CO2 flux is neglected
+PSFCO2(:) = 0.
+!
+!* evaporation
+PEVAP_GARDEN(:) = PLE_GARDEN(:) / XLVTT
+!
+!* Friction flux: assumes neutral formulation with roughness length of 0.1m
+PUW_GARDEN(:) = - (XKARMAN/LOG(PZ_LOWCAN(:)/0.1))**2 * PU_LOWCAN(:)**2
+!
+!* Aerodynamical conductance: neglected because used further only for
+!  implicitation of canyon air temperature when the heat flux depends on the
+!  surface temperature
+!
+PAC_GARDEN(:) = 0.
+!
+!* surface saturation humidity
+PQSAT_GARDEN(:) = QSAT(PT_LOWCAN(:),PPS(:))
+!
+!* Surface temperature : set equal to air temperature
+PTS_GARDEN(:) = PT_LOWCAN(:)
+!
+!* aerocynamical conductance for latent heat and surface humidity
+PAC_AGG_GARDEN(:) = 0.    ! neglected (latent flux does not depend on surface humidity)
+PHU_AGG_GARDEN(:) = 0.8   ! surface humidity set to 80%
+!
+!* garden hydrological diagnostics
+PRUNOFF_GARDEN(:) = 0.    ! garden surface runoff
+PDRAIN_GARDEN (:) = 0.    ! garden total (vertical) drainage
+PIRRIG_GARDEN (:) = 0.    ! garden irrigation during time step
+!-------------------------------------------------------------------------------
+!
+!
+END SUBROUTINE GARDEN