有劳各位帮忙看下我这个图 ，感觉是错的就是不知道错在哪~~

漂洋过海

画的是垂直环流距平场，是五个干旱年合成的。
这是提取资料的gs如下：
'reinit'
'set fwrite H:\drought\analysis\omega\w-dry-years.grd'
'set gxout fwrite'

'sdfopen H:\drought\analysis\omega\data\omega.1965.nc'
'set y 1 73'
'set z 1 12'
tt=152;while(tt<=243)
'set t 'tt
'set x 1'
'define o1=ave(omega,lon=101.5,lon=112.5)'
'd o1'
tt=tt+1;endwhile
'close 1'

'sdfopen H:\drought\analysis\omega\data\omega.1982.nc'
'set y 1 73'
'set z 1 12'
tt=152;while(tt<=243)
'set t 'tt
'set x 1'
'define o2=ave(omega,lon=101.5,lon=112.5)'
'd o2'
tt=tt+1;endwhile
'close 1'

'sdfopen H:\drought\analysis\omega\data\omega.1997.nc'
'set y 1 73'
'set z 1 12'
tt=152;while(tt<=243)
'set t 'tt
'set x 1'
'define o3=ave(omega,lon=101.5,lon=112.5)'
'd o3'
tt=tt+1;endwhile
'close 1'

'disable fwrite'
;


之后我用fortran做了平均

program average
!x=x方向格点数，z=层次，y=年份，d=一年的47天
parameter(x=73,z=12,d=92,y=5)
real o1(x,z,d,y)

real o31(x,z,d),o32(x,z)
character fl1*60,fl2*60,fl3*60
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
fl1="H:\drought\analysis\omega\w-dry-years.grd"

fl3="H:\drought\analysis\omega\w-dry-years-hecheng.grd"
open(1,file=fl1,status='old',form='binary')

open(3,file=fl3,status='new',form='binary')


do i=1,y
   do j=1,d
      do k=1,z
             do m=1,x
                    read(1) o1(m,k,j,i)
         end do
          end do
   end do
end do
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
do i=1,x
   do j=1,z
      do k=1,d
             o31(i,j,k)=0.0
          end do
   end do
end do


do i=1,x
   do j=1,z
      do k=1,d
             do m=1,y
                    o31(i,j,k)=o31(i,j,k)+o1(i,j,k,m)
             end do
          end do
   end do
end do

do i=1,x
   do j=1,z
      do k=1,d
             o31(i,j,k)=o31(i,j,k)/real(y)
          end do
   end do
end do

!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!

do i=1,x
   do j=1,z
      o32(i,j)=0.0
   end do
end do

do i=1,x
   do j=1,z
      do k=1,d
         o32(i,j)=o32(i,j)+o31(i,j,k)
          end do
   end do
end do

do i=1,x
   do j=1,z
      o32(i,j)=o32(i,j)/real(d)
   end do
end do
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!

do i=1,z
   do j=1,x
       write(3) o32(j,i)
           print*, o32(j,i)
   end do

end do


close(1)
close(2)
end

这是画图的gs,
'reinit'
'open H:\drought\analysis\omega\v-dry.ctl'
'open H:\drought\analysis\omega\w-dry.ctl'
'open H:\drought\analysis\omega\v-30y.ctl'
'open H:\drought\analysis\omega\w-30y.ctl'


'enable print H:\drought\analysis\omega\anomaly.gmf'
**********************************************************************
'set parea 2 9 2 7'
'set grads off'
'set grid off'
'set zlog on'
'set xlopts 1 8 0.2'
'set ylopts 1 8 0.2'

'set xlint 5'
'set lat 20 60';
'set x 1';
'set z 1 12';
'set t 1'

'define v=vwnd.1-vwnd.3'
'define o1=omega.2-omega.4'
'define o=-o1*100'

'set gxout shaded'
'set clevs  -5 -3  -2  -1 0 '

'set ccols  8 12 7 10 13  0'
'd o'

'cbarn 0.8 0 5.5 1.3'



'set csmooth on'
'set gxout stream'
'set ccolor 1'
'set strmden 4'
'd v;o'

'draw title 101.5-112.5E dry-years-ave'


'print'
'disable print'
'printim H:\drought\analysis\omega\anomaly.png  white'
;

请教各位，为什么流场那块断线了~~~，十分捉急啊现在这里谢过各位看帖的人啦             

river

楼主这个太麻烦了，看的头疼了

漂洋过海

river 发表于 2013-6-27 11:43

                               
登录/注册后可看大图

楼主这个太麻烦了，看的头疼了

我好像就会这种思路，哎

river

漂洋过海 发表于 2013-6-27 13:13

                               
登录/注册后可看大图

我好像就会这种思路，哎

既然是五个干旱年合成，那为什么提取资料只提取了三年的呢？

whmooncool

会不会是缺测值的问题呢？

漂洋过海

river 发表于 2013-6-27 15:48

                               
登录/注册后可看大图

既然是五个干旱年合成，那为什么提取资料只提取了三年的呢？

我就随手复制了几行····
是不是因为是径向平均距平场，所以断线了啊

漂洋过海

whmooncool 发表于 2013-6-27 16:07

                               
登录/注册后可看大图

会不会是缺测值的问题呢？

也有可能吧~

river

漂洋过海 发表于 2013-6-27 18:23

                               
登录/注册后可看大图

我就随手复制了几行····
是不是因为是径向平均距平场，所以断线了啊

应该不会吧，没碰到过

漂洋过海

river 发表于 2013-6-27 20:31

                               
登录/注册后可看大图

应该不会吧，没碰到过

哎明天了去问问老师吧，谢谢你啊~~

river

漂洋过海 发表于 2013-6-27 20:52

                               
登录/注册后可看大图

哎明天了去问问老师吧，谢谢你啊~~

有了结果回来说说