赤道环流剖面图时出错

海盗船长

学而不厌 发表于 2014-10-16 16:40
'reinit'
'sdfopen d:\omega.mon.mean.nc'
'sdfopen d:\uwnd.mon.mean.nc'

不知道哎

海盗船长

学而不厌 发表于 2014-10-16 16:40
'reinit'
'sdfopen d:\omega.mon.mean.nc'
'sdfopen d:\uwnd.mon.mean.nc'

你这是啥意思？貌似也不行哎。难道是数据的问题？

对雪吟歌

学而不厌 发表于 2014-10-16 16:40
'reinit'
'sdfopen d:\omega.mon.mean.nc'
'sdfopen d:\uwnd.mon.mean.nc'

图片太大了？

xiaocaoqiqiao

楼主是在做季风作业吗？资料的年份用1949-1998年的，你试试这个脚本，和课件上的图还挺接近的。
'reinit'
'sdfopen d:\omega.mon.mean.nc'
'sdfopen d:\uwnd.mon.mean.nc'
'enable print  d:\equator.gmf white'
'set lat 0'
'set lon 40 280'
'set lev 1000 100'
'set dfile 1'
'define dec1=ave(omega.1,t=24,t=612,12)'
'define jan1=ave(omega.1,t=25,t=613,12)'
'define feb1=ave(omega.1,t=26,t=614,12)'
'define w=(dec1+jan1+feb1)/3'   
'set dfile 2'
'define dec2=ave(uwnd.2,t=24,t=612,12)'
'define jan2=ave(uwnd.2,t=25,t=613,12)'
'define feb2=ave(uwnd.2,t=26,t=614,12)'
'define u=(dec2+jan2+feb2)/3'
'set grads off'
'set grid off'   
'set gxout shaded'
'set cmin 1'
*w与z符号相反*
'd -w*100'
'set gxout stream'
'set cthick 5'
'set arrowhead 0.09'
'cbar 1.0 0'
'd u;-w*100'
'draw title Zonal circulation along the equator in winter'
'print'
'disable print'
;

海盗船长

xiaocaoqiqiao 发表于 2014-10-18 16:55
楼主是在做季风作业吗？资料的年份用1949-1998年的，你试试这个脚本，和课件上的图还挺接近的。
'reinit'
...

你可以的~就想知道我的到底在哪里出错了。不过还是谢谢你了！

学而不厌

我那个程序可以的啊。你的程序只做了一年的月的平均，而且只有最后一年的三个月的平均保留了下来，所以我改了一下。



图和老师的还不一样应该是数据的问题，我的数据是2.5*2.5的

学而不厌

xiaocaoqiqiao 发表于 2014-10-18 16:55
楼主是在做季风作业吗？资料的年份用1949-1998年的，你试试这个脚本，和课件上的图还挺接近的。
'reinit'
...

想问一下你的数据是2.5*2.5的吗

xiaocaoqiqiao

学而不厌 发表于 2014-10-22 16:37
想问一下你的数据是2.5*2.5的吗

是的呢  在这儿下载的http://www.esrl.noaa.gov/psd/dat ... lysis.pressure.html

pqman3

1、while循环计算的是每年冬季平均值，由于缺少一个累加的语句，导致最后的结果仅是最后一年的平均值。
2、'define om=ave(omega,t='a1',t='a1'+2)'这个语句写错了。应该是'define om=ave(omega,t='a1',t='a1+2')'

海盗船长

pqman3 发表于 2014-10-23 11:22
1、while循环计算的是每年冬季平均值，由于缺少一个累加的语句，导致最后的结果仅是最后一年的平均值。
2 ...

昨天注意到了这个问题。谢谢你了！