code: def area(poly):
      
      poly_xy = []
      num = len(poly)
      for i in range(num):
          poly[i] = poly[i][0:2] + (0,)
          poly_xy.append(poly[i])
      return surface.area(poly)
docs: Calculation of zone area