add diag

Author: scivision

src/aer.f90

@@ -46,7 +46,6 @@
 
 real(wp) :: x0,y0,z0, e,n,u,dx,dy,dz
 
-
 !> Origin of the local system in geocentric coordinates.
 call geodetic2ecef(lat0, lon0, alt0,x0, y0, z0, spheroid,deg)
 !> Convert Local Spherical AER to ENU

src/ecef.f90

@@ -144,7 +144,6 @@
   sinLon = sin(lon)
 endif
 
-
 x = (N + alt) * cosLat * cosLon
 y = (N + alt) * cosLat * sinLon
 z = (N * (ell%SemiminorAxis / ell%SemimajorAxis)**2 + alt) * sinLat

src/enu.f90

@@ -124,7 +124,7 @@
 
 !> radians
 elev = atan2(u, r)
-az = modulo(atan2(e, n), 2._wp * pi)
+az = modulo(atan2(e, n), 2 * pi)
 
 d=.true.
 if (present(deg)) d = deg

src/maptran.F90

@@ -165,14 +165,14 @@ end function haversine
 elemental real(wp) function degrees(rad)
 real(wp), intent(in) :: rad
 
-degrees = 180._wp / pi * rad
+degrees = 180 / pi * rad
 end function degrees
 
 
 elemental real(wp) function radians(deg)
 real(wp), intent(in) :: deg
 
-radians = pi / 180._wp * deg
+radians = pi / 180 * deg
 end function radians
 
 end module maptran

src/tests/test_mod.f90

@@ -416,23 +416,23 @@ subroutine test_aer2ecef(az, el, rng, lat, lon, alt)
 real(wp) :: x,y,z
 
 call aer2ecef(az, el, rng, lat, lon, alt, x, y, z)
-call assert_allclose([x, y, z], [xl,yl,zl])
+call assert_allclose([x, y, z], [xl,yl,zl], err_msg='aer2ecef1')
 
 call aer2ecef(0._wp, -90._wp, 1._wp, 0._wp, 0._wp, -1._wp, x, y, z)
-call assert_allclose([x, y, z], [ea-1._wp, 0._wp, 0._wp], atol=1e-6_wp)
+call assert_allclose([x, y, z], [ea-1._wp, 0._wp, 0._wp], atol=1e-6_wp, err_msg='aer2ecef2')
 
 call aer2ecef(0._wp, -90._wp, 1._wp, 0._wp, 90._wp, -1._wp, x, y, z)
-call assert_allclose([x,y,z],[0._wp, ea-1._wp, 0._wp], atol=1e-6_wp)
+call assert_allclose([x,y,z],[0._wp, ea-1._wp, 0._wp], atol=1e-6_wp, err_msg='aer2ecef3')
 
 call aer2ecef(0._wp, -90._wp, 1._wp, 90._wp, 0._wp, -1._wp, x, y, z)
-call assert_allclose([x,y,z],[0._wp, 0._wp, eb-1._wp], atol=1e-6_wp)
+call assert_allclose([x,y,z],[0._wp, 0._wp, eb-1._wp], atol=1e-6_wp, err_msg='aer2ecef4')
 
 call aer2ecef(0._wp, -90._wp, 1._wp, 0._wp, 90._wp, -1._wp, x, y, z)
-call assert_allclose([x,y,z],[0._wp, ea-1._wp, 0._wp], atol=1e-6_wp)
+call assert_allclose([x,y,z],[0._wp, ea-1._wp, 0._wp], atol=1e-6_wp, err_msg='aer2ecef5')
 
 call aer2ecef(0._wp, -90._wp, 1._wp, 0._wp, 45._wp, -1._wp, x, y, z)
 call assert_allclose([x,y,z], &
-        [(ea-1._wp)/sqrt(2._wp), (ea-1._wp)/sqrt(2._wp), 0._wp], atol=1e-6_wp)
+        [(ea-1._wp)/sqrt(2._wp), (ea-1._wp)/sqrt(2._wp), 0._wp], atol=1e-6_wp, err_msg='aer2ecef6')
 
 end subroutine test_aer2ecef