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
  A rectangular block of dimensions LX * LY * LZ is
  * held fixed at it north, west and low faces
  * subjected to uniform compression at its south face
  * allowed to contract freely along its east and high faces
 
  The distributions of displacement in the xy plane are computed,
  and also the corresponding stresses and strains in all three
  directions.
  A uniform cartesian grid is used.
  The computed results are compared with the analytically-derived
  values.
 
  ENDDIS
  photon use
  p;;;;
 
  set prop off
  con v1 x 9 y 1 9 z 1 9 fi;0.001
  con v1 y 9 x 1 9 z 1 9 fi;0.001
  con v1 z 9 y 1 9 x 1 9 fi;0.001
  pause
  con u1 x 9 y 1 9 z 1 9 fi;0.001
  con u1 y 9 x 1 9 z 1 9 fi;0.001
  con u1 z 9 y 1 9 x 1 9 fi;0.001
  enduse
 ************************************************************
 IRUNN   =         1 ;LIBREF =         303
 ************************************************************
  Group 1. Run Title
 TEXT(3D block in y-direction compression;S303)
 ************************************************************
  Group 2. Transience
 STEADY  =    T
 ************************************************************
  Groups 3, 4, 5  Grid Information
    * Overall number of cells, RSET(M,NX,NY,NZ,tolerance)
 RSET(M,10,10,10,1.000000E-05)
 ************************************************************
  Group 6. Body-Fitted coordinates
 ************************************************************
  Group 7. Variables: STOREd,SOLVEd,NAMEd
 ONEPHS  =    T
    * Non-default variable names
    * Solved variables list
 SOLVE(P1,U1,V1,W1)
    * Stored variables list
 STORE(PRPS,DRH1,VISL,STRX,STRY,STRZ,STXY,STXZ)
 STORE(STYZ,EPSY,EPSX,EPSZ,U1T ,V1T ,W1T ,U1/T)
 STORE(V1/T,W1/T)
    * Additional solver options
 SOLUTN(U1,Y,Y,Y,N,N,Y)
 SOLUTN(V1,Y,Y,Y,N,N,Y)
 SOLUTN(W1,Y,Y,Y,N,N,Y)
 
 ************************************************************
  Echo InForm settings for Group  7
  INFORM7BEGIN
  Declarations and settings
REAL(FY,POISSON,YOUNG,LY)
LY = 90.e-3
FY= 40.0e6  ! H/m^2 = 40 N/mm^2
YOUNG   = 1/0.5E-11
POISSON=0.3
real(CEPSX,CEPSY,CEPSZ,CEPS)
CEPSX = FY/YOUNG*POISSON
CEPSY =  -FY/YOUNG
CEPSZ =  CEPSX
CEPS = CEPSX+CEPSY+CEPSZ
   **** CALCULATE analytical solution ***
(STORED VAR W1T IS :CEPSZ:*ZW)
(STORED VAR U1T IS :CEPSX:*XU)
(STORED VAR V1T IS :CEPSY:*(YV -:LY:))
 
(STORED VAR U1/T IS U1/(U1T+1.e-20))
(STORED VAR V1/T IS V1/(V1T+1.e-20))
(STORED VAR W1/T IS W1/(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(STYZ) =  1.001000E-10
 FIINIT(STXZ) =  1.001000E-10 ;FIINIT(STXY) =  1.001000E-10
 FIINIT(STRZ) =  1.001000E-10 ;FIINIT(STRY) =  1.001000E-10
 FIINIT(STRX) =  1.001000E-10 ;FIINIT(VISL) =  1.001000E-10
 FIINIT(DRH1) =  1.001000E-10 ;FIINIT(PRPS) =  1.600000E+02
 FIINIT(U1) =  0.000000E+00 ;FIINIT(V1) =  0.000000E+00
 FIINIT(W1) =  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 (UP      ,NWALL ,1,0,0,0,0,0,1,1)
 COVAL (UP      ,V1, 1.000000E+00, 0.000000E+00)
 
 PATCH (FORS01  ,SOUTH ,2,0,0,0,0,0,1,1)
 COVAL (FORS01  ,V1, FIXFLU      , 4.000000E+07)
 
 PATCH (AXESXX  ,WWALL ,3,0,0,0,0,0,1,1)
 COVAL (AXESXX  ,U1, 1.000000E+00, 0.000000E+00)
 
 PATCH (AXESZZ  ,LWALL ,4,0,0,0,0,0,1,1)
 COVAL (AXESZZ  ,W1, 1.000000E+00, 0.000000E+00)
 
 EGWF    =    T
 ************************************************************
  Group 14. Downstream Pressure For PARAB
 ************************************************************
  Group 15. Terminate Sweeps
 LSWEEP  =       100
 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   =       100
 ISG50   =         1
 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   =         8 ;IYMON  =         2 ;IZMON  =         8
 NPRMON  =    100000
 NPRMNT  =         1
 TSTSWP  =        -1
 ************************************************************
  Group 23.Field Print-Out & Plot Control
 NPRINT  =    100000
 ISWPRF  =         1 ;ISWPRL =    100000
   No PATCHes used for this Group
 ************************************************************
  Group 24. Dumps For Restarts
 
 GVIEW(P,6.377395E-01,5.325049E-01,5.565315E-01)
 GVIEW(UP,-4.263606E-01,8.457928E-01,-3.207041E-01)
 
> DOM,    SIZE,        9.000000E-02, 9.000000E-02, 1.200000E-01
> DOM,    MONIT,       6.750000E-02, 1.350000E-02, 9.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,     10, 1.000000E+00
  > GRID,   RSET_Y_1,     10, 1.000000E+00
  > GRID,   RSET_Z_1,     10, 1.000000E+00
 
> OBJ,    NAME,        UP
> OBJ,    POSITION,    0.000000E+00, 9.000000E-02, 0.000000E+00
> OBJ,    SIZE,        9.000000E-02, 0.000000E+00, 1.200000E-01
> OBJ,    GEOMETRY,    default
> OBJ,    ROTATION24,        1
> OBJ,    TYPE,        USER_DEFINED
 
> OBJ,    NAME,        FORS01
> OBJ,    POSITION,    0.000000E+00, 0.000000E+00, 0.000000E+00
> OBJ,    SIZE,        9.000000E-02, 0.000000E+00, 1.200000E-01
> OBJ,    GEOMETRY,    default
> OBJ,    ROTATION24,        1
> OBJ,    TYPE,        USER_DEFINED
 
> OBJ,    NAME,        AXESXX
> OBJ,    POSITION,    0.000000E+00, 0.000000E+00, 0.000000E+00
> OBJ,    SIZE,        0.000000E+00, 9.000000E-02, 1.200000E-01
> OBJ,    GEOMETRY,    default
> OBJ,    ROTATION24,        1
> OBJ,    TYPE,        USER_DEFINED
 
> OBJ,    NAME,        AXESZZ
> OBJ,    POSITION,    0.000000E+00, 0.000000E+00, 0.000000E+00
> OBJ,    SIZE,        9.000000E-02, 9.000000E-02, 0.000000E+00
> OBJ,    GEOMETRY,    default
> OBJ,    ROTATION24,        1
> OBJ,    TYPE,        USER_DEFINED
 ************************************************************
STOP