Changeset 2435 in lmdz_wrf for trunk

Timestamp:

Apr 14, 2019, 10:35:05 PM (6 years ago)

Author:

lfita

Message:

Adding:

• `max_coords_poly': Function to provide the extremes of the coordinates of a polygon

File:

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

trunk/tools/geometry_tools.py

@@ -20 +20 @@
 # deg_deci: Function to pass from degrees [deg, minute, sec] to decimal angles [rad]
 # dist_points: Function to provide the distance between two points
+# max_coords_poly: Function to provide the extremes of the coordinates of a polygon
 # rotate_polygon_2D: Function to rotate 2D plain the vertices of a polygon
 # position_sphere: Function to tranform fom a point in lon, lat deg coordinates to 
@@ -218 +219 @@
     return dist
 
+def max_coords_poly(polygon):
+    """ Function to provide the extremes of the coordinates of a polygon
+      polygon: coordinates [Nvertexs, 2] of a polygon
+    >>> square = np.zeros((4,2), dtype=np.float)
+    >>> square[0,:] = [-0.5,-0.5]
+    >>> square[1,:] = [0.5,-0.5]
+    >>> square[2,:] = [0.5,0.5]
+    >>> square[3,:] = [-0.5,0.5]
+    >>> max_coords_poly(square)
+    [-0.5, 0.5], [-0.5, 0.5], [0.5, 0.5]
+    """
+    fname = 'max_coords_poly'
+
+    nx = np.min(polygon[:,1])
+    xx = np.max(polygon[:,1])
+    ny = np.min(polygon[:,0])
+    xy = np.max(polygon[:,0])
+
+    axx = np.max(np.abs(polygon[:,1]))
+    ayx = np.max(np.abs(polygon[:,0]))
+
+    return [nx, xx], [ny, xy], [axx, ayx]
+
 ####### ###### ##### #### ### ## #
 # Shapes/objects
@@ -352 +376 @@
 
 # FROM: http://www.photographers1.com/Sailing/NauticalTerms&Nomenclature.html
-def zsailing_boat(length=10., beam=3., sternbp=0.5):
+def zsailing_boat(length=10., beam=1., lbeam=0.5, sternbp=0.5):
     """ Function to define an schematic boat from the z-plane
-      length: length of the boat
-      beam: beam of the boat
-      sternbp: beam at stern as percentage of beam
+      length: length of the boat (without stern, default 10)
+      beam: beam of the boat (default 1)
+      lbeam: length at beam (as percentage of length, default 0.5)
+      sternbp: beam at stern (as percentage of beam, default 0.5)
     """
     fname = 'zsailing_boat'
 
-    
+    bow = np.array([length, 0.])
+    maxportside = np.array([length*lbeam, -beam])
+    maxstarboardside = np.array([length*lbeam, beam])
+    portside = np.array([0., -beam*sternbp])
+    starboardside = np.array([0., beam*sternbp])
+
+    # forward section
+    fportsaid = circ_sec(bow,maxportside, length*2)
+    fstarboardsaid = circ_sec(maxstarboardside, bow, length*2)
+    # aft section
+    aportsaid = circ_sec(maxportside, portside, length*2)
+    astarboardsaid = circ_sec(starboardside, maxstarboardside, length*2)
+    # stern
+    stern = circ_sec(portside, starboardside, length*2)
+
+    dpts = stern.shape[0]
+    boat = np.zeros((dpts*5,2), dtype=np.float)
+
+    boat[0:dpts,:] = fportsaid
+    boat[dpts:2*dpts,:] = aportsaid
+    boat[2*dpts:3*dpts,:] = stern
+    boat[3*dpts:4*dpts,:] = astarboardsaid
+    boat[4*dpts:5*dpts,:] = fstarboardsaid
 
     return boat