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

Didier R. 08/06/2021: add Leonard terms (Moeng parametrization) instead of Reynolds terms applied to vertical fluxes of r_np and Thl
for implicit version of turbulence scheme