不同密度资料插值gs_re函数

zhi726

在用fnl资料real WRF的过程中，想将WRF的结果与fnl进行对比，但后处理后的WRF资料与fnl资料的密度和层次均不同（其实层次可以修改namelist.ARWpost调成相同，但当时没想起），所以写了一个gs处理下，主要使用re函数（貌似只有grads 2.0以后才可以用），贴上来有需要的可以节约点时间！

GS文件：
'reinit'
'open C:\dcz\20130630\wrfout\wrfout_d01_2013.06.29.12_ssrafs.vert70.ctl'
'set gxout fwrite'
'set fwrite C:\dcz\20130630\wrfout\comp\wrfout\wrfout.ssrafs.d01.redata.v70.grd'

t=1
while(t<=49)
'set t 't''

**********************redata vars*******************
iz=1
while(iz<=21)
'set z 'iz''
'set x 1 615'
'set y 1 358'
'set undef -9999'
'd re(geopt/9.8,141,linear,70,0.5,81,linear,10,0.5,bl)'
if(iz=1)
iz=iz+2
else
iz=iz+1
endif
endwhile

**********************redata end*************

t=t+3
endwhile
'disable fwrite'
'reinit'



CTL文件：
dset C:\dcz\20130630\wrfout\comp\wrfout\wrfout.ssrafs.d01.redata.v70.grd
title ssrafs.redata
undef -9999
xdef 141 linear 70 0.5
ydef 81 linear 10 0.5
tdef 17 linear 12Z29jun2013 3hr
zdef 20 levels
1000 950 925 900 850 800 750 700 650 600 550 500 450 400 350 300 250 200 150 100
vars 1
hgt  20 99 hgtprs_redata(m2/s2/9.8=gpm)
endvars
画出图：