Changeset 2506 in lmdz_wrf for trunk/tools

Timestamp:

May 4, 2019, 4:35:35 PM (6 years ago)

Author:

lfita

Message:

Adding:

• `join_circ_sec': Function to join aa series of points by circular segments

File:

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

trunk/tools/geometry_tools.py

@@ -23 +23 @@
 # deg_deci: Function to pass from degrees [deg, minute, sec] to decimal angles [rad]
 # dist_points: Function to provide the distance between two points
+# join_circ_sec: Function to join aa series of points by circular segments
 # 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
@@ -285 +286 @@
 
     return reflex
+
+def join_circ_sec(points, radfrac=3., N=200):
+    """ Function to join aa series of points by circular segments
+      points: main points of the island (clockwise ordered, to be joined by circular 
+        segments of radii as the radfrac factor of the distance between 
+        consecutive points)
+      radfrac: multiplicative factor of the distance between consecutive points to 
+        draw the circular segment (3., default)
+      N: number of points (200, default)
+    """
+    fname = 'join_circ_sec'
+
+    jcirc_sec = np.ones((N,2), dtype=np.float)
+
+    # main points
+    lpoints = list(points)
+    Npts = len(lpoints)
+    Np = int(N/(Npts+1))
+    for ip in range(Npts-1):
+        dps = dist_points(lpoints[ip], lpoints[ip+1])
+        jcirc_sec[Np*ip:Np*(ip+1),:] = circ_sec(lpoints[ip], lpoints[ip+1],          \
+          dps*radfrac, Np)
+
+    Np2 = N - (Npts-1)*Np
+    dps = dist_points(lpoints[Npts-1], lpoints[0])
+    jcirc_sec[(Npts-1)*Np:N,:] = circ_sec(lpoints[Npts-1], lpoints[0], dps*3., Np2)
+
+    return jcirc_sec
 
 ####### ###### ##### #### ### ## #
@@ -910 +939 @@
     # Coastline
    
-    # main points
-    lmainpts = list(mainpts)
-    Npts = len(lmainpts)
-    Np = int(N/(Npts+1))
-    for ip in range(Npts-1):
-        dps = dist_points(lmainpts[ip], lmainpts[ip+1])
-        island1[Np*ip:Np*(ip+1),:] = circ_sec(lmainpts[ip], mainpts[ip+1],          \
-          dps*radfrac, Np)
-
-    Np2 = N - (Npts-1)*Np
-    dps = dist_points(lmainpts[Npts-1], lmainpts[0])
-    island1[(Npts-1)*Np:N,:] = circ_sec(lmainpts[Npts-1], lmainpts[0], dps*3., Np2)
+    island1 = join_circ_sec(mainpts, radfrac=3., N=200)
 
     islandsecs = ['coastline']