Update gfs_reader.py

modules/util/gfs_reader.py

@@ -225,7 +225,7 @@ def get_vert_profile_s(xr_dataset, fld_name_s, lons, lats):
     return intrp_fld
 
 
-def get_point(xr_dataset, fld_name, lons, lats):
+def get_point(xr_dataset, fld_name, lons, lats, method='linear'):
     if lon360:  # convert -180/+180 to 0,360
         lons = np.where(lons < 0, lons + 360, lons) # convert -180,180 to 0,360
 
@@ -238,7 +238,7 @@ def get_point(xr_dataset, fld_name, lons, lats):
     dim1 = xr.DataArray(lats, dims='k')
     dim2 = xr.DataArray(lons, dims='k')
 
-    intrp_fld = fld.interp(fakeDim1=dim1, fakeDim2=dim2, method='linear')
+    intrp_fld = fld.interp(fakeDim1=dim1, fakeDim2=dim2, method=method)
 
     return intrp_fld