关于gsn_add_polyline函数的一些疑问

longlivehj

黄小仙儿 发表于 2014-4-24 15:32
哇，出来了。请问这个前后顺序也有影响吗？

放在前面的画，会被后面的图覆盖掉！

黄小仙儿

longlivehj 发表于 2014-4-24 15:57
放在前面的画，会被后面的图覆盖掉！

刚刚用wrf_ll_to_ij算出来的格点

(0)     X/Y location is: 69.3021
(1)     X/Y location is: 132.629
(0)     X/Y location is: 141.581
(1)     X/Y location is: 217.404
来画剖面图，我把格点四舍五入弄成整数，弄到 plane里面plane = (/  69,1332,  141,217 /) ;
然后提取rh_plane = wrf_user_intrp3d(rhs,z,"v",plane,0.,opts)画剖面图。但是提示这一行

fatal:The result of the conditional expression yields a missing value. NCL can not determine branch, see ismissing function
fatal:["Execute.c":8567]:Execute: Error occurred at or near line 167 in file $NCARG_ROOT/lib/ncarg/nclscripts/wrf/WRFUserARW.ncl

fatal:["Execute.c":8567]:Execute: Error occurred at or near line 358 in file $NCARG_ROOT/lib/ncarg/nclscripts/wrf/WRFUserARW.ncl

fatal:["Execute.c":8567]:Execute: Error occurred at or near line 87 in file /home/Huanglei/wrf_CrossSectionsum.ncl
说没有值。
检查了很久，最后换了plane里的值，比如 plane = (/  2,2,  98,98 /) ; start x;y & end x;y point却可以出图。
可以确定是格点的问题，但是不懂是什么原因啊。。。

longlivehj

黄小仙儿 发表于 2014-4-24 16:26
刚刚用wrf_ll_to_ij算出来的格点

(0)     X/Y location is: 69.3021

黄小仙儿，你出过好多这样的错误了哦！

plane = (/  69,1332,  141,217 /)
第二个数是133，不是1332！

黄小仙儿

longlivehj 发表于 2014-4-24 16:33
黄小仙儿，你出过好多这样的错误了哦！

plane = (/  69,1332,  141,217 /)

啊呀~我错了。太粗心了。、、、、这种低级错误~~