snapshot...

modules/util/util.py

@@ -231,6 +231,7 @@ def pressure_to_altitude(pres, temp, prof_pres, prof_temp, sfc_pres=None, sfc_te
 
     if sfc_pres is not None:
         if pres > sfc_pres:  # incoming pressure below surface
+            print(pres, sfc_pres)
             return -999.0
 
         prof_pres = np.array(pres_s)
@@ -261,7 +262,7 @@ def pressure_to_altitude(pres, temp, prof_pres, prof_temp, sfc_pres=None, sfc_te
 
     z = thickness_hydrostatic(prof_pres, prof_temp) + sfc_elev
 
-    return z
+    return z.magnitude
 
 
 # http://fourier.eng.hmc.edu/e176/lectures/NM/node25.html