Q1 files for Couette-flow studies of two-fluid turbulence
    model.
    Run no.  1    2    3     3     3   6    7     8     9     10
    BUOYB   100. 200. 500. 1.E3. 1.E4  0. -100. -200. -500. -1.E3
    GROUP 1. Run identifiers and other preliminaries.
BUOYB=100.0;BUOYC=-BUOYB
TEXT(RUN1  COUETTE  BUOY CONST=1.E2    : 975
TITLE
REAL(CHARW,CHARL,CHART,REYNO)
CHARW=1.0;CHARL=0.01;REYNO=1.E5;CHART=0.2*CHARL/CHARW
    GROUP 2. Transience; time-step specification
STEADY=F;GRDPWR(T,1,1.E3*CHART,1.0)
    GROUP 4. y-direction grid specification.
GRDPWR(Y,-30,CHARL,1.5)
    GROUP 7. Variables (including porosities) named, stored & solved
ONEPHS=F;SOLVE(P1,V1,V2,W1,W2,R1,R2,H1,H2)
    SOLVE(C1,C2,C3,C4,C5,C6)
SOLUTN(C1,Y,Y,P,P,P,P);SOLUTN(C2,Y,Y,P,P,P,P)
SOLUTN(C3,Y,Y,P,P,P,P);SOLUTN(C4,Y,Y,P,P,P,P)
SOLUTN(C5,Y,Y,P,P,P,P);SOLUTN(C6,Y,Y,P,P,P,P)
INTMDT=22;LEN1=23;VIST=24;NAME(INTMDT)=MDOT;NAME(LEN1)=LEN
NAME(INTMDT)=MDOT;NAME(LEN1)=LEN;NAME(VIST)=VIS
SOLUTN(MDOT,Y,N,N,N,N,N);SOLUTN(VIST,Y,N,N,N,N,N)
SOLUTN(LEN1,Y,N,N,N,N,N)
    GROUP 8. Terms (in differential equations) and devices.
TERMS(H1,N,P,P,P,P,P);TERMS(H1,N,P,P,P,P,P)
TERMS(W1,N,Y,Y,Y,Y,Y);TERMS(W2,N,Y,Y,Y,N,Y)
TERMS(C3,P,P,P,P,P,N);TERMS(C4,P,P,P,P,P,N)
    GROUP 9. Properties of the medium (or media).
ENUL=CHARW*CHARL/REYNO;EL1=NIKURCH;ENUT=2FLUID;ENUTA=1.0
PRT(R1)=1.0;PRT(R2)=1.0
    GROUP 10. Interphase-transfer processes and properties.
CFIPS=GRND4;CFIPA=0.0;CFIPB=1.0;CFIPD=-1.0;CFIPC=0.05
CMDOT=GRND1;CMDTA=10.0;CMDTB=0.5;CMDTC=0.0
CINT(C1)=10.0;CINT(C2)=10.0;CINT(C5)=0.1;CINT(C6)=0.1
CINT(C3)=0.0;CINT(C4)=0.0
    GROUP 11. Initial fields of variables, porosities, etc
FIINIT(P1)=0.0;FIINIT(W1)=CHARW;FIINIT(W2)=CHARW
FIINIT(R1)=0.5;FIINIT(R2)=0.5; FIINIT(H1)=0.5;FIINIT(H2)=0.5
FIINIT(V1)=0.01*CHARW;FIINIT(V2)=-0.01*CHARW
FIINIT(C1)=0.5;FIINIT(C2)=0.5;FIINIT(C5)=0.5;FIINIT(C6)=0.5
FIINIT(C3)=0.05*CHARL;FIINIT(C4)=0.05*CHARL
    GROUP 13. Boundary conditions, and special sources.
    Pressure relief
PATCH(RELIEF,CELL,1,1,NY/2,NY/2,1,1,1,1)
COVAL(RELIEF,P1,FIXVAL,0.0)
    SOUTH WALL
REAL(COEF);COEF=1.0
PATCH(FIXED,SOUTH,1,1,1,1,1,NZ,1,1)
COVAL(FIXED,W1,FIXVAL,0.0);COVAL(FIXED,H1,1.0,1.0)
COVAL(FIXED,C1,COEF,1.0);COVAL(FIXED,C5,COEF,1.0)
COVAL(FIXED,C3,COEF,0.0);COVAL(FIXED,C4,COEF,0.0)
    NORTH WALL
PATCH(MOVING,NORTH,1,1,NY,NY,1,NZ,1,1)
COVAL(MOVING,W2,FIXVAL,2.0*CHARW);COVAL(MOVING,H2,1.0,0.0)
COVAL(MOVING,C2,COEF,0.0);COVAL(MOVING,C6,COEF,0.0)
COVAL(MOVING,C3,COEF,0.0);COVAL(MOVING,C4,COEF,0.0)
    WHOLE FIELD
PATCH(SHSOURCE,CELL,1,1,1,NY,1,NZ,1,1)
COVAL(SHSOURCE,V1,FIXFLU,GRND5);COVAL(SHSOURCE,V2,FIXFLU,GRND5)
SHSOA=1.E0
    BODY FORCE
PATCH(BUOY,PHASEM,1,1,1,NY,1,NZ,1,1)
COVAL(BUOY,V1,FIXFLU,LINBC)
IBUOYB=14;IBUOYC=15;CSG2=BUOY;BUOYA=0.0;BUOYD=0.025
    LENGTH-SCALE SOURCE
PATCH(LESO,PHASEM,1,1,1,NY,1,NZ,1,1)
COVAL(LESO,C3,FIXFLU,GRND1);COVAL(LESO,C4,FIXFLU,GRND1)
ELSOA=0.025
    GROUP 15. Termination criteria for sweeps and outer iterations.
LSWEEP=200
RESREF(P1)=1.E-8;RESREF(R1)=1.E-8;RESREF(R2)=1.E-8
RESREF(V1)=1.E-8;RESREF(W1)=1.E-8;RESREF(V2)=1.E-8
RESREF(W2)=1.E-8;RESREF(H1)=1.E-8;RESREF(H2)=1.E-8
RESREF(C5)=1.E-8;RESREF(C6)=1.E-8
LITER(V1)=2;LITER(V2)=2;LITER(W1)=2;LITER(W2)=2
    GROUP 17. Under-relaxation and related devices.
RELAX(V1,FALSDT,10.0*CHART);RELAX(V2,FALSDT,10.0*CHART)
RELAX(W1,FALSDT,10.0*CHART);RELAX(W2,FALSDT,10.0*CHART)
RELAX(H1,FALSDT,10.0*CHART);RELAX(H2,FALSDT,10.0*CHART)
RELAX(R1,LINRLX,0.4);RELAX(R2,LINRLX,0.4)
SPEDAT(SET,GXMONI,TRANSIENT,L,F)
    GROUP 21. Frequency and extent of field printout.
OUTPUT(INTMDT,N,Y,Y,Y,Y,Y);OUTPUT(VIS,N,Y,Y,Y,Y,Y)
OUTPUT(R1,N,N,N,N,N,N);OUTPUT(R2,N,N,N,N,N,N)
OUTPUT(P1,N,N,N,N,N,N);OUTPUT(V1,N,P,P,P,P,P)
OUTPUT(V2,N,P,P,P,P,P);OUTPUT(W1,N,P,P,P,P,P)
OUTPUT(W2,N,P,P,P,P,P);OUTPUT(H1,N,P,P,P,P,P)
OUTPUT(H2,N,P,P,P,P,P);OUTPUT(C1,N,P,P,P,P,P)
OUTPUT(C2,N,P,P,P,P,P);OUTPUT(C3,N,P,P,P,P,P)
OUTPUT(C4,N,P,P,P,P,P);OUTPUT(C5,N,P,P,P,P,P)
OUTPUT(C6,N,P,P,P,P,P)
    GROUP 22. Spot-value & residual printout.
TSTSWP=1;IYMON=NY/2
    GROUP 23. Field printout and plot control
NYPRIN=2;IPLTL=LSWEEP;IPROF=3;ORSIZ=0.4
PATCH(OUTPUT1,PROFIL,1,1,1,NY,NZ,NZ,1,1)
COVAL(OUTPUT1,V1,0.0,CHARW);COVAL(OUTPUT1,V2,-CHARW,0.0)
COVAL(OUTPUT1,W1,0.0,2.0*CHARW);COVAL(OUTPUT1,W2,0.0,2.0*CHARW)
COVAL(OUTPUT1,V1,-2.0,-2.0);COVAL(OUTPUT1,V2,-2.0,-2.0)
COVAL(OUTPUT1,W1,-1.0,-1.0);COVAL(OUTPUT1,W2,-1.0,-1.0)
COVAL(OUTPUT1,R1,0.0,1.0);COVAL(OUTPUT1,INTMDT,0.0,0.0)
PATCH(OUTPUT2,PROFIL,1,1,1,NY,NZ,NZ,1,1)
COVAL(OUTPUT2,LEN1,-0.1,-0.1);COVAL(OUTPUT2,VIS,-0.0,-0.0)
COVAL(OUTPUT2,C3,-0.1,-0.1);COVAL(OUTPUT2,C4,-0.1,-0.1)
 
PATCH(OUTPUT3,PROFIL,1,1,1,NY,NZ,NZ,1,1)
COVAL(OUTPUT3,H1,0.0,1.0);COVAL(OUTPUT3,H2,0.0,1.0)
COVAL(OUTPUT3,C1,0.0,1.0);COVAL(OUTPUT3,C2,0.0,1.0)
COVAL(OUTPUT3,C5,0.0,1.0);COVAL(OUTPUT3,C6,0.0,1.0)
NAME(V1)=A.V1;NAME(V2)=B.V2;NAME(W1)=C.W1;NAME(W2)=D.W2
NAME(H1)=E.H1;NAME(H2)=F.H2;NAME(R1)=UP;NAME(R2)=DOWN
NAME(C1)=G.C1;NAME(C2)=H.C2;NAME(C3)=I.C3;NAME(C4)=J.C4
NAME(C5)=K.C5;NAME(C6)=N.C6
    GROUP 24. Preparations for continuation runs.