Changeset 3198 for palm/trunk/UTIL

Timestamp:

Aug 15, 2018 9:23:10 AM (6 years ago)

Author:

sward

Message:

Added MAS end time, used time_since_reference_point, corrected tolerance_dp in nav_mesh

File:

• palm/trunk/UTIL/nav_mesh/nav_mesh.f90 (modified) (4 diffs)

palm/trunk/UTIL/nav_mesh/nav_mesh.f90

@@ -25 +25 @@
 ! -----------------!
 ! $Id$
+! Reduced tolerance_dp to 3 entries, fixed its initialization
+! 
+! 3168 2018-07-25 06:40:29Z sward
 ! Updated NetCDF ororgraphy and building input
 ! 
@@ -95 +98 @@
     REAL(wp)    ::  start            !< variable for CPU time measurement
 
-    REAL(wp), DIMENSION(0:4) ::  tolerance_dp  !< tolerance in Douglas-Poicker algorithm
+    REAL(wp), DIMENSION(0:2) ::  tolerance_dp = 999999.0_wp  !< tolerance in Douglas-Peucker algorithm
 
     REAL(wp), DIMENSION(:,:), ALLOCATABLE ::  obstacle_height  !< height of obstacles
@@ -1338 +1341 @@
 !
 !--    If tolerance_dp was not set, put in default values
-       IF ( tolerance_dp(0) == 999999.0_wp ) tolerance_dp(0) = SQRT(dx*dy)*2.99
-       IF ( tolerance_dp(1) == 999999.0_wp ) tolerance_dp(1) = SQRT(dx*dy)*1.99
-       IF ( tolerance_dp(2) == 999999.0_wp ) tolerance_dp(2) = SQRT(dx*dy)*1.41
-       IF ( tolerance_dp(3) == 999999.0_wp ) tolerance_dp(3) = SQRT(dx*dy)*.7
-       IF ( tolerance_dp(4) == 999999.0_wp ) tolerance_dp(4) = SQRT(dx*dy)*.35
+       DO i = 0, 2
+          IF ( tolerance_dp(i) == 999999.0_wp ) THEN
+             tolerance_dp(i) = SQRT(dx*dy)*1.41/(2**i)
+          ELSE
+             tolerance_dp(i) = tolerance_dp(i)*SQRT(dx*dy)
+          ENDIF
+       ENDDO
+
 !
 !--    Allocate arrays
@@ -2784 +2790 @@
 !--    of vertices would fall below 4 due to this procedure, the tolerance
 !--    for the algorithm is reduced and it is run again. 
-       DO i_sc = 0, 4
+       DO i_sc = 0, 2
           CALL simplify_polygon(1,nov+1,tolerance_dp(i_sc))
           i_cn = 0