23/07/2018  Correction from Arsene Druel for maj average
23/07/2018 Correction to read the 10-days periods in the PGD file (Rafife Nheili)
23/07/2018 Correction for PGD and NAM_DATA (Arsene Druel)
22/08/2018 Bug for parallel PREP & snow (array sizes) (Q. Rodier)
23/08/2018 Correction for LSPLIT_PATCH and LRESET_BUDGET=F (P. Samuelsson)
23/08/2018 Correction to have the same diagnostic output fields names with LSPLIT_PATCH=F
05/09/2018  bug fix from Aaron for MEB - initialisation
14/09/2018 Bugfix ISBA MEB (Aaron Boone)
25/09/2018 Error of dimension in allocation
10/10/2018 Bugfix ISBA MEB (Aaron Boone)
23/10/2018 Correction for upward compatibility for snow albedo
14/05/2019 BugFix: correction of H_TREE calculation by patch (impact only if H_TREE given in namelist, and if NPATCH different of 1 or NVEGTYPE)
21/06/2019 BugFix: correction of ZMASK used to initialise TEB variables when GTEB=FALSE in the extern file read for initialisation
23/08/2019 Bugfix NVEGTYPE_ECOSG / NVEGTYPE in read_nam_pgd_isba_meb (CEN)
17/11/2020 Bugfix: wrong initialisation of LIGNIN_STRUC and SOILCARB (when CRESPSL=CNT)
17/11/2020 Bugfix: correction for the writing of the variable ROOTFRAC in PREP file
29/03/2021 BugFix: YPHFILETYPE initialized twice, and missing initialization for YFERTFILETYPE
Updated to last commit (e826177a)

Quentin 29/06/2021: bugfix SURFEX avoid abort for interpolation of a field which is fully XUNDEF at spawning with ECOCLIMAP-SG. With ECOSG, a field can be fully XUNDEF after PGD.

Philippe 28/06/2021: small optimisation in the call to ice4_rainfr_vert (create arrays with zeros instead of multiplications by zero)

Philippe 28/06/2021: netCDF: use correct integer datatypes for netCDF integers + respect Doctor norm

Philippe 28/06/2021: IO: better management of pre-5.5 files (default values for missing fields instead of crash)

Philippe 28/06/2021: allocate to zero-size arrays in ini_model_n to prevent use of unallocated arrays later
+ allocate some of them to non zero-size only if necessary

M.Goret 28/06/2021:bugfix increase LNOWRITE table and add abort if writing surfex field needs more memory (for long simulations/large number of surf. field writes without restart)

- LIMA sedimentation (activated as for ICE3/C2R2)
- correction of surface flux reduction with respect to the resolution (LGZ) to activate only if needed (grey-zone resolutions)

Quentin 25/06/2021: bugfix introduced by commit 09f31 (OCEAN LES). Missing a CALL for the atmospheric case for computation of the pressure

- Initialization and coupling aerosols from CAMS forecasts/analysis
- Kappa-Kohler theory to compute CCN activation
- Subgrid condensation with LIMA (LSUBG_COND=T)
- Prognostic supersaturation with LIMA (LSPRO=T and LADJ=F)
- Option to relax the adjustement to saturation

TODO : 014_LIMA test case is not working for now with new options LSPRO=T and LADJ=F

Quentin 16/06/2021: bugfix clarification of some points for the merge of Ocean-Atmosphere coupled LES dev

P.Marquet 15/06/2021: correction on the exponent of final BL89 mixing length according to Lemarié et al. 2021