Changeset 2558 in lmdz_wrf

Timestamp:

May 25, 2019, 7:20:05 PM (6 years ago)

Author:

lfita

Message:

Filling up all the cardinal buoys

File:

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

trunk/tools/geometry_tools.py

@@ -1677 +1677 @@
     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], 'halfk': [halfup, '-', 'k', 1.],    \
-      'halfy': [halfdown, '-', '#FFFF00', 1.]}
+      'sign2': [buoy[N2+N32+2:N,:],'-','k',1.5], 'half1': [halfup, '-', 'k', 1.],    \
+      'half2': [halfdown, '-', '#FFFF00', 1.]}
 
     return buoy, buoysecs, buoydic
@@ -1730 +1730 @@
     buoy = ma.masked_equal(buoy, gen.fillValueF)
 
-    buoysecs = ['buoy', 'sign1', 'sign2']
+    buoysecs = ['buoy', 'sign1', 'sign2', 'third1', 'third2', 'third3']
     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],                                     \
-      'third1y': [halfup, '-', '#FFFF00', 1.], 'thirdk': [halfbtw, '-', 'k', 1.],    \
-      'third2y': [halfdown, '-', '#FFFF00', 1.]}
+      'third1': [halfup, '-', 'k', 1.], 'third2': [halfbtw, '-', '#FFFF00', 1.],     \
+      'third3': [halfdown, '-', 'k', 1.]}
 
     return buoy, buoysecs, buoydic
@@ -1778 +1778 @@
     buoy[N2+N32+2:N,:] = trib + [bottsigns+0.9*lsign,0.]
 
+    # painting it
+    Height = np.max(buoy1v[:,0])
+
+    Ncut, halfdown = cut_ypolygon(buoy1v, yval=Height/2., keep='bottom')
+    Ncut, halfup = cut_ypolygon(buoy1v, yval=Height/2., keep='above')
+
     buoy = ma.masked_equal(buoy, gen.fillValueF)
 
-    buoysecs = ['buoy', 'sign1', 'sign2']
+    buoysecs = ['buoy', 'sign1', 'sign2', 'half1', 'half2']
     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]}
+      'sign2': [buoy[N2+N32+2:N,:],'-','k',1.5], 'half1': [halfup, '-', '#FFFF00', 1.], \
+      'half2': [halfdown, '-', 'k', 1.]}
 
     return buoy, buoysecs, buoydic
@@ -1826 +1833 @@
     buoy[N2+N323+2:N,:] = trib + [bottsigns+1.*lsign,0.]
 
+    # painting it
+    Height = np.max(buoy1v[:,0])
+
+    Ncut, halfdown = cut_ypolygon(buoy1v, yval=Height/3., keep='bottom')
+    Ncut, halfbtw1 = cut_between_ypolygon(buoy1v, yval1=Height/3., yval2=Height*2./3.)
+    Ncut, halfup = cut_ypolygon(buoy1v, yval=Height*2./3., keep='above')
+    Ncut, halfbtw2 = cut_between_ypolygon(buoy1v, yval1=Height*2./3., yval2=Height)
+
     buoy = ma.masked_equal(buoy, gen.fillValueF)
 
-    buoysecs = ['buoy', 'sign1', 'sign2']
+    buoysecs = ['buoy', 'sign1', 'sign2', 'third1', 'third2', 'third3']
     buoydic = {'buoy': [buoy[0:N2,:],'-','k',1.5],                                   \
+      'third1': [halfdown, '-', '#FFFF00', 1.], 'third2': [halfbtw1, '-', 'k', 1.],  \
+      'third3': [halfup, '-', '#FFFF00', 1.],                                        \
       'sign1': [buoy[N2+1:N2+N32+1,:],'-','k',1.5],                                  \
-      'sign2': [buoy[N2+N32+2:N,:],'-','k',1.5]}
+      'sign2': [buoy[N2+N32+2:N,:],'-','#FFFF00',1.5]}
 
     return buoy, buoysecs, buoydic