[已解决]grads对fnl资料绘制湿位涡，出现点问题，求高手指点

quentin011099

请问楼主可不可以借鉴下这个GS相对应的fnl的描述文件？先谢谢了

良辰

fnl直接解码的，1*1资料，原gs找不到了也

o_O_小岛

良辰 发表于 2012-9-30 16:35

                               
登录/注册后可看大图

mikado 说的有道理，但是不是主要解决办法，主要是由于插值造成的，先进行范围较大的插值计算，绘图时候选择 ...

请问楼主，我也是用fnl资料出图，不过还是只有一部分图片，请问第二个问题具体是如何解决的?

.雾霭流岚

我用楼主的程序出不来图啊？错误是“cannot contourgrid - all undefined values”

云巅

楼上的问题该如何解决呢？

xiaoyuyly

2013-10-22
.xiaoyuyly: 关于你在“[已解决]grads对fnl资料绘制湿位涡，出现点问题，求高手指点”的帖子
你好，我纠结湿位涡垂直剖面好几天了，一直出不来图，后来找到了你的这段程序，我就copy过去，把时间行去掉，结果还是不出来啊，请问是什么原因？请你帮我看看程序好吗？非常非常非常感谢
'reinit'
'open E:\2012\ncep\20130823\2013082412.ctl'
'set lon 90 120'
'set lat 10 40'
'set lev 1000 100'
'set z 2 20'
'set zlog on'
*假相当位温**
'define tc=tmpprs-273.15'
'define rh=rhprs'
'define prs=lev'  
'define es=(6.112*exp((17.67*tc)/(tc+243.5)))'
'define qs=0.62197*es/(prs-0.378*es)'
'define q=rh*qs/100'
'define e=prs*q/(0.62197+q)+1e-10'
'define tlcl=55.0+2840.0/(3.5*log(tc+273.16)-log(e)-4.805)'
'define theta=(tc+273.16)*pow((1000/prs),(0.2854*(1.0-0.28*q)))'
'define eqt=theta*exp(((3376./tlcl)-2.54)*q*(1.0+0.81*q))'
*湿位涡*
'define u=UGRDprs'
'define v=VGRDprs'
'define vo=hcurl(u,v)'
'define f=2*7.292*sin(lat*3.14159/180.0)*0.00001'
'define g=9.8'
'define dp=100*(lev(z-1)-lev(z+1))'
'define deqt=eqt(z-1)-eqt(z+1)'
'define du=u(z-1)-u(z+1)'
'define dv=u(z-1)-u(z+1)'
'define dx=2.0*6370949.0*cos(lat*3.14159/180.0)*3.14159/180.0'
'define dy=2.0*6370949.0*3.14159/180.0'
'define dtx=cdiff(eqt,x)'
'define dty=cdiff(eqt,y)'
'define pv1=-g*(vo+f)*deqt/dp'  
'define pv2=g*((dv/dp)*(dtx/dx)-(du/dp)*(dty/dy))'
'define pv=pv1+pv2'
'set lon 95.5'
'set lat 18 32'
'set mpdset yunnan'
'set map 1 1 6'
'set grid off'
'set grads off'
'set gxout contour'
'set csmooth on'
'set cthick 7'
*'set clopts -1 -1 0.10'
'd pv1*10e5'
*'d pv2*10e6'
*'d pv*10e5'
'printim e:\2012\ncep\20130823\2013082412swo.gif gif x800 y600 white'
'enable print  e:\2012\ncep\20130823\2013082412swo.gmf'
'print'
'disable priint'
错误提示：
E:\2012\ncep\20130823\QQ图片20131022210849.jpg

良辰

xiaoyuyly 发表于 2013-10-22 21:15

                               
登录/注册后可看大图

2013-10-22
.xiaoyuyly: 关于你在“[已解决]grads对fnl资料绘制湿位涡，出现点问题，求高手指点”的帖子
...

看不到你的错误提示

xiaoyuyly

非常感谢您能看我给您的回复，错误提示我本来是用图片贴上来的，结果居然显示不出来，我现在把错误提示发给您，请你帮我分析解决一下好吗？非常非常感谢
错误提示：
all files closed; all defined objects released;
all GRADS attributes have been reinitialized
cannot contour grid ----all undefined value
hardcopy output file is closed

良辰

xiaoyuyly 发表于 2013-10-23 11:46

                               
登录/注册后可看大图

非常感谢您能看我给您的回复，错误提示我本来是用图片贴上来的，结果居然显示不出来，我现在把错误提示发给 ...

个人觉得不是程序本身的问题，应该是你grads的问题，你换台机器试一下，在设置一下set t 1  ，我对grads了解也不多，就是用用而已。这个程序我画出过图，应该不是grads命令的问题。你再试一试，或找周边的朋友帮你看看。对了，你用的是fnl数据吗，解码方式可能不一样呀，你打开你的ctl看看里面的变量名和这个gs里面的能对应上不

xiaoyuyly

良辰 发表于 2013-10-24 20:08

                               
登录/注册后可看大图

个人觉得不是程序本身的问题，应该是你grads的问题，你换台机器试一下，在设置一下set t 1  ，我对grads了 ...

变量名没问题，我仔细核对过的，用的是fnl数据，解码方式是用grib2解码的，解码应该不会有问题，因为我用这个资料画出来其它的图，如thse的垂直剖面等等，都是正常的，
我试了一次又一次，就是出不来，我都要崩溃了，耽误了我整整一个礼拜，哎！！
你说的有点道理，也许是我的grads本身的问题，我找别人帮我画画看，非常感谢您，不懂得我再请教您。