Changeset 2333 in lmdz_wrf

Timestamp:

Feb 12, 2019, 3:30:26 PM (6 years ago)

Author:

lfita

Message:

Avoding to remov all polygons when applying `clean_polygons'

Location:

trunk/tools

Files:

• module_ForDiagnostics.f90 (modified) (1 diff)
• module_scientific.f90 (modified) (3 diffs)

trunk/tools/module_ForDiagnostics.f90

@@ -1056 +1056 @@
     ! Homogenizing indices of the ranges
     Lranges = iranges /= 0
-    CALL polygons(.FALSE., d1, d2, Lranges, ranges, Nranges)
+    CALL polygons(.TRUE., d1, d2, Lranges, ranges, Nranges)
 
 !    IF (TRIM(face) == 'WE') THEN

trunk/tools/module_scientific.f90

@@ -2597 +2597 @@
 
   ! The mathematical maximum woiuld be dx*dy/4, but let's be optimistic... (sorry Jero)
-  Nppt = dx*dy/100
+  Nppt = dx*dy/10000
 
   IF (ALLOCATED(borders)) DEALLOCATE(borders)
@@ -2618 +2618 @@
   ! Filling with the points of all the space with .TRUE.
   Npts = COUNT(boolmat)
+  PRINT * ,'Lluis Npts:', Npts, 'Nppt:', Nppt
 
   IF (ALLOCATED(points)) DEALLOCATE(points)
@@ -2753 +2754 @@
     END IF
   END DO
+
+  IF (NnotPol == Npols) THEN
+    PRINT *,'  ' // TRIM(fname) // ": avoiding to remove all polygons !!"
+    NnotPol = 0
+  END IF
 
   IF (dbg) THEN