1C====>Begin Module PROB5_4DIM             File t208x.f
2C---->Makedcls Options: All variables
3C
4C     Parameter variables
5C
6      INTEGER L2NG
7      PARAMETER (L2NG = 6)
8      INTEGER NG
9      PARAMETER (NG = 2**L2NG)
10      INTEGER NB
11      PARAMETER (NB = NG+2)
12      INTEGER NGM1
13      PARAMETER (NGM1 = NG-1)
14      INTEGER NGP1
15      PARAMETER (NGP1 = NG+1)
16      INTEGER NGP2
17      PARAMETER (NGP2 = NG+2)
18      INTEGER NBP1
19      PARAMETER (NBP1 = NB+1)
20      INTEGER NR
21      PARAMETER (NR = 32)
22      INTEGER NT
23      PARAMETER (NT = 16)
24      INTEGER NR2
25      PARAMETER (NR2 = NR+2)
26      INTEGER NT1
27      PARAMETER (NT1 = NT+1)
28      INTEGER NW
29      PARAMETER (NW = 2*NT+15)
30      INTEGER MMAX
31      PARAMETER (MMAX = 200)
32C
33C     Local variables
34C
35      CHARACTER*18 CFILE, E10FILE, E11FILE, E12FILE, E13FILE, E14FILE
36      CHARACTER*18 E15FILE, E16FILE, E17FILE, E18FILE, E19FILE, E1FILE
37      CHARACTER*18 E20FILE, E21FILE, E22FILE, E23FILE, E24FILE, E2FILE
38      CHARACTER*18 E3FILE, E4FILE, E5FILE, E6FILE, E7FILE, E8FILE
39      CHARACTER*18 E9FILE, EFILE, FFILE, FMFILE, MAIFILE, MFILE, MMFILE
40      CHARACTER*18 PAFILE, PARFILE, PNFILE
41      CHARACTER*15 RNAME
42      CHARACTER*18 SYMFILE, UFILE, ZFILE
43C
44      DOUBLE PRECISION CMAX, CMIN, F0, FBMAX, FLMAX, PHIAMAX, PHIAMIN
45      DOUBLE PRECISION PHINMAX, PHINMIN, Q, TCOEF, TFREQ, TIME, TPI
46      DOUBLE PRECISION UMAX, UOLD(0:NBP1,0:NBP1,2), USMAX
47      DOUBLE PRECISION ZAGI(0:NBP1,0:NBP1), ZAGMAX, ZAGMIN
48C
49      INTEGER RLEN
50C
51      LOGICAL AVRAGE(06), AXES3D, BCLIP, BOX, BYROWS, CHKDUP, DERIV(06)
52C
53      REAL DUSUB, DVSUB, DXSUB, DYSUB, DZSUB
54C
55C     Namelists
56C
57      NAMELIST /DATA/ AVRAGE, AXES3D, BCLIP, BOX, BYROWS, CHKDUP, DERIV
58     x, DUSUB, DVSUB, DXSUB, DYSUB, DZSUB
59C
60C     Common variables
61C
62      DOUBLE PRECISION C(0:NBP1,0:NBP1), CN(0:NBP1,0:NBP1)
63C
64      COMMON /ADP/ C, CN
65C
66C     Common variables
67C
68      DOUBLE PRECISION ZAG(0:NBP1,0:NBP1), ZAGN(0:NBP1,0:NBP1)
69C
70      COMMON /AGG/ ZAG, ZAGN
71C
72C     Common variables
73C
74      INTEGER NFIL
75C
76      COMMON /EFNUM/ NFIL
77C
78C     Common variables
79C
80      DOUBLE PRECISION XM(MMAX,2)
81C
82      COMMON /FMARKERS/ XM
83C
84C     Common variables
85C
86      DOUBLE PRECISION F(0:NBP1,0:NBP1,2)
87C
88      COMMON /FORCE/ F
89C
90C     Common variables
91C
92      DOUBLE PRECISION FB(0:NBP1,0:NBP1,2)
93C
94      COMMON /FORCEB/ FB
95C
96C     Common variables
97C
98      DOUBLE PRECISION ELF(0:NBP1,0:NBP1,0:NR2,NT1)
99      DOUBLE PRECISION ELFN(0:NBP1,0:NBP1,0:NR2,NT1)
100C
101      COMMON /LINK/ ELF, ELFN
102C
103C     Common variables
104C
105      INTEGER METHOD, MTHLIM
106C
107      COMMON /MTH/ METHOD, MTHLIM
108C
109C     Common variables
110C
111      DOUBLE PRECISION PHIA(0:NBP1,0:NBP1), PHIAN(0:NBP1,0:NBP1)
112      DOUBLE PRECISION PHIN(0:NBP1,0:NBP1), PHINN(0:NBP1,0:NBP1)
113C
114      COMMON /PHI/ PHIN, PHINN, PHIA, PHIAN
115C
116C     Common variables
117C
118      DOUBLE PRECISION DR(NR2), R(0:NR2)
119C
120      INTEGER NK
121C
122      COMMON /RSIZE/ DR, R, NK
123C
124C     Common variables
125C
126      DOUBLE PRECISION DLT, H
127C
128      COMMON /STEPS/ H, DLT
129C
130C     Common variables
131C
132      DOUBLE PRECISION K0
133C
134      COMMON /STIFF/ K0
135C
136C     Common variables
137C
138      DOUBLE PRECISION U(0:NBP1,0:NBP1,2), UX(0:NBP1,0:NBP1)
139      DOUBLE PRECISION UY(0:NBP1,0:NBP1), VX(0:NBP1,0:NBP1)
140      DOUBLE PRECISION VY(0:NBP1,0:NBP1)
141C
142      COMMON /VEL/ U, UY, VY, UX, VX
143C
144C     Common variables
145C
146      DOUBLE PRECISION WSAVE(NW)
147C
148      COMMON /WAVE/ WSAVE
149C
150C====>End Module   PROB5_4DIM             File t208x.f
151C====>Begin Module UINIT                  File t208x.f
152C---->Makedcls Options: All variables
153C
154C     Parameter variables
155C
156      INTEGER L2NG
157      PARAMETER (L2NG = 6)
158      INTEGER NG
159      PARAMETER (NG = 2**L2NG)
160      INTEGER NB
161      PARAMETER (NB = NG+2)
162      INTEGER NGM1
163      PARAMETER (NGM1 = NG-1)
164      INTEGER NGP1
165      PARAMETER (NGP1 = NG+1)
166      INTEGER NGP2
167      PARAMETER (NGP2 = NG+2)
168      INTEGER NBP1
169      PARAMETER (NBP1 = NB+1)
170      INTEGER NR
171      PARAMETER (NR = 32)
172      INTEGER NT
173      PARAMETER (NT = 16)
174      INTEGER NR2
175      PARAMETER (NR2 = NR+2)
176      INTEGER NT1
177      PARAMETER (NT1 = NT+1)
178      INTEGER MMAX
179      PARAMETER (MMAX = 200)
180C
181C     Local variables
182C
183      DOUBLE PRECISION A0, A2, ACHEM, ADP0, ANG, AX, AY, B0, C0, C1, C2
184      DOUBLE PRECISION CCHAR, CMAX, CMIN, CO, CT, CX, CXH, CY, CYH, D
185      DOUBLE PRECISION DC, DN, ELF0, ELFCHAR, MU, PCHAR, PHI0, PI, R0
186      DOUBLE PRECISION RO, S0, SI, SX, SXH, SY, SYH, TCHAR, TH, TPI, U0
187      DOUBLE PRECISION UCHAR, WL, XCHAR, XH, XLOW, YH, YLOW, Z0
188C
189C     Common variables
190C
191      DOUBLE PRECISION C(0:NBP1,0:NBP1), CN(0:NBP1,0:NBP1)
192C
193      COMMON /ADP/ C, CN
194C
195C     Common variables
196C
197      DOUBLE PRECISION ZAG(0:NBP1,0:NBP1), ZAGN(0:NBP1,0:NBP1)
198C
199      COMMON /AGG/ ZAG, ZAGN
200C
201C     Common variables
202C
203      DOUBLE PRECISION FCHAR, PHICHAR, ZCHAR
204C
205      COMMON /CHAR/ FCHAR, ZCHAR, PHICHAR
206C
207C     Common variables
208C
209      DOUBLE PRECISION CND1, CND2, CND3, CND4, CND5
210C
211      COMMON /CND/ CND1, CND2, CND3, CND4, CND5
212C
213C     Common variables
214C
215      DOUBLE PRECISION A, RE, S1, S2, S3
216C
217      COMMON /COEFS/ S1, S2, S3, RE, A
218C
219C     Common variables
220C
221      DOUBLE PRECISION PEC
222C
223      COMMON /DIFFC/ PEC
224C
225C     Common variables
226C
227      DOUBLE PRECISION PEN
228C
229      COMMON /DIFFN/ PEN
230C
231C     Common variables
232C
233      DOUBLE PRECISION XM(MMAX,2)
234C
235      COMMON /FMARKERS/ XM
236C
237C     Common variables
238C
239      DOUBLE PRECISION FB(0:NBP1,0:NBP1,2)
240C
241      COMMON /FORCEB/ FB
242C
243C     Common variables
244C
245      DOUBLE PRECISION X(0:NBP1), Y(0:NBP1)
246C
247      COMMON /GRID/ X, Y
248C
249C     Common variables
250C
251      DOUBLE PRECISION ELF(0:NBP1,0:NBP1,0:NR2,NT1)
252      DOUBLE PRECISION ELFN(0:NBP1,0:NBP1,0:NR2,NT1)
253C
254      COMMON /LINK/ ELF, ELFN
255C
256C     Common variables
257C
258      DOUBLE PRECISION BETA0
259C
260      COMMON /LINKB/ BETA0
261C
262C     Common variables
263C
264      DOUBLE PRECISION ALPHA0
265C
266      COMMON /LINKF/ ALPHA0
267C
268C     Common variables
269C
270      DOUBLE PRECISION PHIA(0:NBP1,0:NBP1), PHIAN(0:NBP1,0:NBP1)
271      DOUBLE PRECISION PHIN(0:NBP1,0:NBP1), PHINN(0:NBP1,0:NBP1)
272C
273      COMMON /PHI/ PHIN, PHINN, PHIA, PHIAN
274C
275C     Common variables
276C
277      DOUBLE PRECISION THETA
278C
279      COMMON /PSTEPS/ THETA
280C
281C     Common variables
282C
283      DOUBLE PRECISION DR(NR2), R(0:NR2)
284C
285      INTEGER NK
286C
287      COMMON /RSIZE/ DR, R, NK
288C
289C     Common variables
290C
291      DOUBLE PRECISION DLT, H
292C
293      COMMON /STEPS/ H, DLT
294C
295C     Common variables
296C
297      DOUBLE PRECISION K0
298C
299      COMMON /STIFF/ K0
300C
301C     Common variables
302C
303      DOUBLE PRECISION U(0:NBP1,0:NBP1,2), UX(0:NBP1,0:NBP1)
304      DOUBLE PRECISION UY(0:NBP1,0:NBP1), VX(0:NBP1,0:NBP1)
305      DOUBLE PRECISION VY(0:NBP1,0:NBP1)
306C
307      COMMON /VEL/ U, UY, VY, UX, VX
308C
309C====>End Module   UINIT                  File t208x.f
310C====>Begin Module XINIT                  File t208x.f
311C---->Makedcls Options: All variables
312C
313C     Parameter variables
314C
315      INTEGER L2NG
316      PARAMETER (L2NG = 6)
317      INTEGER NG
318      PARAMETER (NG = 2**L2NG)
319      INTEGER NB
320      PARAMETER (NB = NG+2)
321      INTEGER NGM1
322      PARAMETER (NGM1 = NG-1)
323      INTEGER NBP1
324      PARAMETER (NBP1 = NB+1)
325C
326C     Local variables
327C
328      DOUBLE PRECISION PI, SI, SJ, TPON
329C
330C     Common variables
331C
332      DOUBLE PRECISION A, RE, S1, S2, S3
333C
334      COMMON /COEFS/ S1, S2, S3, RE, A
335C
336C     Common variables
337C
338      DOUBLE PRECISION BZERO(0:NBP1,0:NBP1), SINSQ(0:NBP1,0:NBP1)
339C
340      COMMON /FFT/ SINSQ, BZERO
341C
342C     Common variables
343C
344      DOUBLE PRECISION Z(0:NBP1,0:NBP1)
345C
346      COMMON /PTDS/ Z
347C
348C     Common variables
349C
350      DOUBLE PRECISION DLT, H
351C
352      COMMON /STEPS/ H, DLT
353C
354C====>End Module   XINIT                  File t208x.f
355C====>Begin Module NAVS2D                 File t208x.f
356C---->Makedcls Options: All variables
357C
358C     Parameter variables
359C
360      INTEGER L2NG
361      PARAMETER (L2NG = 6)
362      INTEGER NG
363      PARAMETER (NG = 2**L2NG)
364      INTEGER NB
365      PARAMETER (NB = NG+2)
366      INTEGER NGM1
367      PARAMETER (NGM1 = NG-1)
368      INTEGER NBP1
369      PARAMETER (NBP1 = NB+1)
370C
371C     Local variables
372C
373      DOUBLE PRECISION B(0:NBP1,0:NBP1), C(0:NBP1,0:NBP1), LAM(NG), PRO
374      DOUBLE PRECISION R(0:NBP1,0:NBP1), XL, YH(0:NBP1,0:NBP1), YL
375      DOUBLE PRECISION YP(0:NBP1,0:NBP1)
376C
377C     Common variables
378C
379      DOUBLE PRECISION A, RE, S1, S2, S3
380C
381      COMMON /COEFS/ S1, S2, S3, RE, A
382C
383C     Common variables
384C
385      DOUBLE PRECISION BZERO(0:NBP1,0:NBP1), SINSQ(0:NBP1,0:NBP1)
386C
387      COMMON /FFT/ SINSQ, BZERO
388C
389C     Common variables
390C
391      DOUBLE PRECISION F(0:NBP1,0:NBP1,2)
392C
393      COMMON /FORCE/ F
394C
395C     Equivalenced common variables
396C
397      DOUBLE PRECISION W1(0:NBP1,0:NBP1,1:2), W2(0:NBP1,0:NBP1,1:2)
398C
399C     Common variables
400C
401      DOUBLE PRECISION P(0:NBP1,0:NBP1)
402C
403      COMMON /PRES/ P
404C
405C     Equivalenced common variables
406C
407      DOUBLE PRECISION D(0:NBP1,0:NBP1)
408C
409C     Common variables
410C
411      DOUBLE PRECISION Z(0:NBP1,0:NBP1)
412C
413      COMMON /PTDS/ Z
414C
415C     Common variables
416C
417      DOUBLE PRECISION DLT, H
418C
419      COMMON /STEPS/ H, DLT
420C
421C     Common variables
422C
423      DOUBLE PRECISION U(0:NBP1,0:NBP1,2), UX(0:NBP1,0:NBP1)
424      DOUBLE PRECISION UY(0:NBP1,0:NBP1), VX(0:NBP1,0:NBP1)
425      DOUBLE PRECISION VY(0:NBP1,0:NBP1)
426C
427      COMMON /VEL/ U, UY, VY, UX, VX
428C
429C====>End Module   NAVS2D                 File t208x.f
430C====>Begin Module TRIDGX                 File t208x.f
431C---->Makedcls Options: All variables
432C
433C     Parameter variables
434C
435      INTEGER L2NG
436      PARAMETER (L2NG = 6)
437      INTEGER NG
438      PARAMETER (NG = 2**L2NG)
439      INTEGER NB
440      PARAMETER (NB = NG+2)
441      INTEGER NGM1
442      PARAMETER (NGM1 = NG-1)
443      INTEGER N
444      PARAMETER (N = NG-1)
445      INTEGER NP1
446      PARAMETER (NP1 = N+1)
447      INTEGER NM1
448      PARAMETER (NM1 = N-1)
449      INTEGER NBP1
450      PARAMETER (NBP1 = NB+1)
451C
452C     Argument variables
453C
454      DOUBLE PRECISION A, B(0:NBP1,0:NBP1), C(0:NBP1,0:NBP1)
455      DOUBLE PRECISION W(0:NBP1,0:NBP1), Y(0:NBP1,0:NBP1)
456C
457C     Common variables
458C
459      DOUBLE PRECISION D(0:NBP1,0:NBP1), L(0:NBP1,0:NBP1)
460      DOUBLE PRECISION R(0:NBP1,0:NBP1), X(0:NBP1,0:NBP1)
461C
462      COMMON /TDSPACE/ X, L, R, D
463C
464C====>End Module   TRIDGX                 File t208x.f
465C====>Begin Module TRIDGY                 File t208x.f
466C---->Makedcls Options: All variables
467C
468C     Parameter variables
469C
470      INTEGER L2NG
471      PARAMETER (L2NG = 6)
472      INTEGER NG
473      PARAMETER (NG = 2**L2NG)
474      INTEGER NB
475      PARAMETER (NB = NG+2)
476      INTEGER NGM1
477      PARAMETER (NGM1 = NG-1)
478      INTEGER N
479      PARAMETER (N = NG-1)
480      INTEGER NP1
481      PARAMETER (NP1 = N+1)
482      INTEGER NM1
483      PARAMETER (NM1 = N-1)
484      INTEGER NBP1
485      PARAMETER (NBP1 = NB+1)
486C
487C     Argument variables
488C
489      DOUBLE PRECISION A, B(0:NBP1,0:NBP1), C(0:NBP1,0:NBP1)
490      DOUBLE PRECISION W(0:NBP1,0:NBP1), Y(0:NBP1,0:NBP1)
491C
492C     Common variables
493C
494      DOUBLE PRECISION D(0:NBP1,0:NBP1), L(0:NBP1,0:NBP1)
495      DOUBLE PRECISION R(0:NBP1,0:NBP1), X(0:NBP1,0:NBP1)
496C
497      COMMON /TDSPACE/ X, L, R, D
498C
499C====>End Module   TRIDGY                 File t208x.f
500C====>Begin Module NEWU                   File t208x.f
501C---->Makedcls Options: All variables
502C
503C     Parameter variables
504C
505      INTEGER L2NG
506      PARAMETER (L2NG = 6)
507      INTEGER NG
508      PARAMETER (NG = 2**L2NG)
509      INTEGER NB
510      PARAMETER (NB = NG+2)
511      INTEGER NGM1
512      PARAMETER (NGM1 = NG-1)
513      INTEGER NBP1
514      PARAMETER (NBP1 = NB+1)
515C
516C     Argument variables
517C
518      DOUBLE PRECISION P(0:NBP1,0:NBP1), W2(0:NBP1,0:NBP1,2)
519C
520C     Local variables
521C
522      DOUBLE PRECISION S
523C
524C     Common variables
525C
526      DOUBLE PRECISION H, K
527C
528      COMMON /STEPS/ H, K
529C
530C     Common variables
531C
532      DOUBLE PRECISION U(0:NBP1,0:NBP1,2), UX(0:NBP1,0:NBP1)
533      DOUBLE PRECISION UY(0:NBP1,0:NBP1), VX(0:NBP1,0:NBP1)
534      DOUBLE PRECISION VY(0:NBP1,0:NBP1)
535C
536      COMMON /VEL/ U, UY, VY, UX, VX
537C
538C====>End Module   NEWU                   File t208x.f
539C====>Begin Module DIV                    File t208x.f
540C---->Makedcls Options: All variables
541C
542C     Parameter variables
543C
544      INTEGER L2NG
545      PARAMETER (L2NG = 6)
546      INTEGER NG
547      PARAMETER (NG = 2**L2NG)
548      INTEGER NB
549      PARAMETER (NB = NG+2)
550      INTEGER NGM1
551      PARAMETER (NGM1 = NG-1)
552      INTEGER NBP1
553      PARAMETER (NBP1 = NB+1)
554C
555C     Argument variables
556C
557      DOUBLE PRECISION D(0:NBP1,0:NBP1), S3, W2(0:NBP1,0:NBP1,2)
558C
559C====>End Module   DIV                    File t208x.f
560C====>Begin Module FFT2D                  File t208x.f
561C---->Makedcls Options: All variables
562C
563C     Parameter variables
564C
565      INTEGER L2NG
566      PARAMETER (L2NG = 6)
567      INTEGER NG
568      PARAMETER (NG = 2**L2NG)
569      INTEGER NB
570      PARAMETER (NB = NG+2)
571      INTEGER NGM1
572      PARAMETER (NGM1 = NG-1)
573      INTEGER N
574      PARAMETER (N = NG)
575      INTEGER M
576      PARAMETER (M = L2NG)
577      INTEGER NBP1
578      PARAMETER (NBP1 = NB+1)
579C
580C     Argument variables
581C
582      DOUBLE PRECISION A(0:NBP1,0:NBP1), B(0:NBP1,0:NBP1)
583C
584      INTEGER ISIGN
585C
586C     Local variables
587C
588      DOUBLE PRECISION ANG, PI, SSIGN, T1(N), T2(N), T3(N), T4(N), TU1
589      DOUBLE PRECISION U1, U2, W1, W2
590C
591C====>End Module   FFT2D                  File t208x.f
592C====>Begin Module PHIA2D                 File t208x.f
593C---->Makedcls Options: All variables
594C
595C     Parameter variables
596C
597      INTEGER L2NG
598      PARAMETER (L2NG = 6)
599      INTEGER NG
600      PARAMETER (NG = 2**L2NG)
601      INTEGER NB
602      PARAMETER (NB = NG+2)
603      INTEGER NGM1
604      PARAMETER (NGM1 = NG-1)
605      INTEGER NGP1
606      PARAMETER (NGP1 = NG+1)
607      INTEGER NGP2
608      PARAMETER (NGP2 = NG+2)
609      INTEGER NBP1
610      PARAMETER (NBP1 = NB+1)
611C
612C     Argument variables
613C
614      INTEGER N
615C
616C     Common variables
617C
618      DOUBLE PRECISION PHIA(0:NBP1,0:NBP1), PHIAN(0:NBP1,0:NBP1)
619      DOUBLE PRECISION PHIN(0:NBP1,0:NBP1), PHINN(0:NBP1,0:NBP1)
620C
621      COMMON /PHI/ PHIN, PHINN, PHIA, PHIAN
622C
623C     Common variables
624C
625      DOUBLE PRECISION H, K
626C
627      COMMON /STEPS/ H, K
628C
629C     Common variables
630C
631      DOUBLE PRECISION U(0:NBP1,0:NBP1,2), UX(0:NBP1,0:NBP1)
632      DOUBLE PRECISION UY(0:NBP1,0:NBP1), VX(0:NBP1,0:NBP1)
633      DOUBLE PRECISION VY(0:NBP1,0:NBP1)
634C
635      COMMON /VEL/ U, UY, VY, UX, VX
636C
637C====>End Module   PHIA2D                 File t208x.f
638C====>Begin Module ZAG2D                  File t208x.f
639C---->Makedcls Options: All variables
640C
641C     Parameter variables
642C
643      INTEGER L2NG
644      PARAMETER (L2NG = 6)
645      INTEGER NG
646      PARAMETER (NG = 2**L2NG)
647      INTEGER NB
648      PARAMETER (NB = NG+2)
649      INTEGER NGM1
650      PARAMETER (NGM1 = NG-1)
651      INTEGER NGP1
652      PARAMETER (NGP1 = NG+1)
653      INTEGER NGP2
654      PARAMETER (NGP2 = NG+2)
655      INTEGER NBP1
656      PARAMETER (NBP1 = NB+1)
657C
658C     Argument variables
659C
660      INTEGER N
661C
662C     Common variables
663C
664      DOUBLE PRECISION ZAG(0:NBP1,0:NBP1), ZAGN(0:NBP1,0:NBP1)
665C
666      COMMON /AGG/ ZAG, ZAGN
667C
668C     Common variables
669C
670      DOUBLE PRECISION H, K
671C
672      COMMON /STEPS/ H, K
673C
674C     Common variables
675C
676      DOUBLE PRECISION U(0:NBP1,0:NBP1,2), UX(0:NBP1,0:NBP1)
677      DOUBLE PRECISION UY(0:NBP1,0:NBP1), VX(0:NBP1,0:NBP1)
678      DOUBLE PRECISION VY(0:NBP1,0:NBP1)
679C
680      COMMON /VEL/ U, UY, VY, UX, VX
681C
682C====>End Module   ZAG2D                  File t208x.f
683C====>Begin Module PLOT                   File t208x.f
684C---->Makedcls Options: All variables
685C
686C     Parameter variables
687C
688      INTEGER L2NG
689      PARAMETER (L2NG = 6)
690      INTEGER NG
691      PARAMETER (NG = 2**L2NG)
692      INTEGER NB
693      PARAMETER (NB = NG+2)
694      INTEGER NGM1
695      PARAMETER (NGM1 = NG-1)
696      INTEGER NGP1
697      PARAMETER (NGP1 = NG+1)
698      INTEGER NGP2
699      PARAMETER (NGP2 = NG+2)
700      INTEGER NBP1
701      PARAMETER (NBP1 = NB+1)
702      INTEGER NR
703      PARAMETER (NR = 32)
704      INTEGER NT
705      PARAMETER (NT = 16)
706      INTEGER NR2
707      PARAMETER (NR2 = NR+2)
708      INTEGER NT1
709      PARAMETER (NT1 = NT+1)
710      INTEGER MMAX
711      PARAMETER (MMAX = 200)
712C
713C     Argument variables
714C
715      DOUBLE PRECISION C(0:NBP1,0:NBP1), ELF(0:NBP1,0:NBP1,0:NR2,NT1)
716      DOUBLE PRECISION F(0:NBP1,0:NBP1,2), PHIA(0:NBP1,0:NBP1)
717      DOUBLE PRECISION PHIN(0:NBP1,0:NBP1), TIME, U(0:NBP1,0:NBP1,2)
718      DOUBLE PRECISION XM(MMAX,2), ZAG(0:NBP1,0:NBP1)
719C
720      INTEGER N, NK
721C
722C     Common variables
723C
724      INTEGER NFIL
725C
726      COMMON /EFNUM/ NFIL
727C
728C====>End Module   PLOT                   File t208x.f
729C====>Begin Module MARKERS                File t208x.f
730C---->Makedcls Options: All variables
731C
732C     Parameter variables
733C
734      INTEGER L2NG
735      PARAMETER (L2NG = 6)
736      INTEGER NG
737      PARAMETER (NG = 2**L2NG)
738      INTEGER NB
739      PARAMETER (NB = NG+2)
740      INTEGER NGM1
741      PARAMETER (NGM1 = NG-1)
742      INTEGER NGP1
743      PARAMETER (NGP1 = NG+1)
744      INTEGER NGP2
745      PARAMETER (NGP2 = NG+2)
746      INTEGER NBP1
747      PARAMETER (NBP1 = NB+1)
748      INTEGER MMAX
749      PARAMETER (MMAX = 200)
750C
751C     Local variables
752C
753      DOUBLE PRECISION AX, AY, UM, UMN, VM, VMN, X, XN, Y, YN
754C
755C     Common variables
756C
757      DOUBLE PRECISION XM(MMAX,2)
758C
759      COMMON /FMARKERS/ XM
760C
761C     Common variables
762C
763      DOUBLE PRECISION DLT, H
764C
765      COMMON /STEPS/ H, DLT
766C
767C     Common variables
768C
769      DOUBLE PRECISION U(0:NBP1,0:NBP1,2), UX(0:NBP1,0:NBP1)
770      DOUBLE PRECISION UY(0:NBP1,0:NBP1), VX(0:NBP1,0:NBP1)
771      DOUBLE PRECISION VY(0:NBP1,0:NBP1)
772C
773      COMMON /VEL/ U, UY, VY, UX, VX
774C
775C====>End Module   MARKERS                File t208x.f
776C====>Begin Module PHIN2D                 File t208x.f
777C---->Makedcls Options: All variables
778C
779C     Parameter variables
780C
781      INTEGER L2NG
782      PARAMETER (L2NG = 6)
783      INTEGER NG
784      PARAMETER (NG = 2**L2NG)
785      INTEGER NB
786      PARAMETER (NB = NG+2)
787      INTEGER NGM1
788      PARAMETER (NGM1 = NG-1)
789      INTEGER NGP1
790      PARAMETER (NGP1 = NG+1)
791      INTEGER NGP2
792      PARAMETER (NGP2 = NG+2)
793      INTEGER NBP1
794      PARAMETER (NBP1 = NB+1)
795      INTEGER MMAX
796      PARAMETER (MMAX = 200)
797C
798C     Argument variables
799C
800      INTEGER N
801C
802C     Common variables
803C
804      DOUBLE PRECISION PEN
805C
806      COMMON /DIFFN/ PEN
807C
808C     Common variables
809C
810      DOUBLE PRECISION PHIA(0:NBP1,0:NBP1), PHIAN(0:NBP1,0:NBP1)
811      DOUBLE PRECISION PHIN(0:NBP1,0:NBP1), PHINN(0:NBP1,0:NBP1)
812C
813      COMMON /PHI/ PHIN, PHINN, PHIA, PHIAN
814C
815C     Common variables
816C
817      DOUBLE PRECISION H, K
818C
819      COMMON /STEPS/ H, K
820C
821C     Common variables
822C
823      DOUBLE PRECISION U(0:NBP1,0:NBP1,2), UX(0:NBP1,0:NBP1)
824      DOUBLE PRECISION UY(0:NBP1,0:NBP1), VX(0:NBP1,0:NBP1)
825      DOUBLE PRECISION VY(0:NBP1,0:NBP1)
826C
827      COMMON /VEL/ U, UY, VY, UX, VX
828C
829C====>End Module   PHIN2D                 File t208x.f
830C====>Begin Module CHEM2D                 File t208x.f
831C---->Makedcls Options: All variables
832C
833C     Parameter variables
834C
835      INTEGER L2NG
836      PARAMETER (L2NG = 6)
837      INTEGER NG
838      PARAMETER (NG = 2**L2NG)
839      INTEGER NB
840      PARAMETER (NB = NG+2)
841      INTEGER NGM1
842      PARAMETER (NGM1 = NG-1)
843      INTEGER NGP1
844      PARAMETER (NGP1 = NG+1)
845      INTEGER NGP2
846      PARAMETER (NGP2 = NG+2)
847      INTEGER NBP1
848      PARAMETER (NBP1 = NB+1)
849      INTEGER MMAX
850      PARAMETER (MMAX = 200)
851C
852C     Argument variables
853C
854      INTEGER N
855C
856C     Common variables
857C
858      DOUBLE PRECISION C(0:NBP1,0:NBP1), CN(0:NBP1,0:NBP1)
859C
860      COMMON /ADP/ C, CN
861C
862C     Common variables
863C
864      DOUBLE PRECISION PEC
865C
866      COMMON /DIFFC/ PEC
867C
868C     Common variables
869C
870      DOUBLE PRECISION H, K
871C
872      COMMON /STEPS/ H, K
873C
874C     Common variables
875C
876      DOUBLE PRECISION U(0:NBP1,0:NBP1,2), UX(0:NBP1,0:NBP1)
877      DOUBLE PRECISION UY(0:NBP1,0:NBP1), VX(0:NBP1,0:NBP1)
878      DOUBLE PRECISION VY(0:NBP1,0:NBP1)
879C
880      COMMON /VEL/ U, UY, VY, UX, VX
881C
882C====>End Module   CHEM2D                 File t208x.f
883C====>Begin Module CHLFAC                 File t208x.f
884C---->Makedcls Options: All variables
885C
886C     Argument variables
887C
888      DOUBLE PRECISION DIAG(NN), SUBD(NN)
889C
890      INTEGER NN
891C
892C====>End Module   CHLFAC                 File t208x.f
893C====>Begin Module CHLSLV                 File t208x.f
894C---->Makedcls Options: All variables
895C
896C     Parameter variables
897C
898      INTEGER L2NG
899      PARAMETER (L2NG = 6)
900      INTEGER NG
901      PARAMETER (NG = 2**L2NG)
902      INTEGER NB
903      PARAMETER (NB = NG+2)
904      INTEGER NGM1
905      PARAMETER (NGM1 = NG-1)
906      INTEGER NGP1
907      PARAMETER (NGP1 = NG+1)
908      INTEGER NGP2
909      PARAMETER (NGP2 = NG+2)
910      INTEGER NBP1
911      PARAMETER (NBP1 = NB+1)
912      INTEGER MAX
913      PARAMETER (MAX = NBP1)
914C
915C     Argument variables
916C
917      DOUBLE PRECISION B(NM1+1), DIAG(NM1), SUB(NM1), X(NM1)
918C
919      INTEGER NM1
920C
921C     Local variables
922C
923      DOUBLE PRECISION Y(MAX)
924C
925C====>End Module   CHLSLV                 File t208x.f
926C====>Begin Module PERIOD                 File t208x.f
927C---->Makedcls Options: All variables
928C
929C     Parameter variables
930C
931      INTEGER L2NG
932      PARAMETER (L2NG = 6)
933      INTEGER NG
934      PARAMETER (NG = 2**L2NG)
935      INTEGER NB
936      PARAMETER (NB = NG+2)
937      INTEGER NGM1
938      PARAMETER (NGM1 = NG-1)
939      INTEGER NGP1
940      PARAMETER (NGP1 = NG+1)
941      INTEGER NGP2
942      PARAMETER (NGP2 = NG+2)
943      INTEGER NBP1
944      PARAMETER (NBP1 = NB+1)
945C
946C     Argument variables
947C
948      DOUBLE PRECISION PH(0:NBP1,0:NBP1)
949C
950C====>End Module   PERIOD                 File t208x.f
951C====>Begin Module SOL                    File t208x.f
952C---->Makedcls Options: All variables
953C
954C     Argument variables
955C
956      DOUBLE PRECISION L, R, V
957C
958C====>End Module   SOL                    File t208x.f
959C====>Begin Module DIFFU                  File t208x.f
960C---->Makedcls Options: All variables
961C
962C     Parameter variables
963C
964      INTEGER L2NG
965      PARAMETER (L2NG = 6)
966      INTEGER NG
967      PARAMETER (NG = 2**L2NG)
968      INTEGER NB
969      PARAMETER (NB = NG+2)
970      INTEGER NGM1
971      PARAMETER (NGM1 = NG-1)
972      INTEGER NGP1
973      PARAMETER (NGP1 = NG+1)
974      INTEGER NGP2
975      PARAMETER (NGP2 = NG+2)
976      INTEGER NBP1
977      PARAMETER (NBP1 = NB+1)
978C
979C     Argument variables
980C
981      DOUBLE PRECISION ALPHA, K1, PH(0:NBP1,0:NBP1), PHN(0:NBP1,0:NBP1)
982C
983C     Local variables
984C
985      DOUBLE PRECISION B(NB), ETA, GDIAG(NB-1), GSUB(NB-1), LAM, Y(NB-1)
986      DOUBLE PRECISION Z(NB-1)
987C
988C     Common variables
989C
990      DOUBLE PRECISION H, K
991C
992      COMMON /STEPS/ H, K
993C
994C====>End Module   DIFFU                  File t208x.f
995C====>Begin Module REACT                  File t208x.f
996C---->Makedcls Options: All variables
997C
998C     Parameter variables
999C
1000      INTEGER L2NG
1001      PARAMETER (L2NG = 6)
1002      INTEGER NG
1003      PARAMETER (NG = 2**L2NG)
1004      INTEGER NB
1005      PARAMETER (NB = NG+2)
1006      INTEGER NGM1
1007      PARAMETER (NGM1 = NG-1)
1008      INTEGER NGP1
1009      PARAMETER (NGP1 = NG+1)
1010      INTEGER NGP2
1011      PARAMETER (NGP2 = NG+2)
1012      INTEGER NBP1
1013      PARAMETER (NBP1 = NB+1)
1014      INTEGER NR
1015      PARAMETER (NR = 32)
1016      INTEGER NT
1017      PARAMETER (NT = 16)
1018      INTEGER NR2
1019      PARAMETER (NR2 = NR+2)
1020      INTEGER NT1
1021      PARAMETER (NT1 = NT+1)
1022C
1023C     Local variables
1024C
1025      DOUBLE PRECISION ALPHA, C4, RC(0:NBP1), TEMP1, TEMP2, TEMP3, TEMP4
1026      DOUBLE PRECISION Y(0:NBP1,4), YN(0:NBP1,4)
1027C
1028C     Common variables
1029C
1030      DOUBLE PRECISION C(0:NBP1,0:NBP1), CN(0:NBP1,0:NBP1)
1031C
1032      COMMON /ADP/ C, CN
1033C
1034C     Common variables
1035C
1036      DOUBLE PRECISION ZAG(0:NBP1,0:NBP1), ZAGN(0:NBP1,0:NBP1)
1037C
1038      COMMON /AGG/ ZAG, ZAGN
1039C
1040C     Common variables
1041C
1042      DOUBLE PRECISION CND1, CND2, CND3, CND4, CND5
1043C
1044      COMMON /CND/ CND1, CND2, CND3, CND4, CND5
1045C
1046C     Common variables
1047C
1048      DOUBLE PRECISION ELF(0:NBP1,0:NBP1,0:NR2,NT1)
1049      DOUBLE PRECISION ELFN(0:NBP1,0:NBP1,0:NR2,NT1)
1050C
1051      COMMON /LINK/ ELF, ELFN
1052C
1053C     Common variables
1054C
1055      DOUBLE PRECISION PHIA(0:NBP1,0:NBP1), PHIAN(0:NBP1,0:NBP1)
1056      DOUBLE PRECISION PHIN(0:NBP1,0:NBP1), PHINN(0:NBP1,0:NBP1)
1057C
1058      COMMON /PHI/ PHIN, PHINN, PHIA, PHIAN
1059C
1060C     Common variables
1061C
1062      DOUBLE PRECISION DR(NR2), R(0:NR2)
1063C
1064      INTEGER NK
1065C
1066      COMMON /RSIZE/ DR, R, NK
1067C
1068C     Common variables
1069C
1070      DOUBLE PRECISION H, K
1071C
1072      COMMON /STEPS/ H, K
1073C
1074C====>End Module   REACT                  File t208x.f
1075