用gs脚本画辐散风场以及势函数

chuanhao

chuanhao 发表于 2014-3-4 10:52

                               
登录/注册后可看大图

前边有人提供了一个grads函数的网站啊，上面有算速度势的函数，也就是作者用的那个函数，还有计算流函数的 ...

下边是我的gs，可能我改的不对，希望你帮我看看，谢谢！

鲸鱼便子

chuanhao 发表于 2014-3-4 10:53

                               
登录/注册后可看大图

'reinit'
'open F:\xin\sqtl\qu.ctl'
'open F:\xin\sqtl\qv.ctl'

我当初是试着画的是水汽通量的，不知道风场的量级大概是多少。
我觉得你这里 是不是得再看看~~~~
lz给的网站里面对旋转分量是这样定义的 ，我记得mvadv和muadv作用不一样的（分量的方向），而且根据公式负号也得注意一下~~~~~
暂时只看出这里的问题了，你先试试看~~~~~

chuanhao

鲸鱼便子 发表于 2014-3-4 17:13

                               
登录/注册后可看大图

我当初是试着画的是水汽通量的，不知道风场的量级大概是多少。
我觉得你这里是不是得再看看~~~~
lz给的 ...

恩，速度势越大，能说明什么？你能给我说一下么，看书我有点看不懂

鲸鱼便子

chuanhao 发表于 2014-3-5 12:06

                               
登录/注册后可看大图

恩，速度势越大，能说明什么？你能给我说一下么，看书我有点看不懂

这个我也一直没能很好的理解~
我觉得势函数对应的是非旋转分量（也就是辐合辐散），流函数对应的是旋转分量

chuanhao

楼主你好，速度势越大代表辐合越强烈么？看书我有点看不懂，麻烦解释一下

Sarah

若是把辐散风和速度势分开两个gs计算，就不会造成边界问题。保证辐散风和速度势都在全球格点下进行计算。同时写在一个gs里面可能会容易出错。

nuistajun

river 发表于 2012-10-18 11:32
真是高手啊，这么快就搞懂了啊。我用你的gs(稍微改动了一下）画了一张。

大神能否给我看一下我的gs，我想x轴是经度，Y轴表示高度。

'reinit'
'open C:\Users\Lenovo\Desktop\wind\positive\u1.ctl'
'open C:\Users\Lenovo\Desktop\wind\positive\v1.ctl'
'open C:\Users\Lenovo\Desktop\wind\positive\w1.ctl'
'open C:\Users\Lenovo\Desktop\wind\positive\u3.ctl'
'open C:\Users\Lenovo\Desktop\wind\positive\v3.ctl'
'open C:\Users\Lenovo\Desktop\wind\positive\w3.ctl'
'set x 1 144'
'set y 1 73 '
'set z 1 12'
'set annot 8 8'
'set mpdraw on'
'set map auto'
'set mpdset cnworld'
'set grads off'
'set grid off'
'define chi=fish_chi(uwnd.1-uwnd.4,vwnd.2-vwnd.5)'
*pull dummy
'chi=chi-aave(chi,global)'
'set gxout contour'
'set cthick 9'
'set lon 50 300'
'set lat 0'
'set z 1 12'
'set zlog on'
'd ave(chi/1e6,lat=-10,lat=20)'
*pull dummy
'set gxout vector'
'set lon 50 300'
'set lat 0'
'set z 1 12'
'set zlog on'
'd ave(uwnd.1-uwnd.4,lat=-10,lat=20); ave(omega.3-omega.6,lat=-10,lat=20)'
*pull dummy
'printim C:\Users\Lenovo\Desktop\wind\positive\grid('a').png white'
;

nuistajun

nuistajun 发表于 2014-3-30 14:36
大神能否给我看一下我的gs，我想x轴是经度，Y轴表示高度。

'reinit'

忘了补充了我这个gs没有画出来，以下是我处理数据的gs

'reinit'
'sdfopen C:\Users\Lenovo\Desktop\wind\uwnd.mon.mean.nc'
'sdfopen C:\Users\Lenovo\Desktop\wind\omega.mon.mean.nc'
'set x 1 144'
'set y 1 73 '
'set z 1 12'
'define a1=(uwnd.1(time=JUN1979)+uwnd.1(time=JUL1979)+uwnd.1(time=AUG1979))/3'
'define b1=(omega.2(time=JUN1979)+omega.2(time=JUL1979)+omega.2(time=AUG1979))/3'
'define a2=(uwnd.1(time=JUN1983)+uwnd.1(time=JUL1983)+uwnd.1(time=AUG1983))/3'
'define b2=(omega.2(time=JUN1983)+omega.2(time=JUL1983)+omega.2(time=AUG1983))/3'
'define a3=(uwnd.1(time=JUN1987)+uwnd.1(time=JUL1987)+uwnd.1(time=AUG1987))/3'
'define b3=(omega.2(time=JUN1987)+omega.2(time=JUL1987)+omega.2(time=AUG1987))/3'
'define a4=(uwnd.1(time=JUN1988)+uwnd.1(time=JUL1988)+uwnd.1(time=AUG1988))/3'
'define b4=(omega.2(time=JUN1988)+omega.2(time=JUL1988)+omega.2(time=AUG1988))/3'
'define a5=(uwnd.1(time=JUN2003)+uwnd.1(time=JUL2003)+uwnd.1(time=AUG2003))/3'
'define b5=(omega.2(time=JUN2003)+omega.2(time=JUL2003)+omega.2(time=AUG2003))/3'
'define a6=(uwnd.1(time=JUN2007)+uwnd.1(time=JUL2007)+uwnd.1(time=AUG2007))/3'
'define b6=(omega.2(time=JUN2007)+omega.2(time=JUL2007)+omega.2(time=AUG2007))/3'
'define uu=(a1+a2+a3+a4+a5+a6)/6'
'define ww=(b1+b2+b3+b4+b5+b6)/6'
'set  fwrite C:\Users\Lenovo\Desktop\wind\u1.grd '
'set gxout fwrite'
nz=1
while(nz<=12)
'set z 'nz''
'set t 1'
'd uu'
nz=nz+1
endwhile
'disable fwrite'
'set  fwrite C:\Users\Lenovo\Desktop\wind\w1.grd '
'set gxout fwrite'
nz=1
while(nz<=12)
'set z 'nz''
'set t 1'
'd -ww*100'
nz=nz+1
endwhile
'disable fwrite'
;

river

nuistajun 发表于 2014-3-30 14:36
大神能否给我看一下我的gs，我想x轴是经度，Y轴表示高度。

'reinit'

这么复杂的东西，我一个人看不过来，楼主还是单独发个帖子提问吧，请尽量把问题描述清楚

nuistajun

river 发表于 2014-3-30 18:42
这么复杂的东西，我一个人看不过来，楼主还是单独发个帖子提问吧，请尽量把问题描述清楚

http://bbs.06climate.com/forum.p ... mp;page=1#pid269973
大神我已经发帖子啦  十分感谢 能否帮我看一下 实在是搞不出来