------------------------------------
#eq     b  10
#eq    cos45  7.071067812E-01
#eq    mcos45 -7.071067812E-01
		
!       ci is ai x cos45
#eq     c1  35.35533906 
#eq     c2  3.464823228E+01
#eq     c3   3.111269837E+01
#eq     c4  3.040559159E+01 
#eq     c5  2.687005769E+01
#eq     c6  2.616295090E+01 
#eq     c7  2.262741700E+01
#eq     c8  2.192031022E+01
#eq     c9  1.838477631E+01 
#eq     c10 1.767766953E+01
#eq     c11 1.414213562E+01
#eq	c12  1.343502884E+01
#eq	c13 1.060660172E+01

#subd block
1 prism  Fe   0 0 0 / 0 0 0 50  b   c1 c1 /2
2 prism  SCIN 0 2 0 / 0 0 0 49  =   c2 c2 / 3
3 prism  Pb   0 0 0 / 0 0 0 44  =   c3 c3 / 4
4 prism  SCIN 0 2 0 / 0 0 0 43  =   c4 c4 /5
5 prism  Pb   0 0 0 / 0 0 0 38  =   c5 c5 / 6
6 prism  SCIN 0 2 0 / 0 0 0 37  =   c6 c6 /7
7 prism  Pb   0 0 0 / 0 0 0 32  =   c7 c7 / 8
8 prism  SCIN 0 2 0 / 0 0 0 31  =   c8 c8 / 9
9 prism  Pb   0 0 0 / 0 0 0 26  =   c9 c9 /10
10 prism  SCIN 0 2 0 / 0 0 0 25 =   c10 c10 /11
11 prism  Pb   0 0 0 / 0 0 0 20 =   c11 c11 /12
12 prism  SCIN 0 0 0 / 0 0 0 19 =   c12 c12 /13
13 prism  sp   0 0 0 / 0 0 0 15 =   c13 c13
#end block

1 block sp  0 0 0 / 0 0 0
#      rotate it; since "b" is not inclined, b= '0 1 0'  
#      can be omitted  in the last part (cos45 0 cos45 0 1 0 is full)
#      
2 block sp  0 0 0 / 0 0 0   cos45 0  cos45 
3 block sp  0 0 0 / 0 0 0   0 0 1     
4 block sp  0 0 0 / 0 0 0   mcos45 0  cos45 
5 block sp  0 0 0 / 0 0 0   -1  0 0  
6 block sp  0 0 0 / 0 0 0   mcos45 0  mcos45
7 block sp  0 0 0 / 0 0 0   0 0  -1  
8 block sp  0 0 0 / 0 0 0   cos45 0  mcos45 

/*
#news   new-2  cone
#news   new-3  cap
1 ecyl Pb 0 0 0 / 0 0 0    4 1 3
2 cone Pb 0 0 0 / 3 3 3    2 1 3 3
3 cap  Pb 0 0 0 / -3 -3 -3 1 1.5 -0.5
4 pipe	Pb 0 0 0 / 1 2 3  1 2 3
5 pipe_y	Pb 0 0 0 / 5 0 0  1 2 3
*/
-------------------------
