snapshot...

modules/util/gfs_reader.py

@@ -4,17 +4,18 @@ import glob
 import numpy as np
 import xarray as xr
 from util.util import value_to_index
+from metpy import *
 # import cartopy
 # from cartopy import *
-from metpy import *
 # import cartopy.crs as ccrs
-import random
 
 # gfs_directory = '/apollo/cloud/Ancil_Data/clavrx_ancil_data/dynamic/gfs/'
 gfs_directory = '/Users/tomrink/data/gfs/'
 gfs_date_format = '%y%m%d'
 
+# force incoming longitude to (0, 360) to match GFS
 lon360 = True
+# GFS half degree resolution
 lat_coords = np.linspace(-90, 90, 361)
 lon_coords = np.linspace(0, 359.5, 720)