1*9371c9d4SSatish Balay integer probnum, ihod, mx, my, neq, ientro, gorder parameter(neq = 3) double precision theta, pi, time, zero double precision damfac double precision dt, dtmin, dtmax, dtgrow, tfinal, tplot, tcscal, 2*9371c9d4SSatish Balay hcscal 3c4762a1bSJed Brown 4*9371c9d4SSatish Balay double precision eigval, 5*9371c9d4SSatish Balay eigvec, rinv, roestt, fl, fr, deltau, alpha, xnumdif, 6*9371c9d4SSatish Balay froe 7c4762a1bSJed Brown 8*9371c9d4SSatish Balay double precision dx, 9*9371c9d4SSatish Balay xl0, kappa0, kappaa, kappab, visc0, 10*9371c9d4SSatish Balay erg0 11c4762a1bSJed Brown 12*9371c9d4SSatish Balay logical debug, 13*9371c9d4SSatish Balay dampit, wilson, dtcon, 14*9371c9d4SSatish Balay pcnew 15c4762a1bSJed Brown 16*9371c9d4SSatish Balay common 17*9371c9d4SSatish Balay / params / mx, 18*9371c9d4SSatish Balay my, probnum, ihod, ientro, 19*9371c9d4SSatish Balay gorder 20c4762a1bSJed Brown 21*9371c9d4SSatish Balay common 22*9371c9d4SSatish Balay / func / pi, 23*9371c9d4SSatish Balay zero, theta, dx, xl0, damfac, kappa0, kappaa, kappab, visc0, 24*9371c9d4SSatish Balay erg0 25c4762a1bSJed Brown 26*9371c9d4SSatish Balay common 27*9371c9d4SSatish Balay / gudnov / eigval(neq), 28*9371c9d4SSatish Balay eigvec(neq, neq), rinv(neq, neq), roestt(neq), fl(neq), fr(neq), deltau(neq), alpha(neq), xnumdif(neq), 29*9371c9d4SSatish Balay froe(neq) 30c4762a1bSJed Brown 31*9371c9d4SSatish Balay common 32*9371c9d4SSatish Balay / flags / debug, 33*9371c9d4SSatish Balay dampit, wilson, dtcon, 34*9371c9d4SSatish Balay pcnew 35c4762a1bSJed Brown 36*9371c9d4SSatish Balay common 37*9371c9d4SSatish Balay / timcnt / time, 38*9371c9d4SSatish Balay dt, dtmin, dtmax, dtgrow, tfinal, tplot, tcscal, hcscal 39