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
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
|
include <error.h>
# T_ASTTIMES -- Print and record astronomical times for the given date.
procedure t_asttimes ()
int list # List of files
int header # Print header?
int year # Year
int month # Month
int day # Day
double zt # Zone time
double zone # Time zone from Greenwich
double longitude # Longitude for LMST
double ut # Universal time (output)
double epoch # Epoch (output)
double jd # Julian day (output)
double lmst # Local mean siderial time (output)
int fd
char file[SZ_FNAME]
pointer obs
int clpopnu(), clplen(), clgfil(), clgeti(), btoi()
int open(), fscan(), nscan()
bool clgetb()
double clgetd(), obsgetd()
pointer obsopen()
begin
# Get parameters other than date.
list = clpopnu ("files")
header = btoi (clgetb ("header"))
call clgstr ("observatory", file, SZ_FNAME)
obs = obsopen (file)
if (header == YES)
call obslog (obs, "ASTTIMES", "timezone longitude", STDOUT)
zone = obsgetd (obs, "timezone")
longitude = obsgetd (obs, "longitude")
# If no files are given then get dates from the CL.
if (clplen (list) == 0) {
# Get and print times.
year = clgeti ("year")
month = clgeti ("month")
day = clgeti ("day")
zt = clgetd ("time")
call ast_times (year, month, day, zt, zone, longitude, ut, epoch,
jd, lmst, header)
# Record results in the parameter file.
call clputd ("ut", ut)
call clputd ("epoch", epoch)
call clputd ("jd", jd)
call clputd ("lmst", lmst)
} else {
# Scan each file in the list.
while (clgfil (list, file, SZ_FNAME) != EOF) {
iferr (fd = open (file, READ_ONLY, TEXT_FILE)) {
call erract (EA_WARN)
next
}
# Get and print times.
while (fscan (fd) != EOF) {
call gargi (year)
call gargi (month)
call gargi (day)
call gargd (zt)
if (nscan() < 4)
next
call ast_times (year, month, day, zt, zone, longitude, ut,
epoch, jd, lmst, header)
}
call close (fd)
}
call clpcls (list)
}
call obsclose (obs)
end
# TIMES -- Print times.
procedure ast_times (year, month, day, zt, zone, longitude, ut, epoch, jd, lmst,
header)
int year # Year
int month # Month (1-12)
int day # Day of month
double zt # Zone time
double zone # Time zone
double longitude # Longitude
double ut # UT
double epoch # Epoch in 365.25 solar mean days
double jd # Julian date
double lmst # Mean Sidereal Time
int header # Print header?
char dow[3]
int d
double ast_date_to_julday(), ast_mst()
begin
# Determine day of the week in zone time.
if (zt < 0.) {
zt = zt + 24
day = day - 1
}
if (zt >= 24.) {
zt = zt - 24
day = day + 1
}
jd = ast_date_to_julday (year, month, day, zt)
call ast_julday_to_date (jd, year, month, day, zt)
call ast_date_to_epoch (year, month, day, zt, epoch)
call ast_day_of_week (jd, d, dow, 3)
# Determine UT, EPOCH, JD, and MST.
for (ut=zone; ut<-12.; ut=ut+24.)
;
for (; ut>=12.; ut=ut-24.)
;
ut = zt + ut
d = day
if (ut < 0.) {
ut = ut + 24
d = d - 1
}
if (ut >= 24.) {
ut = ut - 24
d = d + 1
}
jd = ast_date_to_julday (year, month, d, ut)
call ast_date_to_epoch (year, month, d, ut, epoch)
lmst = ast_mst (epoch, longitude)
# Print Times.
if (header == YES) {
call printf ("##%2s %3s %6s %10s %10s %10s %12s %10s\n")
call pargstr ("YR")
call pargstr ("MON")
call pargstr (" DAY ")
call pargstr ("ZT")
call pargstr ("UT")
call pargstr ("EPOCH")
call pargstr ("JD")
call pargstr ("LMST")
header = NO
}
call printf ("%4d %3d %2d %3s %10.1h %10.1h %10.5f %12.4f %10.1h\n")
call pargi (year)
call pargi (month)
call pargi (day)
call pargstr (dow)
call pargd (nint (zt*36000.0D0)/36000.0D0)
call pargd (nint (ut*36000.0D0)/36000.0D0)
call pargd (epoch)
call pargd (jd)
call pargd (nint (lmst*36000.0D0)/36000.0D0)
end
|
