Changeset 2622 in lmdz_wrf for trunk

Timestamp:

Jun 22, 2019, 8:45:35 PM (5 years ago)

Author:

lfita

Message:

Adding:

• `mod_vec': Function to compute the module of a vector

File:

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

trunk/tools/geometry_tools.py

@@ -33 +33 @@
 # max_coords_poly: Function to provide the extremes of the coordinates of a polygon
 # mirror_polygon: Function to reflex a polygon for a given axis
+# mod_vec: Function to compute the module of a vector
 # position_sphere: Function to tranform fom a point in lon, lat deg coordinates to 
 #   cartesian coordinates over an sphere
@@ -267 +268 @@
 
     return dist
+
+def mod_vec(vec):
+    """ Function to compute the module of a vector
+      vec: vector [y, x]
+    >>> mod_vec([1., 1.])
+    1.41421356237
+    """
+    fname = 'mod_vec'
+
+    mod = np.sqrt(vec[0]*vec[0] + vec[1]*vec[1])
+
+    return mod
+
+print mod_vec([1., 1.])
+quit()
+
+def angle_vectors2D(veca, vecb):
+    """ Angle between two vectors
+      veca: angle A [ya, xa]
+      vecb: angle B [yb, xb]
+      NOTE: angle from A to B
+    >>> angle_vectors2D([1.,0.], [0.,1.])
+    
+    """
+    fname = 'angle_vectors2D'
+
+    moda = mod_vec(veca)
+    modb = mod_vec(vecb)
+
+    return alpha
 
 def max_coords_poly(polygon):
@@ -1523 +1554 @@
     >>> plgs = {'sqra': polya, 'sqrb': polyb}
     >>> pile_polygons(pns, plgs)
-
+    [[-0.75 -0.5 ]
+     [ 0.75 -0.5 ]
+     [ 0.75  0.5 ]
+     [-0.75  0.5 ]]
     """
     fname = 'pile_polygons'
@@ -1537 +1571 @@
         poly = polygons[polyn]
         Npts = poly.shape[0]
-        print ' LLuis shapes ipoly', ipoly.shape, 'poly', poly.shape
         Nint, inti, intp, pts = fsci.module_scientific.crossingpoints_polys(         \
           nvertexa=iNpts, nvertexb=Npts, nvertexab=iNpts*Npts, polya=ipoly,          \
           polyb=poly)
-        print Nint
-        for ip in range(Nint):
-            print ip, ':', inti[ip], intp[ip], '=', pts[ip]
+        iip=0
+        for ip in Nint:
+            newipoly = ipoly[iip:inti]
+        
        
     return pilepolygons
-
-#pns = ['sqra', 'sqrb']
-#polya = np.array([[-0.5, -0.75], [0.5, -0.75], [0.5, 0.75], [-0.5, 0.75]])
-#polyb = np.array([[-0.75, -0.5], [0.75, -0.5], [0.75, 0.5], [-0.75, 0.5]])
-#plgs = {'sqra': polya, 'sqrb': polyb}
-#pile_polygons(pns, plgs)
-
-
 
 ####### ###### ##### #### ### ## #