ncl作风矢量的垂直剖面图，不知道问题出在哪里

irides

xyan88 发表于 2018-1-12 10:39
请问最后怎么解决的啊，可以提供下代码学习吗

我的代码就是根据贴出的那个链接，结合自己要画的改的，你先试试，有什么问题再一起讨论呀

xyan88

irides 发表于 2018-1-12 10:43
我的代码就是根据贴出的那个链接，结合自己要画的改的，你先试试，有什么问题再一起讨论呀

alfa   = atan2((endlat-startlat),(endlon-startlon))    ;计算线段AB与x轴方向的夹角

对应的end 和start是不是写反了啊？我根据您发的风向画反了

另外alfa结果需要*180/pi 把

irides

xyan88 发表于 2018-1-12 19:52
alfa   = atan2((endlat-startlat),(endlon-startlon))    ;计算线段AB与x轴方向的夹角

对应的end 和s ...

就是end-start，大的减小的
alfa不需要*180/pi呀

xyan88

irides 发表于 2018-1-13 14:55
就是end-start，大的减小的
alfa不需要*180/pi呀

好像我弄错了，多谢啦

xyan88

irides 发表于 2018-1-13 14:55
就是end-start，大的减小的
alfa不需要*180/pi呀

好像我弄错了，多谢啦

xyan88

irides 发表于 2018-1-13 14:55
就是end-start，大的减小的
alfa不需要*180/pi呀

好像我弄错了，多谢啦

xingyue

楼主你好，我也要画类似的图，由于不知道如何做斜线的矢量合成，看了你的答案很受启发。另外看了简介发现你的代码当中有几个错误，可能是导致你画图有问题的原因：
wrf_user_intrp3d这个函数中，使用的角度单位应该是度而不是弧度（rad）,而你用atan2计算出来的角度的单位是rad，需要转换一下；
另外wrf_user_intrp3d这个函数中，角度是与Y轴的角度而不是X轴，即S-N是0度，W-E是90度，所以你这里计算与X轴的夹角是不对的。
顺便，如果你后面换成了与Y轴夹角的话，u,v合成是cos和sin也要换一下。
其实plane是4个变量的话，角度也可以直接设为0，只有在plane是2个变量，即穿过这个点的时候，才设置与Y轴角度。
这是我的几个想法，可能也有问题，但是不同的做法可以尝试一下
祝好~

irides

xingyue 发表于 2018-3-18 17:17
楼主你好，我也要画类似的图，由于不知道如何做斜线的矢量合成，看了你的答案很受启发。另外看了简介发现你 ...

是的，后来修改了一下我的脚本，出来的结果和grads也就差不多了。还有一个问题就是，这样插值出来的垂直层数一定是100层吗？

irides

xingyue 发表于 2018-3-18 17:17
楼主你好，我也要画类似的图，由于不知道如何做斜线的矢量合成，看了你的答案很受启发。另外看了简介发现你 ...

而且插值出来的格点场信息，例如对应的经纬度和高度应该要怎么样获取呢？

irides

xyan88 发表于 2018-1-15 15:55
好像我弄错了，多谢啦

ncl里面sin和cos函数就是直接对弧度进行计算的