fortran求指教

╯①路霺笑

cccccccccc  prepare data for svd analysis ccccccccc
cccccccccc   the program is to normalize sea surface temperature(SST) cccccccccc
cccccccccc   mt: the length of time series;
cccccccccc   mo: the month numbers;           my: the year numbers;
cccccccccc  sst: sea surface temperature data;
cccccccccc sst3: the work array;             avf: the average of SST;
cccccccccc   df: the variance of SST;            
      program main
      parameter(mo=12,my=43,nx=18,ny=12,mt=516)
      dimension sst2(nx,ny,mo,my),avf(mo,nx,ny),df(mo,nx,ny)
      dimension sst(nx,ny,mt),sst3(nx,ny,mt),II(NX,NY),
     & jp(nx,ny,mo,my)
      open(1,file='g:sstpx.grd',form='unformatted',
     &     access='direct',recl=nx*ny)
            irec=1
       do 5 it=1,mt
       read(1,rec=irec)((sst(i,j,it),i=1,nx),j=1,ny)
       irec=irec+1
5      continue
       do i=1,nx
            it1=i
         do j=1,ny
         do it=1,mt
         do m=1,mo
         do myear=1,my
            sst2(i,j,m,myear)=sst(i,j,it1)
            it1=it1+12
         end do
         end do
         end do
         end do
         end do
       do i=1,nx
       do j=1,ny
       do m=1,12
       avf(i,j,m)=0.0
       do myear=1,my
       avf(i,j,m)=avf(i,j,m)+sst2(i,j,m,myear)/my   !求平均值
       end do
       end do
       end do
       end do
         do i=1,nx
       do j=1,ny
       do m=1,mo
       do myear=1,my
       jp(i,j,m,myear)=sst2(i,j,m,myear)-avf(i,j,m)  !求距平
       end do
       end do
       end do
       end do
       open(2,file='g:\jp.grd',form='binary')
       do i=1,nx
       do j=1,ny
       do myear=1,my
       do m=1,mo
       kk=(myear-1)*12+m
       sst3(i,j,kk)=jp(i,j,m,myear)
       end do
       end do
       end do
       end do
       do it=1,mt
       write(2)((sst3(i,j,it),i=1,nx),j=1,ny)
       end do         



ccccccccccccccccccccccccccccccccccccccccccc
      close (1)
      close (2)
      stop
      end


这是我编的一段程序，需要将这段程序和另一段程序结合，构成eof分析的程序，没错，但是在电脑运行不出来，想请大神帮忙看看

no_trust

全忘了，不过两个程序和起来不是有，那个什么子程序吗，调用那类的function

黎大页子

把运行错误提示信息也贴出来

飞幻天狼

把运行错误提示信息贴出来，这样知道是哪里出了问题

guizi

ncl  自带 eof  ......  现成的，就两句话，