TALK=f;RUN(1,1)
  DISPLAY
    Stresses are computed for a rotating disk with zero
    axial-direction (z) stress
     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
 ************************************************************
  Group 1. Run Title and Number
 ************************************************************
readq1=t
  READQ1_BEGIN 
 TEXT(Stress in a rotating disk; 1Dy
  Declarations and settings
REAL(DenBody,RO,OMEGA,POISSON,YOUNG) 
DenBody = 7000 
OMEGA =100 * 6.28
RO = 0.1 
YOUNG   = 1/(0.5E-11) ! Young's modulus
POISSON=0.3           ! Poisson's ratio
INTEGER(CASENO,NYCYL) ! nycyl is y-direction cell number in solid
NYCYL=31

real(FVY)
FVY = DenBody*Omega**2


  READQ1_END 
 ************************************************************
  Group 2. Time dependence
 STEADY  =    T
 ************************************************************
  Group 3. X-Direction Grid Spacing
 CARTES  =    F
 NX=1
 XULAST  = 1.000000E-03
 XFRAC (   1) =  1.000000E+00
 ************************************************************
  Group 4. Y-Direction Grid Spacing
 GRDPWR(Y,NYCYL,RO,1)
 
 ************************************************************
  Group 5. Z-Direction Grid Spacing
 NZ=1
 ZWLAST  = 1.000000E+00
 
   **** !!! ***
ISG59 = 1

 
 ************************************************************
  Group 7. Variables: STOREd,SOLVEd,NAMEd
 ONEPHS  =    T
 SOLVE(P1,V1)
 STORE(PRPS,DEN1,ENUL,DVO1,DRH1)
 STORE(STRX,EPSX,STRY,EPSY,EPSZ)
 STORE(V1TH,V1/T)
 ************************************************************
  GROUP 8. ITERATION NUMBERS ETC
 RESREF(V1)=0.0   ! to prevent premature exit
 LITER(V1) = 20   ! from solver
 ************************************************************
  GROUP 9. PROPERTIES
  
 CSG10='Q1'                  ! materials with various POISSON ratios
  MATFLG=T;NMAT=1         
  160    7800.0    0.3       473.0   43.0      1.0e-5   0.5E-11 
  161    7800.0    0.2       473.0   43.0      1.0e-5   0.5E-11 
  162    7800.0    0.1       473.0   43.0      1.0e-5   0.5E-11 
  163    7800.0    0.0       473.0   43.0      1.0e-5   0.5E-11 
  164    7800.0    0.4       473.0   43.0      1.0e-5   0.5E-11 
  165    7800.0    0.05      473.0   43.0      1.0e-5   0.5E-11 
  166    7800.0    0.01      473.0   43.0      1.0e-5   0.5E-11 
 
 ************************************************************
  GROUP 11. INITIAL VALUES
FIINIT(PRPS)=160

 ************************************************************
  GROUP 13. BOUNDARY & SPECIAL SOURCES
 
PATCH(BODY,VOLUME,1,1,1,NY-1,1,1,1,1)
(SOURCE of V1 at BODY is COVAL(FIXFLU,:FVY:*RV))

SPEDAT(BOUNDARY,XCONST,R,1.E20)
spedat(rlxfac,rlxV1d,r,0.5) 
 
 ************************************************************
  GROUP 15. TERMINATE SWEEPS
 LSWEEP  =  200
 ISG21   =   LSWEEP
 ENDIT(V1) =0.0
 RESFAC=1.e-7 
 RESREF(V1)=0.0                                   
 ************************************************************
  GROUP 17. RELAXATION
#CONPROM
 ************************************************************
  GROUP 19. DATA TRANSMITTED TO GROUND
 STRA    =    T
 ************************************************************
  GROUP 23.FIELD PRINT-OUT & PLOT CONTROL
output(den1,n,n,n,n,n,n)  
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(prps,n,n,n,n,n,n)  
NYPRIN=1
TSTSWP=-1  
IYPRF=1
IYPRL=NY

ISG52   =    3   !   probe & res

 
  inform7begin       
  ! Declarations first:
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:*(RV*:CC1:-:CC2:*RV^3)
(STORED VAR V1TH IS :formV:)
(STORED VAR V1/T IS V1/V1TH)
  inform7end
STOP