12) The output SNOW fields can be written 2D (number of points / number of snow layers), by M. Lafaysse (CEN).

4 new namelist keys are introduced:

• LSNOWDIMNC, in NAM_WRITE_DIAG_SURFn:
  In case CTIMESERIES_FILETYPE = "OFFLIN", this key allows to write the snow output fields in 2D (number of points / number of snow layers). It’s applied in the ISBA, TEB, GARDEN and GREENROOF schemes.
  The concerned fields are:
  • the main snow fields written by writesurf_gr_snow
  • SNOWLIQ, SNOWTEMP, SNOWDZ and new the output snow fields depending on the key LPROSNOW (see below)

• LRESETCUMUL, in NAM_WRITE_DIAG_SURFn:
  Possible to activate only in OFFLINE mode, for the ISBA scheme.
  If LRESETCUMUL = .TRUE., the not cumulated fields are not instantaneous but are the cumulated fields averaged over the output writting time step. Consequently, the cumulated fields are reset to 0. at the end of each writting time step.
  The concerned fields are:
  • LEG, LEGI, LEV, LES, LESL, SNDRIFT, LER, LETR, DRAIN, QSB, RUNOFF, HORT, DRIP, RRVEG, MELT (key LSURF_EVAP_BUDGET)
  • DWG, DWGI, DWR, DSWE, WATBUD (key LWATER_BUDGET)
  • RN, H, LE, LEI, GFLUX, EVAP, SUBL, SWD, SWU, LWD, LWU

• LPROSNOW, in NAM_DIAG_ISBAn:
  This key :
  • controls the calculation and writting of new diagnostics for the CROCUS snow scheme. The new subroutine snowcro_diag called in snow3L_isba does the calculation of these new diagnostics. These new diagnostics are:
    • SNOWDEND: dendricity
    • SNOWSPHER: sphericity
    • SNOWSIZE: grain size
    • SNOWSSA: snow specific surface area
    • SNOWTYPEMEPRA: snow grain type (MEPRA)
    • SNOWRAM: snow ram resistance (MEPRA)
    • SNOWSHEAR: snow shear resistance (MEPRA)
    • SNDPT_1/3/5/7DY: fresh snow depth in 1, 3, 5, 7 days
    • SNSWE_1/3/5/7DY: fresh snow water equivalent in 1, 3, 5, 7 days
    • SNRAM_SONDE: penetration depth of the ram resistance sensor (2 DaN)
    • SN_WETTHCKN: Thickness of wet snow at the top of the snowpack
    • SN_REFRZTHCKN: Thickness of refrozen snow at the top of the snowpack
  • infers that the snow mantel and other diagnostic fields are reprojected on the vertical (according to the slope): new routine reproj_diag_isban.F90.
    
  • with CTIMESERIES_FILETYPE="OFFLIN", merges ISBA_PROGNOSTIC.OUT.nc and ISBA_DIAGNOSTICS.OUT.nc in one only file, ISBA_PROGNOSTIC.OUT.nc.

• LVOLUMETRIC_SNOWLIQ, in NAM_DIAG_ISBAn:
  This key implies that the diagnostic field SNOWLIQ (snow liquid water) is written in kg/m3 and not in m:
  SNOWLIQ = SNOWLIQ * RHOLW / SNOWDZ
  Moreover, if LPROSNOW=.TRUE., SNOWLIQ is corrected by the opposite of the cosine of the typical slope in the grid point (even if LVOLUMETRIC_SNOWLIQ = .FALSE.):
  SNOWLIQ = SNOWLIQ * ( 1./ SLOPECOS)
  This diagnostic is calculated for the ISBA-ES and CROCUS snow scheme.