snapshot...

modules/amv/intercompare.py

@@ -4,13 +4,15 @@ from util.setup import home_dir
 from util.lon_lat_grid import earth_to_indexs
 from util.util import haversine_np
 from deeplearning.amv_raob import get_profile_multi
-#import cartopy
-#from cartopy import *
+import cartopy
+from cartopy import *
 import metpy.calc as mpcalc
-#import cartopy.crs as ccrs
+import cartopy.crs as ccrs
 from metpy.units import units
 import h5py
 import time
+from scipy.interpolate import interp1d
+from util.util import minimize_quadratic
 
 
 goes16_directory = '/arcdata/goes/grb/goes16'  # /year/date/abi/L1b/RadC
@@ -1135,9 +1137,9 @@ def best_fit_altitude(amv_spd, amv_dir, amv_alt, amv_lat, amv_lon, fcst_spd, fcs
         sat_wind_best_fit_u = fcst_uwind[lev_below] * (1.0 - prop) + fcst_uwind[lev_above] * prop
         sat_wind_best_fit_v = fcst_vwind[lev_below] * (1.0 - prop) + fcst_vwind[lev_above] * prop
 
-
     return bf_tup
 
+
 def get_press_bin_ranges(lop, hip, bin_size=100):
     bin_ranges = []
     delp = hip - lop