Index: /issm/trunk-jpl/src/m/classes/mask.py
===================================================================
--- /issm/trunk-jpl/src/m/classes/mask.py	(revision 24901)
+++ /issm/trunk-jpl/src/m/classes/mask.py	(revision 24902)
@@ -65,5 +65,5 @@
             self.ice_levelset[index] = -ice
         elif nargin == 1:
-            ice = *args[0]
+            ice = args[0]
             self.ocean_levelset = -ocean
         else:
Index: /issm/trunk-jpl/src/m/coordsystems/flaglatlongradius.py
===================================================================
--- /issm/trunk-jpl/src/m/coordsystems/flaglatlongradius.py	(revision 24901)
+++ /issm/trunk-jpl/src/m/coordsystems/flaglatlongradius.py	(revision 24902)
@@ -25,8 +25,6 @@
         pass
 
-    distance = []
-    for i in range(len(long)):
-        distance[i] = math.sqrt(pow((lat[i] - lat0[i]), 2) + pow((long[i] - long0[i]), 2))
-    indices = (distance <= radius).np.nonzero() 
+    distance    = np.array(math.sqrt(pow((lat[i] - lat0[i]), 2) + pow((long[i] - long0[i]), 2)))
+    indices     = (distance <= radius).nonzero() 
 
     return indices
Index: /issm/trunk-jpl/src/m/coordsystems/flagradiuselements.m
===================================================================
--- /issm/trunk-jpl/src/m/coordsystems/flagradiuselements.m	(revision 24901)
+++ /issm/trunk-jpl/src/m/coordsystems/flagradiuselements.m	(revision 24902)
@@ -6,5 +6,4 @@
 	y0 = R .* cosd(lat0) .* sind(long0);
 	z0 = R .* sind(lat0);
-
 
 	distance=sqrt( (x-x0).^2+ (y-y0).^2 + (z-z0).^2);
Index: /issm/trunk-jpl/src/m/coordsystems/flagradiuselements.py
===================================================================
--- /issm/trunk-jpl/src/m/coordsystems/flagradiuselements.py	(revision 24902)
+++ /issm/trunk-jpl/src/m/coordsystems/flagradiuselements.py	(revision 24902)
@@ -0,0 +1,28 @@
+import math
+
+import numpy as np
+
+from planetradius import *
+
+
+def flagradiuselements(elements, x, y, z, lat0, long0, radius): # {{{
+    # get x0,y0,z0:
+    R   = planetradius('earth')
+    x0  = R * np.cos(np.deg2rad(lat0)) * np.cos(np.deg2rad(long0))
+    y0  = R * np.cos(np.deg2rad(lat0)) * np.sin(np.deg2rad(long0))
+    z0  = R * np.sin(np.deg2rad(lat0))
+
+    distance    = np.array(math.sqrt(pow((x - x0), 2) + pow((y - y0), 2) + pow((z - z0), 2)))
+    indices     = (distance <= radius * 1000).nonzero()
+
+    #now that we know the indices, determine elements with own these indices:
+    outelements = np.zeros(len(elements))
+    for i in range(len(indices)):
+        pos = np.array()
+        pos = (elements == indices[i]).nonzero()
+        for j in range(len(pos)):
+            outelements[pos[j]] = 1
+    outelements = outelements.nonzero()
+
+    return outelements
+# }}}
Index: /issm/trunk-jpl/src/m/geometry/planetradius.m
===================================================================
--- /issm/trunk-jpl/src/m/geometry/planetradius.m	(revision 24901)
+++ /issm/trunk-jpl/src/m/geometry/planetradius.m	(revision 24902)
@@ -6,5 +6,5 @@
 %
 %   Examples:
-%      earthradius  =planetradius('earth');
+%      earthradius=planetradius('earth');
 
 if strcmpi(planet,'earth'),
Index: /issm/trunk-jpl/src/m/geometry/planetradius.py
===================================================================
--- /issm/trunk-jpl/src/m/geometry/planetradius.py	(revision 24902)
+++ /issm/trunk-jpl/src/m/geometry/planetradius.py	(revision 24902)
@@ -0,0 +1,18 @@
+def planetradius(planet): # {{{
+    '''
+    PLANETRADIUS - return planet radius according to planetary body name
+
+        Usage:
+            radius = planetradius(planet)
+
+        Examples:
+            earthradius = planetradius('earth')
+    '''
+
+    if planet == 'earth':
+        radius = 6.371012 * pow(10, 6)
+    else:
+        raise TypeError("planet type %s not supported yet!" % planet)
+
+    return radius
+# }}}