From 083e45690d298bdb4ab61df59b4fc7aaabcaea3d Mon Sep 17 00:00:00 2001
From: tomrink <rink@ssec.wisc.edu>
Date: Sun, 6 Dec 2020 11:26:35 -0600
Subject: [PATCH] minor

---
 modules/util/gfs_reader.py | 8 ++++----
 1 file changed, 4 insertions(+), 4 deletions(-)

diff --git a/modules/util/gfs_reader.py b/modules/util/gfs_reader.py
index e7cea8f3..eb73761d 100644
--- a/modules/util/gfs_reader.py
+++ b/modules/util/gfs_reader.py
@@ -187,7 +187,7 @@ def get_time_interpolated_layer_s(xr_dataset_s, time_s, time, fld_name_s, press,
     return intrp_lyr
 
 
-def get_vert_profile(xr_dataset, fld_name, lons, lats):
+def get_vert_profile(xr_dataset, fld_name, lons, lats, method='linear'):
     if lon360:  # convert -180/+180 to 0,360
         lons = np.where(lons < 0, lons + 360, lons)
 
@@ -198,12 +198,12 @@ def get_vert_profile(xr_dataset, fld_name, lons, lats):
     dim2 = xr.DataArray(lons, dims='k')
     dim1 = xr.DataArray(lats, dims='k')
 
-    intrp_fld = fld.interp(fakeDim1=dim1, fakeDim2=dim2, fakeDim0=plevs, method='linear')
+    intrp_fld = fld.interp(fakeDim1=dim1, fakeDim2=dim2, fakeDim0=plevs, method=method)
 
     return intrp_fld
 
 
-def get_vert_profile_s(xr_dataset, fld_name_s, lons, lats):
+def get_vert_profile_s(xr_dataset, fld_name_s, lons, lats, method='linear'):
     if lon360:  # convert -180,+180 to 0,360
         lons = np.where(lons < 0, lons + 360, lons)
 
@@ -219,7 +219,7 @@ def get_vert_profile_s(xr_dataset, fld_name_s, lons, lats):
     dim2 = xr.DataArray(lons, dims='k')
     dim1 = xr.DataArray(lats, dims='k')
 
-    intrp_fld = fld.interp(fakeDim1=dim1, fakeDim2=dim2, fakeDim0=plevs, method='linear')
+    intrp_fld = fld.interp(fakeDim1=dim1, fakeDim2=dim2, fakeDim0=plevs, method=method)
 
     return intrp_fld
 
-- 
GitLab