求助：gs - 整层水汽通量积分计算 时间设置，dat结果，gr2stn插值

yutuzhi

程序如下，非常希望大家指出错误或者提供帮助：


'reinit'
#数据都是sub-daily的，一天4次数据，这里尝试画出2009年11月10日6:00-12:00地表到700hpa高度的水汽积分图。
'sdfopen F:\NCEP\shum.2009.nc'
'sdfopen F:\NCEP\uwnd.2009.nc'
'sdfopen F:\NCEP\vwnd.2009.nc'
'sdfopen F:\NCEP\pres.sfc.2009.nc'
'enable print F:\20091110surface-700hpa.gmf'

'set lon 0 180'
'set lat 0 90'
'set z 1 4'
#17 Pressure levels (mb): 1000,925,850,700,600,500,400,300,250,200,150,100,70,50,30,20,10
'set t 1254'
#2009年11月10日6:00-12:00

'define uq=shum.1*uwnd.2*1000'
'define vq=shum.1*vwnd.3*1000'
'define uqs=vint(pres.4(lev=0)/100,uq,700)/9.8'
'define vqs=vint(pres.4(lev=0)/100,vq,700)/9.8'
'define a=mag(uqs,vqs)'

'set gxout vector'
'set z 1'
'set t 1'
##？？？这里应该是t=1,还是和之前的设置一样 t=1254啊？不明白这里的设置是什么意思？？求帮助
'd uqs;vqs'
'draw title 2009/11/10 surface-700hpa'
'print'
'c'

'set fwrite F:20091110surface-700hpa.dat'
#得到的这个文件总是打不开，或者是乱码。后续我不知道该怎么用grads做插值，所以想用MATLAB，可是打不开这个数据。请问应该怎么办呢？
'set gxout fwrite'
'd a'
'disable fwrite'

'reinit'

希望大神们能够提供帮助！！非常感激！

yutuzhi

发现了自己的两个错误，'define uqs=vint(pres.4(lev=0)/100,uq,700)/9.8'中，pres.4的单位是10pa，所以应该除以10，然后vint函数中已经除以9.8了，所以这里也不需要再计算了。改正过来应该是：
'define uqs=vint(pres.4(lev=0)/10,uq,700)'
'define vqs=vint(pres.4(lev=0)/10,vq,700)'

nunu18

'set fwrite F:20091110surface-700hpa.dat' 少了斜杠/

yutuzhi

nunu18 发表于 2016-4-7 16:43
'set fwrite F:20091110surface-700hpa.dat' 少了斜杠/

多谢指出，我贴上来的时候出错了。自己运行的时候是ok的。下次会注意。您知道时间t应该怎么设置吗？

nunu18

yutuzhi 发表于 2016-4-7 16:46
多谢指出，我贴上来的时候出错了。自己运行的时候是ok的。下次会注意。您知道时间t应该怎么设置吗？

你要画某一个时次的水汽通量积分，所以前后时间设置都一样，前面'define uqs=vint(pres.4(lev=0)/100,uq,700)/9.8'
'define vqs=vint(pres.4(lev=0)/100,vq,700)/9.8'
'define a=mag(uqs,vqs)'是u向水汽通量和V向水汽通量积分，最后你输出的是水汽通量积分场和水汽通量流线场，个人理解，这样你画出来的图才好看

yutuzhi

nunu18 发表于 2016-4-7 17:01
你要画某一个时次的水汽通量积分，所以前后时间设置都一样，前面'define uqs=vint(pres.4(lev=0)/100,uq, ...

谢谢！所以我后面就不用再设置一次时间了。太棒了，解决了第一个问题。感谢感谢！

yutuzhi

我找到了另一个方法去实现了，虽然有点儿笨。第二个问题可能是由于nc格式的不能直接存成dat文件吧，我输出的结果总是乱码。现在我就直接输出nc格式啦。然后直接用MATLAB打开nc格式的数据处理。那么第三个问题也解决了，不用grads插值了，直接在MATLAB当中插值计算。