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
  Problem: Vertical prismatic (0.3x03x2.5 m) bar & gravity force
           -Rho*g.
           At the upper end - a uniform normal stress equal to
           total weight of the bar Rho*g*LZ.
   1. Comparison - a analytical solution:
      U = -P*Rho*g*x*z/E,
      V = -P*Rho*g*z*y/E,
      W = Rho*g/E/2*[ (z^2 - LZ^2) + P*(x^2+y^2)]
 
  ENDDIS
 ************************************************************
 IRUNN   =         1 ;LIBREF =       701
 ************************************************************
  Group 1. Run Title
 TEXT(3D Prismatic bar & gravity force; s701  )
 ************************************************************
  Group 2. Transience
 STEADY  =    T
 ************************************************************
  Groups 3, 4, 5  Grid Information
    * Overall number of cells, RSET(M,NX,NY,NZ,tolerance)
 RSET(M,6,6,6)
 ************************************************************
  Group 6. Body-Fitted coordinates
 ************************************************************
  Group 7. Variables: STOREd,SOLVEd,NAMEd
 ONEPHS  =    T
    * Non-default variable names
 NAME(136) =W1/T ; NAME(137) =V1/T
 NAME(138) =U1/T ; NAME(139) =W1T
 NAME(140) =V1T  ; NAME(141) =U1T
 NAME(142) =EPSZ ; NAME(143) =EPSX
 NAME(144) =EPSY ; NAME(145) =VISL
 NAME(146) =DRH1 ; NAME(147) =PRPS
 NAME(148) =DISZ ; NAME(149) =DISY
 NAME(150) =DISX
    * Solved variables list
 SOLVE(DISZ,DISY,DISX)
    * Stored variables list
 STORE(PRPS,DRH1,VISL,EPSY,EPSX,EPSZ,U1T ,V1T )
 STORE(W1T ,U1/T,V1/T,W1/T)
 
 ************************************************************
  Echo InForm settings for Group  7
  inform7begin
 REAL(Rho,GF,FZ,LX,LY,LZ,POISSON,YOUNG)
 Rho = 1.e4
 GF=9.8
 FZ= Rho*GF
 LX=0.15
 LZ=2.5
 LY=0.15
 YOUNG   = 1/0.5E-11   ! Young's modulus
 POISSON=0.3           ! Poisson's ratio
 
  ***** analytical solution *****
 real(CC00,CC11)
 CC00 = FZ/YOUNG*POISSON/2
 char(FrZ1,FrZ2)
 CC11 = FZ/YOUNG/2
 FrZ1=:CC11:*(ZG^2-:LZ:^2)+
 FrZ2=:CC00:*((XG-0.5*:LX:)^2+(YG-0.5*:LY:)^2)
 
 real(CC22)
 CC22 = -FZ/YOUNG*POISSON
 char(FrX,FrY)
 FrX=:CC22:*ZG*(XG-0.5*:LX:)
 FrY=:CC22:*ZG*(YG-0.5*:LY:)
 
(STORED VAR W1T IS :FrZ1::FrZ2:)
 
(STORED VAR U1T IS :frX:)
(STORED VAR V1T IS :frY:)
 
(STORED VAR U1/T IS DISX/(U1T+1.e-20))
(STORED VAR V1/T IS DISY/(V1T+1.e-20))
(STORED VAR W1/T IS DISZ/(W1T+1.e-20))
  inform7end
 ************************************************************
  Group 8. Terms & Devices
 ************************************************************
  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(W1/T) =  1.001000E-10 ;FIINIT(V1/T) =  1.001000E-10
 FIINIT(U1/T) =  1.001000E-10 ;FIINIT(W1T ) =  1.001000E-10
 FIINIT(V1T ) =  1.001000E-10 ;FIINIT(U1T ) =  1.001000E-10
 FIINIT(EPSZ) =  1.001000E-10 ;FIINIT(EPSX) =  1.001000E-10
 FIINIT(EPSY) =  1.001000E-10 ;FIINIT(VISL) =  1.001000E-10
 FIINIT(DRH1) =  1.001000E-10 ;FIINIT(PRPS) =  0.000000E+00
 FIINIT(DISZ) =  0.000000E+00 ;FIINIT(DISY) =  0.000000E+00
 FIINIT(DISX) =  0.000000E+00
   No PATCHes used for this Group
 
 
 INIADD  =    F
 ************************************************************
  Group 12. Convection and diffusion adjustments
   No PATCHes used for this Group
 ************************************************************
  Group 13. Boundary & Special Sources
 
 PATCH (FORS01  ,HIGH  ,1,0,0,0,0,0,1,1)
 COVAL (FORS01  ,DISZ, FIXFLU      , 2.450000E+05)
 
 PATCH (ONEUVW  ,VOLUME,2,0,0,0,0,0,1,1)
 
 EGWF    =    T
 ************************************************************
  Echo InForm settings for Group 13
  inform13begin
 (source of DISZ at ONEUVW is COVAL(1.0e5,W1T))
 (source of DISX at ONEUVW is COVAL(1.0e5,U1T))
 (source of DISY at ONEUVW is COVAL(1.0e5,V1T))
 
(source of DISZ at BODY is COVAL(FIXFLU,-9.800000E+04) with VOLU)

  inform13end
 ************************************************************
  Group 14. Downstream Pressure For PARAB
 ************************************************************
  Group 15. Terminate Sweeps
 LSWEEP  =       300
 RESFAC  = 1.000000E-03
 ************************************************************
  Group 16. Terminate Iterations
 ************************************************************
  Group 17. Relaxation
 ************************************************************
  Group 18. Limits
 ************************************************************
  Group 19. EARTH Calls To GROUND Station
 USEGRD  =    T  ;USEGRX =    T
 STRA    =    T
 ASAP    =    T
 PARSOL  =    F
 ISG21   =       300
 ISG52   =         3
 CSG10   ='q1'
 ************************************************************
  Group 20. Preliminary Printout
 ECHO    =    T
 ************************************************************
  Group 21. Print-out of Variables
 OUTPUT(VISL,N,N,N,N,N,N)
 OUTPUT(DRH1,N,N,N,N,N,N)
 OUTPUT(PRPS,N,N,N,N,N,N)
 ************************************************************
  Group 22. Monitor Print-Out
 IXMON   =         4 ;IYMON  =         4 ;IZMON  =         4
 NPRMON  =    100000
 NPRMNT  =         1
 TSTSWP  =        -1
 ************************************************************
  Group 23.Field Print-Out & Plot Control
 NPRINT  =    100000
 NXPRIN  =         1
 NYPRIN  =         1
 NZPRIN  =         1
 ISWPRF  =         1 ;ISWPRL =    100000
   No PATCHes used for this Group
 ************************************************************
  Group 24. Dumps For Restarts
 
 GVIEW(P,3.063364E-01,6.317725E-01,-7.120544E-01)
 GVIEW(UP,-8.002141E-01,5.760422E-01,1.668313E-01)
 
> DOM,    SIZE,        1.500000E-01, 1.500000E-01, 2.500000E+00
> DOM,    MONIT,       8.750000E-02, 8.750000E-02, 1.458330E+00
> DOM,    SCALE,       1.000000E+00, 1.000000E+00, 1.000000E+00
> DOM,    SNAPSIZE,    1.000000E-02
  > GRID,   RSET_X_1,      6, 1.000000E+00
  > GRID,   RSET_Y_1,      6, 1.000000E+00
  > GRID,   RSET_Z_1,      1, 1.000000E+00
  > GRID,   RSET_Z_2,      5, 1.000000E+00
 
> OBJ,    NAME,        FORS01
> OBJ,    POSITION,    0.000000E+00, 0.000000E+00, 2.500000E+00
> OBJ,    SIZE,        1.500000E-01, 1.500000E-01, 0.000000E+00
> OBJ,    GEOMETRY,    default
> OBJ,    ROTATION24,        1
> OBJ,    TYPE,        USER_DEFINED
 
> OBJ,    NAME,        ONEUVW
> OBJ,    POSITION,    0.000000E+00, 0.000000E+00, 0.000000E+00
> OBJ,    SIZE,        1.500000E-01, 1.500000E-01, 4.166667E-01
> OBJ,    GEOMETRY,    default
> OBJ,    ROTATION24,        1
> OBJ,    TYPE,        USER_DEFINED
 
> OBJ,    NAME,        BODY
> OBJ,    POSITION,    0.000000E+00, 0.000000E+00, 0.000000E+00
> OBJ,    SIZE,        1.500000E-01, 1.500000E-01, 2.500000E+00
> OBJ,    GEOMETRY,    cube14
> OBJ,    ROTATION24,        1
> OBJ,    TYPE,        BLOCKAGE
> OBJ,    COLOR-MODE,  USER
> OBJ,    COLOR-VAL,          8
> OBJ,    MATERIAL,    160,NOTSET
 ************************************************************
STOP