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

supersupergirl · 发表于 2019-8-24 11:56:44

大神还有个问题就是zmax_pos是特定词吗为什么我用zmax_pos1或者zmax_pos2就不行呢

supersupergirl · 发表于 2019-8-26 21:00:38

本帖最后由 supersupergirl 于 2019-8-27 09:11 编辑

还有我想问下大神我如果画风场的涡度散度涡度平流等图还需要处理uv嘛（这个不用回答了是我蠢）

supersupergirl · 发表于 2019-8-27 09:13:36

早上好今天也是等解答的一天

wyanghu · 发表于 2019-8-27 17:31:26

请问楼主，在选取剖面算alfa的时候，线段AB与y轴的角度有可能大于90度，这种情况要如何处理呢

wyanghu · 发表于 2019-8-27 17:36:06

换一种提法，如果在算alfa时，atan2的y和x如果有一个是负数或两个都为负数，要咋办呢，还和两者都是正值的情况做相同处理吗

zhengjingjing · 发表于 2019-10-19 23:48:36

您好楼主下载了你脚本后在284行if (abs(zz(zmax_pos,0)-zmax*1000.) .lt. abs(zz(zmax_pos+1,0)-zmax*1000.) ) then，一直报错（An error occurred reading zz）,可以请教一下是怎么回事么？

小其其格 · 发表于 2019-10-20 07:01:38

wyanghu 发表于 2019-8-27 17:36
换一种提法，如果在算alfa时，atan2的y和x如果有一个是负数或两个都为负数，要咋办呢，还和两者都是正值的 ...

同样处理就行了

小其其格 · 发表于 2019-10-20 07:04:43

zhengjingjing 发表于 2019-10-19 23:48
您好楼主下载了你脚本后在284行if (abs(zz(zmax_pos,0)-zmax*1000.) .lt. abs(zz(zmax_pos+1,0)-zmax*1 ...

那说明变量z读取错误，应该下面这句读取有问题呀
zz = wrf_user_intrp3d(z,z,"v",plane,angle,opts)

你确认没修改过代码？

zhengjingjing · 发表于 2019-10-21 10:56:36

小其其格发表于 2019-10-20 07:04
那说明变量z读取错误，应该下面这句读取有问题呀
zz = wrf_user_intrp3d(z,z,"v",plane,angle,opts)

只修改了路径，把前面的print注释掉了，应该不会有啥问题呀，郁闷

zhengjingjing · 发表于 2019-10-21 10:59:24

zhengjingjing 发表于 2019-10-21 10:56
只修改了路径，把前面的print注释掉了，应该不会有啥问题呀，郁闷

z我可以print 出来

		自动登录	找回密码
密码			立即注册