Quentin 29/07/2022: Packing mode_turb_ver_thermo_flux + move OCEAN/OCOUPLES 3D variables use (XSSTFL) as intent(in),
TODO: LHRAD desactivate : to be transcripted for GPU use later (too much hor. gradient)

Quentin 03/05/2022: turb GPU reduce number of caracters by line (max 130) after expand in physical points

This reverts commit 566b9f7e.

Execution time is now equivalent with or without array-syntax / do loops conversion

Quentin 28/03/2022: to be fixed : NIMAX_ll not available for AROME ==> remove part of code that compute ZDIST for LOCEAN outside of PHYEXé

TO BE TESTED IN AROME

This reverts commit fc7d4b25.

Merge turbulence bit-repro sur vxdev64

Cas test AROME non bit-repro

This reverts commit 6050b2b6, reversing
changes made to 71f963d8.

- One part of the code is not activated in AROME (this part was changed in MNH : is the bug corrected or not) ? To be tested
- A part of the fluxes at IKE was changed during the LOCEAN integration in MNH : Jean-Luc Redelsperger was contacted to see if it was a bug correction or not : to be tested in AROME as well

Arome aux:
- modd_dimn.f90
- modd_blowsnow.F90
- modd_gridn.F90
- modd_metricsn.F90
- modd_ref.F90
- modd_turbn.F90
- modd_frc.F90
- modd_oceanh.F90

- mode_gather_ll.F90 with ABORT
- modi_second_mnh.F90

prandtl.F90 is added in existing mode_prandtl.F90

gradient* files are merged and moved in externals
shuman moved in externals
modification of arguments order for shuman and gradients