snapshot...

modules/util/gfs_reader.py

@@ -344,7 +344,7 @@ def get_time_interpolated_point_s(xr_dataset_s, time_s, fld_name_s, time, lons,
     return intrp_pt
 
 
-def get_time_interpolated_voxel(xr_dataset_s, time_s, time, fld_name, lon, lat, press, x_width=5, y_width=5, z_width=3, method='linear'):
+def get_time_interpolated_voxel(xr_dataset_s, time_s, time, fld_name, lon, lat, press, x_width=3, y_width=3, z_width=3, method='linear'):
     vox_s = []
     for ds in xr_dataset_s:
         vox = get_voxel(ds, fld_name, lon, lat, press, x_width=x_width, y_width=y_width, z_width=z_width)
@@ -358,7 +358,7 @@ def get_time_interpolated_voxel(xr_dataset_s, time_s, time, fld_name, lon, lat,
     return intrp_vox
 
 
-def get_voxel(xr_dataset, fld_name, lon, lat, press, x_width=5, y_width=5, z_width=3):
+def get_voxel(xr_dataset, fld_name, lon, lat, press, x_width=3, y_width=3, z_width=3):
     if lon360:
         if lon < 0:
             lon += 360
@@ -393,7 +393,7 @@ def get_voxel(xr_dataset, fld_name, lon, lat, press, x_width=5, y_width=5, z_wid
     return sub_fld
 
 
-def get_time_interpolated_voxel_s(xr_dataset_s, time_s, time, fld_name_s, lon, lat, press, x_width=5, y_width=5, z_width=3, method='linear'):
+def get_time_interpolated_voxel_s(xr_dataset_s, time_s, time, fld_name_s, lon, lat, press, x_width=3, y_width=3, z_width=3, method='linear'):
     vox_s = []
     for ds in xr_dataset_s:
         vox = get_voxel_s(ds, fld_name_s, lon, lat, press, x_width=x_width, y_width=y_width, z_width=z_width)
@@ -407,7 +407,7 @@ def get_time_interpolated_voxel_s(xr_dataset_s, time_s, time, fld_name_s, lon, l
     return intrp_vox
 
 
-def get_voxel_s(xr_dataset, fld_name_s, lon, lat, press, x_width=5, y_width=5, z_width=3):
+def get_voxel_s(xr_dataset, fld_name_s, lon, lat, press, x_width=3, y_width=3, z_width=3):
     if lon360:
         if lon < 0:
             lon += 360