Changes in Namelists

NEW_PREP

none

NEW_PREP_HIRLAM

  • in PREP_ISBA_SNOW :
    • LSWEMAX : logical switch to set an upper limit on initial snow water equivalent (default : false)
    • XSWEMAX : upper limit of initial snow water equivalent. default : 500 mm

=> to be tested

  • in NAM_ASSIM :
    • NPRINTLEV : verbosity (0, 1, 2 ...). Defautl : 0, lower level
    • LASSIM : Reads diagnostics needed for assimilation as "observations". Default false

=> no assimilations tests for the moment. Add this from Trygve.

  • in NAM_PREP_ISBA
    • Extrapolations of fields in initalization where there is a mismatch in input and output land-sea mask. Default false
      • LEXTRAP_TG : Ground temperature
      • LEXTRAP_WG : Ground water
      • LEXTRAP_WGI : Ground ice
      • LEXTRAP_SN : Snow

=> runs with prep BUFFER. Not tested, for the moment.

  • in NAM_DIAG_SURF_ATMn
    • LT2MMW : Alternative weighting of grid average T2M giving more weight to the land tile. Default false

=> to be tested

V8_GMAP

  • In NAM_ISBAN :

LPERTSURF : if .True. modification of surface fluxes for ensemble forecasting (default : .False.)

  • In NAM_SEAFLUX :

LPERTFLUX : if .True. modification of surface fluxes for ensemble forecasting (default : .False.)

=> to be tested

V8_BV : ISBA-TOP

no changes

V8_VM : TEB

  • in NAM_TEB
    LSOLAR_PANEL logical Default : F
    solar panels on roofs
  • in NAM_TEBn
    XDT_RES real Default 0.
    XDT_OFF real Default 0.
    XDT_RES difference of ordered temperature (heating and cooling) when the building (residential) is not occupied : day (K)
    XDT_OFF difference of ordered temperature (heating and cooling) when the building (offices) is not occupied : nights and WE (K)
  • in NAM_DATA_TEB
    XUNIF_RESIDENTIAL real between 0. and 1. Default 1E+20
    XUNIF_FRAC_PANEL real between 0. and 1. Default 1E+20
    XUNIF_EMIS_PANEL real between 0. and 1. Default 1E+20
    XUNIF_ALB_PANEL real between 0. and 1. Default 1E+20
    XUNIF_EFF_PANEL real between 0. and 1. Default 1E+20

XUNIF_RESIDENTIAL :fraction of residential use (used only for solar panels and occupying schedule)
XUNIF_FRAC_PANEL= fraction of solar panels on roofs
XUNIF_EMIS_PANEL= emissivity of solar panels
XUNIF_ALB_PANEL= albedo of solar panels
XUNIF_EFF_PANEL= efficiency of photovoltaic solar panels
Remarks: the efficiency of thermal solar panels (for the sanitary warm water) is not in the namelist and was fixed to 0.60.
If the buildings in the grid point are purely of residential use, the Sanitary Warm Water is favoured, but after photovoltaïcs.
If the buildings are of other use, there is only photovoltaïcs.
If the use is mixed (for example 0.5), the sharing out is proportional.

  • New NAMELIST NAM_DATA_TEB_IRRIG
    watering of gardens and greenroofs.
    Activated in a specific namelist.
    By default, the watering is non-existent.
    The amounts for watering are in litres/day (or kg/m2/day).
    Time slots and months for are specified by user.
    The watering occuring in ISBA, there is a new input variable in isba.F90.
    There is just one additionnal namelist in the OPTIONS.nam file for the irrigation specifications (garden or green roof irrigation, road watering).
    Default values are equal to XUNDEF (no irrigation or watering).
    &NAM_DATA_TEB_IRRIG
    XUNIF_GD_START_MONTH real default 1.E+20
    XUNIF_GD_END_MONTH real default 1.E+20
    XUNIF_GD_START_HOUR real default 1.E+20.
    XUNIF_GD_END_HOUR real default 1.E+20
    XUNIF_GD_24H_IRRIG real default 1.E+20
    XUNIF_GR_START_MONTH real default 1.E+20
    XUNIF_GR_END_MONTH real default 1.E+20
    XUNIF_GR_START_HOUR real default 1.E+20
    XUNIF_GR_END_HOUR real default 1.E+20
    XUNIF_GR_24H_IRRIG real default 1.E+20
    XUNIF_RD_START_MONTH real default 1.E+20
    XUNIF_RD_END_MONTH real default 1.E+20
    XUNIF_RD_START_HOUR real default 1.E+20
    XUNIF_RD_END_HOUR real default 1.E+20
    XUNIF_RD_24H_IRRIG real default 1.E+20
    • XUNIF_GD_START_MONTH start month for irrigation of gardens (month)
    • XUNIF_GD_END_MONTH end month for irrigation of gardens (month)
    • XUNIF_GD_START_HOUR start hour for irrigation of gardens (solar hour)
    • XUNIF_GD_END_HOUR end hour for irrigation of gardens (solar hour)
    • UNIF_GD_24H_IRRIG mean rate of irrigation for gardens per 24h (kg/m2)
    • XUNIF_GR_START_MONTH start month for irrigation of greenroofs (month)
    • XUNIF_GR_END_MONTH end month for irrigation of greenroofs (month)
    • XUNIF_GR_START_HOUR start hour for irrigation of greenroofs (solar hour)
    • XUNIF_GR_END_HOUR end hour for irrigation of greenroofs (solar hour)
    • XUNIF_GR_24H_IRRIG mean rate of irrigation for greenroofs per 24h (kg/m2)
    • XUNIF_RD_START_MONTH start month for irrigation of roads (month)
    • XUNIF_RD_END_MONTH end month for irrigation of roads (month)
    • XUNIF_RD_START_HOUR start hour for irrigation of roads (solar hour)
    • XUNIF_RD_END_HOUR end hour for irrigation of roads (solar hour)
    • XUNIF_RD_24H_IRRIG mean rate of irrigation for roads per 24h (kg/m2)

The start month number can be greater than the end month number, (example: irrigation starting in december and finishing in February).
Similarly, the start hour can be larger than the end hour (then irrigation starts in the evening and ends in the morning the next day, every day).

There are also additionnal output fields :
RUNOFFC_TW : cumulated aggregated runoff for town (kg/m2)
RUNOFFC_RD : cumulated road surface runoff (kg/m2 road)
RUNOFFC_RF : cumulated aggregated roof runoff (kg/m2 roof)
IRRIGC_RD : cumulated road irrigation (kg/m2 road)
RUNOFFC_GD : cumulated garden surface runoff (kg/m2 garden)
DRAINC_GD : cumulated garden surface drainage (kg/m2 garden)
IRRIGC_GD : cumulated garden irrigation (kg/m2 garden)

for solar panels

PHOTC_BLD ! photovoltaics production cumuled on the length of simulation ( J/m2 of building at floor)
THERC_BLD ! sanitary warm water production cumuled on the length of simulation ( J/m2 of building at floor) Solar panels cover (give shadow and modify the incoming infra-red radiation below) roofs, even if they are also fitted with greenroofs.

CHANGE IN DEFAULTS in NAM_TEBn:
CZ0H = "KAND07" old default value was "MASC92"
CCH_BEM = "DOE-2" old default value was " ".
The option corresponding to this previous default is now called "ROW30". The value " " is no more used.

V8_lafaysse : CROCUS

In NAM_ZS, the choice of YZSFILETYPE to read altitude is extended to " NETCDF ".
It is now possible to force the slope from an external NETCDF file with YSLOPE and YFILESLOPETYPE.
It is also possible to compute the explicit slope during PGD (LEXPLICIT_SLOPE=T) (available only on rectangular grids).

Default values:

&NAM_ZS
XUNIF_ZS = 1.0E+20
YZS = "gtopo30"
YZSFILETYPE = "DIRECT"
COROGTYPE = "MAX"
XENV = 0.2
LIMP_ZS = F
YSLOPE = ""
YFILESLOPETYPE = ""
LEXPLICIT_SLOPE = F
/

In NAM_IO_OFFLINE, two new options are available (LSHADOWS_SLOPE, LSHADOWS_OTHER) to account for shadows respectively of the slope itself or of the surrounding mountains. These options can be set to TRUE only on a rectangular domain, and with XIO_FRAC = 1. and YALG_MPI = "LIN". The PGD has to have been computed with the option LEXPLICIT_SLOPE=T.

Default values:

&NAM_IO_OFFLINE
CSURF_FILETYPE = "NC"
CTIMESERIES_FILETYPE = "NETCDF"
CFORCING_FILETYPE = "ASCII"
CPGDFILE = "PGD"
CPREPFILE = "PREP"
CSURFFILE = "SURFOUT"
LPRINT = T
LRESTART = T
LINQUIRE = T
NHALO = 2
NSCAL = 0
XTSTEP_SURF = 300.
XTSTEP_OUTPUT = 10800.
LDIAG_FA_NOCOMPACT = F
LSET_FORC_ZS = F
LWRITE_COORD = T
LOUT_TIMENAME = F
LLIMIT_QAIR = F
LSHADOWS_SLOPE = F (new in V8)
LSHADOWS_OTHER = F (new in V8)
NB_READ_FORC = 0
LLAND_USE = F
NPROMA = 0
NI = 0
NJ = 0
XIO_FRAC = 1.
YALG_MPI = "LIN"
/

In NAM_PREP_ISBA_SNOW, it is now possible to initialize either snow water equivalent (XWSNOW) either snow depth (XZSNOW).
It is also possible to initialize snow with a non-zero liquid water content (XLWCSNOW).

Default values:
&NAM_PREP_ISBA_SNOW
CSNOW = "D95"
NSNOW_LAYER = 1
CFILE_SNOW = ""
CTYPE_SNOW = ""
CFILEPGD_SNOW = ""
CTYPEPGD_SNOW = ""
LSNOW_IDEAL = F
LSNOW_FRAC_TOT = F
XWSNOW = 1.0E+20
XZSNOW = 1.0E+20 (new in V8)
XTSNOW = 1.0E+20
XLWCSNOW = 0 (new in V8)
XRSNOW = 1.0E+20
XASNOW = 1.0E+20
XSG1SNOW = 1.0E+20
XSG2SNOW = 1.0E+20
XHISTSNOW = 1.0E+20
XAGESNOW = 1.0E+20

The old namelist NAM_CROCUSn is dissociated into two new namelists:

– The namelist NAM_SURF_SNOW_CSTS is dedicated to parameters.
&NAM_SURF_SNOW_CSTS
XZ0ICEZ0SNOW = 10.
XRHOTHRESHOLD_ICE = 850.
XALBICE1 = 0.38
XALBICE2 = 0.23
XALBICE3 = 0.08
XVAGING_NOGLACIER = 60.
XVAGING_GLACIER = 900.

– The namelist NAM_ISBA_SNOWn describes physical options:

&NAM_ISBA_SNOWn
CSNOWMETAMO = "B92" ! metamorphism scheme (new in V8)
CSNOWRAD = "B92" ! radiative transfer scheme. (new in V8)
LSNOWDRIFT = T
LSNOWDRIFT_SUBLIM = F
LSNOW_ABS_ZENITH = F ! if true solar absorption depends on zenithal angle (useful in polar regions)(new in V8)

Available physical options in Crocus:
Metamorphism scheme
CSNOWMETAMO=B92 Brun et al 1992
CSNOWMETAMO=C13 Carmagnola et al 2014
CSNOWMETAMO=T07 Taillandier et al 2007
CSNOWMETAMO=F06 Flanner et al 2006

Radiative transfer scheme
CSNOWRAD=B92 Brun et al 1992
CSNOWRAD=TAR TARTES (Libois et al 2013)
CSNOWRAD=TA1 TARTES with constant impurities
CSNOWRAD=TA2 TARTES with constant impurities as function of ageing

V8_GMGEC_MEB

  • NAM_IO_OFFLINE (modn_io_offline.F90)

Add LADAPT_SW (default = .false.) to activate the simple coherence between solar zenithal angle and radiation coded for TEB. The default is FALSE because this coherence should be computed more realistically. In offline.F90, the code is now:

IF(LADAPT_SW)THEN
 !
 ! coherence between solar zenithal angle and radiation
 ! when solar beam close to horizontal -> reduction of direct radiation to
 ! the benefit of scattered radiation
 ! when pi/2 - 0.1 < ZENITH < pi/2 - 0.05 => weight of direct to scattered radiation decreases linearly with zenith
 ! when pi/2 - 0.05 < ZENITH => all the direct radiation is converted to scattered radiation
 ! coherence between solar zenithal angle and radiation
 !
 ZCOEF(:) = (XPI/2. - XZENITH(:) - 0.05)/0.05
 ZCOEF(:) = MAX(MIN(ZCOEF,1.),0.)
 DO JLOOP=1,SIZE(XDIR_SW,2)
   XSCA_SW(:,JLOOP) = XSCA_SW(:,JLOOP) + XDIR_SW(:,JLOOP) * (1 - ZCOEF)
   XDIR_SW(:,JLOOP) = XDIR_SW(:,JLOOP) * ZCOEF(:)
 ENDDO
 !
ENDIF
!
  • NAM_SURF_CSTS (ini_surf_csts.F90)

New variables :

XALBSEAICE = albedo over sea ice (default depends of LREPROD_OPER flag, see after)

XZ0FLOOD = flood z0 (default= 0.0002)

XALBWATSNOW = snow albedo over water bodies or lakes (default depends of LREPROD_OPER flag, see after)

  • New NAM_REPROD_OPER (ini_surf_csts.F90)

CIMPLICIT_WIND = wind implicitation (previously in NAM_SURF_ATM)

New variables :

CDGDIF =

  1. ’ROOT’ : d2 root depth as Soil depth reference for ISBA-DF (recommended)
  2. ’SOIL’ : d3 soil depth as Soil depth reference for ISBA-DF

LREPROD_OPER = flag to activate reproductibility for SURFEX OPER
Allow to put old value/flag for some variables in order to ensure reproducibility (see ini_surf_csts.F90 for more details)

XEVERG_RSMIN = Rsmin for tropical forest

XEVERG_VEG = Veg fraction for tropical forest

CDGAVG =

  1. ’INV’ : compute harmonic average for Soil depth (recommended)
  2. ’ARI’ : old computation of Soil depth using arithmetic average

CQSAT =

  1. ’NEW’ : qsat computation accounting for phase change. If temperature is <0.°C compute it using ice properties instead of using water properties. (recommended)
  2. ’OLD’ : qsat computation only accounting for water properties.

CCHARNOCK =

  1. ’NEW’ : charnock number vary between 0.011 et 0.018 according to Chris Fairall’s data as in coare3.0 (recommended)
  2. ’OLD’ : charnock number = XVCHRNK

The defaults of some variables depend on LREPROD_OPER

Flag/variables LREPROD_OPER = .F. LREPROD_OPER = .T.
XEVERG_RSMIN 175 250
XEVERG_VEG 1 0.99
CDGAVG INV ARI
CQSAT NEW OLD
CCHARNOCK NEW OLD
XALBWAT 0.06 0.135
XALBSEAICE 0.71 0.85
XALBWATICE 0.40 0.85
XALBWATSNOW 0.60 0.85
XEMISWAT 0.96 0.98
XEMISWATICE 0.97 1.0
XEMISSN 0.99 1.0
  • NAM_SURF_ATM (modn_surf_atm.F90)

LRW_PRECIP replaced by LCPL_GCM

New LVSHIFT_LW and LVSHIFT_PRCP flags to activate/deactivate vertical shift for LongWave radiations and Precip (default=.false.). See forcing_vert_shift.F90 for more details.

CIMPLICIT_WIND is transferred in new namelist NAM_REPROD_OPER (see ini_surf_csts.F90)

  • New NAM_ISBA_AGSn (modn_isban.F90)

New namelist NAM_ISBA_AGSn to read LNITRO_DILU key in order to activate/desactivate nitrogen dilution fct of CO2 (Calvet et al. 2008). Default = .FALSE.
These changes have been also added in modn_teb_vegn.F90.

  • NAM_SPINUP_CARBn (modn_isban.F90)

Add XCO2_START and XCO2_END. During soil carbon spinup with ISBA-CC:

  1. Atmospheric CO2 concentration fixed to Pre-industrial CO2 consentration XCO2_START
  2. Atmospheric CO2 concentration rampin up from XCO2_START to XCO2_END
    These quantities must be prescribed in namelits in ppm. For more detail see with Christine Delire (GMGEC/ASTER).
  • NAM_SGH_ISBAn (modn_isban.F90)

CSOC character is replaced by LSOC logical (default .false.).

CTOPREG (never used), LTRIP and LFLOOD have been deleted due to the new coupling with TRIP via OASIS-MCT. These changes have been also added in modn_teb_greenroofn.F90

  • NAM_FLAKEn (modn_flaken.F90)
    – 
    XICHCE, LPRECIP, LPWEBB have been deleted because ecume can not be call for lake.

These changes have been also added in modn_teb_greenroofn.F90

Add LSKINTEMP flag to use or not the skin temperature computation (default = .false)

NAM_CH_FLAKEn was added to put CCH_DRY_DEP key.

  • NAM_PREP_SEAFLUX (modn_prep_seaflux.F90)

CSEAICE_SCHEME = name of sea ice scheme (NONE (default) or GELATO)

XSSS_UNIF = uniform sea surface salinity concentration used with ECUME6 or Gelato.

XSIC_UNIF = uniform sea ice covert fraction

  • NAM_SEAFLUXn (modn_seafluxn.F90)
    – 
    New option in CSEA_FLUX = ECUME6 to activate new ecumev6 (see with Sophie Belamari GMGEC/NEMO)

Add CINTERPOL_SSS to perform quadratic interpolation of annual or monthly SSS (sea surface salinity) used by ECUME6 and/or Gelato

Add NZ0 = to choose PZ0SEA formulation in ECUME6 (default = 0)

  1. NZ0 = 0 (ARPEGE formulation)
  2. NZ0 = 1 (Smith (1988) formulation)
  3. NZ0 = 2 (Direct computation using the stability functions)
  • New NAM_SEAICEn (modn_seafluxn.F90)
    – 
    New namelist to activate/deactivate sea ice sheme as GELATO. For more detail see Gelato in Surfex documentation.

LDIAG_SEAICE = flag for SeaIce diag variables

CINTERPOL_SIC= Quadratic interpolation of annual or monthly Sea Ice Cover

CINTERPOL_SIT= Quadratic interpolation of annual or monthly Sea Ice Thicknesses

XFREEZING_SST= Value marking frozen sea in SST data (default = -1.8°C)

XSIC_EFOLDING_TIME= For damping of SIC (days) (default = 0.0)

XSIT_EFOLDING_TIME= For damping of SIT (days) (default = 0.0)

XSEAICE_TSTEP= Sea ice model time step

XCD_ICE_CST= Fixed Turbulent exchanges coefficients for sea ice. If equal to 0.0 (default) then the normal (or aladin) formulation is activated to compute Cd and Ch over sea ice.

XSI_FLX_DRV= Derivative of fluxes on seaice w.r.t to the temperature (W m-2 K-1) (default -20)

  • New NAM_SFX_LAND_CPL, NAM_SFX_LAKE_CPL, NAM_SFX_SEA_CPL

These new namelists are used to control the SURFEX/OASIS-MCT coupling. For more details see SURFEX/OASIS-MCT documentation and modn_sfx_oasis.F90.

NAM_SFX_LAND_CPL

XTSTEP_CPL_LAND = Coupling time step for land

CRUNOFF = Name of Surface runoff variable from SFX to TRIP

CDRAIN = Name of Deep drainage variable from SFX to TRIP

CCALVING = Name of Calving flux variable from SFX to TRIP

CRECHARGE = Name of groundwater recharge variable from SFX to TRIP

CPFLOOD = Name of Floodplains recipitation interception variable from SFX to TRIP

CEFLOOD = Name of Floodplains evaporation variable from SFX to TRIP

CIFLOOD = Name of Floodplains infiltration variable from SFX to TRIP

CWTD = Name of water table depth variable from TRIP to SFX

CFWTD = Name of grid-cell fraction of water table rise variable from TRIP to SFX

CFFLOOD = Name of Floodplains fraction variable from TRIP to SFX

CPIFLOOD = Name of Flood potential infiltration variable from TRIP to SFX

For NAM_SFX_LAKE_CPL, NAM_SFX_SEA_CPL, please see modn_sfx_oasis.F90

  • New NAM_MEB_ISBA (in read_nam_pgd_isba_meb.F90)

This new namelist allow to activate MEB.

LMEB_PATCH = vector that activates MEB over selected patches. For example if you use 12 patches: LMEB_PATCH = F,F,F,T,T,T,F,F,F,F,F,F (= only forest patches)

LFORC_MEASURE =

  1. False (default)
  2. True (Forcing data from observations)
  • OASIS: add of coupling namelists NAM_SFX_LAND_CPL, NAM_SFX_SEA_CPL, NAM_SFX_LAKE_CPL
  • GELATO:

Add of CSEAICE_SCHEME in NAM_PREP_SEAFLUX .

Add of new namelist NAM_SEAICE (cf user’s guide).

V8_MESONH

No changes.

V8_ASSIM

  • add of NHALO_PREP in NAM_PREP_SURF_ATM to define the size of the halo where points will be taken for extrapolation. Default: 2
  • YFILESLOPETYPE becomes YSLOPEFILETYPE to harmonize, in NAM_ZS
  • XZSNOW_GR and XLWCSNOW_GR are added in NAM_PREP_GREENROOF_SNOW , because they are used in ES case.
  • CTIMESERIES_FILETYPE = ’NC’ in NAM_IO_OFFLINE (new option)

Also in this section