以下是我用fnl资料计算位温所用的文件输出gs程序，很好用，现共享给大家！

lihuijin

'reinit'
'open H:\fnl201308\fnl1.ctl'
'set lat 0 90'
'set lon 30 150'
*'set lev 100'
*'set t 64'
*'define qse='the(TMPprs,lev,RHprs)
'set fwrite H:\fnl201308\gs\qse.dat'
'set gxout fwrite'
t=1
while(t<=64)
'set t 't''
z=1
while(z<=21)
'set z 'z''
'define qse='the(TMPprs,lev,RHprs)
'd qse'
z=z+1
endwhile
t=t+1
endwhile

function the(tgrid,pgrid,rh)
*------------------------------------------
'define es=6.112*exp(17.67*('tgrid'-273.15)/('tgrid'-29.65))'
'define q1='rh'*(0.62197*es/('pgrid'-0.378*es))/100.'
'define e='pgrid'*q1/(0.62197+q1)+1e-10'
'define qse='tgrid'*pow((1000./'pgrid'),(0.2854*(1.0-0.28*q1)))*exp(((3376./(55.0+2840.0/(3.5*log('tgrid')-log(e)-4.805)))-2.54)*q1*(1.0+0.81*q1))'
return(qse)

qxtlyf

谢谢版主了！！！！！！！！！！！！1

小胡00

谢谢版主的无私共享

Уēｓ、尜玉

谢谢楼主分享

Уēｓ、尜玉

谢谢楼主分享

洗雨

楼主，这个GS输出的qse.dat的ctl是不是这样书写？
dset g:/gs/qse.dat
title qse at fnl20110623
undef 9.999E+20
xdef 21 linear 105.000000 1.000000
ydef 13 linear 34.000000 1.000000
tdef 4 linear 00Z23jun2011 6hr
zdef 21 levels 100000 97500 95000 92500 90000 85000 80000 75000 70000 65000 60000 55000 50000 45000 40000 35000 30000 25000 20000 15000 10000
vars 1
qse 21 99 qse level
endvars

ˇBelieve_

谢谢楼主分享！！！！！！！！

刘大壮壮壮

{:lxm_24:}

拖延症

谢谢楼主分享

hxyj

感谢楼主的分享。