L(796)
 
    CASE 8 As for case 7, but time-dependent
    GROUP 1. Run title
TEXT(Case 8; As For Case 7 But Trans  
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
PORCON=1.E-4
 
    GROUP 2. Transience; time-step specification
STEADY=F;GRDPWR(T,INTVLS,TIME,1.0)
 
    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
 
    RESTRT(P1,U1,V1,W1)
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)
 
LSWEEP=50;NTPRIN=LSTEP
OUTPUT(VPOR,N,N,N,N,N,N)
 
PATCH(TIMEPLOT,PROFIL,NX,NX,2,2,NZ,NZ,1,LSTEP)
PLOT(TIMEPLOT,1STH,T2IN,T1IN);PLOT(TIMEPLOT,2NDH,T2IN,T1IN)