return -99 rather than 1e-99 in safe logs

Author: evbauer

math/public/math_lib_crmath.f90

@@ -155,7 +155,7 @@ elemental function safe_log_ (x) result (log_x)
 
     if (.NOT. IEEE_IS_FINITE(x)) then
 
-       log_x = 1.E-99_dp
+       log_x = -99._dp
 
     else
 
@@ -174,7 +174,7 @@ elemental function safe_log10_ (x) result (log10_x)
 
     if (.NOT. IEEE_IS_FINITE(x)) then
 
-       log10_x = 1.E-99_dp
+       log10_x = -99._dp
 
     else
 

math/public/math_lib_intrinsic.f90

@@ -173,7 +173,7 @@ elemental function safe_log_ (x) result (log_x)
 
     if (.NOT. IEEE_IS_FINITE(x)) then
 
-       log_x = 1.E-99_dp
+       log_x = -99._dp
 
     else
 
@@ -192,7 +192,7 @@ elemental function safe_log10_ (x) result (log10_x)
 
     if (.NOT. IEEE_IS_FINITE(x)) then
 
-       log10_x = 1.E-99_dp
+       log10_x = -99._dp
 
     else