Properly extracts the staggered geopotential height, and removes the 1000 scaling in the Fortran routine. Fixes #75

fortran/calc_uh.f90

@@ -154,8 +154,6 @@ SUBROUTINE DCALCUH(nx, ny, nz, nzp1, zp, mapfct, dx, dy, uhmnhgt, uhmxhgt, us, &
 
     !$OMP END PARALLEL
 
-    uh = uh*1000.   ! Scale according to Kain et al. (2008)
-
     RETURN
 
 END SUBROUTINE DCALCUH

src/wrf/g_helicity.py

@@ -197,7 +197,7 @@ def get_uh(wrfin, timeidx=0, method="cat", squeeze=True,
                           meta=False, _key=_key)
     v = destagger(v_vars[varname], -2) 
     
-    zp = ph + phb
+    zp = (ph + phb) / Constants.G
     
     uh = _udhel(zp, mapfct, u, v, wstag, dx, dy, bottom, top)