snapshot...

modules/util/romio.py

@@ -3,6 +3,7 @@ import numpy as np
 from util.geos_nav import GEOSNavigation
 from netCDF4 import Dataset
 import pickle
+from util.util import haversine_np
 
 
 def romio_to_fgf(filename, goes_e_w='EAST'):
@@ -24,6 +25,11 @@ def romio_to_fgf(filename, goes_e_w='EAST'):
         nav = GEOSNavigation(sub_lon=-137.0)
 
     cc, ll = nav.earth_to_lc_s(lon_array.flatten(), lat_array.flatten())
+    tlons, tlats = nav.lc_to_earth(cc, ll)
+    dist = haversine_np(lon_array.flatten(), lat_array.flatten(), tlons, tlats)
+    ok = np.invert(np.isnan(dist))
+    print(np.average(dist[ok]))
+    print(np.histogram(dist[ok]))
 
     cc = cc.reshape(lon_array.shape)
     ll = ll.reshape(lon_array.shape)