```
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
************************************************************
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

************************************************************
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
```