diff --git a/modules/util/gfs_reader.py b/modules/util/gfs_reader.py
index 3b000a0a0d0e4a264ff9702e368135c1ecccb057..49c67d67a2a15b75e61c5f52a459779cd9504c85 100644
--- a/modules/util/gfs_reader.py
+++ b/modules/util/gfs_reader.py
@@ -171,7 +171,7 @@ def get_time_interpolated_layer(xr_dataset_s, time_s, time, fld_name, press, lon
 def get_time_interpolated_layer_s(xr_dataset_s, time_s, time, fld_name_s, press, lon_range=None, lat_range=None, method='linear'):
     layer_s = []
     for ds in xr_dataset_s:
-        lyr = get_horz_layer_s(ds, fld_name_s, press, lon_range=lon_range, lat_range=lat_range, method=method)
+        lyr = get_horz_layer_s(ds, fld_name_s, press, lon_range=lon_range, lat_range=lat_range)
         layer_s.append(lyr)
 
     lyr = xr.concat(layer_s, 'time')