2 I(MEND1,NEND1 ,JEND1,MNWAV,IMAX,JMAX ,IMX ,JMAXHF,KMAX ,
3 I PNMGC,DPNMGC,GU ,GV ,ER ,SINCLT,IFAX,TRIGS,
7 DIMENSION QROT (2,KMAX,MNWAV),QDIV (2,KMAX,MNWAV)
8 DIMENSION DPNMGC(MNWAV,JMAXHF),PNMGC (MNWAV,JMAXHF),SINCLT(JMAX)
9 DIMENSION GU (IMAX,JMAX,KMAX),GV (IMAX,JMAX,KMAX)
10 DIMENSION GWRK(IMX,JMAX,KMAX),IFAX(10),TRIGS(IMAX)
14 CALL FFT991(GU,TRIGS,IFAX,1,IMX,IMAX,LOT,-1)
15 CALL FFT991(GV,TRIGS,IFAX,1,IMX,IMAX,LOT,-1)
16 C CALL FFT991(GU ,GWRK,TRIGS,IFAX,1,IMX,IMAX,LOT,-1)
17 C CALL FFT991(GV ,GWRK,TRIGS,IFAX,1,IMX,IMAX,LOT,-1)
19 CALL RESET(QROT,2*KMAX*MNWAV)
20 CALL RESET(QDIV,2*KMAX*MNWAV)
28 IF(MOD(JMAX,2).EQ.1.AND.J.EQ.JMAXHF) THEN
30 GWRK(M,JM,1)=ASINCL*(GU(M,JM,K)+GU(M,JP,K))
31 GWRK(M,JM,2)=ASINCL*(GV(M,JM,K)+GV(M,JP,K))
35 GWRK(M,JM,1)=ASINCL*(GU(M,JM,K)+GU(M,JP,K))
36 GWRK(M,JP,1)=ASINCL*(GU(M,JM,K)-GU(M,JP,K))
37 GWRK(M,JM,2)=ASINCL*(GV(M,JM,K)+GV(M,JP,K))
38 GWRK(M,JP,2)=ASINCL*(GV(M,JM,K)-GV(M,JP,K))
46 NMAX=MIN(JEND1+1-M,NEND1)
48 IF(MOD(N-1,2).EQ.0) THEN
52 QROT(1,K,L+N)=QROT(1,K,L+N)-QM*GWRK(2*M ,JM,2)* PNMGC(L+N,J)
53 1 - GWRK(2*M-1,JP,1)*DPNMGC(L+N,J)
54 QROT(2,K,L+N)=QROT(2,K,L+N)+QM*GWRK(2*M-1,JM,2)* PNMGC(L+N,J)
55 1 - GWRK(2*M ,JP,1)*DPNMGC(L+N,J)
56 QDIV(1,K,L+N)=QDIV(1,K,L+N)-QM*GWRK(2*M ,JM,1)* PNMGC(L+N,J)
57 1 + GWRK(2*M-1,JP,2)*DPNMGC(L+N,J)
58 QDIV(2,K,L+N)=QDIV(2,K,L+N)+QM*GWRK(2*M-1,JM,1)* PNMGC(L+N,J)
59 1 + GWRK(2*M ,JP,2)*DPNMGC(L+N,J)
65 QROT(1,K,L+N)=QROT(1,K,L+N)-QM*GWRK(2*M ,JP,2)* PNMGC(L+N,J)
66 1 - GWRK(2*M-1,JM,1)*DPNMGC(L+N,J)
67 QROT(2,K,L+N)=QROT(2,K,L+N)+QM*GWRK(2*M-1,JP,2)* PNMGC(L+N,J)
68 1 - GWRK(2*M ,JM,1)*DPNMGC(L+N,J)
69 QDIV(1,K,L+N)=QDIV(1,K,L+N)-QM*GWRK(2*M ,JP,1)* PNMGC(L+N,J)
70 1 + GWRK(2*M-1,JM,2)*DPNMGC(L+N,J)
71 QDIV(2,K,L+N)=QDIV(2,K,L+N)+QM*GWRK(2*M-1,JP,1)* PNMGC(L+N,J)
72 1 + GWRK(2*M ,JM,2)*DPNMGC(L+N,J)