ce_cent2llh_ellips.f File Reference

Go to the source code of this file.

Functions/Subroutines
subroutine	cecent2llh (xyz, llh)

Function/Subroutine Documentation

cecent2llh()

subroutine cecent2llh	(	type(coord)	xyz,
		type(coord)	llh
	)

Definition at line 22 of file ce_cent2llh_ellips.f.

References coord_types(), d, and d0.

Referenced by ctranscoord2().

!       xyz:  /coord/ structue. Input. The coordinate system is
!            such that the origin is at  the center of the earth
!            x-axis is directed to (0, 0) latitude and longitude.
!            y-axis is directed to (0, 90) latitude and longitude.
!            z-axsi is directed to the north pole.
!       llh:  /coord/ structure. Output. to contain data in latitude,
!                                       longitude, height.
!    ***  note  ***
!           llh can be the same as xyz.  time component is unchanged
!
      implicit none
c----      include '../../Zglobalc.h'
#include  "Zglobalc.h"
c----      include '../Zcoord.h'
#include  "Zcoord.h"
c----      include 'Zearth.h'
#include  "Zearth.h"
      type(coord)::llh,  xyz
      type(coord):: temp
!
      real*8 n, cosphi, sinphi, cosl, sinl, hold
      integer i
!        start with ecentricity 0 approximation
         temp%h = 0.
         n = eradius
         do i = 1, 6
            hold = temp%h
            temp%h = sqrt( xyz%r(1)**2 + xyz%r(2)**2 +(xyz%r(3)/
     *      (1.d0 - n*ecen2/(n+temp%h)))**2 ) - n
             cosphi = sqrt( (xyz%r(1)/(n+temp%h))**2 +
     *                      (xyz%r(2)/(n+temp%h))**2 ) 
             sinphi = xyz%r(3)/(n*one_ecen2 + temp%h)
!      
             cosl = xyz%r(1)/(n+temp%h)/cosphi
             sinl = xyz%r(2)/(n+temp%h)/cosphi
             n = eradius /sqrt(1.d0 - ecen2*sinphi**2)
            if(abs(temp%h -hold) .lt. 1.d-2) goto 100   ! abs error < 1 mm 
         enddo 
 100     continue
!        write(*, *) ' i=', i  !  <i> = 4 
!
         temp%lat = atan2(sinphi, cosphi)*todeg
         temp%long = atan2(sinl, cosl)*todeg
         llh%lat = temp%lat
         llh%long = temp%long
         llh%h = temp%h
         llh%sys = coord_types(2)

Here is the call graph for this function:

Here is the caller graph for this function: