青高感热距平的程序，图画的不对，可能程序错了帮忙看看，我想尽快毕业！

bellatf

这是1979-2012年的求青高夏季感热距平的，路径错误请忽略……


program main
implicit none
integer i,j,m,y
integer,parameter::yr=34,mo=12,nx=18,ny=7
real::shtfl(ny,nx,mo,yr),shtfl1(ny,nx,yr),shtfl3(yr),shtfl2(ny,nx,yr),shtfl4(ny,nx,yr)
open(1,file='d:\data\sh-19792012qinggao.grd',form='binary')
do y=1,yr
do m=1,mo
do i=1,nx
do j=1,ny
   read(1) shtfl(j,i,m,y)
   enddo
  enddo
enddo
enddo
print *,'1'
!!!!!!!!!!!!!!!!!!以上是读入数据!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!


open(3,file='d:\data\sh-19792011qinggaosum.txt')!以下输出结果!!!!!
do i=1,nx
do j=1,ny
  !shtfl3(j,i)=0.
  do y=1,yr
   shtfl2(j,i,y)=0.
   do m=3,5
   shtfl2(j,i,y)=shtfl2(j,i,y)+shtfl(j,i,m,y)/3
   enddo
   !shtfl1(j,i,y)=(shtfl2(j,i,y)+shtfl(j,i,12,y))/3
   !shtfl3(j,i)=shtfl3(j,i)+shtfl1(j,i,y)/33
  enddo
enddo
enddo






do y=1,yr
shtfl3(y)=0.
do i=1,nx
   do j=1,ny
    shtfl3(y)=shtfl3(y)+shtfl2(j,i,y)/(nx*ny)
        !write(3,*) shtfl4(j,i,y)
   enddo
  enddo
  write(3,*) shtfl3(y)
enddo


!shtfl2(j,i)=0.
!do y=1,yr-1
  !shtfl2(j,i)=shtfl2(j,i)+shtfl(j,i,12,y)/33
  !enddo
  !enddo
!enddo
  
!do i=1,nx
!do j=1,ny
  
  !shtfl3(j,i)=shtfl1(j,i)+shtfl2(j,i)
  !write(3) shtfl3(j,i)
  !enddo
!enddo
!enddo


end

bellatf

请大神们花一分钟时间帮我看看吧

bellatf

bellatf 发表于 2015-4-29 11:00
请大神们花一分钟时间帮我看看吧

grd  应该没错，  青藏高原的范围北纬23度-北纬39度 东经73度-104度 这个应该不会出错吧……

hare

没错啊，我下了

bellatf

hare 发表于 2015-4-29 11:36
没错啊，我下了

恩，没错，可是画出来的图 不对…… 没有异常的