****master
***geometry
 **point point0
  *position 0.000000000000000e+00 0.000000000000000e+00 0.000000000000000e+00
 **point point1
  *position 5.000000000000000e+01 0.000000000000000e+00 0.000000000000000e+00
 **point point2
  *position 7.000000000000000e+01 0.000000000000000e+00 0.000000000000000e+00
 **point point3
  *position 1.000000000000000e+02 0.000000000000000e+00 0.000000000000000e+00
 **point point4
  *position 0.000000000000000e+00 1.000000000000000e+02 0.000000000000000e+00
 **point point5
  *position 0.000000000000000e+00 7.000000000000000e+01 0.000000000000000e+00
 **point point6
  *position 0.000000000000000e+00 5.000000000000000e+01 0.000000000000000e+00
 **line line0
  *p1 point1
  *p2 point2
  *cuts 20
 **line line1
  *p1 point2
  *p2 point3
  *cuts 20
 **line line2
  *p1 point4
  *p2 point5
  *cuts 20
 **line line3
  *p1 point5
  *p2 point6
  *cuts 20
 **arc arc0
  *p1 point1
  *p2 point6
  *cuts 40
  *center point0
  *data   5.000000000000000e+01 0.000000000000000e+00 1.570796326794897e+00
 **arc arc1
  *p1 point2
  *p2 point5
  *cuts 40
  *center point0
  *data   7.000000000000000e+01 0.000000000000000e+00 1.570796326794897e+00
 **arc arc2
  *p1 point3
  *p2 point4
  *cuts 40
  *center point0
  *data   1.000000000000000e+02 0.000000000000000e+00 1.570796326794897e+00
 **domain domain1 Ruled5
  *name R0
  *element cax8
  *method 1
  *b1 arc0 
  *b2 line0 
  *b3 arc1 
  *b4 line3 
 **domain domain1 Ruled6
  *name R1
  *element cax8
  *method 1
  *b1 arc1 
  *b2 line1 
  *b3 arc2 
  *b4 line2 
***mesher
  **sweep 
    *elset ALL_ELEMENT 
    *new_elset ALL_ELEMENT 
    *prog 1.00000 
    *angle 90.0000 
    *num 10 
    *axis ( 0.00000 1.00000 0.00000 )
    *center ( 0.00000 0.00000 0.00000 )
    *fusion 0.00100000 
  **bset 
    *type cartesian
    *name right 
    *axes 0 1 2 
    *limit 0.00100000 
    *func x<0.001;
    *use_dimension 3 
    *req_number -1 
  **bset 
    *type cartesian
    *name inner_wall 
    *axes 0 1 2 
    *limit 0.00100000 
    *func x*x+y*y+z*z<50.*50.+0.001;
    *use_dimension 3 
    *req_number -1 
  **bset 
    *type cartesian
    *name outer_wall 
    *axes 0 1 2 
    *limit 0.00100000 
    *func x*x+y*y+z*z>100.*100.-0.001;
    *use_dimension 3 
    *req_number -1 
  **bset 
    *type cartesian
    *name left 
    *axes 0 1 2 
    *limit 0.00100000 
    *func z<0.001;
    *use_dimension 3 
    *req_number -1 
  **bset 
    *type cartesian
    *name bottom 
    *axes 0 1 2 
    *limit 0.00100000 
    *func y<0.001;
    *use_dimension 3 
    *req_number -1 
***plotting
  **default_options 1 0 1
****return