TALK=f;RUN(1,1)
  DISPLAY
  
  
  ENDDIS
  PHOTON USE
  p;;;;
  
  

  set prop off
  cl
  msg PRPS field
  gr ou z 1
  con prps z 1  fi;1.0
  
  pause
  cl
  msg x-displacement field
  gr ou z 1
  con U1 z 1  fi;.0001
  vec z 1  col 0
  
  pause
  cl
  msg y-displacement field
  gr ou z 1
  con V1 z 1 fi;.0001
  
  pause
  cl
  msg epsx field
  gr ou z 1
  con epsx z 1 fi;.0001
  
  pause
  cl
  msg epsy field
  gr ou z 1
  con epsy z 1 fi;.0001

  pause
  cl
  msg strx field
  gr ou z 1
  con strx z 1 fi;.0001
  
  pause
  cl
  msg stry field
  gr ou z 1
  con stry z 1 fi;.0001

  cl
  AUTOPLOT
  file                                                                            
  phi 5 

  data 1                                                                          
  V1 X 5 Y 1 29                                                                   
  col1 1          
  scale
  msg V1(Y)
  plot 1                                                                      
  pause 





  ENDUSE
  
 ************************************************************
  Group 1. Run Title and Number
 ************************************************************
READQ1=T
  begin_readq1 
 TEXT(Model of turbine rotor with big Blades
 
  Declarations and settings
  
REAL(RROT,RO,OMEGA,OMEGA2) 
RROT=0.2
RO  =0.3
OMEGA = 100*6.28
OMEGA2= OMEGA**2
ny =30 
nx =30
nz =1  
integer(IYROT,IYBLAD)
IYROT = (NY*2/3)
IYBLAD = (NY*16/30)

  **** !!! ***
ISG59 = 1
 
 ************************************************************
 
 ************************************************************
  Group 2. Time dependence
 STEADY  =    T
 ************************************************************
  Group 3. X-Direction Grid Spacing
 CARTES  =    F
 XULAST  = 60 /180 * 3.14 
 GRDPWR(X,NX,XULAST,1.0)
 ************************************************************
  Group 4. Y-Direction Grid Spacing
 GRDPWR(Y,NY,RO,1.0)

 ************************************************************
  Group 5. Z-Direction Grid Spacing
 ZWLAST  = 1.000000E+00
 GRDPWR(Z,NZ,1.0,1.0)
 
 ************************************************************
  Group 7. Variables: STOREd,SOLVEd,NAMEd
 ONEPHS  =    T
 SOLVE(P1,V1,U1)
 STORE(PRPS,DEN1,ENUL,DVO1,DRH1)
 STORE(STRX,EPSX,STRY,EPSY,EPSZ)

 ************************************************************
  GROUP 8. ITERATION NUMBERS ETC
 RESFAC=1.e-7 

 RESREF(V1)=0.0  ! to prevent premature exit
 LITER(V1) = 40 ! from solver
 RESREF(u1)=0.0; LITER(u1) = 40
 ************************************************************
  GROUP 9. PROPERTIES
  
 CSG10='Q1'                  ! materials with various POISSON ratios
  MATFLG=T;NMAT=2         
  110    7800.0    0.3       473.0   43.0      1.0e-5   0.5E-11     
  190    2700.0    0.35      896.0   204.0     2.35e-5  1.47e-11

  160 - STEEL
  161 - ALUMINIUM

 ************************************************************
  GROUP 11. INITIAL VALUES
FIINIT(PRPS)=0
PATCH(ROTOR,INIVAL,1,NX,1,IYBLAD-1,1,1,1,1)
COVAL(ROTOR,PRPS,FIXVAL,110)

PATCH(BLADES,INIVAL,1,NX,IYBLAD,NY,1,1,1,1)
COVAL(BLADES,PRPS,FIXVAL,190)

 ************************************************************
  GROUP 13. BOUNDARY & SPECIAL SOURCES
 
char(FV)
FV=:OMEGA2:*RV*Den1 
 
PATCH(BODYP,CELL,1,NX,1,NY,1,1,1,1)
COVAL(BODYP,P1,FIXVAL,0)

PATCH(ROTORV,VOLUME,1,NX,1,NY-1,1,1,1,1)
(SOURCE of V1 at ROTORV is COVAL(FIXFLU,:FV:))


patch(lowx,wwall,1,1,1,NY,1,1,1,1)
coval(lowx,u1,fIXval,0.0)
patch(highx,ewall,nx-1,nx-1,1,NY,1,1,1,1)
coval(highx,u1,fixval,0.0)

SPEDAT(BOUNDARY,ZCONST,R,1.E20)
 

 ************************************************************
  GROUP 15. TERMINATE SWEEPS
 LSWEEP  =  300
 ISG21   =  LSWEEP
 ENDIT(V1) =0.0
 RESREF(V1)=0.0                                   
 ENDIT(u1) =0.0
 RESREF(u1)=0.0                                   
 ************************************************************
  GROUP 17. RELAXATION
#CONPROM
 ************************************************************
  GROUP 19. DATA TRANSMITTED TO GROUND
 STRA    =    T
 ************************************************************
  GROUP 23.FIELD PRINT-OUT & PLOT CONTROL
TSTSWP = - 1   ! graphic-mode
     ISG52   =    3   !   probe & res
#MAXMIN
     
output(drh1,n,n,n,n,n,n)  
output(dvo1,n,n,n,n,n,n)  
output(enul,n,n,n,n,n,n)  
     NYPRIN=1
     NXPRIN=1
IXMON = NX/2
IYMON = NY-1
IZMON = 1
   spedat(rlxfac,rlxu1d,r,0.5)  
   spedat(rlxfac,rlxv1d,r,0.5)  
STOP