PHOTON USE
  p
  n4
  1 1 200
 
 
 
  msg(    film formation and vector field at t= 8s
  gr ou x 1
  con vfol x 1 sh
  0.45 0.55 40
  vec x 1
  pause
  p
  n8
  1 1 200
 
 
 
  msg(    film formation and vector field at t=16s
  gr ou x 1
  con vfol x 1 sh
  0.45 0.55 40
  vec x 1
  pause
  p
  n12
  1 1 200
 
 
 
  msg(    film formation and vector field at t=24s
  gr ou x 1
  con vfol x 1 sh
  0.45 0.55 40
  vec x 1
  pause
  p
  n16
  1 1 200
 
 
 
  msg(    film formation and vector field at t=32s
  gr ou x 1
  con vfol x 1 sh
  0.45 0.55 40
  vec x 1
  pause
  p
  n20
  1 1 200
 
 
 
  msg(    film formation and vector field at t=40s
  gr ou x 1
  con vfol x 1 sh
  0.45 0.55 40
  vec x 1
  msg            -
  msg Press e to END
  enduse
 
    GROUP 1. Run title and other preliminaries
TEXT(FILM FORMATION BY HOL :P101
TITLE
  DISPLAY
  FREE SURFACE FLOW -  Film Formation by HOL method (dip coating)
 
  2-dimensional (y-z), Cartesian, transient, elliptic simulation
 
  One wall of the domain is moving upwards, away from the liquid,
  and draws a film of liquid with it.
 
      moving  ^  |          |
        wall  |  |-         |
                 |- <- film |
                 |-- - - -  |     y ^
                 | - liquid |       |
                 |  - -- -  |       +----->
                 +----------+            z
  enddis
 
    GROUP 2. Transience; time-step specification
STEADY=F;LSTEP=20;TFRAC(1)=-LSTEP; TFRAC(2)=2.0
    GROUP 4. Y-direction grid specification
SUBGRD(Y, 1,20,1.000,1.0)
    GROUP 5. Z-direction grid specification
SUBGRD(Z,1,16,0.001,1.5);SUBGRD(Z,17,20,0.002,1.0)
    GROUP 7. Variables stored, solved & named
STORE(DEN1,PRPS);SOLVE(VFOL)
SOLUTN(P1,Y,Y,Y,N,N,N);SOLUTN(V1,Y,Y,N,N,N,N)
SOLUTN(W1,Y,Y,N,N,N,N)
    GROUP 8. Terms (in differential equations) & devices
GALA=T;TERMS(VFOL,N,N,N,N,P,P)
    GROUP 9. Properties of the medium (or media)

    GROUP 11. Initialization of variable or porosity fields
FIINIT(P1)=0.0;FIINIT(V1)=0.0;FIINIT(W1)=0.0;FIINIT(VFOL)=0.0
FIINIT(DEN1)=1.189
INIADD=F
PATCH(LIQUID,INIVAL,1,NX,1,6,1,NZ,1,1)
COVAL(LIQUID,VFOL,ZERO,1.000);COVAL(LIQUID,DEN1,ZERO,1000.5)
    ** create a large liquid bath
CONPOR(POOL,10000,CELL,1,NX,1,NY,7,NZ)
    GROUP 13. Boundary conditions and special sources
    ** open top
PATCH(REFP,CELL,1,NX,NY,NY, 1,NZ,1,LSTEP);COVAL(REFP,P1,FIXP,ZERO)
    ** moving plate
PATCH(MOVE,LWALL,1,NX,1,NY, 1, 1,1,LSTEP);COVAL(MOVE,V1,1.00,0.01)
    ** stationary bottom wall
PATCH(FIXS,SWALL,1,NX,1, 1, 1,NZ,1,LSTEP);COVAL(FIXS,W1,1.00,0.00)
    ** stationary side wall
PATCH(FIXH,HWALL,1,NX,1,NY,NZ,NZ,1,LSTEP);COVAL(FIXH,V1,1.00,0.00)
    ** apply gravity
PATCH(GRAV,PHASEM,1,NX,1,NY,1,NZ,1,LSTEP)
COVAL(GRAV,V1,FIXFLU,-9.81)
    GROUP 15. Termination of sweeps
LSWEEP=11
    GROUP 16. Termination of iterations
RESREF(P1)=1.E-6;RESREF(V1)=1.E-6;RESREF(W1)=1.E-6
    GROUP 17. Under-relaxation devices
RELAX(V1,FALSDT,1.E-3);RELAX(W1,FALSDT,1.E-3)
    GROUP 19. Data communicated by satellite to GROUND
IDISPB=1;IDISPC=LSTEP;IDISPA=4;CSG1=N;IPRPSA=67;IPRPSB=0
HOL=T
    GROUP 22. Spot-value print-out
TSTSWP=-1;NTPRIN=10;NPRMON=LSWEEP;IYMON=5;IZMON=1
SPEDAT(SET,GXMONI,TRANSIENT,L,F)
    GROUP 23. Field print-out and plot control
OUTPUT(P1,Y,Y,Y,Y,Y,Y);OUTPUT(V1,Y,Y,Y,Y,Y,Y)
OUTPUT(W1,Y,Y,Y,Y,Y,Y);OUTPUT(VFOL,Y,N,y,N,N,N)
OUTPUT(DEN1,Y,N,y,N,N,N);IHOLA=3
    GROUP 24. Dumps for restarts