对OLR的原始逐日资料进行10-90天的滤波

stefan_scofield

我要对OLR的原始逐日资料进行10-90天的滤波，向老师要了一份滤波的fortran程序，我琢磨了很久实在不懂。我就直接修改了部分主程序运行了，可是输出的结果全是NaN，我找不出哪里错了。希望会的老师同学帮我看看程序哪里错了，或是能够提供滤波的原理也行（我在书里没有找到相关内容），我定当感激不尽，人在南信大我请你吃饭。
下面是主程序（需要看全部程序的，去往：http://bbs.06climate.com/forum.p ... mp;page=1#pid266148）：
program mian
      parameter(nx=72,ny=32,nt=5202,freqL=1/90.,freqU=1/10.)

      dimension olr(nx,ny,nt)
        dimension r2(nx,ny,nt)
        dimension mid(nx,ny,nt)
      dimension olrf(nx,ny,nt)
!      complex r1(nt),work(nt)
      real r1(nt),work(nt)

      OPEN(11,FILE='olr.79-12.6-10.grd',FORM='BINARY')
      PRINT*,'NOW READING ...'
      DO IT=1,NT
      DO J=1,NY
      DO I=1,NX
        READ(11) olr(I,J,IT)
      ENDDO
      ENDDO
      enddo
      close(11)


        do it=1,nt
        do j=1,ny
        do i=1,nx
                mid(i,j,it)=olr(i,j,it)
        enddo
        enddo
        enddo

      do 111 i=1,nx
      do 111 j=1,ny
        do 50 k=1,nt
         r1(k)=mid(i,j,k)
   50   continue
        call bandps(nt,r1,1.,freqL,freqU,0,work)
        do 60 k=1,nt
        r2(i,j,k)=real(r1(k))
   60   continue
  111   continue
      do k=1,nt
      do j=1,ny
      do i=1,nx
         olrf(i,j,k)=r2(i,j,k)
      enddo
      enddo
      enddo



      OPEN(12,FILE='bp-filtered-olr.grd',FORM='BINARY')
      PRINT*,'NOW WRITING ...'
      do k=1,nt
      do j=1,ny
      do i=1,nx
        if(abs(olrf(i,j,k))>5000.0) then
         olrf(i,j,k)=-9999
           write(12) olrf(i,j,k)
          else
           write(12) olrf(i,j,k)
          endif
      enddo
      enddo
      enddo
      close(12)
      print*,'NOW COMPLETED!'

      stop
      end

xuexiaotai

楼主你好！我也需要OLR逐日资料但是没找到！能发给我一份么？497710575@qq.com  谢谢了

xuexiaotai

楼主你好！我也需要OLR逐日资料但是没找到！能发给我一份么？497710575@qq.com  谢谢了

xiaocaoqiqiao

xuexiaotai 发表于 2014-4-13 19:13
楼主你好！我也需要OLR逐日资料但是没找到！能发给我一份么？  谢谢了

http://www.esrl.noaa.gov/psd/data/gridded/data.interp_OLR.html
这里有，olr.day.mean.nc就是逐日的