Changeset 2626 in lmdz_wrf

Timestamp:

Jun 23, 2019, 8:04:17 PM (5 years ago)

Author:

lfita

Message:

generalizing `dist_points'

File:

• trunk/tools/geometry_tools.py (modified) (2 diffs)

trunk/tools/geometry_tools.py

@@ -279 +279 @@
 
     v = np.array(vec, dtype=np.float)
-    mod = np.sqrt(v[0]*v[0] + v[1]*v[1])
+    vv = v*v
+    mod = np.sqrt(np.sum(vv[:]))
 
     return mod
@@ -1591 +1592 @@
         intp = intp-1
 
-        # There is a need to sort the cutting points from the perspective of the
-        #   below polygon
+        # Re-constructing below polygon looking for respective crossings
+        linti = list(inti)
+        for ip in range(iNpts):
+            Nc = linti.count(ip)
+            dists = np.zeros((Nc), dtype=np.float)
+            if Nc > 0:
+                iic = gen.multi_index_vec(inti,ip)
+                # Sortering from distance respect the vertex ip
+                for ic in range(Nc):
+                    dists[ic] = dist_points(ipoly[ip,:], pts[inti[ic],:])
+            else:
+                newipoly.append(ipoly[ip,:])
+
+
 
         Npi=0