三个chl8天合成文件，grads求平均求教

血嗜子木

附件三个文件chl8天合成文件，我想利用grads做三个文件chl的融合，及是三个文件数据平均画图，但是利用'd (chlor_a.1+chlor_a.2+chlor_a.3)/3'，最后总是显示和缺失数据最多的那个图类似的图，一直都不知道怎么解决，想请问下各位大神，怎么做呢？前三个是每个图的，后面的是利用'd (chlor_a.1+chlor_a.2+chlor_a.3)/3'画出的，哪位老师可以告诉我下怎么做呢？？万分感谢！
下面是代码：
'reinit'
'sdfopen D:\3\5.nc'
'sdfopen D:\3\6.nc'
'sdfopen D:\3\7.nc'
'set lon 110 122'
'set lat 19 33'
'set xaxis 110 122 1'
'set yaxis 19 33 1'

*左上角第一张
'set font -1'
'set xlopts 1 4 0.12'
'set ylopts 1 4 0.12'
'set parea 1 5.5 4.5 8'
'set grads on'
'set grid on'
'set mpdset hires'
'set gxout shaded'
'set cint 5'
'set clevs 0.08 0.1 0.2 0.3 0.4 0.5 0.6 0.7 1 2.5 10 30'
'set xlpos -20'                                                                           
'set ylpos -20'
'd chlor_a.1'
'run axis.gs -type b  -position m -label on -interval 2 -langle 45'
'run axis.gs -type L  -position o -start 19 -end 33 -interval 2 -suffix `3.'
*右上角第一张
'set parea 6 8.5 4.5 8'
'set grads on'
'set grid on'
'set mpdset hires'
'set gxout shaded'
'set cint 5'
'set clevs 0.08 0.1 0.2 0.3 0.4 0.5 0.6 0.7 1 2.5 10 30'
'set xlpos -20'                                                                           
'set ylpos -20'
'd chlor_a.2'
'run axis.gs -type b  -position m -label on -interval 2 -langle 45'
'run axis.gs -type L  -position o -start 19 -end 33 -interval 2 -suffix `3.'
'cbarn 0.5 1'
*左下角第一张
'set parea 1 5.5 0.5 4'
'set grads on'
'set grid on'
'set mpdset hires'
'set gxout shaded'
'set cint 5'
'set clevs 0.08 0.1 0.2 0.3 0.4 0.5 0.6 0.7 1 2.5 10 30'
'set xlpos -20'                                                                           
'set ylpos -20'
'd chlor_a.3'
'run axis.gs -type b  -position m -label on -interval 2 -langle 45'
'run axis.gs -type L  -position o -start 19 -end 33 -interval 2 -suffix `3.'
*右下角第一张
'set parea 6 8.5 0.5 4'
'set grads on'
'set grid on'
'set mpdset hires'
'set gxout shaded'
'set cint 5'
'set clevs 0.08 0.1 0.2 0.3 0.4 0.5 0.6 0.7 1 2.5 10 30'
'set xlpos -20'                                                                           
'set ylpos -20'
'd (chlor_a.1+chlor_a.2+chlor_a.3)/3'
'run axis.gs -type b  -position m -label on -interval 2 -langle 45'
'run axis.gs -type L  -position o -start 19 -end 33 -interval 2 -suffix `3.'
*画图
'printim d:\3\chl44.png png white'
;

river

求平均那个你把set clevs 先注释掉看看是什么样

血嗜子木

river 发表于 2015-7-18 21:44
求平均那个你把set clevs 先注释掉看看是什么样

D:\3\chl44

river

血嗜子木 发表于 2015-7-18 22:01

脚本问题不大，我建议你把 cint 改小一些试试

血嗜子木

river 发表于 2015-7-19 06:17
脚本问题不大，我建议你把 cint 改小一些试试

grads在使用这种平均的时会出现缺测数据导致平均图缺测，请问您有什么好的方法解决吗？？万分感谢！

river

血嗜子木 发表于 2015-7-19 08:09
grads在使用这种平均的时会出现缺测数据导致平均图缺测，请问您有什么好的方法解决吗？？万分感谢！

我没懂你的意思，资料里原本的缺测值在求平均的时候是不参加运算的······

血嗜子木

river 发表于 2015-7-19 15:02
我没懂你的意思，资料里原本的缺测值在求平均的时候是不参加运算的······

图中的数据平均最后一个图和第一个比较接近，因为第一个图却测数据过多，平均时影响了最后的结果

MeteoInfo

仔细看看这个帖子你更明白这个问题：http://bbs.06climate.com/forum.php?mod=viewthread&tid=10604

但是GrADS里可能没有相应的解决方案。

血嗜子木

MeteoInfo 发表于 2015-7-19 17:28
仔细看看这个帖子你更明白这个问题：http://bbs.06climate.com/forum.php?mod=viewthread&tid=10604

但 ...

谢谢您，好像是的呢，找了很久，没有适合的方式，您这是另外一个软件呀~~还没学呢

naoki1122

刷金币，下资料！！！