Changeset 2625 in lmdz_wrf

Timestamp:

Jun 23, 2019, 1:50:52 AM (5 years ago)

Author:

lfita

Message:

Advancing into the pile_polygons

File:

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

trunk/tools/geometry_tools.py

@@ -1587 +1587 @@
           nvertexa=iNpts, nvertexb=Npts, nvertexab=iNpts*Npts, polya=ipoly,          \
           polyb=poly)
+        # We're in C-mode !
+        inti = inti-1
+        intp = intp-1
+
+        # There is a need to sort the cutting points from the perspective of the
+        #   below polygon
+
         Npi=0
         for ip in range(Nint):
@@ -1601 +1608 @@
             if pip1 > Npts-1: pip1=0
 
+            print ip, ': Lluis iip, iip1', iip, iip1, 'pip, pip1', pip, pip1
+
             iface = ipoly[iip,:]-poly[pip,:]
             pface = poly[pip1,:]-poly[pip,:]
@@ -1613 +1622 @@
                 for iv in range(Npi,Npi+iip): newipoly.append(ipoly[iv,:])
                 Npi = Npi + iip
-                newipoly.append([gen.fillvalueF, gen.fillvalueF])
+                newipoly.append(pts[ip])
+                newipoly.append([gen.fillValueF, gen.fillValueF])
 
             elif efangle < 0.:
                 # Ending point of the cuted face to the left of the overlying polygon
-                newipoly.append([gen.fillvalueF, gen.fillvalueF])
+                newipoly.append([gen.fillValueF, gen.fillValueF])
+                newipoly.append(pts[ip])
                 newipoly.append(ipoly[iip1,:])
                 Npi = iip1
@@ -1623 +1634 @@
                 print errormsg
                 print '  ' + fname + ': No left point found !!'
-                print '    for pile face'
+                print '    crossing point number', ip, 'below polygon vertex:',      \
+                  inti[ip], 'above', intp[ip]
+                print '    for pile face', poly[pip,:], '->', poly[pip1,:],          \
+                  'and below face', ipoly[iip,:], '->', ipoly[iip1,:]
     
         pilepolygons[polyns[Npolys-1]] = newipoly
-       
+    
+    for polyn in polyns:
+        poly = pilepolygons[polyn]
+        poly = ma.masked_equal(poly, gen.fillValueF)
+        pilepolygons[polyn] = poly
+
     return pilepolygons