From 5a433efc5d56ac612e4734e2237c473ea100797c Mon Sep 17 00:00:00 2001
From: tomrink <rink@ssec.wisc.edu>
Date: Fri, 12 May 2023 10:55:31 -0500
Subject: [PATCH] snapshot...
---
modules/util/romio.py | 6 ++++++
1 file changed, 6 insertions(+)
diff --git a/modules/util/romio.py b/modules/util/romio.py
index 1133ca49..ca6b58ca 100644
--- a/modules/util/romio.py
+++ b/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)
--
GitLab