snapshot...

modules/util/gfs_reader.py

@@ -187,7 +187,6 @@ def get_vert_profile(xr_dataset, fld_name, lons, lats):
     if lon360:  # convert -180/+180 to 0,360
         lons = np.where(lons < 0, lons + 360, lons)
 
-    plevs = xr_dataset['pressure levels']
 
     fld = xr_dataset[fld_name]
     fld = fld.assign_coords(fakeDim2=lon_coords, fakeDim1=lat_coords, fakeDim0=plevs)