Changeset 2576 in lmdz_wrf

Timestamp:

Jun 1, 2019, 8:12:33 PM (5 years ago)

Author:

lfita

Message:

Adding:

• `special_buoy1': Function to draw an special mark buoy using buoy1

File:

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

trunk/tools/geometry_tools.py

@@ -71 +71 @@
 # red_buoy1: Function to draw a red mark buoy using buoy1
 # safewater_buoy1: Function to draw a safe water mark buoy using buoy1
+# special_buoy1: Function to draw an special mark buoy using buoy1
 # surface_sphere: Function to provide an sphere as matrix of x,y,z coordinates
 # z_boat: Function to define an schematic boat from the z-plane
@@ -2633 +2634 @@
     return buoy, buoysecs, buoydic
 
+def special_buoy1(height=5., width=10., bradii=1.75, bfrac=0.8, hsigns=0.5, N=300):
+    """ Function to draw an special mark buoy using buoy1
+      height: height of the prism (5., default) 
+      width: width of the prism (10., default)
+      bradii: radii of the ball (1.75, default)
+      bfrac: fraction of the ball above the prism (0.8, default)
+      hisgns: height of the signs [as reg. triangle] as percentage of the height 
+        (0.5, default)
+      N: total number of points of the buoy (300, default)
+    """
+    fname = 'special_buoy1'
+
+    buoy = np.ones((N,2), dtype=np.float)*gen.fillValueF
+
+    # buoy
+    N2 = int(N/2)
+    buoy1v, buoy1vsecs, buoy1vdic = buoy1(height=5., width=10., bradii=1.75,         \
+      bfrac=0.8, N=N2)
+    buoy[0:N2,:] = buoy1v
+
+    Height = np.max(buoy1v[:,0])
+
+    # sign
+    N3 = N - N2 - 1
+    
+    bottsigns = 2.*bradii+height
+    lsign = height*hsigns
+    # up
+    cross, crosssecs, crossdic = p_cross_width(lsign, width=0.2*lsign, Narms=2, N=N3)
+    cross = rotate_polygon_2D(cross, 40.05)
+    buoy[N2+1:N,:] = cross + [Height+1.1*lsign,0.]
+
+    # painting it
+    buoy = ma.masked_equal(buoy, gen.fillValueF)
+
+    buoysecs = ['buoy', 'sign']
+    buoydic = {'buoy': [buoy[0:N2,:],'-','#FFFF00',1.5],                             \
+      'sign': [buoy[N2+1:N,:],'-','#FFFF00',1.5]}
+
+    return buoy, buoysecs, buoydic
+
 ####### ####### ##### #### ### ## #
 # Plotting