earth_radius_mod.F90

Go to the documentation of this file.

 
 
module earth_radius_mod
  implicit none
 
  private
  public :: earth_radius
 
contains
 
  SUBROUTINE earth_radius(IGDTMPL, IGDTLEN, RADIUS, ECCEN_SQUARED)
    IMPLICIT NONE
 
    INTEGER,                INTENT(IN   ) :: igdtlen
    INTEGER,                INTENT(IN   ) :: igdtmpl(igdtlen)
 
    REAL,                   INTENT(  OUT) :: eccen_squared
    REAL,                   INTENT(  OUT) :: radius
 
    REAL                                  :: flat
    REAL                                  :: major_axis, minor_axis
 
    SELECT CASE (igdtmpl(1))
    CASE (0)
       radius        = 6367470.0
       eccen_squared = 0.0
    CASE (1)  ! USER SPECIFIED SPHERICAL
       radius        = float(igdtmpl(3))/float(10**igdtmpl(2))
       eccen_squared = 0.0
    CASE (2)  ! IAU 1965
       radius        = 6378160.0      ! SEMI MAJOR AXIS
       flat          = 1.0/297.0      ! FLATTENING
       eccen_squared = (2.0*flat) - (flat**2)
    CASE (3)  ! USER SPECIFIED ELLIPTICAL (KM)
       major_axis    = float(igdtmpl(5))/float(10**igdtmpl(4))
       major_axis    = major_axis * 1000.0
       minor_axis    = float(igdtmpl(7))/float(10**igdtmpl(6))
       minor_axis    = minor_axis * 1000.0
       eccen_squared = 1.0 - (minor_axis**2 / major_axis**2)
       radius        = major_axis
    CASE (4)  ! IAG-GRS80 MODEL
       radius        = 6378137.0      ! SEMI MAJOR AXIS
       flat          = 1.0/298.2572   ! FLATTENING
       eccen_squared = (2.0*flat) - (flat**2)
    CASE (5)  ! WGS84 DATUM
       radius        = 6378137.0      ! SEMI MAJOR AXIS
       eccen_squared = 0.00669437999013
    CASE (6)
       radius        = 6371229.0
       eccen_squared = 0.0
    CASE (7)  ! USER SPECIFIED ELLIPTICAL (M)
       major_axis    = float(igdtmpl(5))/float(10**igdtmpl(4))
       minor_axis    = float(igdtmpl(7))/float(10**igdtmpl(6))
       eccen_squared = 1.0 - (minor_axis**2 / major_axis**2)
       radius        = major_axis
    CASE (8)
       radius        = 6371200.0
       eccen_squared = 0.0
    CASE DEFAULT
       radius        = -9999.
       eccen_squared = -9999.
    END SELECT
    !
    RETURN
    !
  END SUBROUTINE earth_radius
end module earth_radius_mod