```
TALK=T;RUN(1,1)
DISPLAY

This case is like 369, but use the z-direction instead of x.

An InForm-made box object moves uniformly along a duct, the
cross-section of which it wholly fills.
The duct contains a porous medium causing resistance to flow
The ends of the duct are partially open.
The fluid is compressible.

Parameters are provided which enable the effects on convergence
and accuracy to be noted of:
* box size divided by cell size;
* time to traverse a cell divided by time step;
* the compressibility factor;
* the resistance coefficient;
* pressure-drop coefficients at the ends of the duct.

Line-printer profile plots are used for assisting interpretation

ENDDIS
************************************************************
Group 1. Run Title and Number
************************************************************
************************************************************

TEXT(MOFOR - In-Form: BOX linear motion in Z )

************************************************************
************************************************************

IRUNN = 1 ;LIBREF = 0
************************************************************
Group 2. Time dependence
* Set overall time and no. of steps
TFIRST =0. ;TLAST =1.
FSTEP = 1 ;LSTEP = 50
TFRAC(1)=0.02 ;TFRAC(11)=0.22
TFRAC(21)=0.42 ;TFRAC(31)=0.62
TFRAC(41)=0.82
************************************************************
Group 3. X-Direction Grid Spacing
CARTES = T
NX = 1
XULAST =0.1
XFRAC(1)=1.
************************************************************
Group 4. Y-Direction Grid Spacing
NY = 1
YVLAST =2.
YFRAC(1)=1.
************************************************************
Group 5. Z-Direction Grid Spacing
PARAB = F
NZ = 100
ZWLAST =1.
ZFRAC(1)=4.999996E-04 ;ZFRAC(2)=9.999991E-04
ZFRAC(3)=1.499999E-03 ;ZFRAC(4)=1.999998E-03
ZFRAC(5)=2.499998E-03 ;ZFRAC(6)=2.999997E-03
ZFRAC(7)=3.499997E-03 ;ZFRAC(8)=3.999996E-03
ZFRAC(9)=4.499996E-03 ;ZFRAC(10)=4.999995E-03
ZFRAC(11)=0.016056 ;ZFRAC(12)=0.027111
ZFRAC(13)=0.038167 ;ZFRAC(14)=0.049222
ZFRAC(15)=0.060278 ;ZFRAC(16)=0.071333
ZFRAC(17)=0.082389 ;ZFRAC(18)=0.093444
ZFRAC(19)=0.1045 ;ZFRAC(20)=0.115556
ZFRAC(21)=0.126611 ;ZFRAC(22)=0.137667
ZFRAC(23)=0.148722 ;ZFRAC(24)=0.159778
ZFRAC(25)=0.170833 ;ZFRAC(26)=0.181889
ZFRAC(27)=0.192944 ;ZFRAC(28)=0.204
ZFRAC(29)=0.215056 ;ZFRAC(30)=0.226111
ZFRAC(31)=0.237167 ;ZFRAC(32)=0.248222
ZFRAC(33)=0.259278 ;ZFRAC(34)=0.270333
ZFRAC(35)=0.281389 ;ZFRAC(36)=0.292444
ZFRAC(37)=0.3035 ;ZFRAC(38)=0.314556
ZFRAC(39)=0.325611 ;ZFRAC(40)=0.336667
ZFRAC(41)=0.347722 ;ZFRAC(42)=0.358778
ZFRAC(43)=0.369833 ;ZFRAC(44)=0.380889
ZFRAC(45)=0.391944 ;ZFRAC(46)=0.403
ZFRAC(47)=0.414056 ;ZFRAC(48)=0.425111
ZFRAC(49)=0.436167 ;ZFRAC(50)=0.447222
ZFRAC(51)=0.458278 ;ZFRAC(52)=0.469333
ZFRAC(53)=0.480389 ;ZFRAC(54)=0.491444
ZFRAC(55)=0.5025 ;ZFRAC(56)=0.513556
ZFRAC(57)=0.524611 ;ZFRAC(58)=0.535667
ZFRAC(59)=0.546722 ;ZFRAC(60)=0.557778
ZFRAC(61)=0.568833 ;ZFRAC(62)=0.579889
ZFRAC(63)=0.590944 ;ZFRAC(64)=0.602
ZFRAC(65)=0.613056 ;ZFRAC(66)=0.624111
ZFRAC(67)=0.635167 ;ZFRAC(68)=0.646222
ZFRAC(69)=0.657278 ;ZFRAC(70)=0.668333
ZFRAC(71)=0.679389 ;ZFRAC(72)=0.690444
ZFRAC(73)=0.7015 ;ZFRAC(74)=0.712556
ZFRAC(75)=0.723611 ;ZFRAC(76)=0.734667
ZFRAC(77)=0.745722 ;ZFRAC(78)=0.756778
ZFRAC(79)=0.767833 ;ZFRAC(80)=0.778889
ZFRAC(81)=0.789944 ;ZFRAC(82)=0.801
ZFRAC(83)=0.812056 ;ZFRAC(84)=0.823111
ZFRAC(85)=0.834167 ;ZFRAC(86)=0.845222
ZFRAC(87)=0.856278 ;ZFRAC(88)=0.867333
ZFRAC(89)=0.878389 ;ZFRAC(90)=0.889444
ZFRAC(91)=0.9005 ;ZFRAC(92)=0.911556
ZFRAC(93)=0.922611 ;ZFRAC(94)=0.933667
ZFRAC(95)=0.944722 ;ZFRAC(96)=0.955778
ZFRAC(97)=0.966833 ;ZFRAC(98)=0.977889
ZFRAC(99)=0.988944 ;ZFRAC(100)=1.
************************************************************
Group 6. Body-Fitted Coordinates
************************************************************
Group 7. Variables: STOREd,SOLVEd,NAMEd
ONEPHS = T
NAME(1)=P1 ;NAME(7)=W1
NAME(16)=C1 ;NAME(143)=FLOZ
NAME(144)=RSW1 ;NAME(145)=AHP1
NAME(146)=APW1 ;NAME(147)=SOW1
NAME(148)=RHO1 ;NAME(149)=PRPS
NAME(150)=OBID
* Y in SOLUTN argument list denotes:
* 1-stored 2-solved 3-whole-field
* 4-point-by-point 5-explicit 6-harmonic averaging
SOLUTN(P1,Y,Y,N,N,N,Y)
SOLUTN(W1,Y,Y,N,N,N,Y)
SOLUTN(C1,Y,Y,N,N,N,Y)
SOLUTN(FLOZ,Y,N,N,N,N,Y)
SOLUTN(RSW1,Y,N,N,N,N,Y)
SOLUTN(AHP1,Y,N,N,N,N,Y)
SOLUTN(APW1,Y,N,N,N,N,Y)
SOLUTN(SOW1,Y,N,N,N,N,Y)
SOLUTN(RHO1,Y,N,N,N,N,Y)
SOLUTN(PRPS,Y,N,N,N,N,Y)
SOLUTN(OBID,Y,N,N,N,N,Y)
DEN1 = 148
PRPS = 149
************************************************************
Group 8. Terms & Devices
* Y in TERMS argument list denotes:
* 1-built-in source 2-convection 3-diffusion 4-transient
* 5-first phase variable 6-interphase transport
TERMS(P1,Y,Y,Y,N,Y,Y)
TERMS(W1,Y,Y,Y,Y,Y,Y)
TERMS(C1,N,Y,N,Y,Y,Y)
DIFCUT =0.5 ;ZDIFAC =1.
GALA = F ;ADDDIF = F
ISOLX = -1 ;ISOLY = -1 ;ISOLZ = -1
DENPCO = T
************************************************************
Group 9. Properties used if PRPS is not
stored, and where PRPS = -1.0 if it is!
RHO1 =1. ;TMP1 =0.
EL1 =0.
TSURR =0. ;TEMP0 =0.
PRESS0 =0.
DVO1DT =0. ;DRH1DP =0.
EMISS =0. ;SCATT =0.
ENUL =1.0E-05 ;ENUT =0.
PRNDTL(W1)=1.
PRT(W1)=1.
CP1 =1. ;CP2 =1.
************************************************************
Group 10.Inter-Phase Transfer Processes
************************************************************
Group 11.Initial field variables (PHIs)
FIINIT(P1)=1.0E-10 ;FIINIT(W1)=1.0E-10
FIINIT(C1)=1.0E-10 ;FIINIT(FLOZ)=1.0E-10
FIINIT(RSW1)=1.0E-10 ;FIINIT(AHP1)=1.0E-10
FIINIT(APW1)=1.0E-10 ;FIINIT(SOW1)=1.0E-10
FIINIT(RHO1)=1.0E-10 ;FIINIT(PRPS)=-1.
FIINIT(OBID)=1.0E-10

PATCH(INIC1 ,INIVAL, 1, 1, 1, 1, 1, 20, 1, 1)
INIT(INIC1 ,C1 ,0. ,1. )

PATCH(^OB3 ,INIVAL, 3, 0, 0, 0, 0, 0, 1, 1)
INIT(^OB3 ,PRPS,0. ,198. )
FSWEEP = 1
NAMFI =CHAM
************************************************************
Group 12. Patchwise adjustment of terms
Patches for this group are printed with those
for Group 13.
Their names begin either with GP12 or &
************************************************************
Group 13. Boundary & Special Sources

PATCH(RESIST ,VOLUME, 1, 1, 1, 1, 1, 100, 1, 50)
COVAL(RESIST ,W1 ,10. ,0. )

PATCH(INLET ,CELL , 1, 1, 1, 1, 1, 1, 1, 50)
COVAL(INLET ,P1 ,1. ,0. )
COVAL(INLET ,C1 ,0. ,1. )

PATCH(OUTLET ,CELL , 1, 1, 1, 1, 100, 100, 1, 50)
COVAL(OUTLET ,P1 ,1. ,0. )
XCYCLE = F
EGWF = T
WALLCO = GRND2
************************************************************
Group 14. Downstream Pressure For PARAB
************************************************************
Group 15. Terminate Sweeps
LSWEEP = 60 ;ISWC1 = 1
LITHYD = 1 ;LITFLX = 1 ;LITC = 1 ;ITHC1 = 1
SELREF = T
RESFAC =1.0E-02
************************************************************
Group 16. Terminate Iterations
LITER(P1)=20 ;LITER(W1)=10
LITER(C1)=20
ENDIT(P1)=1.0E-03 ;ENDIT(W1)=1.0E-03
ENDIT(C1)=1.0E-03
************************************************************
Group 17. Relaxation
RELAX(P1,LINRLX,0.1)
RELAX(W1,FALSDT,1.)
RELAX(C1,FALSDT,1.0E+09)
RELAX(FLOZ,LINRLX,1.)
RELAX(RSW1,LINRLX,1.)
RELAX(AHP1,LINRLX,1.)
RELAX(APW1,LINRLX,1.)
RELAX(SOW1,LINRLX,1.)
RELAX(RHO1,LINRLX,0.1)
RELAX(PRPS,LINRLX,1.)
RELAX(OBID,LINRLX,1.)
OVRRLX =0.
EXPERT = F ;NNORSL = F
************************************************************
Group 18. Limits
VARMAX(P1)=1.0E+10 ;VARMIN(P1)=-1.0E+10
VARMAX(W1)=1.0E+06 ;VARMIN(W1)=-1.0E+06
VARMAX(C1)=1.0E+10 ;VARMIN(C1)=-1.0E+10
VARMAX(FLOZ)=1.0E+10 ;VARMIN(FLOZ)=-1.0E+10
VARMAX(RSW1)=1.0E+10 ;VARMIN(RSW1)=-1.0E+10
VARMAX(AHP1)=1.0E+10 ;VARMIN(AHP1)=-1.0E+10
VARMAX(APW1)=1.0E+10 ;VARMIN(APW1)=-1.0E+10
VARMAX(SOW1)=1.0E+10 ;VARMIN(SOW1)=-1.0E+10
VARMAX(RHO1)=1.0E+10 ;VARMIN(RHO1)=-1.0E+10
VARMAX(PRPS)=1.0E+10 ;VARMIN(PRPS)=-1.0E+10
VARMAX(OBID)=1.0E+10 ;VARMIN(OBID)=-1.0E+10
************************************************************
Group 19. Data transmitted to GROUND
CONWIZ = T
ISG52 = 2
ISG62 = 0
SPEDAT(SET,MOFOR,MOFFILE,C,NOTSET)
SPEDAT(SET,GXMONI,TRANSIENT,L,F)
SPEDAT(SET,PROPERTY,RHO1,C,=1.0+1.0E-04*P1)
SPEDAT(SET,PROPERTY,DRH1DP,C,=1.0E-04)
SPEDAT(SET,MOVOB,BOX,C,=POS(0.0&0.0&0.15*ZWLAST+0.666667*TIM&0&0&\$)
SPEDAT(SET,MOVOB,BOX,C,0&0&0&0&0))
SPEDAT(SET,STORED,SOW1,C,=SORC(W1)!SOURCS)
SPEDAT(SET,STORED,APW1,C,=APCO(W1)!EQCOEF)
SPEDAT(SET,STORED,AHP1,C,=AHCO(P1)!EQCOEF)
SPEDAT(SET,STORED,RSW1,C,=RESI(W1)!RESIDU)
SPEDAT(SET,STORED,FLOZ,C,=RHO1*W1!ZSLFIN)
SPEDAT(SET,OBJNAM,^OB3,C,BOX)
SPEDAT(SET,OBJTYP,^OB3,C,BLOCKAGE)
SPEDAT(SET,BOX,MATERIAL,R,198.)
SPEDAT(SET,FACETDAT,NUMOBJ,I,3)
SPEDAT(SET,MATERIAL,198,L,T)
************************************************************
Group 20. Preliminary Printout
DISTIL = T ;NULLPR = F
NDST = 0
DSTTOL =1.0E-02
EX(P1)=0.01829 ;EX(W1)=0.13
EX(C1)=0.9001 ;EX(FLOZ)=0.2128
EX(RSW1)=6.341E-06 ;EX(AHP1)=0.1449
EX(APW1)=0.2856 ;EX(SOW1)=1.118E-04
EX(RHO1)=0.9996 ;EX(PRPS)=0.
EX(OBID)=0.
************************************************************
Group 21. Print-out of Variables
INIFLD = F ;SUBWGR = F
* Y in OUTPUT argument list denotes:
* 1-field 2-correction-eq. monitor 3-selective dumping
* 4-whole-field residual 5-spot-value table 6-residual table
OUTPUT(P1,Y,N,Y,Y,Y,Y)
OUTPUT(W1,Y,N,Y,Y,Y,Y)
OUTPUT(C1,Y,N,Y,Y,Y,Y)
OUTPUT(FLOZ,Y,N,Y,N,N,N)
OUTPUT(RSW1,Y,N,Y,N,N,N)
OUTPUT(AHP1,Y,N,Y,N,N,N)
OUTPUT(APW1,Y,N,Y,N,N,N)
OUTPUT(SOW1,Y,N,Y,N,N,N)
OUTPUT(RHO1,Y,N,Y,N,N,N)
OUTPUT(PRPS,Y,N,Y,N,N,N)
OUTPUT(OBID,Y,N,Y,N,N,N)
************************************************************
Group 22. Monitor Print-Out
IXMON = 1 ;IYMON = 1 ;IZMON = 1
NPRMON = 60 ;NPRMNT = 1 ;TSTSWP = -1
UWATCH = T ;USTEER = T
HIGHLO = F
************************************************************
Group 23.Field Print-Out & Plot Control
NPRINT = 100000 ;NUMCLS = 5
NTPRIN = 10 ;ISTPRF = 1 ;ISTPRL = 100000
NZPRIN = -1 ;IZPRF = 1 ;IZPRL = 10000
XZPR = F ;YZPR = F
IPLTF = 1 ;IPLTL = -1 ;NPLT = -1
ISWPRF = 1 ;ISWPRL = 100000
ITABL = 3 ;IPROF = 1
ABSIZ =0.5 ;ORSIZ =0.4
NTZPRF = 1 ;NCOLPF = 50
ICHR = 2 ;NCOLCO = 45 ;NROWCO = 20

PATCH(ZPLOT10 ,PROFIL, 1, 1, 1, 1, 1, 99, 10, 10)
PLOT(ZPLOT10 ,P1 ,0. ,0. )
PLOT(ZPLOT10 ,W1 ,0. ,0. )
PLOT(ZPLOT10 ,C1 ,0. ,2. )
PLOT(ZPLOT10 ,FLOZ,0. ,0. )
PLOT(ZPLOT10 ,RHO1,0.8 ,1.2 )

PATCH(ZPLOT20 ,PROFIL, 1, 1, 1, 1, 1, 99, 20, 20)
PLOT(ZPLOT20 ,P1 ,0. ,0. )
PLOT(ZPLOT20 ,W1 ,0. ,0. )
PLOT(ZPLOT20 ,C1 ,0. ,2. )
PLOT(ZPLOT20 ,FLOZ,0. ,0. )
PLOT(ZPLOT20 ,RHO1,0.8 ,1.2 )

PATCH(ZPLOT30 ,PROFIL, 1, 1, 1, 1, 1, 99, 30, 30)
PLOT(ZPLOT30 ,P1 ,0. ,0. )
PLOT(ZPLOT30 ,W1 ,0. ,0. )
PLOT(ZPLOT30 ,C1 ,0. ,2. )
PLOT(ZPLOT30 ,FLOZ,0. ,0. )
PLOT(ZPLOT30 ,RHO1,0.8 ,1.2 )

PATCH(ZPLOT40 ,PROFIL, 1, 1, 1, 1, 1, 99, 40, 40)
PLOT(ZPLOT40 ,P1 ,0. ,0. )
PLOT(ZPLOT40 ,W1 ,0. ,0. )
PLOT(ZPLOT40 ,C1 ,0. ,2. )
PLOT(ZPLOT40 ,FLOZ,0. ,0. )
PLOT(ZPLOT40 ,RHO1,0.8 ,1.2 )

PATCH(ZPLOT50 ,PROFIL, 1, 1, 1, 1, 1, 99, 50, 50)
PLOT(ZPLOT50 ,P1 ,0. ,0. )
PLOT(ZPLOT50 ,W1 ,0. ,0. )
PLOT(ZPLOT50 ,C1 ,0. ,2. )
PLOT(ZPLOT50 ,FLOZ,0. ,0. )
PLOT(ZPLOT50 ,RHO1,0.8 ,1.2 )
************************************************************
Group 24. Dumps For Restarts
SAVE = T ;NOWIPE = F
NSAVE =CHAM
STOP
```