diff --git a/modules/amv/intercompare.py b/modules/amv/intercompare.py index 395cda0480962f0e87c303ec158180792ce99748..53bf931d32fc2ae78d3ff7facae77ab3a1b0bb9a 100644 --- a/modules/amv/intercompare.py +++ b/modules/amv/intercompare.py @@ -1,29 +1,25 @@ import numpy as np import xarray as xr -from util.setup import home_dir +# from util.setup import home_dir from util.lon_lat_grid import earth_to_indexs from util.util import haversine_np from util.gfs_reader import get_vert_profile_s -import cartopy +# import cartopy from cartopy import * import metpy.calc as mpcalc import cartopy.crs as ccrs from metpy.units import units import h5py import time -from scipy.interpolate import interp1d +# from scipy.interpolate import interp1d from util.util import minimize_quadratic from netCDF4 import Dataset import datetime from datetime import timezone from util.util import LatLonTuple -#-- AMV intercompare stuff ------------------------------------------ - -class MyGenericException(Exception): - def __init__(self, message): - self.message = message - +# -- AMV intercompare stuff ------------------------------------------ +# -------------------------------------------------------------------- amv_hdr_list = ['lat', 'lon', 'tbox', 'sbox', 'spd', 'dir', 'pres', 'lowl', 'mspd', 'mdir', 'alb', 'corr', 'tmet', 'perr', 'qi', 'cqi', 'qif'] @@ -45,6 +41,7 @@ amv_qif_idx = amv_hdr_list.index('qif') amv_centers_list = ['EUM', 'BRZ', 'JMA', 'KMA', 'NOA', 'NWC'] + def get_amv_nd(filename): header = None data = []