求教用GrADS处理多个NC文件。

|▎花落°莫相离

我想要求CCN3的冬夏差异，其中6，7，8求平均作为夏季，12，1，2求平均作为冬季。改用循环写，结果说是无效的预算对象：
资料CTL：


我的GS：
'reinit'
'sdfopen F:\data\cam5_con\ccn_cam5_con_atm_06_1993-2007.nc'
'sdfopen F:\data\cam5_con\ccn_cam5_con_atm_07_1993-2007.nc'
'sdfopen F:\data\cam5_con\ccn_cam5_con_atm_08_1993-2007.nc'
'sdfopen F:\data\cam5_aero\ccn_cam5_aero_atm_06_1993-2007.nc'
'sdfopen F:\data\cam5_aero\ccn_cam5_aero_atm_07_1993-2007.nc'
'sdfopen F:\data\cam5_aero\ccn_cam5_aero_atm_08_1993-2007.nc'
'enable print F:\data\ccn_xia.gmf'
'set lat -90 90'
'set lon 0 360'
'set z 1'
aa1=0
aa2=0
aa3=0
aa4=0
aa5=0
aa6=0
it=1
         while(it<=15)
         'set t 'it''
      'define aa1=aa1+ccn3.1(t=it)'
      'define aa2=aa2+ccn3.2(t=it)'
      'define aa3=aa3+ccn3.3(t=it)'
      'define aa4=aa4+ccn3.4(t=it)'
      'define aa5=aa5+ccn3.5(t=it)'
      'define aa6=aa6+ccn3.6(t=it)'
        it=it+1
       endwhile
'define xia1=(aa1+aa2+aa3)/45'
'define xia2=(aa4+aa5+aa6)/45'
'define xia3=xia1-xia2'
'd xia3'
'draw title 'CCN cencentration at S=0.1％''
'print'
'disable print'
;

程序报错：

我不知道哪里出错了，求各位大神指点，轻拍。。。

传说中的谁

你家的纬度是从0-360的么？

mofangbao

传说中的谁 发表于 2013-4-24 23:08

                               
登录/注册后可看大图

你家的纬度是从0-360的么？

这个有点绕啊

river

你家的经度是从-90—90的么？

|▎花落°莫相离

传说中的谁 发表于 2013-4-24 23:08

                               
登录/注册后可看大图

你家的纬度是从0-360的么？

如果是画全球的，上面那样程序报错的原因是什么呀？求大神指点

|▎花落°莫相离

river 发表于 2013-4-25 09:05

                               
登录/注册后可看大图

你家的经度是从-90—90的么？

如果是画全球的，上面那样程序报错的原因是什么呀？求大神指点

|▎花落°莫相离

mofangbao 发表于 2013-4-24 23:17

                               
登录/注册后可看大图

这个有点绕啊

如果是画全球的，上面那样程序报错的原因是什么呀？求大神指点

river

|▎花落°莫相离 发表于 2013-4-26 17:48

                               
登录/注册后可看大图

如果是画全球的，上面那样程序报错的原因是什么呀？求大神指点

你先把经纬度设置对了，再考虑下面那个循环

小晶123

学习了。谢谢