TALK=F;RUN( 1, 1)
 
 ************************************************************
   Q1 created by VDI menu, Version 2009, Date 17/02/09
 CPVNAM=VDI;SPPNAM=Core
 ************************************************************
  Echo DISPLAY / USE settings
  DISPLAY
    Stresses are considered in rotating disk with StrZ = 0:
     Omega   - Angular velocity of the rotation,
     DenBody - density,
     RO      - radius of disk.
    Exact solution (Timoshenko, paragraph 32):
 
     V = DenBody*Omega^2*RO^3 * (1-P)/8/E * { R1*(3+P) -
         R1^3*(1+P) },
     R1 = R/RO
 
  ENDDIS
 ************************************************************
 IRUNN = 1 ;LIBREF = 508
 ************************************************************
  Group 1. Run Title
 TEXT(USP: Stress in a 2D rotating disk; S508    )
 ************************************************************
  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,1,3.22581E-05)
 ************************************************************
  Group 6. Body-Fitted coordinates
 ************************************************************
  Group 7. Variables: STOREd,SOLVEd,NAMEd
    * Non-default variable names
 NAME(134)=TET7 ;NAME(135)=DISR
 NAME(136)=R7 ;NAME(137)=V1/T
 NAME(138)=V1TH ;NAME(139)=EPSZ
 NAME(140)=EPSY ;NAME(141)=EPSX
 NAME(142)=STRY ;NAME(143)=STRX
 NAME(144)=DRH1 ;NAME(145)=DVO1
 NAME(146)=ENUL ;NAME(147)=DEN1
 NAME(148)=PRPS ;NAME(149)=DISY
 NAME(150) =DISX
    * Solved variables list
 SOLVE(DISY,DISX)
    * Stored variables list
 STORE(PRPS,DEN1,ENUL,DVO1,DRH1,STRX,STRY,EPSX)
 STORE(EPSY,EPSZ,V1TH,V1/T,R7,DISR,TET7)
 
 ************************************************************
  Echo save-block settings for Group  7
  SAVE7BEGIN
 
REAL(DenBody,RO,OMEGA,POISSON,YOUNG)
DenBody = 7800
OMEGA =100 * 6.28
RO = 0.1
YOUNG   = 1/(0.5E-11) ! Young's modulus
POISSON=0.3           ! Poisson's ratio
 
real(FVY)
FVY = DenBody*Omega**2
 
REAL(CC0, CC1, CC2)
CC0 = FVY*RO**3*(1-POISSON)/YOUNG/8
CC1 = (3+POISSON)/RO
CC2 = (1+POISSON)/RO**3
char(formV)
formV=:CC0:*(R7*:CC1:-:CC2:*R7^3)
(STORED VAR V1TH IS :formV:)
(STORED VAR V1/T IS DISR/V1TH)
(STORED VAR DISR IS SQRT(DISX^2+DISY^2))
(STORED VAR R7 IS SQRT(XG^2+YG^2))
(STORED VAR TET7 IS ATAN(YG/(XG+1.e-10)))
 
char(frmSX,frmSY)
frmSX=:FVY:*R7*COS(TET7)*VOL
frmSY=:FVY:*R7*SIN(TET7)*VOL
 
  SAVE7END
 ************************************************************
  Group 8. Terms & Devices
 USP = T ; USPDBG= F
 UTCPLT = T ; USPIMB= F
 UAUTO = T ; UGTGRD= F
 USPVTK = T ; USPTER= F
 USCANX = T ; USCANY= T
 USCANZ = T
 MXLEV = 0 ;MYLEV = 0 ;MZLEV = 0
 DOMAT = -1
 CELLST = 10 ; FACEST= 10
 MINPRP = -1 ; MAXPRP= 250
 UMXLEV = 4 ; UMXCLS= 10000
 UNOLAY = 2 ; UIOLEV= 1
 UWALEV = 1 ; UGRITR= 2
 UMXCRS = 16
 USPREL =0.7 ; UMNCEL=1.
 UDVREL =1.
 ************************************************************
  Group 9. Properties
    * Domain material index is 198 signifying:
    * Solid with smooth-wall friction
 SETPRPS(1,198,Q1)
 ENUT =0.
   *  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(V1/T)=1.001E-10 ;FIINIT(V1TH)=1.001E-10
 FIINIT(EPSZ)=1.001E-10 ;FIINIT(EPSY)=1.001E-10
 FIINIT(EPSX)=1.001E-10 ;FIINIT(STRY)=1.001E-10
 FIINIT(STRX)=1.001E-10 ;FIINIT(DRH1)=1.001E-10
 FIINIT(DVO1)=1.001E-10 ;FIINIT(ENUL)=1.001E-10
 FIINIT(DEN1)=1.001E-10 ;FIINIT(PRPS)=198.
   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
   No PATCHes used for this Group
 
 EGWF = T
 ************************************************************
  Echo save-block settings for Group 13
  SAVE13BEGIN
(SOURCE of DISX at LEFT is COVAL(1,0)!VOLU)
(SOURCE of DISY at BOT is COVAL(1,0)!VOLU)
(SOURCE of DISX at BODY is COVAL(FIXFLU,:frmSX:)!VOLU)
(SOURCE of DISY at BODY is COVAL(FIXFLU,:frmSY:)!VOLU)
 
  SAVE13END
 ************************************************************
  Group 14. Downstream Pressure For PARAB
 ************************************************************
  Group 15. Terminate Sweeps
 LSWEEP = 400
 RESFAC =1.0E-07
 ************************************************************
  Group 16. Terminate Iterations
 ************************************************************
  Group 17. Relaxation
 ************************************************************
  Group 18. Limits
 ************************************************************
  Group 19. EARTH Calls To GROUND Station
 STRA = T
 CONWIZ = T
 ISG21 = 100
 ISG52 = 3
 ISG59 = 1
 CSG10   ='Q1'
 ************************************************************
  Group 20. Preliminary Printout
 ************************************************************
  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)
 OUTPUT(PRPS,N,N,N,N,N,N)
 ************************************************************
  Group 22. Monitor Print-Out
 IXMON = 5 ;IYMON = 1 ;IZMON = 1
 NPRMON = 100000
 NPRMNT = 1
 TSTSWP = -1
 ************************************************************
  Group 23.Field Print-Out & Plot Control
 NPRINT = 100000
 IYPRF = 1 ;IYPRL = 31
 ISWPRF = 1 ;ISWPRL = 100000
   No PATCHes used for this Group
 ************************************************************
  Group 24. Dumps For Restarts
 
 GVIEW(P,-0.111061,-0.02338,0.993539)
 GVIEW(UP,-0.093718,0.995515,0.01295)
 
> DOM,    SIZE,        1.200000E-01, 1.200000E-01, 2.000000E-02
> DOM,    MONIT,       5.400000E-02, 1.935484E-03, 1.000000E-02
> DOM,    SCALE,       1.000000E+00, 1.000000E+00, 1.000000E+00
> DOM,    INCREMENT,   1.000000E-02, 1.000000E-04, 1.000000E-04
  > GRID,   RSET_X_1,     10, 1.000000E+00
  > GRID,   RSET_Y_1,     10, 1.000000E+00
  > GRID,   RSET_Z_1,      1, 1.000000E+00
 
> OBJ,    NAME,        LEFT
> OBJ,    POSITION,    0.000000E+00, 0.000000E+00, 0.000000E+00
> OBJ,    SIZE,        0.000000E+00, 1.200000E-01, 2.000000E-02
> OBJ,    GEOMETRY,    cube11
> OBJ,    TYPE,        PLATE
 
> OBJ,    NAME,        BOT
> OBJ,    POSITION,    0.000000E+00, 0.000000E+00, 0.000000E+00
> OBJ,    SIZE,        1.200000E-01, 0.000000E+00, 2.000000E-02
> OBJ,    GEOMETRY,    cube11
> OBJ,    TYPE,        PLATE
 
> OBJ,    NAME,        BODY
> OBJ,    POSITION,    0.000000E+00, 0.000000E+00, 0.000000E+00
> OBJ,    SIZE,        1.000000E-01, 1.000000E-01, 2.000000E-02
> OBJ,    GEOMETRY,    quarter-cylinder
> OBJ,    GRID,        NO
> OBJ,    TYPE,        BLOCKAGE
> OBJ,    MATERIAL,    160,NOTSET
 ************************************************************
STOP