COSMOS v7.655  COSMOSv7655
(AirShowerMC)
Functions/Subroutines
cllh2e_cent_ellips.f File Reference

Go to the source code of this file.

Functions/Subroutines

subroutine cllh2ecent (llh, xyz)
 

Function/Subroutine Documentation

◆ cllh2ecent()

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

Definition at line 17 of file cllh2e_cent_ellips.f.

References coord_types(), and cos.

Referenced by ctranscoord2().

17 ! llh: /coord/ structure. input. containing data in latitude,
18 ! longitude, height.
19 ! xyz: /coord/ structue. output. The coordinate system is
20 ! such that the origin is at the center of the earth
21 ! x-axis is directed to (0, 0) latitude and longitude.
22 ! y-axis is directed to (0, 90) latitude and longitude.
23 ! z-axsi is directed to the north pole.
24 !
25 ! xyz.r(1) is x coordinate value in m
26 ! xyz.r(2) is y
27 ! xyz.r(3) is z
28 ! *** note ***
29 ! xyz can be the same as llh. time component is unchanged
30 !
31  implicit none
32 c---- include '../../Zglobalc.h'
33 #include "Zglobalc.h"
34 c---- include '../../Tracking/Zcoord.h'
35 #include "Zcoord.h"
36 c---- include 'Zearth.h'
37 #include "Zearth.h"
38  type(coord)::llh, xyz
39  type(coord):: temp
40 !
41  real*8 n
42 !
43  n = eradius/sqrt(1. -eecen2 * sin(llh%lat*torad)**2)
44  temp%x = (n + llh%h) *cos(llh%lat*torad)* cos(llh%long*torad)
45  temp%y = (n + llh%h) *cos(llh%lat*torad)* sin(llh%long*torad)
46  xyz%r(3) = (n* eone_ecen2 + llh%h)*sin(llh%lat * torad)
47  xyz%r(1) = temp%x
48  xyz%r(2) = temp%y
49  xyz%sys = coord_types(1)
cos
! for length to thickness conversion or v v ! integer maxnodes real Hinf ! This is used when making table for dim simulation ! The slant length for vertical height to km is ! divided by LenStep m steps ! It can cover the slant length of about km for cos
Definition: Zatmos.h:8
xyz
latitude latitude this system is used *****************************************************************! type coord sequence union map real z z in m endmap xyz map real ! latitude in deg is to the north ! longitude in deg is to the east *h ! height in m endmap llh map real ! polar angle ! azimuthal angle *radius ! radial distance endmap sph endunion character *sys ! which system xyz
Definition: Zcoord.h:25
llh
latitude latitude this system is used *****************************************************************! type coord sequence union map real z z in m endmap xyz map real ! latitude in deg is to the north ! longitude in deg is to the east *h ! height in m endmap llh map real ! polar angle ! azimuthal angle *radius ! radial distance endmap sph endunion character *sys ! which system llh
Definition: Zcoord.h:25
coord_types
integer max_coord_types * coord_types(2)/'llh'/
coord
Definition: Zcoord.h:43
n
integer n
Definition: Zcinippxc.h:1
Here is the call graph for this function:
Here is the caller graph for this function: