NCL绘制wrfout风场垂直剖面并叠加地形

小其其格

NCl菜鸟一枚，原来论文绘图基本是用Grads，但是Grads代码没有编程逻辑可言。

用了NCL发现这代码编写更人性化。参照官网例子绘制了wrfout的风场和云水混合比的垂直剖面图，

并叠加地形。

参照例子为：http://www2.mmm.ucar.edu/wrf/OnL ... smooth_terrain4.ncl

需要特别注意的是，函数atan2计算出来的角度是以正北为0度的，因此要算出直线AB与X轴的交角，就需要换算一下，具体我代码里有标注。

另外wrf_user_ll_to_ij函数把经纬度转换为屏幕坐标点后，是按一下规则来对应经纬度最大最小值的
        ;--------------------------------------------------------------------------
        ;经纬度转为网格点后对应的最小最大值
        ;--------------------------------------------------------------------------
        ; Xstart = loc(0,0)         ;经度最小
        ; Xend   = loc(0,1)         ;经度最大
        ; Ystart = loc(1,0)         ;纬度最小
        ; Yend   = loc(1,1)         ;纬度最大


最后自己设定了坐标qcloud_uw,并用gsn_define_colormap函数读取。

经过版友@Bozhero 的认真修订，楼主发现原来绘图脚本中有些错误，具体为：1、atan2函数求角度问题
2、求取uv问题
现在作了修正，具体脚本以及出图效果如下：



脚本：

随缘

{:eb502:}{:eb502:}

随缘

vector = wrf_vector(a,wks,u_plane(0:zmax_pos,:),w_plane(0:zmax_pos,:)*20.,vcres)
这个乘以20是啥意思？

小其其格

随缘 发表于 2018-6-7 10:21
vector = wrf_vector(a,wks,u_plane(0:zmax_pos,:),w_plane(0:zmax_pos,:)*20.,vcres)
这个乘以20是啥意思 ...

因为实际情况来说，w比u和v小一个量级，因此在做u-w剖面图时，一般要放大10倍或者20倍

m337872010

{:eb502:}{:eb502:}

_﹎__╱/●勄

楼主，请教一下，我想把剖面的纬度缩小一点也就是X轴范围变小，应该怎么办，用哪个命令？

Bozhero

谢谢楼主分享

stcya

谢谢楼主分享！

雨萌萌

感谢楼主分享！我试试看！

雨萌萌

楼主您好！我用了你的程序试着画了一下，可是地形与物理量对不上，如下图C:\Users\kristina\Desktop\plt_CrossSection4.000003是怎么回事呢