gvc: provide hvel after restart

src/3d/variables_3d.F90

@@ -194,7 +194,7 @@
 #include "dynamic_allocations_3d.h"
 #endif
 
-   hn = _ZERO_ ; hvel = _ZERO_ ; hun = _ZERO_ ; hvn = _ZERO_
+   ho = _ZERO_ ; hn = _ZERO_ ; hvel = _ZERO_ ; hun = _ZERO_ ; hvn = _ZERO_
    uu = _ZERO_ ; vv = _ZERO_ ; ww = _ZERO_
 #ifdef _MOMENTUM_TERMS_
    tdv_u = _ZERO_ ; adv_u = _ZERO_ ; vsd_u = _ZERO_ ; hsd_u = _ZERO_

src/getm/initialise.F90

@@ -56,7 +56,8 @@
    use turbulence, only: init_turbulence
    use mtridiagonal, only: init_tridiagonal
    use rivers, only: init_rivers
-   use variables_3d, only: avmback,avhback
+   use variables_3d, only: ho,hn,hvel,avmback,avhback
+   use vertical_coordinates, only: restart_with_ho,restart_with_hn
 #ifdef SPM
    use suspended_matter, only: init_spm
 #endif
@@ -271,6 +272,12 @@
 
 #ifndef NO_3D
       if (runtype .ge. 2) then
+         if ( restart_with_ho .and. restart_with_hn ) then
+            hvel = _HALF_ * ( ho + hn )
+         else
+!           throw a warning here, if ho or hn are needed before they
+!           are calculated in postinit_3d()
+         end if
          Uint=_ZERO_
          Vint=_ZERO_
       end if