【求助】grads 如何显示出计算均方根误差的最后结果

shangbn305

以下是我的grads程序，目的是计算d02-28800.ctl和sta.ctl的风速的均方根误差，sta.ctl是站点数据，因此先对sta.ctl做平均。现在的问题是，程序中可以计算出跟均方误差，但是不知道该如何显示出这个值，使用print但是不好用。变量cha是最后需要显示的值'reinit'
'open C:\GrADS20\win32\grid.ctl'
'open C:\GrADS20\win32\sta.ctl'
'open C:\GrADS20\win32\d02-28800.ctl'
'set t 1'
'define a=oacres(g,big.2)'
'define aa=aave(a,x=1,x=537,y=1,y=347)'
'set t 4'
'set z 1'
i=1
j=1
'define cha1=0'
while(i<538)
'set x ' i
while(j<348)
'set y 'j
'define quan=mag(u.3,v.3)'
'define cha1=(aa-quan)*(aa-quan)+cha1'
j=j+1
endwhile
i=i+1
endwhile
'define cha2=cha1/(537*347)'
'define cha=sqrt(cha2)';

river

算出来就是一个数字？那直接d cha就显示在屏幕上了

做个霸气的木头

你可以在运行完gs，在界面d cha，就能在下一行看到结果了。或者使用write函数写入到一个txt文件中

shangbn305

做个霸气的木头 发表于 2013-3-24 19:17

                               
登录/注册后可看大图

你可以在运行完gs，在界面d cha，就能在下一行看到结果了。或者使用write函数写入到一个txt文件中

您好，比如我想定义wrf中的变量，当u.3=1e+30时，令u.3=0，是用set还是define呢？如下：
if(u.3=1.e30)
'define u.3=0' 还是
'set u.3=0'
endif

做个霸气的木头

shangbn305 发表于 2013-3-24 19:41

                               
登录/注册后可看大图

您好，比如我想定义wrf中的变量，当u.3=1e+30时，令u.3=0，是用set还是define呢？如下：
if(u.3=1.e30)
...

define那种

shangbn305

做个霸气的木头 发表于 2013-3-24 21:06

                               
登录/注册后可看大图

define那种

非常感谢

3120101008

楼主，你好，我现在在做wrf与观测资料的对比，均方根误差，是根据你的脚本修改的，算出来的值7.22511e+07，感觉不对啊，请教一下，这是什么问题，正常吗？还是我的数据有问题？