Vertical gradients and shuman

[SebastienRietteMTO]

In Meso-NH, arrays are overdimensioned along the vertical axis.
In AROME, we can use "normal" arrays for all parametrisation except the turbulence one.

This issue list the steps to realise in order to suppress the unphysical points for AROME:

1. Update vertical halo before calling the parametrizations in Meso-NH: At least when calling the PMMC09 scheme, the values in the unphysical level under ground is not equal to the values on the first physical level. Could it be possible to update the vertical halo before calling each parametrisation?
2. Try to suppress computation on unphysical points in shuman_mf, keeping the results unchanged
3. If the preceding rewritting is successful, EDKF can be rewritten to suppress computation on unphysical points
4. Again, if the shuman_mf is successful, the turb scheme could use these routines (all the physics could use the same set of routines to compute gradients and shuman)
5. The last step would be suppress the computations on unphysical points in the turbulence scheme, giving the possibility to suppress the overdimensioning in apl_arome