diff --git a/modules/util/gfs_reader.py b/modules/util/gfs_reader.py
index cafa3794bcfc7a6c8d1d0b793e3de4d0b76578be..3abc8c90a93f144433e6f5dc4de1c88e4a26fc4f 100644
--- a/modules/util/gfs_reader.py
+++ b/modules/util/gfs_reader.py
@@ -156,7 +156,7 @@ def get_time_interpolated_layer(xr_dataset_s, time_s, time, fld_name, press, lon
     lyr = xr.concat(layer_s, 'time')
     lyr = lyr.assign_coords(time=time_s)
 
-    intrp_lyr = lyr.interp(time=time, method='linear')
+    intrp_lyr = lyr.interp(time=time, method=method)
     intrp_lyr = intrp_lyr.expand_dims('channel')
     intrp_lyr = intrp_lyr.assign_coords(channel=[fld_name])