SURFEX

Documentation provided by developers

Minvielle — 2026-01-27T10:41:35Z

Preparation of the SURFEX Version 8

2016-05-03T14:02:58Z

This page is devoted to the building of the surfex version 8:

A reminder of the general recommendations to mergers /Updated 6 Feb. 2014
Description of the changes /Updated 5 Jun. 2014
Tests with the test database /Updated 5 Jun. 2014
Changes in namelists /Updated 5 Jun. 2014
The draft of the schedule for the building of SURFEX V8

Scheduled contributions (autumn 2013):

The different contributions concern:

branch NEW_PREP from MOSAYC

SODA from HIRLAM (branch new_soda)
OI_MAIN + PREP from GMAP
Ensemble prediction in AROME from MICADO
VARASSIM from VEGEO (branch DEV_VARASSIM)

CROCUS from CEN (branch lafaysse_dev_for_test)
TOPMODEL from MICADO (branch isbatop_df_Npatches)
TEB from TURBAU (branch TEB_DEV_1514)
rotated grid from P. Samuelsson

a set of modifications from ALADIN

SURFEX in MESONH V4_10_0 from LA (branch BUG)
technical modifications from HIRLAM
suppression of global variables from CERFACS (branch GABRIEL_CERFACS)

When all contributions will be defined, we will set:

Changes in Namelists

2014-02-14T15:22:05Z

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 =

'ROOT' : d2 root depth as Soil depth reference for ISBA-DF (recommended)
'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 =

'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 ISBA-CC:

Atmospheric CO2 concentration fixed to Pre-industrial 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 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)

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 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 =

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)

Results of tests and analyses

2014-02-06T11:14:05Z

This section summarizes the test results (from the database, based originally on NEW_PREP).

NEW_PREP_HIRLAM

From the database :
NO_ECO_TEB_OPT_FILES Emax/=0: 4
The fix in vegetation_update.F90 had a side effect on the update of TEB_garden (with no scientific practical effect). The TEB_garden update will be corrected by V. Masson.

PREP_EXTE_FILE_ASCII_ISBA_DIF83L_NEW Emax/=0:112
Some differences appeared when old version was run with lfi and new version in netcdf

PREP_EXTE_FILE_ASCII_ISBA_SNOW_3LD95_NEW Emax/=0: 104
PREP_EXTE_FILE_ASCII_ISBA_SNOW_3LEBA_NEW Emax/=0: 102
Bug present in both version, corrected by rev1835 in the trunk and corrected in V8_GMAP branch

V8_GMAP
tests made by adding the 3 files of the trunk revision 1835 to NEW_PREP_HIRLAM. Tests not done in parallel version.

– > no changes

introduction of bug update corresponding to revision 1835 in the trunk

revision 1835: corrections of mode_soil for W33_FUNC (field capacity at -0.33bar for ISBA-DF only) (B. Decharme) and correction for snow fields in prep (when passing from a multilayer to a single layer model and corection of the vertical shift, E. Martin). Files : mode_soil.F90, prep_snow_extern.F90, prep_ver_snow.F90

– > changes in the ISBA-DF and ISBA PREP SNOW tests cases

TOPMODEL - V8_BV

The script database has been corrected for TOP MODEL (and updated in the trunk)

change in the script script_cas_parts.sh
two new scripts in the "TESTS" directory script_exec_topd.sh and script_exec_topd_restart.sh
new files in the directory "TESTS/CAS_PART/TOPD/files/"

With these corrections :

– > no changes

TEB - V8_VM

Summary :
Limited changes are observed on scientific results due to :
a correction of a bug in sunpos.F90 ( 1E-12)
The implicitation of the coupling of the snow and the roof (important changes at the beginning of the simulation because of an important difference between Tsnow and troof at the beginning of the simulation

other differences comme from bug corrections :
NO_ECO_SOME_TEB_OPT
greenroof (diagnostic values where no greenroof are present)

****************************
* PHASE 1 de la VALIDATION *
****************************

Tests sans implicitation du couplage neige/toit et sans correction de sunpos.F90 :

– 4 diagnostiques sont toujours différents (HU2MMIN/MAX T2MMIN/MAX) : car avant ils n'étaient pas réellement calculés dans TEB (seulement la valeur instantanée était prise), ni dans Flake.

– 2 diagnostiques differents (T2MMAX_ISBA et HU2MMAX_ISBA) qui étaient égaux à XUNDEF dans l'ancienne version quand on utilise canopy pour ISBA, et qui sont corrects maintenant (correction de bug dans average_diag_isban.F90).

– jusqu' à 40 diagnostiques diff&ecaute;rents quand il y a greenroof (correction de bug: quand GREENROOF est activ&ecaute; mais que la fraction de GREENROOF vaut zero sur certains point (ce qui est le cas dans la base de tests, où elle vaut zero pour tous les cas tests de TEB), les diagnostiques sont mis àÃƒƒ‚Ãƒ‚Ã‚Â XUNDEF pour ces points).

– cas NO_ECO_SOME_TEB_OPT (185 différences) : correction d'un bug, déjàÃƒƒ‚Ãƒ‚Ã‚Â pointé par Stéphanie, qui faisait que l'on ne pouvait pas utiliser ecoclimap pour GARDEN et NAM_DATA_ISBA pour initialiser ISBA pour le tile NATURE. En effet, les paramètres d'ISBA-NATURE étaient alors affectés aussi à ISBA-GARDEN. Et ceci ne pouvait marcher que si les nombres de points des deux tilesétaient identiques (cas de ce cas test). Sinon, c'est le 'segmentation fault' assuré. Maintenant, dans la routine convert_patch_isba, utilisée aussi par garden, seul le tile isba peut utiliser les paramètres issus de NAM_DATA_ISBA. Les paramètres pour GARDEN sont alors calculés avec ecoclimap. Ceci induit que les résultats du cas test sont très différents, puisque GARDEN est maintenant (correctement) initialisé par ecoclimap et non en recopiant les champs d'ISBA du tile NATURE.

– cas NO_ECO_TEB_OPT_FILES : 4 champs PGD qui ne sont plus updatés pour des points non pertinents (pas de jardins) dans la vegetation_update.F90. Les résultats de simulation ne sont pas affectés.

****************************
* PHASE 2 de la VALIDATION *
****************************

Tests AVEC implicitation du couplage neige/toit et sans correction de sunpos.F90 :

L'implicitation entre neige et toit de l'equation de conduction de la chaleur dans TEB modifie les variables pronostiques de TEB, et donc un nombre important de variables:
– toutes la variables liées au toit quand l'on n'a ni BEM (pas de couplage entre le toit et le mur via le bâtiment) ni CANOPY (pas de couplage via l'air dans la canyon).
– toutes les variables de TEB si l'on a soit BEM, soit CANOPY, soit les deux.

L'amplitude de modification maximum a lieu pour la température de la couche supérieure du toit (TROOF1). Dans les simulations de la base de test, les conditions initiales provoquent un choc important sur les températures de neige et toit, du fait des 3 conditions suivantes:
– il y a beaucoup de neige
– la température initiale de la neige est de 271.52K
– la température initiale de la couche superficielle du toit, en contact immédiat avec la neige, est de 283.38K

Ceci induit un refroidissement très rapide de la température du toit. Au premier pas de temps, la différence entre les deux versions du code est d'environ 1°C et converge au cours du temps. Plusieurs cas se présentent :
– soit on n'a pas de toiture végétalisée (fraction nulle), ce qui est le cas par exemple dans les tests sur TEB, alors en 3h (pas de temps des diagnostiques de sortie dans la base de test), et alors la température des deux simulations a convergé et elle reste inférieure àÃƒƒ‚Ãƒ‚Ã‚Â 0.1°C quelle que soit la configuration de simulation. Cf figure evol_TROOF1_sans_GREENROOF.png.
– soit on a 50% de toits végétalisés (cas avec ECO_ALL par exemple), et alors la température du toit ne peut pas décroître aussi vite (du fait que la moitié est occupée par la toiture végétalisée, et que le flux de conduction depuis ce toit empêche le refroidissement). Alors l'écart entre les deux version vaut 0.7°C (valeur à 3h de simulation), puis décroit au cours du temps jusqu'àÃƒƒ‚Ãƒ‚Ã‚Â converger. Cf figure evol_TROOF1_avec_moitie_GREENROOF.png.
– soit on a 100% de toit végétalisé (cas hors base de test), et les deux versions du code sont bien identiques (pas de neige directement sur le toit).

****************************
* PHASE 3 de la VALIDATION *
****************************

Tests de la nouvelle version (même version en _old et _new dans la base de tests) en multiprocesseurs.
8 cas posent problèmes:

6 viennent sans doute d'un problème de code. Cependant, seulement 6 parmi tous les tests ISBA posent problème:
ISBA_2L_ALL4_SGH_OMP4 Emax/=0: 112
ISBA_2L_ALL4_SGH_MPI2OMP2 Emax/=0: 112
ISBA_2L_ALL4_SGH_CANOPY_DRAG_MPI2OMP2 Emax/=0: 150
ISBA_3L_SN3L_ALL4_SGH_OMP4 Emax/=0: 140
ISBA_3L_SN3L_ALL4_SGH_MPI2OMP2 Emax/=0: 140
ISBA_3L_SN3L_ALL4_SGH_CANOPY_DRAG_MPI2OMP2 Emax/=0: 178

Ces erreurs semblent être dues à la parallélisation OPEN-MP pour l'option ALL4_SGH d'ISBA.
Ces erreurs sont déjà présentes dans la version que j'ai récupérée de Béatrice Vincendon.
==> Il faudra itérer à nouveau avec elle pour voir d'où cela peux venir.

Les deux suivants, sont complètement isolés dans leur série dans la base de test. Je n'ai pas de piste pour comprendre le problème.
TEB_BEM_TWO_WALLS_DXCOIL_FINCAP_AUTOSIZE_OMP4 Emax/=0: 2
2 nouveaux diagnostiques qui sont différents (durées cumulées dans deux seuils de stress thermique donnés différentes d'un pas de temps entre le run en OPEN-MP 4 et les autres runs). Le calcul de ces variables se fait dans la routine utcic_stress.F90, où il y a un where.
Ce n'est pas reproductible. J'ai refait une simulation de ce cas sans problème.

SNAP_NC_NETCDF Emax/=0: 2
Les champs E_CO et E_O3 sont différents. Alors que ce sont exactement les mêmes exécutables qui ont tourné.

****************************
* PHASE 4 de la VALIDATION *
****************************

Tests AVEC implicitation du couplage neige/toit et AVEC correction de sunpos.F90 :

Diagnostics are only sligthly affected (difference of about 1E-12).

TEB - V8_lafaysse

The branch V8_lafaysse is compared with the branch V8_VM_BUGS (r2037).

Results are exactly identical for all tests which do not use Crocus.

Differences in tests with Crocus mainly come from the correction of solar radiation absorption (point 1). They can also come from the new wind formulation (point 4) and from some rewriting of the code (threshold effects sensitive to numerical approximations).
( to be completed)

*************
* REMARQUES *
*************

==> changement dans la base de tests du fichier de donnees architecturales de bati (fichier csv) par une version plus recente pour tenir compte des nouveaux champs dans TEB (panneaux solaires). A l'avenir, la lecture de nouveaux champs pourrait être mise sous clé (mais pour l'instant ce n'est pas une option utilisée ailleurs qu'à TURBAU, et on a mis à niveau ce fichier dans nos dernières simulations en 7_3).

==> modifier le fichier script_teb_run.sh dans la base de tests: changement de "" en "ROW30" pour CCH_BEM dans la namelist.

==> Attention, j'ai changé les numeros de VERSION et Bugfix à 8 et 0 respectivement. Faites attention que vos divers tests de lecture de variable dans un fichier soient bien écrits (et pas uniquement en fonction du numero de bugfix par exemple).

==> ATTENTION, amélioration à prévoir: les champs T2MMIN/MAX HU2MMIN/MAX sont remis à XUNDEF ou -XUNDEF dès qu'ils sont écrits. Donc dans uns simulation offline, ils sont corrects dans les fichiers de diagnostiques, mais ils sont faux dans le fichier de restart à la fin !! (ça n'a pas changé entre les 2 versions, sauf qu'il y a les champs pour TEB et Flake en plus).

==> NAM_LANDUSE_PREP_NC: OLD plante systématiquement, et ce quelque soit la version du code en _old (même si c'est le nouveau, comme lorsque j'ai fait le test de parallélisation). C'est donc probablement un problème lié à la base de test.

Note concerning the test database (version of January 2014) :

The following tests are suspicious for all versions (systematic crash)

NAM_LANDUSE_PREP_NC: OLD

=> because of script_to_old.sh, that needs to be adapted to use NC PGD and PREP file for OLD version of the code (currently, it's LFI, due to the first version of STRATO that compared NEW_PREP and trunk_for_test)

SURF_FILETYPE_FA_RESTART_OLD: OLD

=> works with NEW_PREP. Intermediate versions have to be tested.

SEABATHY_ETOPO2: NEW

=> works with NEW_PREP. Intermediate versions have to be tested.

SEABATHY_ETOPO2: OLD

=> works with NEW_PREP. Intermediate versions have to be tested.

And that the non reproductibility diagnostic appears systematically for all versions :

SNAP_NC_NETCDF Emax/=0: 2

V8_lafaysse - V8_for_assimilation revision 2246

Tests were realized on PC and BEAUFIX in O2 and DEBUG modes and give exactly the same results, except for the assimilation tests that where corrected in V8_FOR_UNIFIED_ASSIMILATION.

V8_for_assimilation revision 2246 - V8_GMGEC_MEB

Tests are done with options LREPROD_OPER so that new defaults aren't taken into accound during the tests.

Given the amount of modifications in V8_GMGEC_MEB, many simulations give different results but nothing so big or unexpectable we can think there is a problem.

V8_GMGEC_MEB - V8_MESONH

The only differences are numerical and linked to the little changes of Juan.

V8_MESONH - V8_ASSIM

Ok, no changes.

V8_MESONH - V8_RJTECH

V8 RJTECH report

(R. Jasinskas)

Main changes in SURFEX v8

2014-02-06T10:32:37Z

Branch NEW_PREP (MOSAYC)

NB: some default values for namelist options changed in V8.

First, see:

http://www.cnrm.meteo.fr/surfex-lab/spip.php?article357&var_mode=calcul

Then, default for LNOSOF and LVERTSHIFT from NAM_SURF_ATM are inverted:

http://www.cnrm.meteo.fr/surfex-lab/spip.php?article77&var_mode=calcul

Also, values for GMES_ST (used with CPHOTO = NIT,AST,LST,NCB) in ini_data_param.F90 changed :

GMES_ST	V7_3	V8
NO	0.003	0.003
ROCK	0.003	0.003
SNOW	0.003	0.003
TEBD	0.003	0.0018
TRBD	0.003	0.0012
TEBE	0.003	0.0019
BOBD	0.003	0.0018
SHRB	0.003	0.0016
BONE	0.002	0.0019
TENE	0.002	0.0019
BOND	0.002	0.0012
TRBE	0.002	0.0012
C3	0.001	0.001 / 0.00175
C4	0.001 / 0.009	0.006 / 0.0098
IRR	0.001 / 0.009	0.006 / 0.0098
GRAS	0.001	0.0024
BOGR	0.001	0.0024
TROG	0.006	0.017
PARK	0.001	0.0024

For crops types, "/" indicates LAGRI_TO_GRASS=T/LAGRI_TO_GRASS=F.

Contains:

Short-term optimizations of PREP step proposed by ALADIN in March 2013 (Tayfun Dalkilic, Daan Degrauwe):
- PATCH AVERAGING:
  - the loop on covers is reduced to a loop on effective present covers on the area.
  - open-mp directives are added to still improve this loop.
  - generalization by reducing the size of array XCOVER to the present covers in the area, in PREP but also in PGD step (leads to use less memory)
- BILIN INTERPOLATION:
  - hor_interpol is no longer systematically called for all 12 vegtypes.
  - the fortran function count that takes much resource is replaced by an optimized code.
  - open-mp directives are added to improve the more expensive loops.
    
    Short-term optimizations of PREP
    
    Tayfun Dalkilic, Daan Degrauwe
    March 11, 2013
optimizations of PGD in GAUSSIAN grid:
- changes in the writing of the loops
- adds of open-mp directives to improve the more expensive loops during the index calculation.
PGD and PREP files can be written and read in NETCDF format:
- new specific I/O routines
Library XRD38 is replaced by XRD39:
- lead to include USE MODI of FA routines in SURFEX code, what is not totally satisfying (incompatible with the GMKPACK compilation) and should be revised.
A PREP file can enter the namelist only for the snow part:
- the entry for the PGD corresponding file was missing
Corrections in the chemical part:
- in collaboration with CAIAC (M. Michou, P. Hamer, D. Saint Martin)
NEW_PREP was used to test and validate the new base of automatic tests. This lead to 3 categories of changes in the code:
- bugs were detected thanks to the base of tests and corrected in NEW_PREP, notably for the parallelization with OPEN-MP,
- for an optimal use of the base of tests, names of namelists entries were changed so that they are specific to each block (a same namelist entry can't appear in two different blocks),
- new diagnostics were added in NETCDF format, especially for the TEB part.
CROCUS files are re-written following SURFEX rules for names of variables and typography, and trying to make common subroutines when the same code is used in different places:
- this lead to include a part of the CEN contribution

Remarks:

Some adds for PGD & PREP files in NETCDF format are not set yet, for example the possibility to use NETCDF format in case of a LIMP_ZS, LIMP_COVER or other LIMP_... in running PGD.

Branch NEW_PREP_HIRLAM (MET.no, SMHI, R. Hamdi) built on top of NEW_PREP

Reminder:

SODA aims at gathering and externalizing developments around data assimilation in SURFEX.
For the moment, VARASSIM from Alina Barbu is not merged with SODA. So, a question remains if VARASSIM has to be put in SURFEX V8, and if yes, if in or out of SODA. It has been agreed to merge VARASSIM developments by Alina and SODA by HIRLAM, but it has not been agreed at which time this should happen. SURFEX V8 might come too soon for this.

Branch contains:

Separate IO in SODA from the assimilation part. This allows us to use the main IO infrastructure of surfex and to be able to call SODA from CANARI. Assimilation related code moved to src/ASSIM.
EKF initialization updates
Extrapolation of variables. Allow optional extrapolation in PREP step through namelist variables LEXTRAP_TG,LEXTRAP_WG,LEXTRAP_WGI,LEXTRAP_SN. Disabled by default.
Set an upper limit for extrapolation distance
Allow search over a subdomain by providing NDIM2 in the call to OI_HOR_EXTRAPOL_SURF.F90
Fix to prevent use of uninitialized values in solar zenith calculations when running soda.
Optional weighting of grid average T2M giving more weight to the land tiles. Controlled by LT2MMW and switched off by default
Add exchange coefficient PCD to Get_FLUX_n interface ( Rafiq Hamdi (ALADIN), but needed when running ALARO for coupled mode)
Control of upper limit on initial snow water equivalent through LSWEMAX and XSWEMAX
Optional removal of small snow values through LAESNM
Fix for vegetation_update.F90 to be OpenMP safe for ECOCLIMAP-II and use the correct 10 day periods for months with 31 days.

Branch V8_GMAP, built on top of NEW_PREP_HIRLAM

This branch correspond to the code that was entered in Cy40t1 and is also known as surfexv7.3+

1) Optimisations PREP (openMP) and bfs for oi_control ( from Mohanmed Jidane and Francoise Taillefer):

surfex/ASSIM/oi_control.F90 (correction for 10m winds when it is above forcing level)
surfex/SURFEX/hor_extrapol_surf.F90 ( OpenMP for PREP)
surfex/SURFEX/hor_interpol_gauss.F90
surfex/SURFEX/prep_grid_gauss.F90 (debugging PREP, lecture lfi global gauss grid)

2) For stochastic flux perturbations (F. Bouttier) : (under logical switches LPERTFLUX for sea in NAM_SEAFLUXn, LPERTSURF for isba in NAM_ISBAn)

Purpose: implement stochastic surface perturbations for ensemble prediction. This code will be used in the MF AROME ensemble prediction system.
Description:
- initial perturbation of surface fields are prepared before the forecast by an external program that manipulates the Surfex initial conditions file to
  - (1) alter some prognostic fields, and
  - (2) append some perturbation patterns for physiographic fields.
- The code modifications perform two tasks:
  - to read and apply perturbations to fields that cannot be easily perturbed directly in the initial file.
  - to preserve initial perturbations of quasi-constant fields (ISBA veg, lai, cv,albedo, Z0), even when there is a decadal update of physiographic fields from PGD.
- Details:
  - Both are achieved by storing the field perturbations in new modd_isban module fields XPERTVEG, XPERTLAI, XPERTCV, XPERTALB, XPERTZ0.
  - The perturbations are read from the initial file in compute_isba_parameters.
  - The same albedo perturbation pattern is applied to all ISBA albedoes.
  - The same Z0 perturbation pattern is applied to all ISBA roughness parameters, except vegetation Z0.
  - The physiographies for non-ISBA surfaces (ice, snow, town, flood, etc) are not perturbed.
- The surface perturbation code is activated by two namelist switches:
  - NAM_ISBAn/LPERTSURF: perturbations of ISBA (veg,lai,cv,albedo,Z0). Requires fields PERTALB and PERTZ0LAND from initial file.
  - NAM_SEAFLUXn/LPERTFLUX: multiplicative perturbation of Ecume fluxes. Requires field PERTSEAFLUX from initial file.
- When these switches are off, the code is bit-reproducible, and standard initial files can be used.
compute_isba_parameters.F90
hor_extrapol_surf.F90
modn_seafluxn.F90
coupling_isban.F90
hor_interpol_gauss.F90
oi_control.F90
coupling_seafluxn.F90
modd_isban.F90
prep_grid_gauss.F90
ecume_flux.F90
modd_seafluxn.F90
read_seafluxn.F90
ecume_seaflux.F90
modn_isban.F90
vegetation_update.F90

3) Pb compilation on IBM (from Ryad El Khatib and Oldrich Spaniel) : most of the modifications are linked with the renaming of "NINDX2" from MODD_SURFEX_OMP as "NINDX2SFX"

ch_init_snapn.F90
gather_and_write_mpi.F90
gather_and_write_mpi_k4.F90
get_size_fulln.F90
init_io_surf_fan.F90
init_surf_atmn.F90
modd_surfex_omp.F90
pgd_grid.F90
pgd_surf_atm.F90
pgd_teb_veg.F90
prep_isba_netcdf.F90
prep_ocean_ascllv.F90
read_and_send_mpi.F90
read_csvdata_teb.F90
read_isban.F90

Branch V8_BV, built on top of V8_GMAP

ISBA-TOP improvements :

isba-top V8

Branch V8_VM, built on top of V8_BV

TEB improvements

watering:
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.
fraction of residential use:
XUNIF_RESIDENTIAL=0.5, read in namelist NAM_DATA_TEB
used only for solar panels and occupying schedule
solar panels on roofs:
1. in NAM_TEB: LSOLAR_PANEL=T
2. in NAM_DATA_TEB:
  1. XUNIF_FRAC_PANEL= 0.50 ! fraction of solar panels on roofs
  2. XUNIF_EMIS_PANEL=0.90 ! emissivity of solar panels
  3. XUNIF_ALB_PANEL=0.10 ! albedo of solar panels
  4. XUNIF_EFF_PANEL=0.14 ! efficiency of photovoltaic solar panels
3. new diagnostics:
  1. PHOTC_BLD ! photovoltaics production cumuled on the length of simulation ( J/m2 of building at floor)
  2. 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.
4. 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 photovoltaics.
  If the buildings are of other use, there is only photovoltaics.
  If the use is mixed (for example 0.5), the sharing out is proportional.
schedule of occupying:
1. in NAM_TEBn:
  1. XDT_RES = 3. ! difference of ordered temperature (heating and cooling) when the building (residential) is not occupied : day (K)
  2. XDT_OFF = 3. ! difference of ordered temperature (heating and cooling) when the building (offices) is not occupied : nights and WE (K)
2. In grid points where buildings are mainly residential (>0.5), it's admitted that during day, buildings are unoccupied, the ordered heating decreases of XDT_RES degrees Celsius and the ordered cooling grows of XDT_RES degrees Celsius.
  In grid points where builings are mainly not residential (<0.5),> during nights and WE, buildings are unoccupied, the ordered heating decreases of XDT_OFF degrees Celsius and the ordered cooling grows of XDT_OFF degrees Celsius.
  Simulations done until now supposed that the ordered temperature would not be modified even if the buildings were unoccupied (XDT_RES=0. and XDT_OFF=0.). The study can go on with XDT_RES=0. and XDT_OFF=0. anyhow, the sensibility study will so give the influence of the ordered temperature.
Cumulative UTCI stresses:
Values are in seconds.
1. UTCIC_IN_ECS corresponds to the time in seconds spent below the threshold of Extrem Cold Stress inside.
2. UTCIC_IN_VCS corresponds to the time in seconds spent between the thresholds of Very Strong Cold Stress and Extrem Cold Stress inside.
3. UTCIC_IN_NHS corresponds to the time in seconds spent in comfortable index inside.
4. UTCIC_IN_MHS corresponds to the time in seconds spent between the thresholds of Medium Heat Stress and Strong Heat Stress inside.
5. These times are on the total length of simulation (it can make a lot of seconds...)
6. Glossary:
  1. ECS : extreme cold stress
  2. VCS : very strong cold stress
  3. SCS : strong cold stress
  4. MCS : medium cold stress
  5. LCS : light cold stress
  6. NHS : no heat stress (= comfortable)
  7. MHS : medium heat stress (to be noted that there is no Light Heat stress)
  8. SHS : strong heat stress
  9. VHS : very strong heat stress
  10. EHS : extreme heat stress
  11. _IN_ : indoors
  12. _SH_ : outdoors in shadow
  13. _SU_ : outdoors in sun
Mechanical and controlled ventilation:
new option for the type of ventilation in buildings (namelist entry: NATVENT ).
Change of default values for roughtness of TEB walls:
&NAM_TEBn
CZ0H = "KAND07"
CCH_BEM = "DOE-2"
add/correction of diagnostics for min and max temperatures:
add in Flake (didn't exist)
corrections in TEB
Various bugs and diagnostics, for TEB or GARDEN.

Branch V8_lafaysse, built on top of V8_VM and V8_VM_BUG

Contribution for V8 from CEN (M. Lafaysse)

Branch V8_GMGEC_MEB, built on top of V8_FOR_UNIFIED_ASSIMILATION revision 2246

Contribution for V8 from GMGEC & GMME (MEB)

Bertrand Decharme, Aaron Boone

Description for the coupling of SURFEX with OASIS

Bertrand Decharme

Introduction of GELATO in SURFEX

Stéphane Sénési, Bertrand Decharme

Branch V8_MESONH, built on top of V8_GMGEC_MEB revision 2745

Changes from Mesonh for SURFEX V8

Gaëlle Tanguy, Juan Escobar

Branch V8_ASSIM, built on top of V8_MESONH revision 2779

a correction for TREEDRAG from Mesonh (XH_TREE must be initialized if LTREEDRAG=T and CPHOTO='NON'
DR_HOOK is now a library compiled with the rest of the SURFEX code (simpler). For the moment, the DR_HOOK sources are in directory DRHOOK_SIMPLE.
Corrections for TOPMODEL, following running STRATO in Debug mode
To print an error message if the namelist is not present in the run directory (E. Martin)
An add from P. Tulet for the diagnostic for the wind at 10 meters
An optimization for extrapolation during PREP when the interpolation is bilinear, and add of the possibility to write the ALBEDO diagnostics during PREP (L. Raynaud)
Merge of OI_MAIN & SODA in SODA-OI (same results are obtained)
Merge of VARASSIM & SODA in SODA-EKF (same results are obtained)
Add of other CSURF_FILETYPE format than ASCII in SODA
It's now possible to set CTIMESERIES_FILETYPE='NC'

Branch V8_RJTECH, built on top of V8_ASSIM revision 2851

See first: technical modifications from Hirlam about modifications inside the code.

Then, concerning compilation:

1). all scripts where made born shell compatible and point to /bin/sh for portability.
As a general recommendation all other shells should have '#!/usr/bin/env SHELL'.

2). In bin/spll tool explicit .mod and .o moving by makefile rules rules disabled, because often breaks on lustre attached network systems, specially at partial recompile.

Also to avoid using VPATH outside SURFEX compilation directory, that can create lots of unexpected problems. Now handled directly by compiler options with $(REALFC) and $(MODOUT_FLAG) variables.

If compiler type is not known, the ".NOTPARRALLEL:" makefile rule is applied and compilation will be done as with make -j1, by using implicit mv commands for .mod files.

3). configure script is simplified to provide sane defaults that can be overriden by:

 export ARCH=${ARCH} #arch name export VER_MPI=${VER_MPI} #Version of MPI {VER_MPI=NOMPI,MPIVIDE,LAMMPI,MPIAUTO,...} export VER_OMP=${VER_OMP} #Version of OMP {VER_OMP=NOOMP,OMP}, actually handled by user Rules file only export OPTLEVEL=${OPTLEVEL} #OPTIMIZATION LEVEL {OPTLEVEL=DEBUG,O2} export VER_USER=${VER_USER} #USER VERSION export CONF_DOLLAR=${CONF_DOLLAR} # 'Scandollar' Post Processing default variables export VERSION_CDF=${VERSION_CDF} #used only for CDFAUTO case, which redistributable archive to extract export VER_CDF=${VER_CDF} #Version of NETCDF {VER_CDF=CDFAUTO,CDFEXT,...} export VERSION_GRIBAPI=${VERSION_GRIBAPI} #used only for GRIBAPI_AUTO case, which redistributable archive to extractexport VERSION_OASIS=${VERSION_OASIS} export VER_GRIBAPI=${VER_GRIBAPI} #Version of GRIBAPI {VER_GRIBAPI=GRIBAPI_AUTO,GRIBAPI_EXT} export VERSION_OASIS=${VERSION_OASIS} #which redistributable archive to extract export VER_OASIS=${VER_OASIS} #Version of OASIS {VER_OASIS=NOmct,mct} export VER_DRHOOK=${VER_DRHOOK} #Version of DR_HOOK {VER_DRHOOK=SIMPLE,BYPASS,MINI}

4). AUTO gribapi and netcdf redistributable compilation is forced to be done before compiling MASTER.

This simplifies makefiles greatly and ensures that redistributables can be compiled for $(REALFC) compiler.

MPI versions were disabled, since rest of makefile assumes that $(F90) is $(REALFC). Also disabled user F90FLAGS propagation into these redistributables, to avoid using non standard/broken libraries.

Special note about -i8/-fdefault-integer-8 in these libraries, this can lead to bad output files.

5). In bin/spll added EXTERNAL variable that does module dependencies blacklisting currently it includes netcdf and gribapi, since external libraries should be assumed static.

6). In makefiles CPPFLAGS were renamed to FPPFLAGS to separate cpp defines for fortran and C parts.

Mainly done to avoid name collisions in new post XRD39 versions of FA/LFI c part variable names.

7). OASIS is off by default(NOmct), since it requires MPI and for running TESTS it is not needed. Oasis compilation flags are no longer static (make.inc_cnrm) but dynamic from user Rules file through oasis_sfx_make.inc copy after extracting redistributable. This allows to have agnostic setup.

8). 'make clean' was fixed not to loop over and over

9). 'make user' now can be done without having done 'make master' previously. 'make user' program linking list separated from $(PROG_LIST) into $(PROG_LIST_USER).

10). VER_OMP handling should be done by user Rules files, currently I just check if NOOMP is used and null out FCFLAGS_OMP CFLAGS_OMP flags.

11). added splr.pl script that replaces all spll/spn/awk/cpp functionality.

It uses my updated Fortran90_stuff.pm (with some bugfixes after testing on surfex).

Creates much cleaner/correct modi interfaces. Still -fimplicit-none is recommended.

It's capable to :

create .D dependencies (—deps flag)
create spll_modi_* intefaces (—modi flag), optionaly to make generic interfaces
split into separate units based on (—split flag)
perform simple fpp (—fpp) by commenting out cpp directives with '!zFPP ' prefixes
wrap splitted units into modules (—wrap flag)

Created inteface has name pattern spll_modi_NAME.f90.

Created splitted source has name pattern spll_NAME.F90

As a special case splr.pl contains some exceptions for SURFEX:

read_surf/write_surf is not splitted until bug in SODA_OI2 (LEXTRAPOL_ISBA case) is fixed
if file name contains (lonlat_rot|gltools_newice_r|glt_constrain_r) genric interface generation is disabled, since these files are not TKR (Type,Kind,Rank) compatible, passing za(*) to xa(*,*).
$f90s_ALLOW_DIRTY_MOD_REGEX allows to select which non "USE ...,ONLY:" modules get included in interface for dependencies in dummy arguments declarations.
If file name contains 'goto_' then —wrap is disabled, since those subroutines are RECURSIVE

12). .F90 rules were added to makefile implicit compilation rules.

13). Special user selectable flags were added:

$(USE_SPLR)=="YES" selects to use splr.pl for DEP step. Otherwise spll chain.
$(AVOID_CPP)=="YES" used to bypass cpp usage on input source files. Sometimes it is useful to work with partially preprocessed files.
$(USE_SPLR_WRAP)=="YES" to not generate separate spll_modi_ file, but to wrap whole unit into spll_modi_UNITNAME.F90. This implements part No.7 of Hirlam technical modifications document. Also it is great functionality for checking missing USE MODI in code tree.
$(DO_ASM)=="YES" for assembly listing generation (.s), if not useful can be removed. I used it to ensure freeform changes were lexically identical to f77 file logic.

14) Several DRHOOK library variants were added:

LIB/DRHOOK/BYPASS(dummy library), needed for my tools.
LIB/DRHOOK/MINI (minimal functionality required for signal handling, tracebacks, drhook traces)

Still can do memory usage profiling, catch crashes, etc, just without MPL layer that allows to compile freely with/without MPI and/or OpenMP LIB/DRHOOK/MPL (together with MINI forms DRHOOK_SIMPLE), just for testing compatibility -DSFX_MPL.

15). few fixes in XRD39 FA/LFI so that use of FA modi interfaces no longer needed, -DOFF define removed.

16). custom XRD40 FA/LFI library with selective changes between XRD39 and XRD40, but without whole new c layer.

17). I/O specific CPP flags in surfex now prefixed with 'SFX_' to avoid various problems.

18). -DNOMPI inverted to -DSFX_MPI to better explain it's usage in GELATO and to avoid double NOT while reading code (#ifndef NOMPI, #if !defined NOMPI)

Also -DSFX_MPI handled automatically by $(VER_MPI) value if it is not "NOMPI"

19). temporary cpp define '-DRJ_OFIX' was added to fix NOOMP version of OFFLINE execution. After review it should be removed in final version. Same for -DRJ_PFIX.

OI_MAIN + PREP from GMAP

2013-12-03T12:00:16Z

Contains:

modifications to run PREP with input LFI AROME file (M. Jidane)
correction in oi_control (surface analysis, F. Taillefer)

This contribution could be added to a general GMAP contribution containing the SURFEX version that is included in AROME CY40t1 (containing also SODA from HIRLAM and Ensemble prediction in AROME from MICADO).

The question of who will test the whole package or each part of it is still unresolved.

rotated grid from P. Samuelsson

2013-12-02T16:22:32Z

Contains:

a new grid option: rotated lat lon

P. Samuelsson will include this contribution in SURFEX V8.

It's considered with the intermediate developments.

VARASSIM from VEGEO

2013-11-25T11:33:22Z

Contains:
VARASSIM will be in the end part of SODA, but, as it was not done for the moment (lack of tie), and because of the need to have a running version of VARASSIM maintained in SURFEX, the external version is introduced in version 8.
A specific test case will be provided to enter STRATO (base of automatic tests).
This contribution is among those concerning a very specific part of the code.

- Surfex v8_0

Draft of the schedule for the building of SURFEX V8

2013-11-22T17:26:18Z

NEW_PREP is the starting point of SURFEX V8.

January 2014: Bertrand merges and tests CNRM_CM6_DEV (+ ECUME).
February 2014: Bertrand merges MEB_DEV2 and begins to test it
March 2014: Bertrand, Aaron, Eric Brun, StÃ©phane SÃ©nÃ©si, Christine Delire, Roland SÃ©fÃ©rian, Patrick Le Moigne test NEW_PREP+CNRM_CM6_DEV (+ ECUME) +MEB_DEV2

1-15 April 2014: Matthieu Lafaysse merges and tests lafaysse_dev_for_test
1-15 May 2014: Valery Masson merges and tests TEB_DEV_1514
16-31 May 2014: Patrick Samuelsson merges and tests rotated latlon grid

// NB: the 3 following contributions were already merged with NEW_PREP and included in AROME CY40t1. To be seen if the resulting common version couldn't be taken as basis for SURFEX V8. The question of tests remains to be solved. //

1-10 June 2014: Ulf Andrae merges and tests new_soda
11-20 June 2014: FranÃ§ois Bouttier merges and tests stochastic surface perturbations
21-30 June 2014: ?? merges ans tests OI_MAIN + PREP from GMAP

1-10 July 2014: Alina Barbu merges ans tests VARASSIM

11-20 July 2014: Rafiq Hamdi merges and tests ALADIN modifications // could be included in another contribution because concerns only 3 files - to be seen with Rafiq //
contribution included in branch NEW_PREP December,6^th 2013.

21-31 July 2014: GaÃ«lle Tanguy merges and tests MESONH modifications
1-15 August 2014: Rimvydas Jasinskas merges and tests HIRLAM technical modifications
16 August - 15 September 2014: Gabriel Jonville merges and tests the suppression of global variables

branch MEB_DEV2 from MOANA & SMHI

2013-11-22T16:37:22Z

Contains:

An energy budget is explicitly taken into account for the vegetation.

This development is finished and a lot of tests were already done by A. Boone and P. Samuelsson, but for more validation, the tests will continue until February 2014.

MEB will be merged in SURFEX V8 by B. Decharme after its CNRM_CM6_DEV contribution. This work could be finished around the end of March or during April 2014.

March will be dedicated to the ending debugging, and other people will help with it:

Aaron for MEB
Eric Brun for MEB + Snow
Christine & Roland for Landuse and dynamic vegetation
Patrick Le Moigne for FLAKE

New test cases should also be defined for the base of automatic tests (STRATO).