correct conversion to radians

rgknox · Aug 3, 2022 · 8186d50 · 8186d50
1 parent 8b0ed78
commit 8186d50
Show file tree

Hide file tree

Showing 2 changed files with 6 additions and 5 deletions.
diff --git a/main/FatesConstantsMod.F90 b/main/FatesConstantsMod.F90
@@ -233,5 +233,6 @@ module FatesConstantsMod
 
   ! PI
   real(fates_r8), parameter, public :: pi_const = 3.14159265359_fates_r8
+  real(fates_r8), parameter, public :: rad_per_deg = pi_const/180.0_fates_r8
 
 end module FatesConstantsMod
diff --git a/main/FatesUtilsMod.F90 b/main/FatesUtilsMod.F90
@@ -86,7 +86,7 @@ end subroutine check_var_real
   real(r8) function GreatCircleDist(slons,slonf,slats,slatf)
 
      use FatesConstantsMod, only : earth_radius_eq &
-                                 , pi_const 
+                                 , rad_per_deg 
      implicit none
 
      !----- Local variables. ----------------------------------------------------------------!
@@ -107,10 +107,10 @@ real(r8) function GreatCircleDist(slons,slonf,slats,slatf)
      !---------------------------------------------------------------------------------------!
 
      !----- Convert the co-ordinates to double precision and to radians. --------------------!
-     lons = slons * pi_const
-     lonf = slonf * pi_const
-     lats = slats * pi_const
-     latf = slatf * pi_const
+     lons = slons * rad_per_deg
+     lonf = slonf * rad_per_deg
+     lats = slats * rad_per_deg
+     latf = slatf * rad_per_deg
      dlon = lonf - lons
      dlat = latf - lats