求可以画出涡度和散度的gs文件

小霸王和小霸王

大神你们好，求可以画出涡度和散度的gs文件，本人想要分析涡度散度和急流的关系，但是不会编写，总是出问题~

weinihou

右上，搜索

有具体问题问问题，有槽也可以考虑吐吐，说一圈其他到头还是伸手党就没意思了

尽头的尽头

grads自带的求涡度和散度的函数，好好学习一下手册

伽蓝鸟

hdivg, hcurl，考虑了球面性

yuchen

* horizontal divergence unit:/s
'define div=hdivg(ugrdprs,vgrdprs)'
* vertical vorticity unit:/s
'define vor=hcurl(ugrdprs,vgrdprs)'
贴出你的错误大家才好帮你啊。

小霸王和小霸王

yuchen 发表于 2018-3-9 13:58
* horizontal divergence unit:/s
'define div=hdivg(ugrdprs,vgrdprs)'
* vertical vorticity unit:/s
...

'reinit'
'sdfopen e:\lcx\ncdata\uwnd.2011.nc'
'sdfopen e:\lcx\ncdata\vwnd.2011.nc'

'set grads off'
'set grid off'
'set lon 20 180'
'set lat 0 80'
'set z 3'
'set t 30'
'define div=hdivg(uwnd.1,vwnd.2)'
'd div'
'draw title 2011'
'printim e:\lcx\ncdata\sandu.png'
'reinit'
'reinit'
散度


'reinit'
'reinit'
'sdfopen e:\lcx\ncdata\uwnd.2011.nc'
'sdfopen e:\lcx\ncdata\vwnd.2011.nc'

'set grads off'
'set grid off'
'set lon 20 180'
'set lat 0 80'
'set lev 500'
'set t 1'
'define vort=hcurl(uwnd.1,vwnd.2)'
'd vort'
'draw title 2011'
'printim e:\lcx\ncdata\wodu.png'
'reinit'
'reinit'

涡度

这两者都是显示的不是数值

yuchen

d uwnd.1和vwnd.2看看风速显示正常不。

小霸王和小霸王

yuchen 发表于 2018-3-10 18:08
d uwnd.1和vwnd.2看看风速显示正常不。

风速显示的是正常的，只是显示的不是数值，是数字和字母的组合~

小霸王和小霸王

yuchen 发表于 2018-3-10 18:08
d uwnd.1和vwnd.2看看风速显示正常不。

风速显示的是正常的，只是显示的不是数值，是数字和字母的组合~

小霸王和小霸王

yuchen 发表于 2018-3-10 18:08
d uwnd.1和vwnd.2看看风速显示正常不。

风速显示的是正常的，只是显示的不是数值，是数字和字母的组合~