PHOTON USE
  p
  phi
 
 
 
  msg grid
  gr x 1;gr x m
  gr y 1;gr y m
  msg Press  to continue
  pause
  cl
  vi x
  up z
  scale
  cl
  msg 1st-phase (light phase) velocity vectors on plane ix=3
  set vec ref;0.1
  vec x 3 sh
  gr ou x 3
  msg Press  to continue
  pause
  cl
  msg 2nd-phase (dense phase) velocity vectors on plane ix=3
  set vec ref;0.004
  set vec comp
  - v2 w2
  vec x 3 sh
  gr ou x 3
  msg Press  to continue
  pause
  cl
  msg contour of pressure on plane ix=3
  con p1 x 3 fi;0.002
  gr ou x 3
  msg Press  to continue
  pause
  cl
  msg contour of 1st-phase (vapour) volume fraction on plane ix=3
  con gas x 3 fi;0.002
  gr ou x 3
  msg Press  to continue
  pause
  cl
  msg contour of shell-side vapour enthalpy on plane ix=3
  con 1sth x 3 fi;0.002
  gr ou x 3
  msg Press  to continue
  pause
  cl
  msg contour of shell-side liquid enthalpy on plane ix=3
  con 2ndh x 3 fi;0.002
  gr ou x 3
  msg Press  to continue
  pause
  cl
  msg contour of tube-side fluid enthalpy on plane ix=3
  con tubh x 3 fi;0.002
  gr ou x 3
  msg Press e to END
  enduse
 
    CASE W800 As for case W799, but time-dependent
L(W799)
    GROUP 1. Run title
TEXT(AS FOR CASE W799 BUT TRANS. :  W800
TITLE
    User-defined variables:
 
    TIME = duration of process simulated
    INTVLS = number of time intervals
REAL(TIME,DELT,CAP1,CAP2,PORCON);INTEGER(INTVLS)
TIME=20.0;INTVLS=10;DELT=2.0;CAP1=0.5;CAP2=0.3
TIME=10.0;INTVLS=5;DELT=2.0;CAP1=0.5;CAP2=0.3
PORCON=1.E-5
    GROUP 2. Transience; time-step specification
STEADY=F;GRDPWR(T,INTVLS,TIME,1.0)
SPEDAT(SET,GXMONI,TRANSIENT,L,F) 
    A small volume porosity is used to nullify the effect of
    the built-in transient terms, which are replaced
    by the use of the special patch named NEOLD below.
 
STORE(VPOR);FIINIT(VPOR)=PORCON
RELAX(U1,FALSDT,1.E-1);RELAX(V1,FALSDT,1.E-1)
RELAX(W1,FALSDT,1.E-1)
PATCH(NEOLD,VOLUME,1,NX,1,NY,1,NZ,1,LSTEP)
COVAL(NEOLD,1STH,CAP1/DELT,LOCNE)
COVAL(NEOLD,2NDH,CAP2/DELT,LOCNE)
NTPRIN=LSTEP
OUTPUT(VPOR,N,N,N,N,N,N)
lsweep=20
PATCH(TIMEPLOT,PROFIL,NX,NX,2,2,NZ,NZ,1,LSTEP)
PLOT(TIMEPLOT,1STH,T2IN,T1IN);PLOT(TIMEPLOT,2NDH,T2IN,T1IN)