TEXT(YX BFC UNIFORM FLOW ACROSS SKEWED BOX
TITLE
  DISPLAY
  This problem concerns uniform vertical flow of a scalar across
  a skewed box. The scalar is convected through the south boundary
  only. The MUSCL scheme is used for the momentum equations, and
  the QUICK scheme for the scalar concentration. This case provides
  a simple test of the higher-order schemes on non-orthogonal
  meshes.
  ENDDIS
  PHOTON USE
  P
 
 
 
   0.20443E+04 0.15633E+04 CR
  GR Z 1;VEC Z 1 SH
  msg velocity vectors & mesh
  msg press  to continue
  pause
  CL;gr ou z 1;CON C1 Z 1 FI;.01
  msg QUICK scheme for scalar concentration
  msg press  to continue
  pause
  ENDUSE
REAL(VTOT,YLEN,RELX);VTOT=1.0;YLEN=1.0
BFC=T
GSET(P,P01,0.0000,0.0000,0.0000)
GSET(P,P02,1.0000,0.0000,0.0000)
GSET(P,P03,1.0000,1.0000,0.0000)
GSET(P,P04,0.0000,1.0000,0.0000)
GSET(L,L01,P01,P02,20,1.0)
GSET(L,L02,P02,P03,20,1.0)
GSET(L,L03,P03,P04,20,1.0)
GSET(L,L04,P04,P01,20,1.0)
GSET(F,F01,P01,-,P02,-,P03,-,P04,-)
GSET(D,20,20,1,1.0000,1.0000,0.1000)
GSET(M,F01,+I+J,1,1,1,TRANS)
GSET(P,P01,0.0000,0.0000,0.0000)
GSET(P,P02,1.0000,0.0000,0.0000)
GSET(P,P03,1.7071,0.7071,0.0000)
GSET(P,P04,0.7071,0.7071,0.0000)
GSET(L,L01,P01,P02,20,1.0)
GSET(L,L02,P02,P03,20,1.0)
GSET(L,L03,P03,P04,20,1.0)
GSET(L,L04,P04,P01,20,1.0)
GSET(F,F01,P01,-,P02,-,P03,-,P04,-)
GSET(D,20,20,1,1.0000,1.0000,0.1000)
GSET(M,F01,+I+J,1,1,1,TRANS)
GSET(C,K2,F,K1,1,20,1,20,+,0,0,1.0000E-02)
GSET(C,K:NZ+1:,F,K1,1,NX,1,NY,+,0.0,0.0,0.1,INC,1.0)
NONORT=T
   GROUP 7.  Variables named, stored & solved.
SOLVE(P1,U1,V1);SOLUTN(P1,Y,Y,Y,N,N,N);SOLVE(C1)
SCHEME(MUSCL,U1,V1);SCHEME(QUICK,C1)
   GROUP 9.  Properties of the medium (or media).
RHO1=1.161;ENUL=1.5890E-05;ENUT=0
   GROUP 11. Initialization of fields of variables, etc.
FIINIT(U1)=0.7071;FIINIT(V1)=0.7071
FIINIT(UCRT)=0.0;FIINIT(VCRT)=1.0
   GROUP 12. Convection and diffusion adjustments
   GROUP 13. Boundary conditions and special sources
   * INLET boundary condition, name INB
INLET(BFCINB,SOUTH,#1,#1,#1,#1,#1,#1,#1,#NREGT)
VALUE(BFCINB,U1,GRND1);VALUE(BFCINB,UCRT,0.0)
VALUE(BFCINB,V1,GRND1);VALUE(BFCINB,VCRT,1.0)
VALUE(BFCINB,P1,GRND1);VALUE(BFCINB,C1,1);BFCA=RHO1
   * INLET boundary condition, name INS
INLET(BFCINS,EAST,#1,#1,#1,#1,#1,#1,#1,#NREGT)
VALUE(BFCINS,U1,GRND1);VALUE(BFCINS,UCRT,0.0)
VALUE(BFCINS,V1,GRND1);VALUE(BFCINS,VCRT,1.0)
VALUE(BFCINS,P1,GRND1);VALUE(BFCINS,C1,0)
   * OUTLET boundary condition, name OUTT
OUTLET(OUTT,NORTH,#1,#1,#1,#1,#1,#1,#1,#NREGT)
VALUE(OUTT,P1,0);VALUE(OUTT,C1,SAME)
   * OUTLET boundary condition, name OUTS
OUTLET(OUTS,WEST,#1,#1,#1,#1,#1,#1,#1,#NREGT)
VALUE(OUTS,P1,0);VALUE(OUTS,C1,SAME)
   GROUP 15. Termination criteria for sweeps and outer iterations.
LSWEEP=140
   GROUP 16. Termination criteria for inner iterations.
   GROUP 17. Under-relaxation and related devices.
RELX=YLEN/VTOT/NX;RELAX(P1,LINRLX,0.3)
RELAX(U1,FALSDT,RELX);RELAX(V1,FALSDT,RELX);RELAX(C1,FALSDT,RELX)
   GROUP 22. Location of spot-value & residual printout.
IXMON=11;IYMON=11;TSTSWP=-1
   GROUP 23. Variable-by-variable field printout and plots
ITABL=3;NPLT=10