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 landsea mask. Default false
 LEXTRAP_TG : Ground temperature
 LEXTRAP_WG : Ground water
 LEXTRAP_WGI : Ground ice
 LEXTRAP_SN : Snow
 Extrapolations of fields in initalization where there is a mismatch in input and output landsea mask. Default false
=> 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 : ISBATOP
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 nonexistent.
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 infrared 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 = "DOE2" 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 nonzero 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 =
 ’ROOT’ : d2 root depth as Soil depth reference for ISBADF (recommended)
 ’SOIL’ : d3 soil depth as Soil depth reference for ISBADF
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 =
 ’INV’ : compute harmonic average for Soil depth (recommended)
 ’ARI’ : old computation of Soil depth using arithmetic average
CQSAT =
 ’NEW’ : qsat computation accounting for phase change. If temperature is <0.Â°C compute it using ice properties instead of using water properties. (recommended)
 ’OLD’ : qsat computation only accounting for water properties.
CCHARNOCK =
 ’NEW’ : charnock number vary between 0.011 et 0.018 according to Chris Fairall’s data as in coare3.0 (recommended)
 ’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 ISBACC:
 Atmospheric CO2 concentration fixed to Preindustrial CO2 consentration XCO2_START
 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 OASISMCT. 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)
 NZ0 = 0 (ARPEGE formulation)
 NZ0 = 1 (Smith (1988) formulation)
 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 m2 K1) (default 20)
 New NAM_SFX_LAND_CPL, NAM_SFX_LAKE_CPL, NAM_SFX_SEA_CPL
These new namelists are used to control the SURFEX/OASISMCT coupling. For more details see SURFEX/OASISMCT 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 gridcell 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 =
 False (default)
 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)