; This script takes downloaded GEFS ensemble grib files, selects only the variables of interest, 
; and outputs a netCDF file.

begin

prefix = "gefs"
times = ispan(1,40,1)

suffix = times

inf = addfiles(prefix+suffix+".grb", "r")

ListSetType(inf,"join")
data := inf[:]->UGRD_P1_L100_GLL0
data!0 = "time"
data&time = times

ncout = addfile("gefsens"+ens+"_"+monday+".nc","c")
ncout->UGRD_P1_L100_GLL0 = data

end