2 ! { dg
-options
"-O3 -fwhole-file" }
3 CALL UVSET
(NX
,NY
,NZ
,HVAR
,ZET
,NP
,DZ
,DKM
,UM
,VM
,UG
,VG
,TM
,DCDX
, ! { dg
-warning
"Rank mismatch" }
5 CALL DCTDX
(NX
,NY
,NX1
,NFILT
,C
(MLAG
),DCDX
(MLAG
),HELP
,HELPA
,
6 * HELP
,HELPA
,FY
,FYC
,SAVEY
)
8 SUBROUTINE PADEC
(DKS
,DKDS
,HVAR
,WM
,WG
,FN
,NS
,AN
,BN
,CN
,IT
)
9 COMPLEX*16 WM
(*),WG
(*),FN
(*),AN
(*),BN
(*),CN
(*)
12 FN
(J
)=(AS
+F4
-GAMMA*F2
)*H2
+(F4
-AS
-GAMMA*F2
)*H0
+
13 * H1*
(F3
-GAMMA
/3.D0
)+GAMMA*WG
(J
)-CONST
15 SUBROUTINE UVSET
(NX
,NY
,NZ
,HVAR
,ZET
,NP
,DZ
,DKM
,UM
,VM
,UG
,VG
,TM
,
17 DIMENSION HVAR
(*),ZET
(*),TM
(*),DKM
(*),UM
(*),VM
(*),UG
(*),VG
(*),
20 CALL PADEC
(DKM
,VM
,HVAR
,WORK
(LWM
),WORK
(LWG
), ! { dg
-warning
"Rank mismatch" }
21 * WORK
(LF
),NZ
,WORK
(LA
),WORK
(LB
),WORK
(LC
),ITY
)