Changeset 2553 in lmdz_wrf

Timestamp:

May 23, 2019, 4:28:00 PM (6 years ago)

Author:

lfita

Message:

Adding looking for polygon vertices within grid in `grid_within_polygon'

File:

• trunk/tools/module_scientific.f90 (modified) (1 diff)

trunk/tools/module_scientific.f90

@@ -5056 +5056 @@
             END IF
           END IF
+
+          ! Inverting process, look polygon within grid
+          IF (.NOT.within(ix,iy)) THEN
+            DO iv = 1, NvertexA
+              vertex = polygonA(iv,:)
+              IF (point_inside(vertex, Nvertexmax, xBvals(ix,iy,:))) THEN
+                Ngridsin = Ngridsin + 1
+                ! Getting coordinates
+                gridsin(Ngridsin,1) = ix
+                gridsin(Ngridsin,2) = iy
+                within(ix,iy) = .TRUE.
+                CYCLE
+              END IF
+            END DO
+          END IF
         END IF
       END DO