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

begin

prefix = "gfsanl"
times = ispan(1,85,1)

suffix = times

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

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

ncout = addfile("gfsanl_"+monday+".nc","c")
ncout->UGRD_P0_L100_GLL0 = data

end