****mesher
     ***mesh sphere.geof
      **open sphere.geof

      **elset mat0         % elset for material 0
       *func x*x+y*y+z*z<70.*70.+0.001;

      **elset mat1         % elset for material 1 
       *func x*x+y*y+z*z>70.*70.-0.001;

    ****return