var/convertor/kma_wave2grid/G2W.inc

   1       SUBROUTINE G2W
   2      I(MEND1,NEND1 ,JEND1,MNWAV,IMAX,JMAX  ,IMX,JMAXHF,KMAX ,
   3      I PNM  ,GDAT ,IFAX  ,TRIGS,GW  ,
   4      O QDAT ,
   5      W GWRK)
   6 C
   7       DIMENSION QDAT(2,KMAX,MNWAV)  ,PNM(MNWAV,JMAXHF)
   8       DIMENSION GDAT(IMAX,JMAX,KMAX),GWRK(IMX,JMAX,KMAX)
   9       DIMENSION IFAX(10),TRIGS(IMAX),GW(JMAX)
  10 C
  11       LOT =JMAX*KMAX
  12        print*,' rizvi in G2W input Gauss grid data '
  13        do j=100,JMAX,50
  14        print*,j,' Lat (Gauss) grid val ',(GDAT(i,j,KMAX),i=200,210)
  15        end do
  16       print*,' in G2W calling FFT991 with LOT= ',LOT
  17       CALL FFT991(GDAT,TRIGS,IFAX,1,IMX,IMAX,LOT,-1)
  18 C     CALL FFT991(GDAT,GWRK,TRIGS,IFAX,1,IMX,IMAX,LOT,-1)
  19 C
  20       CALL RESET(QDAT,2*KMAX*MNWAV)
  21 C
  22       CALL WEIHT2(IMAX,MEND1,JMAX/2,JMAX,KMAX,   0,   0,2*KMAX,
  23      1             GW,GDAT,GWRK)
  24
  25       CALL G2WPP (MEND1,NEND1,JEND1,1,2*KMAX,2*KMAX,MNWAV,
  26      1             JMAX/2,JMAX,PNM,QDAT,GWRK)
  27 C
  28       RETURN
  29       END SUBROUTINE G2W
  30