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 36 CHARACTER*18 e13file, e14file, e15file, e16file 37 CHARACTER*18 e17file, e18file, e19file, e1file 38 CHARACTER*18 e20file, e21file, e22file, e23file 39 CHARACTER*18 e24file, e2file, e3file, e4file 40 CHARACTER*18 e5file, e6file, e7file, e8file 41 CHARACTER*18 e9file, efile, ffile, fmfile 42 CHARACTER*18 maifile, mfile, mmfile, pafile 43 CHARACTER*18 parfile, pnfile 44 CHARACTER*15 rname 45 CHARACTER*18 symfile, ufile, zfile 46C 47 DOUBLE PRECISION cmax, cmin, f0, fbmax 48 DOUBLE PRECISION flmax, phiamax, phiamin, phinmax 49 DOUBLE PRECISION phinmin, q, tcoef, tfreq 50 DOUBLE PRECISION time, tpi, umax 51 DOUBLE PRECISION uold(0:nbp1,0:nbp1,2), usmax 52 DOUBLE PRECISION zagi(0:nbp1,0:nbp1), zagmax, zagmin 53C 54 INTEGER rlen 55C 56 LOGICAL avrage(06), axes3d, bclip, box 57 LOGICAL byrows, chkdup, deriv(06) 58C 59 REAL dusub, dvsub, dxsub, dysub 60 REAL dzsub 61C 62C Namelists 63C 64 NAMELIST / data / avrage, axes3d, bclip, box 65 x, byrows, chkdup, deriv, dusub 66 x, dvsub, dxsub, dysub, dzsub 67C 68C Common variables 69C 70 DOUBLE PRECISION c(0:nbp1,0:nbp1), cn(0:nbp1,0:nbp1) 71C 72 COMMON / adp / c, cn 73C 74C Common variables 75C 76 DOUBLE PRECISION zag(0:nbp1,0:nbp1), zagn(0:nbp1,0:nbp1) 77C 78 COMMON / agg / zag, zagn 79C 80C Common variables 81C 82 INTEGER nfil 83C 84 COMMON / efnum / nfil 85C 86C Common variables 87C 88 DOUBLE PRECISION xm(mmax,2) 89C 90 COMMON / fmarkers/ xm 91C 92C Common variables 93C 94 DOUBLE PRECISION f(0:nbp1,0:nbp1,2) 95C 96 COMMON / force / f 97C 98C Common variables 99C 100 DOUBLE PRECISION fb(0:nbp1,0:nbp1,2) 101C 102 COMMON / forceb / fb 103C 104C Common variables 105C 106 DOUBLE PRECISION elf(0:nbp1,0:nbp1,0:nr2,nt1) 107 DOUBLE PRECISION elfn(0:nbp1,0:nbp1,0:nr2,nt1) 108C 109 COMMON / link / elf, elfn 110C 111C Common variables 112C 113 INTEGER method, mthlim 114C 115 COMMON / mth / method, mthlim 116C 117C Common variables 118C 119 DOUBLE PRECISION phia(0:nbp1,0:nbp1), phian(0:nbp1,0:nbp1) 120 DOUBLE PRECISION phin(0:nbp1,0:nbp1), phinn(0:nbp1,0:nbp1) 121C 122 COMMON / phi / phin, phinn, phia, phian 123C 124C Common variables 125C 126 DOUBLE PRECISION dr(nr2), r(0:nr2) 127C 128 INTEGER nk 129C 130 COMMON / rsize / dr, r, nk 131C 132C Common variables 133C 134 DOUBLE PRECISION dlt, h 135C 136 COMMON / steps / h, dlt 137C 138C Common variables 139C 140 DOUBLE PRECISION k0 141C 142 COMMON / stiff / k0 143C 144C Common variables 145C 146 DOUBLE PRECISION u(0:nbp1,0:nbp1,2), ux(0:nbp1,0:nbp1) 147 DOUBLE PRECISION uy(0:nbp1,0:nbp1), vx(0:nbp1,0:nbp1) 148 DOUBLE PRECISION vy(0:nbp1,0:nbp1) 149C 150 COMMON / vel / u, uy, vy, ux 151 COMMON / vel / vx 152C 153C Common variables 154C 155 DOUBLE PRECISION wsave(nw) 156C 157 COMMON / wave / wsave 158C 159C====>End Module PROB5_4DIM File t208x.f 160C====>Begin Module UINIT File t208x.f 161C---->Makedcls Options: All variables 162C 163C Parameter variables 164C 165 INTEGER l2ng 166 PARAMETER (l2ng = 6) 167 INTEGER ng 168 PARAMETER (ng = 2**l2ng) 169 INTEGER nb 170 PARAMETER (nb = ng+2) 171 INTEGER ngm1 172 PARAMETER (ngm1 = ng-1) 173 INTEGER ngp1 174 PARAMETER (ngp1 = ng+1) 175 INTEGER ngp2 176 PARAMETER (ngp2 = ng+2) 177 INTEGER nbp1 178 PARAMETER (nbp1 = nb+1) 179 INTEGER nr 180 PARAMETER (nr = 32) 181 INTEGER nt 182 PARAMETER (nt = 16) 183 INTEGER nr2 184 PARAMETER (nr2 = nr+2) 185 INTEGER nt1 186 PARAMETER (nt1 = nt+1) 187 INTEGER mmax 188 PARAMETER (mmax = 200) 189C 190C Local variables 191C 192 DOUBLE PRECISION a0, a2, achem, adp0 193 DOUBLE PRECISION ang, ax, ay, b0 194 DOUBLE PRECISION c0, c1, c2, cchar 195 DOUBLE PRECISION cmax, cmin, co, ct 196 DOUBLE PRECISION cx, cxh, cy, cyh 197 DOUBLE PRECISION d, dc, dn, elf0 198 DOUBLE PRECISION elfchar, mu, pchar, phi0 199 DOUBLE PRECISION pi, r0, ro, s0 200 DOUBLE PRECISION si, sx, sxh, sy 201 DOUBLE PRECISION syh, tchar, th, tpi 202 DOUBLE PRECISION u0, uchar, wl, xchar 203 DOUBLE PRECISION xh, xlow, yh, ylow 204 DOUBLE PRECISION z0 205C 206C Common variables 207C 208 DOUBLE PRECISION c(0:nbp1,0:nbp1), cn(0:nbp1,0:nbp1) 209C 210 COMMON / adp / c, cn 211C 212C Common variables 213C 214 DOUBLE PRECISION zag(0:nbp1,0:nbp1), zagn(0:nbp1,0:nbp1) 215C 216 COMMON / agg / zag, zagn 217C 218C Common variables 219C 220 DOUBLE PRECISION fchar, phichar, zchar 221C 222 COMMON / char / fchar, zchar, phichar 223C 224C Common variables 225C 226 DOUBLE PRECISION cnd1, cnd2, cnd3, cnd4 227 DOUBLE PRECISION cnd5 228C 229 COMMON / cnd / cnd1, cnd2, cnd3, cnd4 230 COMMON / cnd / cnd5 231C 232C Common variables 233C 234 DOUBLE PRECISION a, re, s1, s2 235 DOUBLE PRECISION s3 236C 237 COMMON / coefs / s1, s2, s3, re 238 COMMON / coefs / a 239C 240C Common variables 241C 242 DOUBLE PRECISION pec 243C 244 COMMON / diffc / pec 245C 246C Common variables 247C 248 DOUBLE PRECISION pen 249C 250 COMMON / diffn / pen 251C 252C Common variables 253C 254 DOUBLE PRECISION xm(mmax,2) 255C 256 COMMON / fmarkers/ xm 257C 258C Common variables 259C 260 DOUBLE PRECISION fb(0:nbp1,0:nbp1,2) 261C 262 COMMON / forceb / fb 263C 264C Common variables 265C 266 DOUBLE PRECISION x(0:nbp1), y(0:nbp1) 267C 268 COMMON / grid / x, y 269C 270C Common variables 271C 272 DOUBLE PRECISION elf(0:nbp1,0:nbp1,0:nr2,nt1) 273 DOUBLE PRECISION elfn(0:nbp1,0:nbp1,0:nr2,nt1) 274C 275 COMMON / link / elf, elfn 276C 277C Common variables 278C 279 DOUBLE PRECISION beta0 280C 281 COMMON / linkb / beta0 282C 283C Common variables 284C 285 DOUBLE PRECISION alpha0 286C 287 COMMON / linkf / alpha0 288C 289C Common variables 290C 291 DOUBLE PRECISION phia(0:nbp1,0:nbp1), phian(0:nbp1,0:nbp1) 292 DOUBLE PRECISION phin(0:nbp1,0:nbp1), phinn(0:nbp1,0:nbp1) 293C 294 COMMON / phi / phin, phinn, phia, phian 295C 296C Common variables 297C 298 DOUBLE PRECISION theta 299C 300 COMMON / psteps / theta 301C 302C Common variables 303C 304 DOUBLE PRECISION dr(nr2), r(0:nr2) 305C 306 INTEGER nk 307C 308 COMMON / rsize / dr, r, nk 309C 310C Common variables 311C 312 DOUBLE PRECISION dlt, h 313C 314 COMMON / steps / h, dlt 315C 316C Common variables 317C 318 DOUBLE PRECISION k0 319C 320 COMMON / stiff / k0 321C 322C Common variables 323C 324 DOUBLE PRECISION u(0:nbp1,0:nbp1,2), ux(0:nbp1,0:nbp1) 325 DOUBLE PRECISION uy(0:nbp1,0:nbp1), vx(0:nbp1,0:nbp1) 326 DOUBLE PRECISION vy(0:nbp1,0:nbp1) 327C 328 COMMON / vel / u, uy, vy, ux 329 COMMON / vel / vx 330C 331C====>End Module UINIT File t208x.f 332C====>Begin Module XINIT File t208x.f 333C---->Makedcls Options: All variables 334C 335C Parameter variables 336C 337 INTEGER l2ng 338 PARAMETER (l2ng = 6) 339 INTEGER ng 340 PARAMETER (ng = 2**l2ng) 341 INTEGER nb 342 PARAMETER (nb = ng+2) 343 INTEGER ngm1 344 PARAMETER (ngm1 = ng-1) 345 INTEGER nbp1 346 PARAMETER (nbp1 = nb+1) 347C 348C Local variables 349C 350 DOUBLE PRECISION pi, si, sj, tpon 351C 352C Common variables 353C 354 DOUBLE PRECISION a, re, s1, s2 355 DOUBLE PRECISION s3 356C 357 COMMON / coefs / s1, s2, s3, re 358 COMMON / coefs / a 359C 360C Common variables 361C 362 DOUBLE PRECISION bzero(0:nbp1,0:nbp1), sinsq(0:nbp1,0:nbp1) 363C 364 COMMON / fft / sinsq, bzero 365C 366C Common variables 367C 368 DOUBLE PRECISION z(0:nbp1,0:nbp1) 369C 370 COMMON / ptds / z 371C 372C Common variables 373C 374 DOUBLE PRECISION dlt, h 375C 376 COMMON / steps / h, dlt 377C 378C====>End Module XINIT File t208x.f 379C====>Begin Module NAVS2D File t208x.f 380C---->Makedcls Options: All variables 381C 382C Parameter variables 383C 384 INTEGER l2ng 385 PARAMETER (l2ng = 6) 386 INTEGER ng 387 PARAMETER (ng = 2**l2ng) 388 INTEGER nb 389 PARAMETER (nb = ng+2) 390 INTEGER ngm1 391 PARAMETER (ngm1 = ng-1) 392 INTEGER nbp1 393 PARAMETER (nbp1 = nb+1) 394C 395C Local variables 396C 397 DOUBLE PRECISION b(0:nbp1,0:nbp1), c(0:nbp1,0:nbp1) 398 DOUBLE PRECISION lam(ng), pro, r(0:nbp1,0:nbp1) 399 DOUBLE PRECISION xl, yh(0:nbp1,0:nbp1), yl 400 DOUBLE PRECISION yp(0:nbp1,0:nbp1) 401C 402C Common variables 403C 404 DOUBLE PRECISION a, re, s1, s2 405 DOUBLE PRECISION s3 406C 407 COMMON / coefs / s1, s2, s3, re 408 COMMON / coefs / a 409C 410C Common variables 411C 412 DOUBLE PRECISION bzero(0:nbp1,0:nbp1), sinsq(0:nbp1,0:nbp1) 413C 414 COMMON / fft / sinsq, bzero 415C 416C Common variables 417C 418 DOUBLE PRECISION f(0:nbp1,0:nbp1,2) 419C 420 COMMON / force / f 421C 422C Equivalenced common variables 423C 424 DOUBLE PRECISION w1(0:nbp1,0:nbp1,1:2) 425 DOUBLE PRECISION w2(0:nbp1,0:nbp1,1:2) 426C 427C Common variables 428C 429 DOUBLE PRECISION p(0:nbp1,0:nbp1) 430C 431 COMMON / pres / p 432C 433C Equivalenced common variables 434C 435 DOUBLE PRECISION d(0:nbp1,0:nbp1) 436C 437C Common variables 438C 439 DOUBLE PRECISION z(0:nbp1,0:nbp1) 440C 441 COMMON / ptds / z 442C 443C Common variables 444C 445 DOUBLE PRECISION dlt, h 446C 447 COMMON / steps / h, dlt 448C 449C Common variables 450C 451 DOUBLE PRECISION u(0:nbp1,0:nbp1,2), ux(0:nbp1,0:nbp1) 452 DOUBLE PRECISION uy(0:nbp1,0:nbp1), vx(0:nbp1,0:nbp1) 453 DOUBLE PRECISION vy(0:nbp1,0:nbp1) 454C 455 COMMON / vel / u, uy, vy, ux 456 COMMON / vel / vx 457C 458C====>End Module NAVS2D File t208x.f 459C====>Begin Module TRIDGX File t208x.f 460C---->Makedcls Options: All variables 461C 462C Parameter variables 463C 464 INTEGER l2ng 465 PARAMETER (l2ng = 6) 466 INTEGER ng 467 PARAMETER (ng = 2**l2ng) 468 INTEGER nb 469 PARAMETER (nb = ng+2) 470 INTEGER ngm1 471 PARAMETER (ngm1 = ng-1) 472 INTEGER n 473 PARAMETER (n = ng-1) 474 INTEGER np1 475 PARAMETER (np1 = n+1) 476 INTEGER nm1 477 PARAMETER (nm1 = n-1) 478 INTEGER nbp1 479 PARAMETER (nbp1 = nb+1) 480C 481C Argument variables 482C 483 DOUBLE PRECISION a, b(0:nbp1,0:nbp1) 484 DOUBLE PRECISION c(0:nbp1,0:nbp1), w(0:nbp1,0:nbp1) 485 DOUBLE PRECISION y(0:nbp1,0:nbp1) 486C 487C Common variables 488C 489 DOUBLE PRECISION d(0:nbp1,0:nbp1), l(0:nbp1,0:nbp1) 490 DOUBLE PRECISION r(0:nbp1,0:nbp1), x(0:nbp1,0:nbp1) 491C 492 COMMON / tdspace/ x, l, r, d 493C 494C====>End Module TRIDGX File t208x.f 495C====>Begin Module TRIDGY File t208x.f 496C---->Makedcls Options: All variables 497C 498C Parameter variables 499C 500 INTEGER l2ng 501 PARAMETER (l2ng = 6) 502 INTEGER ng 503 PARAMETER (ng = 2**l2ng) 504 INTEGER nb 505 PARAMETER (nb = ng+2) 506 INTEGER ngm1 507 PARAMETER (ngm1 = ng-1) 508 INTEGER n 509 PARAMETER (n = ng-1) 510 INTEGER np1 511 PARAMETER (np1 = n+1) 512 INTEGER nm1 513 PARAMETER (nm1 = n-1) 514 INTEGER nbp1 515 PARAMETER (nbp1 = nb+1) 516C 517C Argument variables 518C 519 DOUBLE PRECISION a, b(0:nbp1,0:nbp1) 520 DOUBLE PRECISION c(0:nbp1,0:nbp1), w(0:nbp1,0:nbp1) 521 DOUBLE PRECISION y(0:nbp1,0:nbp1) 522C 523C Common variables 524C 525 DOUBLE PRECISION d(0:nbp1,0:nbp1), l(0:nbp1,0:nbp1) 526 DOUBLE PRECISION r(0:nbp1,0:nbp1), x(0:nbp1,0:nbp1) 527C 528 COMMON / tdspace/ x, l, r, d 529C 530C====>End Module TRIDGY File t208x.f 531C====>Begin Module NEWU File t208x.f 532C---->Makedcls Options: All variables 533C 534C Parameter variables 535C 536 INTEGER l2ng 537 PARAMETER (l2ng = 6) 538 INTEGER ng 539 PARAMETER (ng = 2**l2ng) 540 INTEGER nb 541 PARAMETER (nb = ng+2) 542 INTEGER ngm1 543 PARAMETER (ngm1 = ng-1) 544 INTEGER nbp1 545 PARAMETER (nbp1 = nb+1) 546C 547C Argument variables 548C 549 DOUBLE PRECISION p(0:nbp1,0:nbp1), w2(0:nbp1,0:nbp1,2) 550C 551C Local variables 552C 553 DOUBLE PRECISION s 554C 555C Common variables 556C 557 DOUBLE PRECISION h, k 558C 559 COMMON / steps / h, k 560C 561C Common variables 562C 563 DOUBLE PRECISION u(0:nbp1,0:nbp1,2), ux(0:nbp1,0:nbp1) 564 DOUBLE PRECISION uy(0:nbp1,0:nbp1), vx(0:nbp1,0:nbp1) 565 DOUBLE PRECISION vy(0:nbp1,0:nbp1) 566C 567 COMMON / vel / u, uy, vy, ux 568 COMMON / vel / vx 569C 570C====>End Module NEWU File t208x.f 571C====>Begin Module DIV File t208x.f 572C---->Makedcls Options: All variables 573C 574C Parameter variables 575C 576 INTEGER l2ng 577 PARAMETER (l2ng = 6) 578 INTEGER ng 579 PARAMETER (ng = 2**l2ng) 580 INTEGER nb 581 PARAMETER (nb = ng+2) 582 INTEGER ngm1 583 PARAMETER (ngm1 = ng-1) 584 INTEGER nbp1 585 PARAMETER (nbp1 = nb+1) 586C 587C Argument variables 588C 589 DOUBLE PRECISION d(0:nbp1,0:nbp1), s3 590 DOUBLE PRECISION w2(0:nbp1,0:nbp1,2) 591C 592C====>End Module DIV File t208x.f 593C====>Begin Module FFT2D File t208x.f 594C---->Makedcls Options: All variables 595C 596C Parameter variables 597C 598 INTEGER l2ng 599 PARAMETER (l2ng = 6) 600 INTEGER ng 601 PARAMETER (ng = 2**l2ng) 602 INTEGER nb 603 PARAMETER (nb = ng+2) 604 INTEGER ngm1 605 PARAMETER (ngm1 = ng-1) 606 INTEGER n 607 PARAMETER (n = ng) 608 INTEGER m 609 PARAMETER (m = l2ng) 610 INTEGER nbp1 611 PARAMETER (nbp1 = nb+1) 612C 613C Argument variables 614C 615 DOUBLE PRECISION a(0:nbp1,0:nbp1), b(0:nbp1,0:nbp1) 616C 617 INTEGER isign 618C 619C Local variables 620C 621 DOUBLE PRECISION ang, pi, ssign, t1(n) 622 DOUBLE PRECISION t2(n), t3(n), t4(n), tu1 623 DOUBLE PRECISION u1, u2, w1, w2 624C 625C====>End Module FFT2D File t208x.f 626C====>Begin Module PHIA2D File t208x.f 627C---->Makedcls Options: All variables 628C 629C Parameter variables 630C 631 INTEGER l2ng 632 PARAMETER (l2ng = 6) 633 INTEGER ng 634 PARAMETER (ng = 2**l2ng) 635 INTEGER nb 636 PARAMETER (nb = ng+2) 637 INTEGER ngm1 638 PARAMETER (ngm1 = ng-1) 639 INTEGER ngp1 640 PARAMETER (ngp1 = ng+1) 641 INTEGER ngp2 642 PARAMETER (ngp2 = ng+2) 643 INTEGER nbp1 644 PARAMETER (nbp1 = nb+1) 645C 646C Argument variables 647C 648 INTEGER n 649C 650C Common variables 651C 652 DOUBLE PRECISION phia(0:nbp1,0:nbp1), phian(0:nbp1,0:nbp1) 653 DOUBLE PRECISION phin(0:nbp1,0:nbp1), phinn(0:nbp1,0:nbp1) 654C 655 COMMON / phi / phin, phinn, phia, phian 656C 657C Common variables 658C 659 DOUBLE PRECISION h, k 660C 661 COMMON / steps / h, k 662C 663C Common variables 664C 665 DOUBLE PRECISION u(0:nbp1,0:nbp1,2), ux(0:nbp1,0:nbp1) 666 DOUBLE PRECISION uy(0:nbp1,0:nbp1), vx(0:nbp1,0:nbp1) 667 DOUBLE PRECISION vy(0:nbp1,0:nbp1) 668C 669 COMMON / vel / u, uy, vy, ux 670 COMMON / vel / vx 671C 672C====>End Module PHIA2D File t208x.f 673C====>Begin Module ZAG2D File t208x.f 674C---->Makedcls Options: All variables 675C 676C Parameter variables 677C 678 INTEGER l2ng 679 PARAMETER (l2ng = 6) 680 INTEGER ng 681 PARAMETER (ng = 2**l2ng) 682 INTEGER nb 683 PARAMETER (nb = ng+2) 684 INTEGER ngm1 685 PARAMETER (ngm1 = ng-1) 686 INTEGER ngp1 687 PARAMETER (ngp1 = ng+1) 688 INTEGER ngp2 689 PARAMETER (ngp2 = ng+2) 690 INTEGER nbp1 691 PARAMETER (nbp1 = nb+1) 692C 693C Argument variables 694C 695 INTEGER n 696C 697C Common variables 698C 699 DOUBLE PRECISION zag(0:nbp1,0:nbp1), zagn(0:nbp1,0:nbp1) 700C 701 COMMON / agg / zag, zagn 702C 703C Common variables 704C 705 DOUBLE PRECISION h, k 706C 707 COMMON / steps / h, k 708C 709C Common variables 710C 711 DOUBLE PRECISION u(0:nbp1,0:nbp1,2), ux(0:nbp1,0:nbp1) 712 DOUBLE PRECISION uy(0:nbp1,0:nbp1), vx(0:nbp1,0:nbp1) 713 DOUBLE PRECISION vy(0:nbp1,0:nbp1) 714C 715 COMMON / vel / u, uy, vy, ux 716 COMMON / vel / vx 717C 718C====>End Module ZAG2D File t208x.f 719C====>Begin Module PLOT File t208x.f 720C---->Makedcls Options: All variables 721C 722C Parameter variables 723C 724 INTEGER l2ng 725 PARAMETER (l2ng = 6) 726 INTEGER ng 727 PARAMETER (ng = 2**l2ng) 728 INTEGER nb 729 PARAMETER (nb = ng+2) 730 INTEGER ngm1 731 PARAMETER (ngm1 = ng-1) 732 INTEGER ngp1 733 PARAMETER (ngp1 = ng+1) 734 INTEGER ngp2 735 PARAMETER (ngp2 = ng+2) 736 INTEGER nbp1 737 PARAMETER (nbp1 = nb+1) 738 INTEGER nr 739 PARAMETER (nr = 32) 740 INTEGER nt 741 PARAMETER (nt = 16) 742 INTEGER nr2 743 PARAMETER (nr2 = nr+2) 744 INTEGER nt1 745 PARAMETER (nt1 = nt+1) 746 INTEGER mmax 747 PARAMETER (mmax = 200) 748C 749C Argument variables 750C 751 DOUBLE PRECISION c(0:nbp1,0:nbp1) 752 DOUBLE PRECISION elf(0:nbp1,0:nbp1,0:nr2,nt1) 753 DOUBLE PRECISION f(0:nbp1,0:nbp1,2), phia(0:nbp1,0:nbp1) 754 DOUBLE PRECISION phin(0:nbp1,0:nbp1), time 755 DOUBLE PRECISION u(0:nbp1,0:nbp1,2), xm(mmax,2) 756 DOUBLE PRECISION zag(0:nbp1,0:nbp1) 757C 758 INTEGER n, nk 759C 760C Common variables 761C 762 INTEGER nfil 763C 764 COMMON / efnum / nfil 765C 766C====>End Module PLOT File t208x.f 767C====>Begin Module MARKERS File t208x.f 768C---->Makedcls Options: All variables 769C 770C Parameter variables 771C 772 INTEGER l2ng 773 PARAMETER (l2ng = 6) 774 INTEGER ng 775 PARAMETER (ng = 2**l2ng) 776 INTEGER nb 777 PARAMETER (nb = ng+2) 778 INTEGER ngm1 779 PARAMETER (ngm1 = ng-1) 780 INTEGER ngp1 781 PARAMETER (ngp1 = ng+1) 782 INTEGER ngp2 783 PARAMETER (ngp2 = ng+2) 784 INTEGER nbp1 785 PARAMETER (nbp1 = nb+1) 786 INTEGER mmax 787 PARAMETER (mmax = 200) 788C 789C Local variables 790C 791 DOUBLE PRECISION ax, ay, um, umn 792 DOUBLE PRECISION vm, vmn, x, xn 793 DOUBLE PRECISION y, yn 794C 795C Common variables 796C 797 DOUBLE PRECISION xm(mmax,2) 798C 799 COMMON / fmarkers/ xm 800C 801C Common variables 802C 803 DOUBLE PRECISION dlt, h 804C 805 COMMON / steps / h, dlt 806C 807C Common variables 808C 809 DOUBLE PRECISION u(0:nbp1,0:nbp1,2), ux(0:nbp1,0:nbp1) 810 DOUBLE PRECISION uy(0:nbp1,0:nbp1), vx(0:nbp1,0:nbp1) 811 DOUBLE PRECISION vy(0:nbp1,0:nbp1) 812C 813 COMMON / vel / u, uy, vy, ux 814 COMMON / vel / vx 815C 816C====>End Module MARKERS File t208x.f 817C====>Begin Module PHIN2D File t208x.f 818C---->Makedcls Options: All variables 819C 820C Parameter variables 821C 822 INTEGER l2ng 823 PARAMETER (l2ng = 6) 824 INTEGER ng 825 PARAMETER (ng = 2**l2ng) 826 INTEGER nb 827 PARAMETER (nb = ng+2) 828 INTEGER ngm1 829 PARAMETER (ngm1 = ng-1) 830 INTEGER ngp1 831 PARAMETER (ngp1 = ng+1) 832 INTEGER ngp2 833 PARAMETER (ngp2 = ng+2) 834 INTEGER nbp1 835 PARAMETER (nbp1 = nb+1) 836 INTEGER mmax 837 PARAMETER (mmax = 200) 838C 839C Argument variables 840C 841 INTEGER n 842C 843C Common variables 844C 845 DOUBLE PRECISION pen 846C 847 COMMON / diffn / pen 848C 849C Common variables 850C 851 DOUBLE PRECISION phia(0:nbp1,0:nbp1), phian(0:nbp1,0:nbp1) 852 DOUBLE PRECISION phin(0:nbp1,0:nbp1), phinn(0:nbp1,0:nbp1) 853C 854 COMMON / phi / phin, phinn, phia, phian 855C 856C Common variables 857C 858 DOUBLE PRECISION h, k 859C 860 COMMON / steps / h, k 861C 862C Common variables 863C 864 DOUBLE PRECISION u(0:nbp1,0:nbp1,2), ux(0:nbp1,0:nbp1) 865 DOUBLE PRECISION uy(0:nbp1,0:nbp1), vx(0:nbp1,0:nbp1) 866 DOUBLE PRECISION vy(0:nbp1,0:nbp1) 867C 868 COMMON / vel / u, uy, vy, ux 869 COMMON / vel / vx 870C 871C====>End Module PHIN2D File t208x.f 872C====>Begin Module CHEM2D File t208x.f 873C---->Makedcls Options: All variables 874C 875C Parameter variables 876C 877 INTEGER l2ng 878 PARAMETER (l2ng = 6) 879 INTEGER ng 880 PARAMETER (ng = 2**l2ng) 881 INTEGER nb 882 PARAMETER (nb = ng+2) 883 INTEGER ngm1 884 PARAMETER (ngm1 = ng-1) 885 INTEGER ngp1 886 PARAMETER (ngp1 = ng+1) 887 INTEGER ngp2 888 PARAMETER (ngp2 = ng+2) 889 INTEGER nbp1 890 PARAMETER (nbp1 = nb+1) 891 INTEGER mmax 892 PARAMETER (mmax = 200) 893C 894C Argument variables 895C 896 INTEGER n 897C 898C Common variables 899C 900 DOUBLE PRECISION c(0:nbp1,0:nbp1), cn(0:nbp1,0:nbp1) 901C 902 COMMON / adp / c, cn 903C 904C Common variables 905C 906 DOUBLE PRECISION pec 907C 908 COMMON / diffc / pec 909C 910C Common variables 911C 912 DOUBLE PRECISION h, k 913C 914 COMMON / steps / h, k 915C 916C Common variables 917C 918 DOUBLE PRECISION u(0:nbp1,0:nbp1,2), ux(0:nbp1,0:nbp1) 919 DOUBLE PRECISION uy(0:nbp1,0:nbp1), vx(0:nbp1,0:nbp1) 920 DOUBLE PRECISION vy(0:nbp1,0:nbp1) 921C 922 COMMON / vel / u, uy, vy, ux 923 COMMON / vel / vx 924C 925C====>End Module CHEM2D File t208x.f 926C====>Begin Module CHLFAC File t208x.f 927C---->Makedcls Options: All variables 928C 929C Argument variables 930C 931 DOUBLE PRECISION diag(nn), subd(nn) 932C 933 INTEGER nn 934C 935C====>End Module CHLFAC File t208x.f 936C====>Begin Module CHLSLV File t208x.f 937C---->Makedcls Options: All variables 938C 939C Parameter variables 940C 941 INTEGER l2ng 942 PARAMETER (l2ng = 6) 943 INTEGER ng 944 PARAMETER (ng = 2**l2ng) 945 INTEGER nb 946 PARAMETER (nb = ng+2) 947 INTEGER ngm1 948 PARAMETER (ngm1 = ng-1) 949 INTEGER ngp1 950 PARAMETER (ngp1 = ng+1) 951 INTEGER ngp2 952 PARAMETER (ngp2 = ng+2) 953 INTEGER nbp1 954 PARAMETER (nbp1 = nb+1) 955 INTEGER max 956 PARAMETER (max = nbp1) 957C 958C Argument variables 959C 960 DOUBLE PRECISION b(nm1+1), diag(nm1), sub(nm1), x(nm1) 961C 962 INTEGER nm1 963C 964C Local variables 965C 966 DOUBLE PRECISION y(max) 967C 968C====>End Module CHLSLV File t208x.f 969C====>Begin Module PERIOD File t208x.f 970C---->Makedcls Options: All variables 971C 972C Parameter variables 973C 974 INTEGER l2ng 975 PARAMETER (l2ng = 6) 976 INTEGER ng 977 PARAMETER (ng = 2**l2ng) 978 INTEGER nb 979 PARAMETER (nb = ng+2) 980 INTEGER ngm1 981 PARAMETER (ngm1 = ng-1) 982 INTEGER ngp1 983 PARAMETER (ngp1 = ng+1) 984 INTEGER ngp2 985 PARAMETER (ngp2 = ng+2) 986 INTEGER nbp1 987 PARAMETER (nbp1 = nb+1) 988C 989C Argument variables 990C 991 DOUBLE PRECISION ph(0:nbp1,0:nbp1) 992C 993C====>End Module PERIOD File t208x.f 994C====>Begin Module SOL File t208x.f 995C---->Makedcls Options: All variables 996C 997C Argument variables 998C 999 DOUBLE PRECISION l, r, v 1000C 1001C====>End Module SOL File t208x.f 1002C====>Begin Module DIFFU File t208x.f 1003C---->Makedcls Options: All variables 1004C 1005C Parameter variables 1006C 1007 INTEGER l2ng 1008 PARAMETER (l2ng = 6) 1009 INTEGER ng 1010 PARAMETER (ng = 2**l2ng) 1011 INTEGER nb 1012 PARAMETER (nb = ng+2) 1013 INTEGER ngm1 1014 PARAMETER (ngm1 = ng-1) 1015 INTEGER ngp1 1016 PARAMETER (ngp1 = ng+1) 1017 INTEGER ngp2 1018 PARAMETER (ngp2 = ng+2) 1019 INTEGER nbp1 1020 PARAMETER (nbp1 = nb+1) 1021C 1022C Argument variables 1023C 1024 DOUBLE PRECISION alpha, k1, ph(0:nbp1,0:nbp1) 1025 DOUBLE PRECISION phn(0:nbp1,0:nbp1) 1026C 1027C Local variables 1028C 1029 DOUBLE PRECISION b(nb), eta, gdiag(nb-1) 1030 DOUBLE PRECISION gsub(nb-1), lam, y(nb-1), z(nb-1) 1031C 1032C Common variables 1033C 1034 DOUBLE PRECISION h, k 1035C 1036 COMMON / steps / h, k 1037C 1038C====>End Module DIFFU File t208x.f 1039C====>Begin Module REACT File t208x.f 1040C---->Makedcls Options: All variables 1041C 1042C Parameter variables 1043C 1044 INTEGER l2ng 1045 PARAMETER (l2ng = 6) 1046 INTEGER ng 1047 PARAMETER (ng = 2**l2ng) 1048 INTEGER nb 1049 PARAMETER (nb = ng+2) 1050 INTEGER ngm1 1051 PARAMETER (ngm1 = ng-1) 1052 INTEGER ngp1 1053 PARAMETER (ngp1 = ng+1) 1054 INTEGER ngp2 1055 PARAMETER (ngp2 = ng+2) 1056 INTEGER nbp1 1057 PARAMETER (nbp1 = nb+1) 1058 INTEGER nr 1059 PARAMETER (nr = 32) 1060 INTEGER nt 1061 PARAMETER (nt = 16) 1062 INTEGER nr2 1063 PARAMETER (nr2 = nr+2) 1064 INTEGER nt1 1065 PARAMETER (nt1 = nt+1) 1066C 1067C Local variables 1068C 1069 DOUBLE PRECISION alpha, c4, rc(0:nbp1), temp1 1070 DOUBLE PRECISION temp2, temp3, temp4 1071 DOUBLE PRECISION y(0:nbp1,4), yn(0:nbp1,4) 1072C 1073C Common variables 1074C 1075 DOUBLE PRECISION c(0:nbp1,0:nbp1), cn(0:nbp1,0:nbp1) 1076C 1077 COMMON / adp / c, cn 1078C 1079C Common variables 1080C 1081 DOUBLE PRECISION zag(0:nbp1,0:nbp1), zagn(0:nbp1,0:nbp1) 1082C 1083 COMMON / agg / zag, zagn 1084C 1085C Common variables 1086C 1087 DOUBLE PRECISION cnd1, cnd2, cnd3, cnd4 1088 DOUBLE PRECISION cnd5 1089C 1090 COMMON / cnd / cnd1, cnd2, cnd3, cnd4 1091 COMMON / cnd / cnd5 1092C 1093C Common variables 1094C 1095 DOUBLE PRECISION elf(0:nbp1,0:nbp1,0:nr2,nt1) 1096 DOUBLE PRECISION elfn(0:nbp1,0:nbp1,0:nr2,nt1) 1097C 1098 COMMON / link / elf, elfn 1099C 1100C Common variables 1101C 1102 DOUBLE PRECISION phia(0:nbp1,0:nbp1), phian(0:nbp1,0:nbp1) 1103 DOUBLE PRECISION phin(0:nbp1,0:nbp1), phinn(0:nbp1,0:nbp1) 1104C 1105 COMMON / phi / phin, phinn, phia, phian 1106C 1107C Common variables 1108C 1109 DOUBLE PRECISION dr(nr2), r(0:nr2) 1110C 1111 INTEGER nk 1112C 1113 COMMON / rsize / dr, r, nk 1114C 1115C Common variables 1116C 1117 DOUBLE PRECISION h, k 1118C 1119 COMMON / steps / h, k 1120C 1121C====>End Module REACT File t208x.f 1122