blob: 5c200090eddcbcd43d6d3fc52170fc5e3cc803fa (
plain) (
blame)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
|
SUBROUTINE sla_DCS2C (A, B, V)
*+
* - - - - - -
* D C S 2 C
* - - - - - -
*
* Spherical coordinates to direction cosines (double precision)
*
* Given:
* A,B dp spherical coordinates in radians
* (RA,Dec), (Long,Lat) etc
*
* Returned:
* V dp(3) x,y,z unit vector
*
* The spherical coordinates are longitude (+ve anticlockwise
* looking from the +ve latitude pole) and latitude. The
* Cartesian coordinates are right handed, with the x axis
* at zero longitude and latitude, and the z axis at the
* +ve latitude pole.
*
* P.T.Wallace Starlink October 1984
*
* Copyright (C) 1995 Rutherford Appleton Laboratory
*-
IMPLICIT NONE
DOUBLE PRECISION A,B,V(3)
DOUBLE PRECISION COSB
COSB=COS(B)
V(1)=COS(A)*COSB
V(2)=SIN(A)*COSB
V(3)=SIN(B)
END
|
