L(249) 
  PHOTON USE
  p;;;;;
 
  msg the grid. press return
  gr z 1;pause;gr off;red;gr ou z 1;vec z 1 sh
  msg Velocity vectors. press return
  pause;vec off;red;con v1 z 1 fi;0.001
  msg vertical-velocity contours. press return
  pause;con off;red;con u1 z 1 fi;0.001
  msg horizontal-velocity contours. press return
  pause;con off;red;con h1 z 1 fi;0.001
  msg temperature contours. Press e to END
  enduse
 
  *************************************************
    GROUP 1. Run title and other preliminaries
TEXT(Expert Test-The Driven-Cavity Flow
TITLE
mesg(PC486/50 time last reported as appx. 1. min
iswc1=1
tstswp=-10;LSWEEP=500;ipltl=lsweep;nplt=1
ITABL=1;ipltf=2
  Begin interactive session
  DISPLAY
 
 
 
  The "driven-cavity" problem; PHOENICS library case 249.
  The flow is laminar, two-dimensional and steady
  The Navier-Stokes equations are solved by PHOENICS
  The grid is coarse (10*10) to save time
  The under-relaxation of velocity has been made very heavy,
  as though by an inexperienced user
 
 
  ***********************************************************
  ENDDIS
real(dtf0,reyno)
dtf0=1.0e-4
mesgm(initial dtfals for u1 & v1 is :dtf0: If not OK, insert value
readvdu(dtf0,real,dtf0)
dtfals(3)=dtf0;dtfals(5)=dtf0
dtfals(1)=-1.0
 
reyno=1.e3
mesgm(Reynolds number is :reyno: If not OK, insert value
readvdu(reyno,real,reyno)
enul=rho1*xulast*1.0/reyno
mesgm(Reynolds number = :reyno:
 
mesga(Ask EXPERT to find optimum relaxation parameters? (y/n)
do ii=1,10
 mesg(
enddo
readvdu(ans,char,n)
if(:ans:.eq.Y) then
 mesg(
 mesg(EXPERT has been switched ON
#exprt
endif
 
resfac=0.01
selref=t
mesga(Computations will be terminated when all residuals fall
mesgb(below representative values computed by EARTH times  :resfac:
mesg(OK? If not, insert desired value
readvdu(resfac,real,resfac)
resfac
delay(100)
STOP