************************************************************
  Group 1. Run Title
 TEXT(INCLINED JET-PROP FLOW OVER FLAT BED :P113
 TITLE
 ************************************************************
  Groups 3, 4, 5  Grid Information
    * Overall number of cells, RSET(M,NX,NY,NZ,tolerance)
 GRDPWR(X,1,6.283E-02,1.0)
 GRDPWR(Y,12,6.0,1.0)
 NREGZ=2
 IREGZ=1;GRDPWR(Z,6,4.0,1.0)
 IREGZ=2;GRDPWR(Z,1,1.0,1.0)
    * Cylindrical-polar grid
    **Not really correct, but at least represents spreading.
 CARTES=F
 ************************************************************
  Group 7. Variables: STOREd,SOLVEd,NAMEd
 ONEPHS  =    T
    * Non-default variable names
 NAME(C11)=SPT1 ; NAME(C21)=WPT1
char(setaslp); setaslp=$P001
#setaslp
    * Stored variables list
 STORE(DEN1,VISL,IMB1,PCOR)
    **Include special variables for sand model.
    **(Remember to add EPOR for 3-d cases).
 STORE(SPT1,PRPS)
    *Store particle vertical velocities at cell "low" faces,
      derived in the local version of GXASLP.
 STORE(WPT1,VPOR)
 SOLUTN(P1  ,Y,Y,Y,N,N,Y)
 ************************************************************
  Group 8. Terms & Devices
    **Include terms for ASM.
 NEWENT  =    T
 ************************************************************
  Group 9. Properties
    **Include special density terms for ASM.
 RHO2 = 1.026E+03
 PRESS0  = 1.000E+06 ; TEMP0   = 2.730E+02
    **Set density information for ASM.
 PHINT(PT0) = RHO2
 PHINT(PT1) = 2.400E+03
    **Set particle diameters for ASM.
 CINT(PT0) = 1.000E-10
 CINT(PT1) = 7.000E-05
     **Set viscosities for ASM.
 PRNDTL(PT0) = 1.350E-06
 PRNDTL(PT1) = 1.350E-04
 ************************************************************
  Group 11.Initialise Var/Porosity Fields
 FIINIT(DEN1) = 1.026E+03 ; FIINIT(VISL) = 1.350E-06
 FIINIT(PRPS) = 0.000E+00 ; FIINIT(SPT1) = 0.000E+00
 FIINIT(PT0)  = 1.000E+00 ; FIINIT(PT1)  = 0.000E+00
 
 INIADD  =    F
 PATCH (SEA     ,INIVAL,1,1,1,NY,1,NZ-1,1,1)
    *ASM needs some particles present to start with.
 COVAL (SEA     ,PT0 ,FIXVAL,1.0-1.0E-06)
 COVAL (SEA     ,PT1 ,FIXVAL,1.000E-06)
 
 PATCH (BED     ,INIVAL,1,1,1,NY,NZ,NZ,1,1)
 COVAL (BED     ,SPT1,FIXVAL,1.000E-01)
 COVAL (BED     ,PRPS,FIXVAL,1.990E+02)
 WALLA = 1.000E-03
 ************************************************************
  Group 13. Boundary & Special Sources
 
 INLET (ASM_PROP,LOW   ,1,1,3,4,1,1,1,1)
 VALUE (ASM_PROP,P1  , RHO2*2.828)
 VALUE (ASM_PROP,W1  , 2.828E+00)
 VALUE (ASM_PROP,V1  , 2.828E+00)
 VALUE (ASM_PROP,PT0 , 1.000E+00)
 VALUE (ASM_PROP,PT1 , 0.000E+00)
 VALUE (ASM_PROP,VFOL, 1.0/RHO2)
 
 PATCH (ASM_TOP ,LOW   ,1,1,1,2,1,1,1,1)
 COVAL (ASM_TOP ,P1  , 1.00E+01, 0.000E+00)
 COVAL (ASM_TOP ,V1  , ONLYMS  , 0.000E+00)
 COVAL (ASM_TOP ,W1  , ONLYMS  , 0.000E+00)
 COVAL (ASM_TOP ,PT0 , ONLYMS  , 1.000E+00)
 COVAL (ASM_TOP ,PT1 , ONLYMS  , 0.000E+00)
 COVAL (ASM_TOP ,VFOL, ONLYMS  , 1.0/RHO2)
 
 PATCH (ASM_TOP2 ,LOW   ,1,1,5,NY,1,1,1,1)
 COVAL (ASM_TOP2 ,P1  , 1.00E+01, 0.000E+00)
 COVAL (ASM_TOP2 ,V1  , ONLYMS  , 0.000E+00)
 COVAL (ASM_TOP2 ,W1  , ONLYMS  , 0.000E+00)
 COVAL (ASM_TOP2 ,PT0 , ONLYMS  , 1.000E+00)
 COVAL (ASM_TOP2 ,PT1 , ONLYMS  , 0.000E+00)
 COVAL (ASM_TOP2 ,VFOL, ONLYMS  , 1.0/RHO2)
 
 PATCH (ASM_EDG2,NORTH ,1,1,NY,NY,1,NZ-1,1,1)
 COVAL (ASM_EDG2,P1  , 1.00E+03, 0.000E+00)
 COVAL (ASM_EDG2,W1  , ONLYMS  , 0.000E+00)
 COVAL (ASM_EDG2,V1  , ONLYMS  , 0.000E+00)
 COVAL (ASM_EDG2,PT0 , ONLYMS  , 1.000E+00)
 COVAL (ASM_EDG2,PT1 , ONLYMS  , 0.000E+00)
 COVAL (ASM_EDG2,VFOL, ONLYMS  , 1.0/RHO2)
 ************************************************************
  Group 15. Terminate Sweeps
 LSWEEP  =    100 ; SELREF  =    T ; RESFAC  = 1.000E-06
 ************************************************************
  Group 16. Terminate Iterations
 ENDIT (P1  ) =  1.000E-03 ;ENDIT (V1  ) =  1.000E-03
 ENDIT (W1  ) =  1.000E-03
   **Special limits for ASM.
 LITER(PT0)=5              ; LITER(PT1)=10
 ENDIT(PT0)=1E-8           ; ENDIT(PT1)=1E-6
 ************************************************************
  Group 17. Relaxation
 RELAX(P1  ,LINRLX, 6.000E-01)
 RELAX(V1  ,FALSDT, 1.000E+01)
 RELAX(W1  ,FALSDT, 5.000E+00)
 RELAX(PT0 ,LINRLX, 5.000E-01)
 RELAX(PT1 ,LINRLX, 3.000E-01)
 RELAX(VFOL,LINRLX, 0.5)
 ************************************************************
  Group 18. Limits
 VARMIN(P1  )  = -1.000E+06 ; VARMAX(P1  )  =  1.000E+07
 ************************************************************
  Group 19. EARTH Calls To GROUND Station
 GENK    =    T
   **Activate ASM
 LASLPB  =    T ; LASLPA  =    T
 ************************************************************
  Group 22. Monitor Print-Out
 IXMON   =       1 ;IYMON  =      2 ;IZMON  =      5
 TSTSWP  =      -1
 ************************************************************
  Group 23.Field Print-Out & Plot Control
   No PATCHes used for this Group
 NZPRIN  =      1   ; IZPRF =      1 ; IZPRL =   NZ
 NYPRIN  =      1   ; IYPRF =      1 ; IYPRL =   NY
 ************************************************************
  Group 24. Dumps For Restarts