A sphere with three ocean waves in differing shades of blue. Next to the sphere, there is the bolded text "Ocean Color". Under that, there are two acronyms, separated by a vertical pipe: OB.DAAC (Ocean Biology Distributed Active Archive Center) and OBPG (Ocean Biology Processing Group).
Lorem
ipsum
dolor
sit
amet
Toggle navigation
Ocean Color Science Software
Jump to content
ocssw
V2022
web
ocssw
ocssw_src
oel_hdf4
libnav
setsun.f
Go to the documentation of this file.
1
SUBROUTINE
setsun
(T,ES)
2
C VERSION 1/31/87
3
C PURPOSE
4
C SETS UP SOME SOLAR ORBITAL ELEMENTS OF THE BUILT-IN LUNI-SOLAR
5
C EPHEMERIDES
6
C INPUT ARGUMENTS
7
C T = CURRENT JULIAN TIME (SEC)
8
C OUTPUT ARGUMENTS
9
C ES = MEAN ORBITAL ELEMENTS OF THE SUN IN EARTH MEAN EQUATOR
10
C AND EQUINOX OF DATE
11
C (1) = A, SEMI-MAJOR AXIS (KM)
12
C (2) = E, ECCENTRICITY
13
C (3) = I, INCLINATION (RAD)
14
C (4) = CAPW, LONGITUDE OF ASCENDING NODE, =0 BY DEFINITION
15
C (5) = W, ARGUMENT OF PERIAPSIS (RAD)
16
C CALL SUBROUTINES
17
C NONE
18
C REFERENCES
19
C JPL EM 312/86-153, 20 APRIL 1987
20
C EXPLANATORY SUPPLEMENT TO THE ASTRONOMICAL EPHEMERIS AND THE
21
C AMERICAN EPHEMERIS AND NAUTICAL ALMANAC
22
C ANALYSIS
23
C JOHNNY H. KWOK - JPL
24
C PROGRAMMER
25
C JOHNNY H. KWOK - JPL
26
C MODIFICATIONS
27
C NONE
28
C COMMENTS
29
C NONE
30
C
31
IMPLICIT DOUBLE PRECISION
(a-h,o-z)
32
dimension es(7)
33
DATA
rjd/2.41502d6/
34
DATA
dts/8.64d4/
35
dt=t/dts-rjd
36
dt1=dt*1.d-4
37
dt2=dt1*dt1
38
dt3=dt1*dt2
39
es(1)=1.496d8
40
es(2)=1.675104d-2-1.1444d-5*dt1-9.4d-9*dt2
41
es(3)=.4093197474d0-6.217910d-5*dt1-2.1468d-9*dt2+1.7977d-10*dt3
42
es(4)=0.d0
43
es(5)=4.908229653d0+8.2149855d-7*dt+5.9167d-7*dt2+1.22d-9*dt3
44
RETURN
45
END
setsun
subroutine setsun(T, ES)
Definition:
setsun.f:2
