code: def run(xmin, ymin, xmax, ymax, step, range_, range_x, range_y, t):
      pt = [[0] * range_y for _ in range(range_x)]
      "omp parallel for"
      for i in range(range_x):
          for j in range(range_y):
              s = 0
              for k in t:
                  tmp = 6368.0 * math.acos(
                      math.cos(xmin + step * i)
                      * math.cos(k[0])
                      * math.cos((ymin + step * j) - k[1])
                      + math.sin(xmin + step * i) * math.sin(k[0])
                  )
                  if tmp < range_:
                      s += k[2] / (1 + tmp)
              pt[i][j] = s
      return pt
  
label_desc: Use image node_id 1 for calling a global function with example usage: run(xmin, ymin, xmax, ymax, step, range_, range_x, range_y, t) and returns: pt