sep = ','
output = '/home/subdomains/forecasts/public_html/tmp/|outputfilename|'

'open /home/subdomains/forecasts/public_html/grib/grib2/|rundate||runhour|/gfs.|runhour|z.ctl'
'set t plottime'
'set lat |latreplace|'
'set lon |lonreplace|'

'd lftxsfc'
li=subwrd(result,4)

'tc=(tmpprs-273.16)'
'd tc-( (14.55+0.114*tc)*(1-0.01*rhprs) + pow((2.5+0.007*tc)*(1-0.01*rhprs),3) + (15.9+0.117*tc)*pow((1-0.01*rhprs),14) )'
dp=subwrd(result,4)

'set lev 850'
'd tmpprs'
t850=subwrd(result,4)
t850=t850-273.16

'd ugrdprs'
v1 = subwrd(result,4)

values = |rundate||runhour| sep li sep dp sep t850
write (output,values,append)

'quit'