Interpolation was not correct for U and V (did not use correct position indices)

Quentin 16/07/2021: stations : move LDIAG_RESULTS to LDIAG_SURFRAD and apply it (developpement was missing on station and write_station)

Quentin 16/07/2021: station change default time step to 60s and remove TYPE variable of csv reading.
TODO: remove the STATION%TYPE which is not used

F.Couvreux 30/06/2021: add LRELAX_UVMEAN_FRC forcing. Instead of LRELAX_UV_FRC which applies the u,v relaxation read in namelist on each atmospheric column,
LRELAX_UVMEAN_FRC applies the relaxation on the horizontal mean of u and v at each level. To be used in LES mode.

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

- 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

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

F. Auguste, T. Nagel, R. Schoetter, V. Masson et al. (IBM and Recycling turb.) + E. Jezequel (station in namelists)
TODO: check I/O restart and CALL of write_diachro in write_stationn

Quentin 27/04/2021: NAM_BLANK and DUMMY becomes NAM_BLANK_n for grid_nesting use. The values of the BLANK DUMMY variables can be different in n models. Cleaning of non-used CALL MODD_BLANK

Quentin 28/10/2020: default value of LRED is now TRUE (LRES = reduction of time dependency of ICE3/ICE4 microphysics scheme)

(cherry picked from commit 47004fb1)

Christine 27/11/2019: correction in the drag formula and application to building in addition to tree

RODIER Quentin authored 5 years ago and WAUTELET Philippe committed 5 years ago

(cherry picked from commit 254aa18e)

RODIER Quentin authored 5 years ago and WAUTELET Philippe committed 5 years ago

(cherry picked from commit cd621ccc)

Soline Bielli : 01/03/2019: Sea salt update: the significant sea wave height now influences salt emission.
The sea salt is either set to a default value of 2meters or read in ECMWF Grib file (if the variable exists).
Now 5 modes are available (3 before) and the order in output is changed.
The large-scale sea wave height from ECMWF grib is 2D. HALO2 update of the ls variable in modeln is not yet adapted to 2D field (it is switch off in this version for this variable).

Quentin 23/01/2019: modification of initialization of PM09 convection scheme at hectometric scale.The mass-flux at the surface is reduced and converge to EDKF value when resolution is meso-scale. Flag to turn on : LGZ = TRUE by default. Modifs by Rachel Honnert