爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

新浪微博登陆

只需一步, 快速开始

搜索
楼主: hustlh

[分享资料] 求1948-2011年的每月平均水汽通量,每月平均水汽通量图一样,不知道是哪儿错了?

[复制链接]

新浪微博达人勋

发表于 2012-8-15 17:08:53 | 显示全部楼层
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2012-8-15 17:11:48 | 显示全部楼层
冰点 发表于 2012-8-15 17:08
太惭愧了,多谢版主才是!赔钱赔贡献的……

呵呵,有愧于大家啊!算是共同进步了1
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2012-8-15 17:13:25 | 显示全部楼层
做个霸气的木头 发表于 2012-8-15 17:00
总算理解对了,循环的东西你试试吧,不难的!

恩,好的,这个我会。要还有不懂,还望赐教啊!

密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2012-8-15 17:23:05 | 显示全部楼层
冰点 发表于 2012-8-15 17:08
太惭愧了,多谢版主才是!赔钱赔贡献的……

说是呢,际遇太差了
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2012-8-15 17:23:52 | 显示全部楼层
hustlh 发表于 2012-8-15 17:13
恩,好的,这个我会。要还有不懂,还望赐教啊!

共同提高,加油吧~
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2012-8-20 10:46:03 | 显示全部楼层
做个霸气的木头 发表于 2012-8-15 17:00
总算理解对了,循环的东西你试试吧,不难的!

版主,你好,上次你说的最后一段加时间t循环,我加了一下,确实是每个月循环,但是所有图的都覆盖在一起了,最后只出了一张图。最近一直在网上看帖子,还是没想到解决办法,我估计可能还是最后一段的命令写错了,也不知道怎样改。
怎样可以输出12张图,就是让它们分别一张一张存储起来,不要叠加在一张图上了???

还得麻烦你帮我看看啊。

最后一段的部分GS:
i=1
while(i<=12)
'set t 'i''
'set gxout shaded'
'd w1'
'set grads off'
'set gxout contour'
'set cint 40'
'd w1'
'set grads off'
'set gxout vector'
'd w2;w3'
'run cbarn'
'printim g:\whole.tif x1000 y800 white'
'print'
i=i+1
endwhile
'disable print'
;


wh.jpg
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2012-8-20 12:35:42 | 显示全部楼层
hustlh 发表于 2012-8-20 10:46
版主,你好,上次你说的最后一段加时间t循环,我加了一下,确实是每个月循环,但是所有图的都覆盖在一起了 ...

grads存图有两种格式啦,一种是printim的,一种是enable print的,你两种都混着来的,这是基础知识,你选择一种吧。想不覆盖,在i=i+1前面加个'c'就行了,意思是每次画完图清下屏。
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2012-8-20 15:07:15 | 显示全部楼层
做个霸气的木头 发表于 2012-8-20 12:35
grads存图有两种格式啦,一种是printim的,一种是enable print的,你两种都混着来的,这是基础知识,你选 ...

恩,解决了,多谢斑竹了!
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2012-8-20 16:10:17 | 显示全部楼层
做个霸气的木头 发表于 2012-8-20 12:35
grads存图有两种格式啦,一种是printim的,一种是enable print的,你两种都混着来的,这是基础知识,你选 ...

惭愧啊,,我还有个问题想请教哈你,还望斑竹不吝赐教啊。

关于求64年(1948-2011)年平均水汽通量,我初步思路是:将求出的1-12月的一共12个月的月平均水气通量求和然后再除以12,自己改了一下脚本命令,但最后求出的图有问题,还得你帮我看看是哪儿出了问题??? 是红色标注的命令有问题吗???

'reinit'
'set grid off'
'set grads off'
'set poli on'
'set parea 1 10 1 8'
'sdfopen f:\data\monthly\shum.mon.mean.nc'  
'sdfopen f:\data\monthly\uwnd.mon.mean.nc'
'sdfopen f:\data\monthly\vwnd.mon.mean.nc'
'sdfopen f:\data\monthly\pres.mon.mean.nc'
'set lat 4 53'
'set lon 73 135'
'set ccols 9 14 4 11 5 13 3 10 7 12 8 2 6'
'set mpdset cnworld hires'
'set map 1 1 1'
'draw map'
'set ylopts 1 4 0.12'
'set xlopts 1 4 0.12'
'set cthick 4'
'set t 1 768'
'set z 1 8'
'define uq=shum.1*uwnd.2'
'define vq=shum.1*vwnd.3'
'set z 1'
'set t 1 768'
'define uqs=vint(pres.4(z=1),uq,300)/1000'
'define vqs=vint(pres.4(z=1),vq,300)/1000'
'define sq=mag(uqs,vqs)'
'set t 1 12'
'define w1=ave(sq,t+0,t=768,12)'
'define w2=ave(uqs,t+0,t=768,12)'
'define w3=ave(vqs,t+0,t=768,12)'
'define m1=sum(w1,t=1,t=12)'
'define m2=sum(w2,t=1,t=12)'
'define m3=sum(w3,t=1,t=12)'

'set gxout shaded'
'd m1'
'set grads off'
'set gxout contour'
'd m1'
'set grads off'
'set gxout vector'
'd m2;m3'
'run cbarn'
'draw title 1948-2011 year average vopor flux china'
'printim g:\year.tif x1000 y800 white'
;



1111.jpg
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2012-8-21 09:34:34 | 显示全部楼层
hustlh 发表于 2012-8-20 16:10
惭愧啊,,我还有个问题想请教哈你,还望斑竹不吝赐教啊。

关于求64年(1948-2011)年平均水 ...

还是上次那样的设定时间问题,还有按照你的思路,相当于将768个月总的平均,你的是分两步求,第二部步骤也是求平均的,你只是求了个和
密码修改失败请联系微信:mofangbao
您需要登录后才可以回帖 登录 | 立即注册 新浪微博登陆

本版积分规则

Copyright ©2011-2014 bbs.06climate.com All Rights Reserved.  Powered by Discuz! (京ICP-10201084)

本站信息均由会员发表,不代表气象家园立场,禁止在本站发表与国家法律相抵触言论

快速回复 返回顶部 返回列表