Group 1. Run Title
INTEGER(IRR,IMCD,IBIN,ITHM)
MESG(
MESG(Select test case- 1  BATCH
MESG(                  2  WILKE
MESG(                  3  FICK
MESG(                  4  STEFMAX
MESG(                  5  WILKEREF
MESG(                  6  SORCJLJ
MESG(                  7  SORCJRS
MESG(                  8  SOREXLJ
MESG(                  9  SOREXRS
READVDU(IRR,INT,1)
CASE (IRR) OF
WHEN 1
TEXT(SIMPLE FLOW TUBE TEST CASE - BATCH :D201
TITLE
WHEN 2
TEXT(SIMPLE FLOW TUBE TEST CASE - WILKE :D201
TITLE
WHEN 3
TEXT(SIMPLE FLOW TUBE TEST CASE - FICK  :D201
TITLE
WHEN 4
TEXT(SIMPLE FLOW TUBE TEST CASE- STEFMAX:D201
TITLE
WHEN 5
TEXT(SIMPLE FLOW TUBE TEST CASE-WILKEREF:D201
TITLE
WHEN 6
TEXT(SIMPLE FLOW TUBE TEST CASE- SORCJLJ:D201
TITLE
WHEN 7
TEXT(SIMPLE FLOW TUBE TEST CASE- SORCJRS:D201
TITLE
WHEN 8
TEXT(SIMPLE FLOW TUBE TEST CASE- SOREXLJ:D201
TITLE
WHEN 9
TEXT(SIMPLE FLOW TUBE TEST CASE- SOREXRS:D201
TITLE
ENDCASE
  Groups 3, 4, 5  Grid Information
CARTES=F
NREGX=1;REGEXT(X,1)
IREGX=1;GRDPWR(X,1,0.01,1)
NREGY=2;REGEXT(Y,1)
IREGY=1;GRDPWR(Y,4,0.05,-2)
IREGY=2;GRDPWR(Y,2,0.0005,1)
NREGZ=3;REGEXT(Z,1)
IREGZ=1;GRDPWR(Z,4,0.07,1)
IREGZ=2;GRDPWR(Z,3,0.02,1)
IREGZ=3;GRDPWR(Z,2,0.02,1)
  Group 7. Variables: STOREd,SOLVEd,NAMEd
ONEPHS=T
NAME(C1)=S177;NAME(C2)=S80 ;NAME(C3)=S97
NAME(96)=BLOK;NAME(97)=TEM1;NAME(98)=DEN1;NAME(99)=PRPS
SOLVE(P1,V1,W1,TEM1,S80,S97)
STORE(S177,PRPS,DEN1,BLOK,CNDT,SPHT,DEPO,ENUL)
SOLUTN(P1,Y,Y,Y,N,N,Y) ; SOLUTN(TEM1,Y,Y,Y,N,N,Y)
SOLUTN(S80,Y,Y,Y,N,N,Y); SOLUTN(S97,Y,Y,Y,N,N,Y)
IVARBK=-1;ISOLBK=1
  Group 8. Terms & Devices
NEWRH1=T ; NEWENL=T
ISOLX=0  ; ISOLY=0  ;ISOLZ=0
USOURC=T ; UDIFNE=T ;DIFCUT=0.0
  Group 9. Properties
 ;  ;  ; 
PRNDTL(S177)=-GRND8; PRNDTL(S80)=-GRND8; PRNDTL(S97)=-GRND8
PRESS0=133.0
IF(IRR.EQ.5) THEN
+ TMP1A=273.0
ENDIF
CSG10='q1'
  MATFLG=T;IMAT=2
    70      GRND8   GRND8   GRND8   GRND8   1.000
    0.0 0.0 0.0
    0.0 0.0 0.0
    0.0 0.0 0.0
    0.0 0.0 0.0
    154    1.926E4   1.0   134.0   178.0   0.0
  Group 11.Initialise Var/Porosity Fields
FIINIT(P1)=0.0     ; FIINIT(W1)=0.0   ; FIINIT(V1)=0.0
FIINIT(TEM1)=290   ; FIINIT(S177)=0.94; FIINIT(S80)=0.06
FIINIT(S97)=1.0E-10; FIINIT(BLOK)=1.0 ; FIINIT(PRPS)=70
CONPOR(WAF,-1.00,CELL,-#1,-#1,-#2,-#2,-#2,-#2)
INIT(WAF,BLOK,0.000E+00,2.000E+00)
INIT(WAF,PRPS,0.000E+00,154)
RSTPRP=F;INIADD=F
  Group 13. Boundary & Special Sources
REAL(RHOIN,UAV);RHOIN=1.7E-3;UAV=1.345*3.
INLET(IN,LOW,#1,#1,#1,#2,#1,#1,#1,#1)
VALUE(IN,P1,UAV*RHOIN);VALUE(IN,W1,UAV)
VALUE(IN,TEM1,290);VALUE(IN,S177,0.94)
VALUE(IN,S80,0.06);VALUE(IN,S97,0.000)
OUTLET(OUT,HIGH,#1,#1,#1,#2,#3,#3,#1,#1)
VALUE(OUT,V1,SAME);VALUE(OUT,W1,SAME)
VALUE(OUT,P1,0.0) ;VALUE(OUT,TEM1,SAME)
VALUE(OUT,S177,SAME);VALUE(OUT,S80,SAME)
VALUE(OUT,S97,SAME)
PATCH(WALL,NWALL,#1,#1,#2,#2,#1,#3,#1,#1)
COVAL(WALL,W1,1.0,0.0);COVAL(WALL,TEM1,1.0,290)
PATCH(WAFER,CELL,#1,#1,#2,#2,#2,#2,#1,#1)
COVAL(WAFER,TEM1,FIXVAL,733)
PATCH(SURF1,NORTH,#1,#1,#1,#1,#2,#2,1,LSTEP)
COVAL(SURF1,S80,FIXFLU,GRND1)
COVAL(SURF1,S97,FIXFLU,GRND1)
IF(IRR.EQ.1) THEN
+ COVAL(SURF1,P1,4.000,GRND1)
ELSE
+ COVAL(SURF1,P1,1.000,GRND1)
ENDIF
COVAL(SURF1,TEM1,FIXFLU,GRND1)
PATCH(RELT,PHASEM,1,NX,1,NY,1,NZ,1,LSTEP)
COVAL(RELT,S80,GRND1,SAME)
COVAL(RELT,S97,GRND1,SAME)
EGWF=T
  Group 15. Terminate Sweeps
LSWEEP=50 ; SELREF=T ; RESFAC=1.000E-03
  Group 16. Terminate Iterations
ENDIT(P1)=1.000E-03 ; ENDIT(V1)=1.000E-03
ENDIT(W1)=1.000E-03 ; ENDIT(TEM1)=1.000E-03
  Group 17. Relaxation
RELAX(P1,LINRLX,1.000E+00) ; RELAX(V1,FALSDT,1.042E-01)
RELAX(W1,FALSDT,1.042E-01) ; RELAX(TEM1,FALSDT,1.042E+01)
REAL(TREL) ; TREL=1.0E2
RELAX(S80,FALSDT,TREL)     ; RELAX(S97,FALSDT,TREL)
  Group 18. Limits
VARMIN(S177)=1.0E-20 ; VARMAX(S177)=1.0
VARMIN(S80)=1.0E-20  ; VARMAX(S80)=1.0
VARMIN(S97)=1.0E-20  ; VARMAX(S97)=1.0
VARMIN(TEM1)=200.0   ; VARMAX(TEM1)=1000.0
  Group 19. EARTH Calls To GROUND Station
GENK=T
IF(IRR.EQ.5) THEN
+ IBIN=2
ELSE
+ IBIN=4
ENDIF
IF(IRR.EQ.3) THEN
+ IMCD=1
ELSE
+ IF(IRR.EQ.4) THEN
+   IMCD=3
+ ELSE
+   IMCD=2
+ ENDIF
ENDIF
SPEDAT(SET,CVD,MCDOPT,I,IMCD)
SPEDAT(SET,CVD,BINOPT,I,IBIN)
IF(IRR.GT.5) THEN
CASE (IRR) OF
WHEN 6
+ ITHM=2
WHEN 7
+ ITHM=1
WHEN 8
+ ITHM=4
WHEN 9
+ ITHM=3
ENDCASE
SPEDAT(SET,CVD,THMDIF,L,T)
SPEDAT(SET,CVD,THMOPT,I,ITHM)
SPEDAT(SET,CVD,THMFRQ,I,3)
SPEDAT(SET,CVD,THMRLX,R,1.0)
ENDIF
SPEDAT(SET,CVD,MCPROP,I,3)
SPEDAT(SET,CVD,CHMRLX,R,0.5)
SPEDAT(SET,CVD,NGREAC,I,0)
SPEDAT(SET,CVD,NSREAC,I,1)
SPEDAT(SET,CVD,SREAC(1),I,17)
NAMGRD=CVD
  Group 21. Print-out of Variables
OUTPUT(BLOK,Y,N,Y,N,N,N)
OUTPUT(DEN1,Y,N,Y,N,N,N)
OUTPUT(PRPS,Y,N,Y,N,N,N)
  Group 22. Monitor Print-Out
IXMON=1 ; IYMON=4 ; IZMON=5
TSTSWP=-1
  Group 24. Dumps For Restarts