FINISH
/CLEAR

Ri=0.15
R0=0.23
Hacc=0.06
Pi=2e6
Pe=6e6
omega=20
ndiv=5

/PREP7
  
ET,1,PLANE182   
!*  
KEYOPT,1,1,0
KEYOPT,1,3,1
KEYOPT,1,6,0
!*  

MP,EX,1,200e9  
MP,PRXY,1,0.3  
MP,DENS,1,7800 

K,1,Ri,0,0
K,2,R0,0,0
K,3,R0,Hacc/2,0
K,4,Ri,Hacc/2,0

LSTR,       1,       2  
LSTR,       2,       3  
LSTR,       3,       4  
LSTR,       4,       1  

AL,1,2,3,4 


LESIZE,2, , ,ndiv,1
LESIZE,4, , ,ndiv,1
LESIZE,1, , ,50,-5
LESIZE,3, , ,50,-5

MSHKEY,1
AMESH,1   



SFL,4,PRES,Pi,  
SFL,2,PRES,Pe,
OMEGA,0,omega,0,   
 
DL,1, ,UY,0
  
FINISH  
/SOL

SOLVE   

FINISH  
/POST1  

PLNSOL, S,EQV, 0,1.0

/WAIT,5

PATH,Pathrad,2,30,20,   
PPATH,1,,Ri,0,0
PPATH,2,,R0,0,0

PDEF,SigmaR,S,X,AVG 
PDEF,SigmaC,S,Z,AVG 
PDEF,SigmaVM,S,EQV,AVG  

PLPATH,SIGMAR,SIGMAC,SIGMAVM