TALK=T;RUN(1,1)
  DISPLAY
  This case is like 369, but use the y-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.

  PHOTON USE commands are included to assist display of results.
  ENDDIS
  photon use
  p
  parphi


  con v1 x 1 z 2 m fi;0.01
  pause
  con floy x 1  fi;0.01
  pause
  con obid x 1 fi;0.01
  pause
  con p1 x 1 z 2 m fi;0.01
  pause
  con anp1 x 1 fi ;0.01
  pause
  con apv1 x 1 fi;0.01
  pause
  enduse
 ************************************************************
  Group 1. Run Title and Number
 ************************************************************
 ************************************************************
 
 TEXT(MOFOR - In-Form: BOX linear motion in Y )
 
 ************************************************************
 ************************************************************
 
 IRUNN = 1 ;LIBREF = 0
 ************************************************************
  Group 2. Time dependence
 STEADY = F
    * Set overall time and no. of steps
 TFIRST =0. ;TLAST =1.
 FSTEP = 1 ;LSTEP = 50
 TFRAC(1)=0.02 ;TFRAC(2)=0.04
 TFRAC(3)=0.06 ;TFRAC(4)=0.08
 TFRAC(5)=0.1 ;TFRAC(6)=0.12
 TFRAC(7)=0.14 ;TFRAC(8)=0.16
 TFRAC(9)=0.18 ;TFRAC(10)=0.2
 TFRAC(11)=0.22 ;TFRAC(12)=0.24
 TFRAC(13)=0.26 ;TFRAC(14)=0.28
 TFRAC(15)=0.3 ;TFRAC(16)=0.32
 TFRAC(17)=0.34 ;TFRAC(18)=0.36
 TFRAC(19)=0.38 ;TFRAC(20)=0.4
 TFRAC(21)=0.42 ;TFRAC(22)=0.44
 TFRAC(23)=0.46 ;TFRAC(24)=0.48
 TFRAC(25)=0.5 ;TFRAC(26)=0.52
 TFRAC(27)=0.54 ;TFRAC(28)=0.56
 TFRAC(29)=0.58 ;TFRAC(30)=0.6
 TFRAC(31)=0.62 ;TFRAC(32)=0.64
 TFRAC(33)=0.66 ;TFRAC(34)=0.68
 TFRAC(35)=0.7 ;TFRAC(36)=0.72
 TFRAC(37)=0.74 ;TFRAC(38)=0.76
 TFRAC(39)=0.78 ;TFRAC(40)=0.8
 TFRAC(41)=0.82 ;TFRAC(42)=0.84
 TFRAC(43)=0.86 ;TFRAC(44)=0.88
 TFRAC(45)=0.9 ;TFRAC(46)=0.92
 TFRAC(47)=0.94 ;TFRAC(48)=0.96
 TFRAC(49)=0.98 ;TFRAC(50)=1.
 ************************************************************
  Group 3. X-Direction Grid Spacing
 CARTES = T
 NX = 1
 XULAST =2.
 XFRAC(1)=1.
 ************************************************************
  Group 4. Y-Direction Grid Spacing
 NY = 100
 YVLAST =1.
 YFRAC(1)=4.999995E-03 ;YFRAC(2)=0.015051
 YFRAC(3)=0.025101 ;YFRAC(4)=0.035152
 YFRAC(5)=0.045202 ;YFRAC(6)=0.055253
 YFRAC(7)=0.065303 ;YFRAC(8)=0.075354
 YFRAC(9)=0.085404 ;YFRAC(10)=0.095455
 YFRAC(11)=0.105505 ;YFRAC(12)=0.115556
 YFRAC(13)=0.125606 ;YFRAC(14)=0.135657
 YFRAC(15)=0.145707 ;YFRAC(16)=0.155758
 YFRAC(17)=0.165808 ;YFRAC(18)=0.175859
 YFRAC(19)=0.185909 ;YFRAC(20)=0.19596
 YFRAC(21)=0.20601 ;YFRAC(22)=0.216061
 YFRAC(23)=0.226111 ;YFRAC(24)=0.236162
 YFRAC(25)=0.246212 ;YFRAC(26)=0.256263
 YFRAC(27)=0.266313 ;YFRAC(28)=0.276364
 YFRAC(29)=0.286414 ;YFRAC(30)=0.296465
 YFRAC(31)=0.306515 ;YFRAC(32)=0.316566
 YFRAC(33)=0.326616 ;YFRAC(34)=0.336667
 YFRAC(35)=0.346717 ;YFRAC(36)=0.356768
 YFRAC(37)=0.366818 ;YFRAC(38)=0.376869
 YFRAC(39)=0.386919 ;YFRAC(40)=0.39697
 YFRAC(41)=0.40702 ;YFRAC(42)=0.417071
 YFRAC(43)=0.427121 ;YFRAC(44)=0.437172
 YFRAC(45)=0.447222 ;YFRAC(46)=0.457273
 YFRAC(47)=0.467323 ;YFRAC(48)=0.477374
 YFRAC(49)=0.487424 ;YFRAC(50)=0.497475
 YFRAC(51)=0.507525 ;YFRAC(52)=0.517576
 YFRAC(53)=0.527626 ;YFRAC(54)=0.537677
 YFRAC(55)=0.547727 ;YFRAC(56)=0.557778
 YFRAC(57)=0.567828 ;YFRAC(58)=0.577879
 YFRAC(59)=0.587929 ;YFRAC(60)=0.59798
 YFRAC(61)=0.60803 ;YFRAC(62)=0.618081
 YFRAC(63)=0.628131 ;YFRAC(64)=0.638182
 YFRAC(65)=0.648232 ;YFRAC(66)=0.658283
 YFRAC(67)=0.668333 ;YFRAC(68)=0.678384
 YFRAC(69)=0.688434 ;YFRAC(70)=0.698485
 YFRAC(71)=0.708535 ;YFRAC(72)=0.718586
 YFRAC(73)=0.728636 ;YFRAC(74)=0.738687
 YFRAC(75)=0.748737 ;YFRAC(76)=0.758788
 YFRAC(77)=0.768838 ;YFRAC(78)=0.778889
 YFRAC(79)=0.788939 ;YFRAC(80)=0.79899
 YFRAC(81)=0.80904 ;YFRAC(82)=0.819091
 YFRAC(83)=0.829141 ;YFRAC(84)=0.839192
 YFRAC(85)=0.849242 ;YFRAC(86)=0.859293
 YFRAC(87)=0.869343 ;YFRAC(88)=0.879394
 YFRAC(89)=0.889444 ;YFRAC(90)=0.899495
 YFRAC(91)=0.909545 ;YFRAC(92)=0.919596
 YFRAC(93)=0.929646 ;YFRAC(94)=0.939697
 YFRAC(95)=0.949748 ;YFRAC(96)=0.959798
 YFRAC(97)=0.969849 ;YFRAC(98)=0.979899
 YFRAC(99)=0.989949 ;YFRAC(100)=1.
 ************************************************************
  Group 5. Z-Direction Grid Spacing
 PARAB = F
 NZ = 1
 ZWLAST =0.1
 ZFRAC(1)=1.
 ************************************************************
  Group 6. Body-Fitted Coordinates
 ************************************************************
  Group 7. Variables: STOREd,SOLVEd,NAMEd
 ONEPHS = T
 NAME(1)=P1 ;NAME(5)=V1
 NAME(143)=FLOY ;NAME(144)=RSV1
 NAME(145)=ANP1 ;NAME(146)=APV1
 NAME(147)=SOV1 ;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(V1,Y,Y,N,N,N,Y)
 SOLUTN(FLOY,Y,N,N,N,N,Y)
 SOLUTN(RSV1,Y,N,N,N,N,Y)
 SOLUTN(ANP1,Y,N,N,N,N,Y)
 SOLUTN(APV1,Y,N,N,N,N,Y)
 SOLUTN(SOV1,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(V1,Y,Y,Y,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.
 RADIA =0. ;RADIB =0.
 ENUL =1.0E-05 ;ENUT =0.
 PRNDTL(V1)=1.
 PRT(V1)=1.
 CP1 =1. ;CP2 =1.
 ************************************************************
  Group 10.Inter-Phase Transfer Processes
 ************************************************************
  Group 11.Initial field variables (PHIs)
 FIINIT(P1)=1.0E-10 ;FIINIT(V1)=1.0E-10
 FIINIT(FLOY)=1.0E-10 ;FIINIT(RSV1)=1.0E-10
 FIINIT(ANP1)=1.0E-10 ;FIINIT(APV1)=1.0E-10
 FIINIT(SOV1)=1.0E-10 ;FIINIT(RHO1)=1.0E-10
 FIINIT(PRPS)=-1. ;FIINIT(OBID)=1.0E-10
 
 PATCH(^OB3 ,INIVAL, 3, 0, 0, 0, 0, 0, 1, 1)
 INIT(^OB3 ,PRPS,0. ,198. )
 INIADD = F
 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, 100, 1, 1, 1, 50)
 COVAL(RESIST ,V1 ,10. ,0. )
 
 PATCH(INLET ,CELL , 1, 1, 1, 1, 1, 1, 1, 50)
 COVAL(INLET ,P1 ,1. ,0. )
 
 PATCH(OUTLET ,CELL , 1, 1, 100, 100, 1, 1, 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(V1)=10
 ENDIT(P1)=1.0E-03 ;ENDIT(V1)=1.0E-03
 ************************************************************
  Group 17. Relaxation
 RELAX(P1,LINRLX,1.)
 RELAX(V1,FALSDT,1.)
 RELAX(FLOY,LINRLX,1.)
 RELAX(RSV1,LINRLX,1.)
 RELAX(ANP1,LINRLX,1.)
 RELAX(APV1,LINRLX,1.)
 RELAX(SOV1,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(V1)=1.0E+06 ;VARMIN(V1)=-1.0E+06
 VARMAX(FLOY)=1.0E+10 ;VARMIN(FLOY)=-1.0E+10
 VARMAX(RSV1)=1.0E+10 ;VARMIN(RSV1)=-1.0E+10
 VARMAX(ANP1)=1.0E+10 ;VARMIN(ANP1)=-1.0E+10
 VARMAX(APV1)=1.0E+10 ;VARMIN(APV1)=-1.0E+10
 VARMAX(SOV1)=1.0E+10 ;VARMIN(SOV1)=-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.15*YVLAST+0.666667*TIM&0.0&0&0&$)
 SPEDAT(SET,MOVOB,BOX,C,0&0&0&0&0))
 SPEDAT(SET,STORED,SOV1,C,=SORC(V1)!SOURCS)
 SPEDAT(SET,STORED,APV1,C,=APCO(V1)!EQCOEF)
 SPEDAT(SET,STORED,ANP1,C,=ANCO(P1)!EQCOEF)
 SPEDAT(SET,STORED,RSV1,C,=RESI(V1)!RESIDU)
 SPEDAT(SET,STORED,FLOY,C,=RHO1*V1!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)=2.807 ;EX(V1)=0.6666
 EX(FLOY)=0.6598 ;EX(RSV1)=7.599E-06
 EX(ANP1)=0.1535 ;EX(APV1)=1000.
 EX(SOV1)=9.851001E-04 ;EX(RHO1)=0.9999
 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(V1,Y,N,Y,Y,Y,Y)
 OUTPUT(FLOY,Y,N,Y,N,N,N)
 OUTPUT(RSV1,Y,N,Y,N,N,N)
 OUTPUT(ANP1,Y,N,Y,N,N,N)
 OUTPUT(APV1,Y,N,Y,N,N,N)
 OUTPUT(SOV1,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 = 26 ;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 = 1 ;ISTPRF = 1 ;ISTPRL = 100000
 NYPRIN = 1 ;IYPRF = 1 ;IYPRL = 10000
 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
   No PATCHes yet used for this Group
 ************************************************************
  Group 24. Dumps For Restarts
 SAVE = T ;NOWIPE = F
 NSAVE =CHAM
 IDISPA = 1 ;IDISPB = 0 ;IDISPC = 0
STOP