Changeset 2659 in lmdz_wrf for trunk

Timestamp:

Jul 3, 2019, 5:01:12 AM (5 years ago)

Author:

lfita

Message:

Adding: arc' and side' into 'join_circ_sec'

File:

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

trunk/tools/geometry_tools.py

@@ -403 +403 @@
     return reflex
 
-def join_circ_sec(points, radfrac=3., N=200):
+def join_circ_sec(points, radfrac=3., arc='short', side='left', 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 
@@ -410 +410 @@
       radfrac: multiplicative factor of the distance between consecutive points to 
         draw the circular segment (3., default)
+      arc: type of arc ('short', default)
+      pos: position of arc ('left', default)
       N: number of points (200, default)
     """
@@ -424 +426 @@
         p2 = lpoints[ip+1]
         dps = dist_points(p1, p2)
-        jcirc_sec[Np*ip:Np*(ip+1),:] = circ_sec(p1, p2, dps*radfrac, 'short', Np)
+        jcirc_sec[Np*ip:Np*(ip+1),:] = circ_sec(p1,p2,dps*radfrac, arc, side, Np)
 
     Np2 = N - (Npts-1)*Np
@@ -430 +432 @@
     p2 = lpoints[0]
     dps = dist_points(p1, p2)
-    jcirc_sec[(Npts-1)*Np:N,:] = circ_sec(p1, p2, dps*3., 'short', Np2)
+    jcirc_sec[(Npts-1)*Np:N,:] = circ_sec(p1, p2, dps*3., arc, side, Np2)
 
     return jcirc_sec