PHOTON USE
  p
  n4
  1 1 200
 
 
 
  msg(    film formation and vector field at t= 8s
  gr ou x 1
  con den1 x 1 fi
  0.001
  vec x 1
  pause
  p
  n8
  1 1 200
 
 
 
  msg(    film formation and vector field at t=16s
  gr ou x 1
  con den1 x 1 fi
  0.001
  vec x 1
  pause
  p
  n12
  1 1 200
 
 
 
  msg(    film formation and vector field at t=24s
  gr ou x 1
  con den1 x 1 fi
  0.001
  vec x 1
  pause
  p
  n16
  1 1 200
 
 
 
  msg(    film formation and vector field at t=32s
  gr ou x 1
  con den1 x 1 fi
  0.001
  vec x 1
  pause
  p
  n20
  1 1 200
 
 
 
  msg(    film formation and vector field at t=40s
  gr ou x 1
  con den1 x 1 fi
  0.001
  vec x 1
  msg            -
  msg Press e to END
  enduse
 
    GROUP 1. Run title and other preliminaries
TEXT(FILM FORMATION BY SEM METHOD:P105
TITLE
  DISPLAY
  FREE SURFACE FLOW -  Film Formation by SEM 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
solve(surn,vfol);store(den1,prps)
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(surn,n,n,n,n,p,p)
    GROUP 9. Properties of the medium (or media)
enul=1.e-3
    GROUP 11. Initialization of variable or porosity fields
fiinit( p1 )=0.0;fiinit( v1 )=0.0
fiinit( w1 )=0.0;fiinit(surn)=0.0
fiinit(den1)=1.161;iniadd=f
patch(liquid,inival,1,nx,1,6,1,nz,1,1)
coval(liquid,surn,zero,1.000)
coval(liquid,den1,zero,1000.5)
    ** create a pool of liquid
conpor(pool,10000,cell,1,nx,1,ny,7,nz)
    GROUP 13. Boundary conditions and special sources
    ** open top of liquid pool
patch(refp,cell,1,nx,ny,ny,1,nz,1,lstep)
coval(refp,p1,fixp,zero)
    **
patch(tcon,cell,1,nx,1,ny,1,nz,1,lstep)
coval(tcon,surn,grnd,grnd)
    ** moving plate on which a film of liquid will form
patch(move,lwall,1,nx,1,ny, 1, 1,1,lstep)
coval(move,v1,1.00,0.01)
    ** stationay bottom wall
patch(fixs,swall,1,nx,1, 1, 1,nz,1,lstep)
coval(fixs,w1,1.00,0.00)
    ** stationay 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;liter(surn)=1
    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,0.01);relax(w1,falsdt,0.1)
    GROUP 19. Data communicated by satellite to GROUND
iprpsa=67;iprpsb=0;idispa=4;csg1=n
SPEDAT(SET,GXMONI,TRANSIENT,L,F)
SURF=T ; rlolim=0.3 ; ruplim=0.7
varmin(surn)=0.0; varmax(surn)=1.0
    GROUP 22. Spot-value print-out
tstswp=-1;ntprin=10
nprmon=lsweep;iymon=5;izmon=1
    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(surn,y,n,y,n,n,n)
    GROUP 24. Dumps for restarts