Quentin 01/12/2022: merge ice4_sedimentation_stat with Ryads opti for new SNOW sedim computation from Mesonh 5.5 to 5.6

For GPU, ice4_nucleation cannot be put in an include file
as included subroutine must not contain local variables.

mesonh/micro/rain_ice.f90 -> mesonh/micro/rain_ice_old.f90
mesonh/micro/rain_ice_red.f90 -> mesonh/micro/rain_ice.f90
bugfix for the new condensation interface
interface harmonisation betwwen mesonh/micro/rain_ice_red.f90 and arome/micro/rain_ice.F90

PCOMPUTE/ZCOMPUTE replaced by logical array
Last where statement decorated with mnh_expand

Hail and graupel (with ICE4 option) budgets was wrong in the initial commit and different,
but still wrong, in the previsous commit.
Bug correction in fast_rh protected by a REPRO48 macro

Hail and graupel (with ICE4 option) budgets was wrong in the initial commit and different,
but still wrong, in the previsous commit.
Bug correction in fast_rh protected by a REPRO48 macro

All modd_* replaceed by structures in ICE3

nucleation optimised
docs/TODO updated

RIETTE Sébastien authored 3 years ago and RODIER Quentin committed 3 years ago

OCND2 included in common code (from AROME version)
OCND2 specific code moved in common/micro or arome/dead_code directories
High/Low cloud included for condensation (from Méso-NH version)
Harmonie SPP code moved outside of physics
condensation and ice_adjust optimised by Ryad
compute_frac_ice as an elemental procedure

Adaptation of mode_ice4_compute_pdf
Revisit of nucleation

Modules transformed into structures (rain_ice_param, rain_ice_descr, param_ice)
Both sedimentations are done
rain_ice main routine is done
Nucleation routines are done (using an elemental function)

OCND2 included in common code (from AROME version)
OCND2 specific code moved in common/micro or arome/dead_code directories
High/Low cloud included for condensation (from Méso-NH version)
Harmonie SPP code moved outside of physics
condensation and ice_adjust optimised by Ryad
compute_frac_ice as an elemental procedure

The 3 routines are merged into mode_compute_frac_ice
using optimisation (with pointers) developped by Ryad.

Merge arome rain_ice with:
- Ryad's versio: optimisation, KPROMA blocs
- Meso-NH version: subgrid precipiation
Arome aux:
- mode_ll to provide get_indice_ll
- mode_budget

Known issues:
- Bug around KPROMA blocs: results are different when used
- precipitation fraction cannot be used as long as condensation is not merged
- mode_ice4_sedimentation_momentum is still not plugged

Modifications prepared by Ryad are stored in arome/modset_Ryad
To help the merging of ice4_fast_ri, the array syntax was replaced
by a DO loop