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