G2W.inc
References to this file elsewhere.
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