TALK=T;RUN(1,1)
  DISPLAY
      A piston with a flat crown having an open cavity with a
  central cusp moves inside a flat topped cylinder. The solution
  domain is represented by a 2-part moving grid, the first
  part of which expands and contracts in the Z-direction.
  The second part of the grid moves with the piston, but does
  not change its size. This part of the grid includes the
  piston cavity, which is represented by porosities. Since the
  model simulates only the compression and expansion parts of the
  four-stroke cycle, there are no inflow or outflow boundaries.
  The flow is represented in two-dimensional axi-symmetric
  polar coordinates in the Y-Z plane but with swirl velocity
  in the X-direction. An initial field of U-velocity is set up
  to represent 'solid body' rotation of the fluid. The density
  varies isentropically with pressure. This is based on the
  PHOENICS-81 PDR/CHAM UK/6.
  ENDDIS
 ************************************************************
  Group 1. Run Title and Number
 ************************************************************
 ************************************************************
 
 TEXT(Case 1 AXI-Symmetric Piston Bowl        )
 
 ************************************************************
 ************************************************************
 
 IRUNN = 1 ;LIBREF = 14
 ************************************************************
  Group 2. Time dependence
 STEADY = F
    * Set overall time and no. of steps
 TFIRST =1.157407E-04 ;TLAST =0.033681
 FSTEP = 1 ;LSTEP = 10
   Method of pairs used for grid setting.
 TFRAC(1)=-7. ;TFRAC(2)=0.034483
 TFRAC(3)=18. ;TFRAC(4)=0.017241
 TFRAC(5)=40. ;TFRAC(6)=3.448276E-03
 TFRAC(7)=4. ;TFRAC(8)=0.017241
 TFRAC(9)=7. ;TFRAC(10)=0.034483
 ************************************************************
  Group 3. X-Direction Grid Spacing
 CARTES = F
 NX = 1
 XULAST =0.1
 XFRAC(1)=1.
 ************************************************************
  Group 4. Y-Direction Grid Spacing
 NY = 11
 YVLAST =0.0689
 YFRAC(1)=0.073 ;YFRAC(2)=0.145
 YFRAC(3)=0.218 ;YFRAC(4)=0.29
 YFRAC(5)=0.363 ;YFRAC(6)=0.421
 YFRAC(7)=0.458 ;YFRAC(8)=0.508
 YFRAC(9)=0.653 ;YFRAC(10)=0.753
 YFRAC(11)=0.88632
 ************************************************************
  Group 5. Z-Direction Grid Spacing
 PARAB = F
 NZ = 11
 ZWLAST =0.1707
 ZFRAC(1)=0.227 ;ZFRAC(2)=0.518
 ZFRAC(3)=0.673 ;ZFRAC(4)=0.776
 ZFRAC(5)=0.825 ;ZFRAC(6)=0.853
 ZFRAC(7)=0.882 ;ZFRAC(8)=0.926
 ZFRAC(9)=0.953 ;ZFRAC(10)=0.982
 ZFRAC(11)=1.
 ************************************************************
  Group 6. Body-Fitted Coordinates
 ************************************************************
  Group 7. Variables: STOREd,SOLVEd,NAMEd
 ONEPHS = T
 NAME(1)=P1 ;NAME(3)=U1
 NAME(5)=V1 ;NAME(7)=W1
 NAME(149)=PRPS
    * 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,Y,N,N,N)
 SOLUTN(U1,Y,Y,N,N,N,Y)
 SOLUTN(V1,Y,Y,N,N,N,Y)
 SOLUTN(W1,Y,Y,N,N,N,Y)
 SOLUTN(PRPS,Y,N,N,N,N,N)
 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(U1,Y,Y,Y,Y,Y,Y)
 TERMS(V1,Y,Y,Y,Y,Y,Y)
 TERMS(W1,Y,Y,Y,Y,Y,Y)
 DIFCUT =0.5 ;ZDIFAC =1.
 GALA = F ;ADDDIF = F
 NEWRH1 = T
 U1AD =0. ;V1AD =0. ;W1AD = ZMOVE
 ISOLX = -1 ;ISOLY = -1 ;ISOLZ = -1
 ************************************************************
  Group 9. Properties used if PRPS is not
  stored, and where PRPS = -1.0 if it is!
 RHO1 = GRND3 ;TMP1 =0. ;EL1 =0.
 TSURR =0. ;TEMP0 =0. ;PRESS0 =1.03E+05
 DVO1DT =0. ;DRH1DP = GRND3
 RHO1A =1.4839E-04 ;RHO1B =0.765 ;RHO1C =0.
 EMISS =0. ;SCATT =0.
 RADIA =0. ;RADIB =0.
 ENUL =1.0E-05 ;ENUT =0.1
 PRNDTL(U1)=1. ;PRNDTL(V1)=1.
 PRNDTL(W1)=1.
 PRT(U1)=1. ;PRT(V1)=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(U1)=407.130005
 FIINIT(V1)=1.0E-10 ;FIINIT(W1)=1.0E-10
 FIINIT(PRPS)=-1.
 
 PATCH(CMP0 ,INIVAL, 1, 1, 8, 11, 6, 11, 1, 1)
 INIT(CMP0 ,PRPS,0. ,198. )
 INIADD = F
 FSWEEP = 1
 IURINI = -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(HEAD ,LWALL , 1, 1, 1, 11, 1, 1, 1, 10)
 COVAL(HEAD ,U1 , GRND2 ,0. )
 COVAL(HEAD ,V1 , GRND2 ,0. )
 
 PATCH(WALL ,NWALL , 1, 1, 11, 11, 1, 5, 1, 10)
 COVAL(WALL ,U1 , GRND2 ,0. )
 COVAL(WALL ,W1 , GRND2 ,0. )
 XCYCLE = F
 EGWF = T
 WALLCO = GRND2
 ************************************************************
  Group 14. Downstream Pressure For PARAB
 ************************************************************
  Group 15. Terminate Sweeps
 LSWEEP = 20 ;ISWC1 = 1
 LITHYD = 1 ;LITFLX = 1 ;LITC = 1 ;ITHC1 = 1
 SELREF = T
 RESFAC =1.0E-05
 ************************************************************
  Group 16. Terminate Iterations
 LITER(P1)=10 ;LITER(U1)=10
 LITER(V1)=10 ;LITER(W1)=10
 ENDIT(P1)=1.0E-03 ;ENDIT(U1)=1.0E-03
 ENDIT(V1)=1.0E-03 ;ENDIT(W1)=1.0E-03
 ************************************************************
  Group 17. Relaxation
 RELAX(P1,LINRLX,1.)
 RELAX(U1,FALSDT,1.)
 RELAX(V1,FALSDT,0.1)
 RELAX(W1,FALSDT,0.1)
 RELAX(PRPS,LINRLX,1.)
 OVRRLX =0.
 EXPERT = F ;NNORSL = F
 ************************************************************
  Group 18. Limits
 VARMAX(P1)=1.0E+10 ;VARMIN(P1)=5.0E+04
 VARMAX(U1)=1.0E+06 ;VARMIN(U1)=-1.0E+06
 VARMAX(V1)=1.0E+06 ;VARMIN(V1)=-1.0E+06
 VARMAX(W1)=1.0E+06 ;VARMIN(W1)=-1.0E+06
 VARMAX(PRPS)=1.0E+10 ;VARMIN(PRPS)=-1.0E+10
 ************************************************************
  Group 19. Data transmitted to GROUND
 PARSOL = F
 IZW1 = 5
 AZW1 =150.796005 ;BZW1 =0.06985
 CZW1 =3.7265 ;DZW1 =0.
 ISG62 = 1
 SPEDAT(SET,GXMONI,TRANSIENT,L,F)
 SPEDAT(SET,GXMONI,PLOTALL,L,T)
 SPEDAT(SET,MATERIAL,198,L,T)
 ************************************************************
  Group 20. Preliminary Printout
 ************************************************************
  Group 21. Print-out of Variables
 INIFLD = F ;SUBWGR = T
    * 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(U1,Y,N,Y,Y,Y,Y)
 OUTPUT(V1,Y,N,Y,Y,Y,Y)
 OUTPUT(W1,Y,N,Y,Y,Y,Y)
 OUTPUT(PRPS,Y,N,Y,N,N,N)
 ************************************************************
  Group 22. Monitor Print-Out
 IXMON = 1 ;IYMON = 1 ;IZMON = 5
 NPRMON = 100000 ;NPRMNT = 1 ;TSTSWP = 1
 UWATCH = T ;USTEER = T
 HIGHLO = F
 ************************************************************
  Group 23.Field Print-Out & Plot Control
 NPRINT = 100000 ;NUMCLS = 5
 NTPRIN = 100000 ;ISTPRF = 1 ;ISTPRL = 100000
 NYPRIN = -1 ;IYPRF = 1 ;IYPRL = 10000
 NZPRIN = -1 ;IZPRF = 1 ;IZPRL = 10000
 XZPR = F ;YZPR = F
 IPLTF = 1 ;IPLTL = -1 ;NPLT = -1
 ISWPRF = 1 ;ISWPRL = 100000
 ITABL = 2 ;IPROF = 1
 ABSIZ =0.5 ;ORSIZ =0.4
 NTZPRF = 1 ;NCOLPF = 50
 ICHR = 2 ;NCOLCO = 45 ;NROWCO = 20
 
 PATCH(PROF ,PROFIL, 1, 1, 1, 11, 3, 3, 10, 10)
 PLOT(PROF ,U1 ,0. ,0. )
 PLOT(PROF ,V1 ,0. ,0. )
 PLOT(PROF ,W1 ,0. ,0. )
 
 PATCH(CONT ,CONTUR, 1, 1, 1, 11, 1, 11, 10, 10)
 PLOT(CONT ,U1 ,0. ,20. )
 PLOT(CONT ,W1 ,0. ,20. )
 ************************************************************
  Group 24. Dumps For Restarts
 SAVE = T ;NOWIPE = F
 NSAVE =CHAM
STOP