Changeset 2570 in lmdz_wrf for trunk

Timestamp:

May 29, 2019, 11:06:08 PM (6 years ago)

Author:

lfita

Message:

Adding:

• `isolateddanger_buoy1': Function to draw an isolated danger buoy using buoy1

File:

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

trunk/tools/geometry_tools.py

@@ -51 +51 @@
 # ellipse_polar: Function to determine an ellipse from its center and polar coordinates
 # green_buoy1: Function to draw a green mark buoy using buoy1
+# isolateddanger_buoy1: Function to draw an isolated danger buoy using buoy1
 # p_angle_triangle: Function to draw a triangle by an initial point and two 
 #   consecutive angles and the first length of face. The third angle and 2 and 3rd 
@@ -2409 +2410 @@
     return buoy, buoysecs, buoydic
 
+def isolateddanger_buoy1(height=5., width=10., bradii=1.75, bfrac=0.8, hsigns=0.5,   \
+  N=300):
+    """ Function to draw an isolated danger 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 = 'isolateddanger_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
+
+    # signs
+    N3 = N - N2 - 2
+    
+    bottsigns = 2.*bradii+height
+    lsign = height*hsigns
+    # up
+    N32 = int(N3/2) 
+    circle = p_circle(lsign/2., N=N32)
+    trib = circle + [0.,0.] 
+
+    buoy[N2+1:N2+1+N32,:] = trib + [bottsigns+3.2*lsign,0.]
+
+    # up
+    N323 = N - N32 - N2 - 2
+    trid = p_circle(lsign/2., N=N32)
+    trib = circle + [0.,0.] 
+    buoy[N2+N32+2:N,:] = trib + [bottsigns+2.*lsign,0.]
+
+    # painting it
+    Height = np.max(buoy1v[:,0])
+
+    Ncut, third1 = cut_ypolygon(buoy1v, yval=Height/3., keep='bottom')
+    Ncut, third2 = cut_between_ypolygon(buoy1v, yval1=Height/3., yval2=Height*2./3.)
+    Ncut, third3 = cut_ypolygon(buoy1v, yval=Height*2./3., keep='above')
+
+    buoy = ma.masked_equal(buoy, gen.fillValueF)
+
+    buoysecs = ['buoy', 'sign1', 'sign2', 'halfk', 'halfy']
+    buoydic = {'buoy': [buoy[0:N2,:],'-','k',1.5],                                   \
+      'sign1': [buoy[N2+1:N2+N32+1,:],'-','k',1.5],                                  \
+      'sign2': [buoy[N2+N32+2:N,:],'-','k',1.5], 'third1': [third1, '-', 'k', 1.],   \
+      'third2': [third2, '-', 'r', 1.], 'third3': [third3, '-', 'k', 1.]}
+
+    return buoy, buoysecs, buoydic
+
 ####### ####### ##### #### ### ## #
 # Plotting