Changeset 4568 for palm/trunk/UTIL/inifor/src/inifor_util.f90

Timestamp:

Jun 19, 2020 11:56:30 AM (4 years ago)

Author:

eckhard

Message:

Handle COSMO soil data with and without additional surface temperature

File:

• palm/trunk/UTIL/inifor/src/inifor_util.f90 (modified) (2 diffs)

palm/trunk/UTIL/inifor/src/inifor_util.f90

@@ -26 +26 @@
 ! -----------------
 ! $Id$
+! Added function for checking floating point equality
+! 
+! 
+! 4523 2020-05-07 15:58:16Z eckhard
 ! respect integer working precision (iwp) specified in inifor_defs.f90
 ! 
@@ -462 +466 @@
  END SUBROUTINE
 
+
+!------------------------------------------------------------------------------!
+! Description:
+! ------------
+!> Check for exact or near equality of two floating point numbers. Inspired by
+!> https://randomascii.wordpress.com/2012/02/25/comparing-floating-point-numbers-2012-edition/
+!------------------------------------------------------------------------------!
+ LOGICAL FUNCTION nearly_equal(a, b, max_rel_diff)
+
+    REAL(wp), INTENT(IN) ::  a, b, max_rel_diff
+    REAL(wp)             ::  diff, mag
+ 
+    diff = ABS( a - b ) 
+    mag = MAX( ABS(a), ABS(b) )
+    nearly_equal = ( diff .LE. mag * max_rel_diff )
+
+ END FUNCTION nearly_equal
+
  END MODULE inifor_util