W2GPXY.inc

1       SUBROUTINE W2GPXY
2      I(MEND1,NEND1,JEND1,MNWAV,IMAX,JMAX,IMX ,JMAXHF,
3      I IFAX ,TRIGS,PNM  ,DPNM ,IPXY,IPSL,QPAI,
4      O GDAT ,
5      W GWRK)
6 C
7       DIMENSION QPAI(2,MNWAV)
8 
9 c990408
10 c      DIMENSION IFAX(10),TRIGS(500)
11       DIMENSION IFAX(10),TRIGS(IMAX)
12 C...WARNING: THE DIMENSION OF GDAT,GWRK SHOULD BE .GE. IMX*JMAX*3
13 C            IMX=IMAX+2
14       DIMENSION GDAT(IMAX,JMAX,3) ,GWRK(IMAX,JMAX,3)
15       DIMENSION PNM (MNWAV,JMAXHF),DPNM(MNWAV,JMAXHF)
16 C...IPXY=1; OUTPUT OF THIS SUBR. IS  PS,D(PS)/DX,D(PS)/DY
17 C...     0;                          PS
18 C
19       CALL LGNW2G
20      I(MEND1,NEND1,JEND1,MNWAV,IMAX,JMAX,JMAXHF,   1,   1, PNM,QPAI,
21      O GDAT ,
22      W GWRK)
23       LOT=JMAX
24       IF(IPXY.EQ.1) THEN
25       CALL LGNW2G
26      I(MEND1,NEND1,JEND1,MNWAV,IMAX,JMAX,JMAXHF,   1,   0,DPNM,QPAI,
27      O GDAT(1,1,3),
28      W GWRK)
29       DO 100 J=1,JMAX
30       DO 100 M=1,MEND1
31       GDAT(2*M-1,J,2)=-FLOAT(M-1)*GDAT(2*M  ,J,1)
32       GDAT(2*M  ,J,2)=+FLOAT(M-1)*GDAT(2*M-1,J,1)
33   100 CONTINUE
34       DO 120 J=1,JMAX
35       DO 120 I=2*MEND1+1,IMAX
36       GDAT(I,J,2)=0.0
37       GDAT(I,J,3)=0.0
38   120 CONTINUE
39       LOT=JMAX*3
40       END IF
41 C
42       CALL FFT991(GDAT,TRIGS,IFAX,1,IMX,IMAX,LOT,1)
43 C     CALL FFT991(GDAT,GWRK,TRIGS,IFAX,1,IMX,IMAX,LOT,1)
44 C
45       IF(IPSL.EQ.1) THEN
46       DO 200 J=1,JMAX
47       DO 200 I=1,IMAX
48       X          =        EXP(GDAT(I,J,1))
49       GDAT(I,J,1)=X
50       GDAT(I,J,2)=GDAT(I,J,2)*X
51       GDAT(I,J,3)=GDAT(I,J,3)*X
52   200 CONTINUE
53       END IF
54 C
55       RETURN
56       END SUBROUTINE W2GPXY
57