C     Specification statements for streamf3
      INTEGER JMAX,JM1,JM2,JMAX2,NSTEP
      PARAMETER(JMAX=41, JM1=JMAX-1, JM2=JMAX-2, JMAX2=2*JM2*JM2
     *          ,JM3=JMAX-3, NSTEP=4)
      INTEGER i,j,k,IA,NN,icase,ibc,nframes,nframe
      DOUBLE PRECISION
     *      a(JMAX,JMAX,2,2),b(JMAX,JMAX,2,2),c(JMAX,JMAX,2,2)
     *     ,d(JMAX,JMAX,2,2),e(JMAX,JMAX,2,2),f(JMAX,JMAX,2)
     *     ,s(JMAX,JMAX),w(JMAX,JMAX),u(JMAX,JMAX),v(JMAX,JMAX)
     *     ,p(JMAX,JMAX),tau(JMAX,JMAX),r(JMAX)
     *     ,sbx(2,JMAX),sby(JMAX,2),dw(JMAX,JMAX)
     *     ,vtbx(2,JMAX),vtby(JMAX,2)
     *     ,vnbx(2,JMAX),vnby(JMAX,2)
     *     ,del,alf,del2,alf2,Re,t,dt,dtinit,dwmax,dwm,tmax,tim(40),dtold,dtmax
     *     ,uim1,uip1,vjm1,vjp1
     *     ,eps,beta,gamma,gam2
C     Specification statements for routine linbcg below
C     bA & xA are the right side & unknown of linear system
C     If you change NMAX here, you must also change it in the subroutines
      INTEGER NP,NP1,NMAX,ITOL,ITMAX
      DOUBLE PRECISION TOL
      PARAMETER(NP=JMAX2,NP1=NP/2,NMAX=50000,ITOL=3,TOL=1.d-9,ITMAX=300)
      INTEGER iter,ija,it
      DOUBLE PRECISION bA(NP),xA(NP),bcmp(NP),sa,err
      COMMON /mat/ sa(NMAX),ija(NMAX)
      COMMON i,j,k,IA,NN,icase,ibc,iter,it,nframes,nframe
     * ,a,b,c,d,e,f,s,w,u,v,p,tau,r,sbx,sby,dw,vtbx,vtby,vnbx,vnby
     * ,del,alf,del2,alf2,Re,t,dt,dtinit,dwmax,dwm,tmax,tim,dtmax
     * ,uim1,uip1,vjm1,vjp1
     * ,eps,beta,gamma,gam2