[Dart-dev] [6177] DART/branches/development/models/PBL_1d/src/module_map_utils.F90: Change calls from ALOG and ALOG10 to plain LOG and LOG10.

[nancy at ucar.edu]

Revision: 6177
Author:   nancy
Date:     2013-05-29 10:02:21 -0600 (Wed, 29 May 2013)
Log Message:
-----------
Change calls from ALOG and ALOG10 to plain LOG and LOG10.
The compiler will take care of calling either the single or
double precision versions of these functions.

Modified Paths:
--------------
    DART/branches/development/models/PBL_1d/src/module_map_utils.F90

-------------- next part --------------
Modified: DART/branches/development/models/PBL_1d/src/module_map_utils.F90
===================================================================
--- DART/branches/development/models/PBL_1d/src/module_map_utils.F90	2013-05-29 16:01:22 UTC (rev 6176)
+++ DART/branches/development/models/PBL_1d/src/module_map_utils.F90	2013-05-29 16:02:21 UTC (rev 6177)
@@ -631,10 +631,10 @@
     ! intersects the Earth's surface at each of the distinctly different
     ! latitudes
     IF (ABS(truelat1-truelat2) .GT. 0.1_r8) THEN
-      cone = ALOG10(COS(truelat1*rad_per_deg)) - &
-             ALOG10(COS(truelat2*rad_per_deg))
-      cone = cone /(ALOG10(TAN((45.0_r8 - ABS(truelat1)/2.0_r8) * rad_per_deg)) - &
-             ALOG10(TAN((45.0_r8 - ABS(truelat2)/2.0_r8) * rad_per_deg)))        
+      cone = LOG10(COS(truelat1*rad_per_deg)) - &
+             LOG10(COS(truelat2*rad_per_deg))
+      cone = cone /(LOG10(TAN((45.0_r8 - ABS(truelat1)/2.0_r8) * rad_per_deg)) - &
+             LOG10(TAN((45.0_r8 - ABS(truelat2)/2.0_r8) * rad_per_deg)))        
     ELSE
        cone = SIN(ABS(truelat1)*rad_per_deg )  
     ENDIF
@@ -795,7 +795,7 @@
 
     proj%rsw = 0.0_r8
     IF (proj%lat1 .NE. 0.0_r8) THEN
-      proj%rsw = (ALOG(TAN(0.5_r8*((proj%lat1+90.0_r8)*rad_per_deg))))/proj%dlon
+      proj%rsw = (LOG(TAN(0.5_r8*((proj%lat1+90.0_r8)*rad_per_deg))))/proj%dlon
     ENDIF
     RETURN
   END SUBROUTINE set_merc
@@ -816,7 +816,7 @@
     IF (deltalon .LT. -180.0_r8) deltalon = deltalon + 360.0_r8
     IF (deltalon .GT. 180.0_r8) deltalon = deltalon - 360.0_r8
     i = 1.0_r8 + (deltalon/(proj%dlon*deg_per_rad))
-    j = 1.0_r8 + (ALOG(TAN(0.5_r8*((lat + 90.0_r8) * rad_per_deg)))) / &
+    j = 1.0_r8 + (LOG(TAN(0.5_r8*((lat + 90.0_r8) * rad_per_deg)))) / &
            proj%dlon - proj%rsw
 
     RETURN