diff --git a/modules/amv/intercompare.py b/modules/amv/intercompare.py index afe68d6ac35c59e3669984069677d3813fb302ff..ea3c9635dba4cc703fc49eb7b75173640abd14d4 100644 --- a/modules/amv/intercompare.py +++ b/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