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 hollow, thick-walled cylinder, immersed in an outer fluid, 
  contains a second fluid having a different pressure.
  Recirculating square-wave cavity is made in internal surface of 
  the cylinder. Cylinder is cooled to temperature -30. 
     No analytical solution exists; but the computed displacements,
  strains and stresses are in close agreement with those computed
  by the finite-element code ELCUT.
 
  ENDDIS
  PHOTON USE
  p;;;;
  
  
  set prop off
  msg
  msg V1 contours
  msg
  cont V1 x 1 y 1 21 z 2 19 fil;.0001
  vec x 1 y 2 21 z 2 19 col 1
  dump
  S225_V
  
  pause
  
  cl
  msg
  msg W1 contours
  msg
  cont W1 x 1 y 2 21 y 1 19 fil;.0001
  vec x 1 y 2 21 z 2 19 col 1
  dump
  S225_W


  pause
 
 
  ENDUSE
 ************************************************************
 IRUNN   =         1 ;LIBREF =         225
 ************************************************************
  Group 1. Run Title
 TEXT(Pressurised cylinder with Hole and not zero TEM)
 ************************************************************
  Group 2. Transience
 STEADY  =    T
 ************************************************************
  Groups 3, 4, 5  Grid Information
    * Overall number of cells, RSET(M,NX,NY,NZ,tolerance)
 RSET(M,1,22,20,8.000016E-04)
    * Cylindrical-polar grid
 CARTES=F
 ************************************************************
  Group 6. Body-Fitted coordinates
 ************************************************************
  Group 7. Variables: STOREd,SOLVEd,NAMEd
 ONEPHS  =    T
    * Solved variables list
 SOLVE(P1  ,V1  ,W1  )
    * Stored variables list
 STORE(P1TH,PRPS,DEN1,ENUL,DVO1,DRH1,STRX,EPSX)
 STORE(STRY,EPSY,STRZ,EPSZ)
 STORE(TEM1,EPST)
 
 ************************************************************
  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=2
   *  Name
   *Ind. Dens.  Viscos.  Spec.heat  Conduct.  Expans.  Compr.
   *  
    160 7800.0 0.3 473.0 43.0 1.0E-5 0.5E-11
   *  
    68 1.0 1.E+20 1000.0 0.0 0 0
 ************************************************************
  Group 10.Inter-Phase Transfer Processes
 ************************************************************
  Group 11.Initialise Var/Porosity Fields
 FIINIT(EPSZ) =  1.001000E-10 ;FIINIT(STRZ) =  1.001000E-10
 FIINIT(EPSY) =  1.001000E-10 ;FIINIT(STRY) =  1.001000E-10
 FIINIT(EPSX) =  1.001000E-10 ;FIINIT(STRX) =  1.001000E-10
 FIINIT(DRH1) =  1.001000E-10 ;FIINIT(DVO1) =  1.001000E-10
 FIINIT(DEN1) =  1.001000E-10 ;FIINIT(PRPS) =  1.600000E+02
 FIINIT(P1TH) =  1.001000E-10
   No PATCHes used for this Group

   inform13begin
real(DTE_R,TE_R)
DTE_R = 0.0
TE_R = -30.0      
(STORED VAR TEM1 IS (:TE_R:+:DTE_R:*(YG-0.2)) with imat>100)
   inform13end
 
 INIADD  =    F
 ************************************************************
  Group 12. Convection and diffusion adjustments
   No PATCHes used for this Group
 ************************************************************
  Group 13. Boundary & Special Sources
 
 PATCH (INNER1  ,CELL,4,0,0,0,0,0,1,1)
 COVAL (INNER1  ,P1  , FIXVAL      , 1.000000E+07)
 COVAL (INNER1  ,W1  , FIXVAL      , 0.000000E+00)
 
 PATCH (INNER2  ,CELL  ,5,1,0,0,0,0,1,1)
 COVAL (INNER2  ,W1  , FIXVAL      , 0.000000E+00)
 
 PATCH (INNER3  ,CELL  ,6,1,0,0,0,0,1,1)
 COVAL (INNER3  ,V1  , FIXVAL      , 0.000000E+00)
 
 PATCH (INNER   ,CELL,7,0,0,0,0,0,1,1)
 COVAL (INNER   ,P1  , FIXVAL      , 1.000000E+07)
 
 PATCH (OUTER   ,CELL,8,0,0,0,0,0,1,1)
 COVAL (OUTER   ,P1  , FIXVAL      , 1.000000E+06)
 COVAL (OUTER   ,W1  , FIXVAL      , 0.000000E+00)
 
 PATCH (Z01     ,LOW   ,9,0,0,0,0,0,1,1)
 COVAL (Z01     ,W1  , FIXVAL      , 0.000000E+00)
 
 PATCH (Z02     ,HIGH  ,10,0,0,0,0,0,1,1)
 COVAL (Z02     ,W1  , FIXVAL      , 0.000000E+00)
 
 EGWF    =    T
 ************************************************************
  Group 14. Downstream Pressure For PARAB
 ************************************************************
  Group 15. Terminate Sweeps
 LSWEEP  =       400
 RESREF(V1  ) =  0.000000E+00 ;RESREF(W1  ) =  0.000000E+00
 RESFAC  = 1.000000E-07
 ************************************************************
  Group 16. Terminate Iterations
 LITER (V1  ) =      100 ;LITER (W1  ) =      100
 ENDIT (V1  ) =  0.000000E+00
 ************************************************************
  Group 17. Relaxation
 RELAX(P1  ,LINRLX, 1.000000E+00)
 ************************************************************
  Group 18. Limits
 VARMAX(V1  ) = 1.000000E+06 ;VARMIN(V1  ) =-1.000000E+06
 VARMAX(W1  ) = 1.000000E+06 ;VARMIN(W1  ) =-1.000000E+06
 ************************************************************
  Group 19. EARTH Calls To GROUND Station
 USEGRD  =    T  ;USEGRX =    T
 STRA    =    T
 ASAP    =    T
 PARSOL  =    F
 CONWIZ  =    T
 ISG21   =       400
 ISG52   =         1
 ISG59   =         1
 CSG10   ='q1'
 SPEDAT(SET,BOUNDARY,XCONST,R,1.00000E+20)
 ************************************************************
  Group 20. Preliminary Printout
 ECHO    =    T
 ************************************************************
  Group 21. Print-out of Variables
 OUTPUT(DRH1,N,N,N,N,N,N)
 OUTPUT(DVO1,N,N,N,N,N,N)
 OUTPUT(ENUL,N,N,N,N,N,N)
 OUTPUT(DEN1,N,N,N,N,N,N)
 ************************************************************
  Group 22. Monitor Print-Out
 IXMON   =         1 ;IYMON  =        21 ;IZMON  =        10
 NPRMON  =    100000
 NPRMNT  =         1
 TSTSWP  =        -1
 ************************************************************
  Group 23.Field Print-Out & Plot Control
 NPRINT  =    100000
 NYPRIN  =         1
 NZPRIN  =         1
 ISWPRF  =         1 ;ISWPRL =    100000
   No PATCHes used for this Group
 ************************************************************
  Group 24. Dumps For Restarts
 
 GVIEW(P,-9.118605E-01,2.358848E-01,3.359592E-01)
 GVIEW(UP,2.395464E-01,9.703851E-01,-3.115272E-02)
 
> DOM,    SIZE,        5.000000E-01, 8.080000E-01, 1.000000E+00
> DOM,    MONIT,       2.500000E-01, 7.850000E-01, 4.750000E-01
> 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,      1, 1.000000E+00
> GRID,   RSET_Y_2,     10, 1.000000E+00
> GRID,   RSET_Y_3,     10, 1.000000E+00
> GRID,   RSET_Y_4,      1, 1.000000E+00
> GRID,   RSET_Z_1,      4, 1.000000E+00
> GRID,   RSET_Z_2,      6, 1.000000E+00
> GRID,   RSET_Z_3,      9, 1.000000E+00
> GRID,   RSET_Z_4,      1, 1.000000E+00
 
> OBJ,    NAME,        POUT
> OBJ,    POSITION,    0.000000E+00, 8.000000E-01, 0.000000E+00
> OBJ,    SIZE,        5.000000E-01, 7.999957E-03, 1.000000E+00
> OBJ,    GEOMETRY,    polcubt2
> OBJ,    ROTATION24,        1
> OBJ,    TYPE,        BLOCKAGE
> OBJ,    MATERIAL,    0,  Air at 20 deg C, 1 atm, treated as i
> OBJ,    SCAL_FIXF,   Per-unit-volume
 
> OBJ,    NAME,        PIN
> OBJ,    POSITION,    0.000000E+00, 0.000000E+00, 0.000000E+00
> OBJ,    SIZE,        5.000000E-01, 2.000000E-01, 1.000000E+00
> OBJ,    GEOMETRY,    polcubt2
> OBJ,    ROTATION24,        1
> OBJ,    TYPE,        BLOCKAGE
> OBJ,    MATERIAL,    0,  Air at 20 deg C, 1 atm, treated as i
> OBJ,    SCAL_FIXF,   Per-unit-volume
 
> OBJ,    NAME,        PHOLE
> OBJ,    POSITION,    0.000000E+00, 2.000000E-01, 2.000000E-01
> OBJ,    SIZE,        5.000000E-01, 3.000000E-01, 3.000000E-01
> OBJ,    GEOMETRY,    polcubt2
> OBJ,    ROTATION24,        1
> OBJ,    TYPE,        BLOCKAGE
> OBJ,    MATERIAL,    0,  Air at 20 deg C, 1 atm, treated as i
> OBJ,    SCAL_FIXF,   Per-unit-volume
 
> OBJ,    NAME,        INNER1
> OBJ,    POSITION,    0.000000E+00, 0.000000E+00, 0.000000E+00
> OBJ,    SIZE,        5.000000E-01, 2.000000E-01, 1.000000E+00
> OBJ,    GEOMETRY,    poldef
> OBJ,    ROTATION24,        1
> OBJ,    TYPE,        USER_DEFINED
 
> OBJ,    NAME,        INNER2
> OBJ,    POSITION,    0.000000E+00, 2.010000E-01, 2.010000E-01
> OBJ,    SIZE,        5.000000E-01, 2.980000E-01, 2.480000E-01
> OBJ,    GEOMETRY,    poldef
> OBJ,    ROTATION24,        1
> OBJ,    TYPE,        CELLTYPE
 
> OBJ,    NAME,        INNER3
> OBJ,    POSITION,    0.000000E+00, 0.000000E+00, 2.010000E-01
> OBJ,    SIZE,        5.000000E-01, 4.700000E-01, 2.980000E-01
> OBJ,    GEOMETRY,    poldef
> OBJ,    ROTATION24,        1
> OBJ,    TYPE,        CELLTYPE
 
> OBJ,    NAME,        INNER
> OBJ,    POSITION,    0.000000E+00, 2.000000E-01, 2.000000E-01
> OBJ,    SIZE,        5.000000E-01, 3.000000E-01, 3.000000E-01
> OBJ,    GEOMETRY,    poldef
> OBJ,    ROTATION24,        1
> OBJ,    TYPE,        USER_DEFINED
 
> OBJ,    NAME,        OUTER
> OBJ,    POSITION,    0.000000E+00, 8.000000E-01, 0.000000E+00
> OBJ,    SIZE,        5.000000E-01, 7.999957E-03, 1.000000E+00
> OBJ,    GEOMETRY,    poldef
> OBJ,    ROTATION24,        1
> OBJ,    TYPE,        USER_DEFINED
 
> OBJ,    NAME,        Z01
> OBJ,    POSITION,    0.000000E+00, 2.000000E-01, 0.000000E+00
> OBJ,    SIZE,        5.000000E-01, 6.000000E-01, 0.000000E+00
> OBJ,    GEOMETRY,    poldef
> OBJ,    ROTATION24,        1
> OBJ,    TYPE,        USER_DEFINED
 
> OBJ,    NAME,        Z02
> OBJ,    POSITION,    0.000000E+00, 2.000000E-01, 9.500000E-01
> OBJ,    SIZE,        5.000000E-01, 6.000000E-01, 0.000000E+00
> OBJ,    GEOMETRY,    poldef
> OBJ,    ROTATION24,        1
> OBJ,    TYPE,        USER_DEFINED
 ************************************************************
STOP