TALK=F;RUN( 1, 1)
 
 ************************************************************
   Q1 created by VDI menu, Version 2007, Date 26/09/07
 CPVNAM=VDI;SPPNAM=Core
 ************************************************************
  Echo DISPLAY / USE settings
  DISPLAY
  The steep propagating 1DY wave.
  ENDDIS
  PHOTON USE
  p;parphi
  1 1 4000
 
  clear
  set prop off
  view x
  gr ou x 1
  con disy x 1 fi;0.0001
  msg DISY(Y,Z=Time)
  pause
 
 
  autoplot
  FILE
  parphi 5
 
  clear
  data 1
  DISY Z 10 Y 1 m
  col1 1
  data 1
  DISY Z 20 Y 1 m
  colA 2
  data 1
  DISY Z 40 Y 1 m
  colF 3
  scale
  plot 1 2 3
  msg DISY(Y) for Time=10(W),20(Y),40(R)
  pause
 
  clear
  data 1
  DISY Y 50 Z 1 m
  col1 1
  scale
  plot 1
  msg DISY(Time) for Y=LY/2
  pause
 
  ENDUSE
 ************************************************************
 IRUNN   =         1 ;LIBREF =       674
 ************************************************************
  Group 1. Run Title
 TEXT(The steep propagating 1DY wave; S674    )
 ************************************************************
  Group 2. Transience
 STEADY=F
    * Set overall time and no. of steps
 RSET(U,0.000000E+00,4.740000E-05,100)
    * Modify regions
 RSET(T,1,100,1.000000E+00)
 ************************************************************
  Groups 3, 4, 5  Grid Information
    * Overall number of cells, RSET(M,NX,NY,NZ,tolerance)
 RSET(M,1,120,1,1.000000E-05)
 ************************************************************
  Group 6. Body-Fitted coordinates
 ************************************************************
  Group 7. Variables: STOREd,SOLVEd,NAMEd
 ONEPHS  =    T
    * Non-default variable names
 NAME(141) =EPSZ ; NAME(142) =EPSX
 NAME(143) =EPSY ; NAME(144) =DRH1
 NAME(145) =DVO1 ; NAME(146) =ENUL
 NAME(147) =DEN1 ; NAME(148) =PRPS
 NAME(149) =VDSY ; NAME(150) =DISY
    * Solved variables list
 SOLVE(DISY)
    * Stored variables list
 STORE(VDSY,PRPS,DEN1,ENUL,DVO1,DRH1,EPSY,EPSX)
 STORE(EPSZ)
 
 ************************************************************
  Echo InForm settings for Group  7
  INFORM7BEGIN
REAL(LY,POISSON,YOUNG)
LY=120.e-3
YOUNG   = 1/0.5E-11
POISSON=0.3
INTEGER(NYBODY)
NYBODY = 120
 
real(D0,SSound,V0,DY00, TimeMax)
D0 = 1.e-6
SSound = (YOUNG/7800.0)**0.5   !  = 5.06E+03
DY00 = LY/NYBODY
V0 = SSound*D0/(2.0*DY00)
TimeMax = 2.0*LY/SSound        ! = 4.74E-05
  INFORM7END
 ************************************************************
  Group 8. Terms & Devices
 TERMS (DISY,N,N,Y,Y,Y,N)
 ************************************************************
  Group 9. Properties
    * Domain material index is 160 signifying:
    * NOTSET
 SETPRPS(1,160,Q1)
 ENUT    = 0.000000E+00
 DRH1DP  = 5.000000E-12
 DVO1DT  = 1.000000E-05
   *  List of user-defined materials to be read by EARTH
    MATFLG=T;IMAT=1
   *  Name
   *Ind. Dens.  Viscos.  Spec.heat  Conduct.  Expans.  Compr.
   *  
    160 7800.0 0.3 473.0 43.0 1.0E-5 0.5E-11
 ************************************************************
  Group 10.Inter-Phase Transfer Processes
 ************************************************************
  Group 11.Initialise Var/Porosity Fields
 FIINIT(EPSZ) =  1.001000E-10 ;FIINIT(EPSX) =  1.001000E-10
 FIINIT(EPSY) =  1.001000E-10 ;FIINIT(DRH1) =  1.001000E-10
 FIINIT(DVO1) =  1.001000E-10 ;FIINIT(ENUL) =  1.001000E-10
 FIINIT(DEN1) =  1.001000E-10 ;FIINIT(PRPS) =  1.600000E+02
 FIINIT(VDSY) =  0.000000E+00 ;FIINIT(DISY) =  0.000000E+00
 
 PATCH (INBL10  ,INIVAL,1,0,0,0,0,0,1,1)
 
 PATCH (INBL11  ,INIVAL,2,0,0,0,0,0,1,1)
 
 PATCH (INBL12  ,INIVAL,3,0,0,0,0,0,1,1)
 
 
 INIADD  =    F
 ************************************************************
  Echo InForm settings for Group 11
  INFORM11BEGIN
(INITIAL of DISY at InBl10 is :D0:)
(INITIAL of DISY at InBl11 is :D0:/2.0)
(INITIAL of VDSY at InBl11 is :V0:)
(INITIAL of VDSY at InBl12 is :V0:)
  INFORM11END
 ************************************************************
  Group 12. Convection and diffusion adjustments
   No PATCHes used for this Group
 ************************************************************
  Group 13. Boundary & Special Sources
 
 PATCH (TOP     ,NWALL ,5,0,0,0,0,0,1,100)
 COVAL (TOP     ,DISY, 1.000000E+00, 0.000000E+00)
 PATCH (BOTTOM     ,SWALL ,4,0,0,0,0,0,1,100)
 
 EGWF    =    T
 ************************************************************
  Echo InForm settings for Group 13
  INFORM13BEGIN
(SOURCE of DISY at Bottom is COVAL(1,:D0:))
  INFORM13END
 ************************************************************
  Group 14. Downstream Pressure For PARAB
 ************************************************************
  Group 15. Terminate Sweeps
 LSWEEP  =        60
 RESFAC  = 1.000000E-03
 ************************************************************
  Group 16. Terminate Iterations
 ************************************************************
  Group 17. Relaxation
 RELAX(DISY,LINRLX, 1.000000E+00)
 SPEDAT(SET,GXMONI,TRANSIENT,L,F)
 ************************************************************
  Group 18. Limits
 ************************************************************
  Group 19. EARTH Calls To GROUND Station
 USEGRD  =    T  ;USEGRX =    T
 STRA    =    T
 ASAP    =    T
 PARSOL  =    F
 ISG21   =        60
 ISG52   =         3
 CSG10   ='q1'
 ************************************************************
  Group 20. Preliminary Printout
 ECHO    =    T
 ************************************************************
  Group 21. Print-out of Variables
 ************************************************************
  Group 22. Monitor Print-Out
 NPRMON  =    100000
 NPRMNT  =         1
 TSTSWP  =        -1
 ************************************************************
  Group 23.Field Print-Out & Plot Control
 NPRINT  =    100000
 NTPRIN  =    100000 ;ISTPRF =         1 ;ISTPRL =    100000
 ISWPRF  =         1 ;ISWPRL =    100000
   No PATCHes used for this Group
 ************************************************************
  Group 24. Dumps For Restarts
 IDISPA  =         1 ;IDISPB =         0 ;IDISPC =         0
 
 GVIEW(P,-9.988757E-01,1.449645E-02,-4.513222E-02)
 GVIEW(UP,1.501839E-02,9.998240E-01,-1.124712E-02)
 
> DOM,    SIZE,        4.000000E-02, 1.200000E-01, 4.000000E-02
> DOM,    MONIT,       0.000000E+00, 5.000000E-04, 2.000000E-02
> DOM,    SCALE,       1.000000E+00, 1.000000E+00, 1.000000E+00
> DOM,    INCREMENT,   1.000000E-02, 1.000000E-02, 1.000000E-02
  > GRID,   RSET_X_1,      1, 1.000000E+00
> GRID,   RSET_Y_1,     10, 1.000000E+00
> GRID,   RSET_Y_2,      1, 1.000000E+00
> GRID,   RSET_Y_3,      1, 1.000000E+00
> GRID,   RSET_Y_4,    108, 1.000000E+00
  > GRID,   RSET_Z_1,      1, 1.000000E+00
 
> OBJ,    NAME,        INBL10
> OBJ,    POSITION,    0.000000E+00, 0.000000E+00, 0.000000E+00
> OBJ,    SIZE,        4.000000E-02, 1.000000E-02, 4.000000E-02
> OBJ,    GEOMETRY,    default
> OBJ,    ROTATION24,        1
> OBJ,    TYPE,        USER_DEFINED
 
> OBJ,    NAME,        INBL11
> OBJ,    POSITION,    0.000000E+00, 1.000000E-02, 0.000000E+00
> OBJ,    SIZE,        4.000000E-02, 1.000001E-03, 4.000000E-02
> OBJ,    GEOMETRY,    default
> OBJ,    ROTATION24,        1
> OBJ,    TYPE,        USER_DEFINED
 
> OBJ,    NAME,        INBL12
> OBJ,    POSITION,    0.000000E+00, 1.100000E-02, 0.000000E+00
> OBJ,    SIZE,        4.000000E-02, 1.000000E-03, 4.000000E-02
> OBJ,    GEOMETRY,    default
> OBJ,    ROTATION24,        1
> OBJ,    TYPE,        USER_DEFINED
 
> OBJ,    NAME,        BOTTOM
> OBJ,    POSITION,    0.000000E+00, 0.000000E+00, 0.000000E+00
> OBJ,    SIZE,        4.000000E-02, 0.000000E+00, 4.000000E-02
> OBJ,    GEOMETRY,    default
> OBJ,    ROTATION24,        1
> OBJ,    TYPE,        USER_DEFINED
 
> OBJ,    NAME,        TOP
> OBJ,    POSITION,    0.000000E+00, 1.200000E-01, 0.000000E+00
> OBJ,    SIZE,        4.000000E-02, 0.000000E+00, 4.000000E-02
> OBJ,    GEOMETRY,    default
> OBJ,    ROTATION24,        1
> OBJ,    TYPE,        USER_DEFINED
 ************************************************************
STOP