Main changes in Surfex v7_3

  • ISBA coupling with the TOPMODEL library (B. Vincendon)
TOPODYN option in Surfex v7.3
  • Chemical part: NO fluxes (C. Delon)
    • if LCH_NO_FLUX (NAM_CH_ISBAn), NO flux is caculated in ISBA part
    • if LCH_NO_FLUX (NAM_CH_ISBAn), NO flux is written in NOFLUX output
    • with LCH_NO_FLUX, PH and FERT fields are given at PGD step (NAM_ISBA)
  • Chemical part: SNAP emissions (S. Queguiner)
    • to use other types of cadastres for the chemical emissions (smaller input files)
  • Contribution from GMGEC (B. Decharme)
Contribution from GMGEC for 7.3
  • New initialization for ALBEDO (D. Carrer & C. Meurey)
    • Like for the LAI, soil and vegetation ALBEDO are derived from satellite products (SAF) and defined for each cover and each vegetation type in ECOCLIMAP.
  • Contribution from CEN (M. Lafaysse)
    • possibility during PREP to initialize TG, WG and WGI with NETCDF files in ISBA part.
    • modification of the management of CROCUS optimal thicknesses for very thick coverings of snow.
    • partial merge of routines SNOWCRORAD and SNOWCROALB to supress the code that was duplicated
    • possibility to simulate glaciers with crocus (albedo forced in the 3 bands and different roughness for the apparent ice)
    • all constants concerning albedo routines are put in modd_snow_par to be declared only once.
    • new constants for ice albedo and roughness can be defined in namelist NAM_SURF_CSTS. 
      Default values :
! Roughness length ratio between ice and snow
XZ0ICEZ0SNOW=10.
! 3 bands spectral albedo for glacier ice (CROCUS)
! Default values from Lejeune et al 2009 (Zongo, Bolivia)
XALBICE1=0.38
XALBICE2=0.23
XALBICE3=0.08
! Density threshold for ice detection kg.m-3
XRHOTHRESHOLD_ICE=850.
    • variable PSNOWAGE is now an age and no longer a date (what allows to do easier spinup when CROCUS is activated)
    • print option for CROCUS debugging, conditioned by the definition of environment variables
User’s guide for CROCUS debugging
  • Modifications for flake fluxes and 1D ocean model (P. Le Moigne & P. Peyrillé)
    • possibility to activate ECUME in Flake (CFLK_FLUX=’ECUME’ in NAM_FLAKEn) + associated options LPRECIP and LPWEBB
    • add of relaxation for the ocean model: temperature, salinity, currents
    • possibility to correct fluxes with a coefficient given in namelist (XCORFLX)
    • possibility to set surface fluxes to 0
    • flag diapycnal : calculation of mixing coefficients at internal waves and wind shearing
  • Contribution for the assimilation part (SODA) (Trygve Aspelien)
    • changes for the snow assimilation
  • New TEB (V. Masson, A. Lemonsu, G. Pigeon, C. De Munck)
    • inclusion of BEM scheme (Building Energy Model)
      BEM documentation
    • add of patches in TEB to describe the roads orientation
    • add of GREENROOF scheme
    • suppression of PATCH argument in GARDEN variables
  • Parallelization of the OFFLINE driver (F. Suzat & S. Faroux)
    • parallelization possible with MPI and / or OPEN-MP
    • possibility to unactivate the MPI implementation (export VER_MPI=NOMPI during configure part of installation)
  • For Méso-NH (M. Tomasini):
    • Add an index in the second dimension of the ISBA array variables for BILIN interpolation routine to not bug in case 2D (this is not the more beautiful method; the BILIN routine should better be adapted)
  • Side developments:
    • all fields names have now 12 characters and no longer 16.
      => In NETCDF forcing files, FORC_TIME_STEP becomes FRC_TIME_STP. To change it easily, do:
      ncrename -v FORC_TIME_STEP,FRC_TIME_STP FORCING.nc

Also in this section