Changeset 2496 in lmdz_wrf for trunk/tools

Timestamp:

May 2, 2019, 1:33:01 AM (6 years ago)

Author:

lfita

Message:

Adding:

• `zisland1': Function to draw an island from z-axis as the union of a series of points by circular segments

File:

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

trunk/tools/geometry_tools.py

@@ -52 +52 @@
 # z_boat: Function to define an schematic boat from the z-plane
 # zsailing_boat: Function to define an schematic sailing boat from the z-plane with sails
-
+# zisland1: Function to draw an island from z-axis as the union of a series of points by
+#   circular segments
 ## Plotting
 # plot_sphere: Function to plot an sphere and determine which standard lines will be 
@@ -889 +890 @@
     return sailingboat, sailingboatsecs, dicsailingboat
 
+def zisland1(mainpts= np.array([[-0.1,0.], [-1.,1.], [-0.8,1.2], [0.1,0.6], [1., 0.9],\
+  [2.8, -0.1], [0.1,-0.6]], dtype=np.float), radfrac=3., N=200):
+    """ Function to draw an island from z-axis as the union of a series of points by
+        circular segments
+      mainpts: 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)
+          * default= np.array([[-0.1,0.], [-1.,1.], [-0.8,1.2], [0.1,0.6], [1., 0.9],
+            [2.8, -0.1], [0.1,-0.6]], dtype=np.float)
+      radfrac: multiplicative factor of the distance between consecutive points to 
+        draw the circular segment (3., default)
+      N: number of points (200, default)
+    """
+    import numpy.ma as ma
+    fname = 'zisland1'
+
+    island1 = np.ones((N,2), dtype=np.float)*gen.fillValueF
+
+    # 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)
+
+    islandsecs = ['coastline']
+    islanddic = {'coastline': [island1, '-', '#161616', 2.]}
+
+    island1 = ma.masked_equal(island1, gen.fillValueF)
+
+    return island1, islandsecs, islanddic
+
 def write_join_poly(polys, flname='join_polygons.dat'):
     """ Function to write an ASCII file with the combination of polygons