求助！处理并输出固定月份nc数据

夏无且

选取的数据是1948年-2010年的u月均值，想求1950-2010年夏季风指数，即6，7，8月份的两个区域平均值相减的结果，刚开始学看了相关一些gs改了几遍还是哪有问题，跪求大家指教！
这是我写的gs


运行结果


来自群组: 成于大气 信达天下

夏无且

经各位大大指正后的运行成功的脚本，供大家参考嘿嘿

f117hqh

目前看到几个问题：
第一，既然求面积平均，就固定x和y，即set x 1;set y 1
第二，set t那一行有语法错误，应该是'set t 'ya''

夏无且

f117hqh 发表于 2016-12-16 13:46
目前看到几个问题：
第一，既然求面积平均，就固定x和y，即set x 1;set y 1
第二，set t那一行有语法错误 ...

谢谢指教，但还是有问题，还有看有些帖子也没固定xy，好像也能输出样。。。

liyf

j<30
while(j<751)
i=0
while(i<3)
ya=j+i
'set t 'ya
'set z 1'
'set lon 100 150'
'set lat  10 20'
'define u1=aave(uwnd,lon=100,lon=150,lat=10,lat=20)'
'define u2=aave(uwnd,lon=100,lon=150,lat=10,lat=20)'.....我不知道你的另一个区域?
'set lon 100'
'set lat 20'
'set z 1'
'd u1-u2'
i=i+1
endwhile
j=j+12
endwhile
......因为我觉得你说的意思有点模糊，所以仅仅是按照你的程序改了一下
。。这样得到的结果，无法写对应ctl.
因此，如果你想要得到三个月的数据，最好分开求，每一个月dat写一个ctl，如果你想要求678平均，那就好办，再求一个平均就好

夏无且

liyf 发表于 2016-12-16 19:05
j

新改过后这样（抱歉没发现没截全图）
'define u1=aave(uwnd,lon=100,lon=150,lat=10,lat=20)'
'define u2=aave(uwnd,lon=100,lon=150,lat=25,lon=35)'
'set lon 100'
'set lat 20'
'set z 3'
'd u1-u2'
之前也试过这样分开算平均，再d，就会出现

因为想提取出每年6，7，8在dat里后，就直接用Fortran读，算滑动平均及和另一量的相关系数，就可以不分开每月写了吧。
想另问下，define u1，u2，后面set x，y的值是在区域内取得任意值么？谢谢谢谢！

liyf

你写错了'define u2=aave(uwnd,lon=100,lon=150,lat=25,lon=35)..lat=35
...如果你仅仅是用来计算，怎么存都没事。
define u1，u2，后面set x，y的值是在区域内取得任意值

夏无且

liyf 发表于 2016-12-19 09:14
你写错了'define u2=aave(uwnd,lon=100,lon=150,lat=25,lon=35)..lat=35
...如果你仅仅是用来计算，怎么存 ...

哈哈哈，蜜汁尴尬。。。谢谢！！！