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
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
218
219
220
221
222
223
224
225
226
227
228
229
230
231
232
233
234
235
236
237
238
239
240
241
242
243
244
245
246
247
248
249
250
251
252
253
254
255
256
257
258
//
// GENERATED FILE
//
use super::*;
use crate::SpiceContext;
use f2rust_std::*;
const NCONST: i32 = 1;
const BEGEL1: i32 = (NCONST + 1);
/// Read SPK record from segment, type 14
///
/// Read a single data record from a type 14 SPK segment.
///
/// # Required Reading
///
/// * [SPK](crate::required_reading::spk)
///
/// # Brief I/O
///
/// ```text
/// VARIABLE I/O DESCRIPTION
/// -------- --- --------------------------------------------------
/// HANDLE I Handle of the open SPK file.
/// DESCR I Descriptor of the segment with the desired record.
/// ET I Epoch used to identify the desired record.
/// RECORD O The desired type 14 SPK record.
/// ```
///
/// # Detailed Input
///
/// ```text
/// HANDLE is the handle of the open SPK file which contains the
/// segment of interest.
///
/// DESCR is the descriptor for a type 14 SPK segment that contains
/// the record of interest.
///
/// ET is the target epoch used to determine the particular
/// record to be obtained from the SPK segment.
/// ```
///
/// # Detailed Output
///
/// ```text
/// RECORD is the record from the specified segment which,
/// when evaluated at epoch ET, will give the state
/// (position and velocity) of some body, relative
/// to some center, in some inertial reference frame.
/// ```
///
/// # Exceptions
///
/// ```text
/// 1) It is assumed that the descriptor and handle supplied are
/// for a properly constructed type 14 segment. No checks are
/// performed to ensure this.
///
/// 2) If the input ET value is not within the range specified
/// in the segment descriptor, the error SPICE(TIMEOUTOFBOUNDS)
/// is signaled.
///
/// 3) If any other error occurs while looking up SPK data, the error
/// is signaled by a routine in the call tree of this routine.
/// ```
///
/// # Files
///
/// ```text
/// See argument HANDLE.
/// ```
///
/// # Particulars
///
/// ```text
/// This subroutine will read a single record from a type 14 SPK
/// segment. The record read will provide the data necessary to
/// compute the state for a some body in some inertial frame at epoch
/// ET.
///
/// See the SPK Required Reading file for a description of the
/// structure of a type 14 SPK segment.
/// ```
///
/// # Examples
///
/// ```text
/// The data returned by the SPKRnn routine is in a raw form, taken
/// directly from the segment. As such, it will be not be directly
/// useful to a user unless they have a complete understanding of the
/// structure of the data type. Given that understanding, however,
/// the SPKRnn routines could be used to "dump" and check segment data
/// for a particular epoch, as in the example which follows.
///
///
/// C
/// C Get a segment applicable to a specified body and epoch.
/// C
/// CALL SPKSFS ( BODY, ET, HANDLE, DESCR, IDENT, FOUND )
///
/// C
/// C Look at parts of the descriptor.
/// C
/// CALL DAFUS ( DESCR, 2, 6, DCD, ICD )
/// CENTER = ICD( 2 )
/// REF = ICD( 3 )
/// TYPE = ICD( 4 )
///
/// IF ( TYPE .EQ. 14 ) THEN
/// CALL SPKR14 ( HANDLE, DESCR, ET, RECORD )
/// .
/// . Look at the RECORD data.
/// .
/// END IF
/// ```
///
/// # Restrictions
///
/// ```text
/// 1) This subroutine should not be called directly by a casual
/// user. It is intended for use by the subroutine SPKPV, and
/// certain tests for error conditions are not performed here, as
/// SPKPV will have already performed them.
/// ```
///
/// # Author and Institution
///
/// ```text
/// J. Diaz del Rio (ODC Space)
/// K.R. Gehringer (JPL)
/// ```
///
/// # Version
///
/// ```text
/// - SPICELIB Version 1.0.1, 14-APR-2021 (JDR)
///
/// Added IMPLICIT NONE statement.
///
/// Edited the header to comply with NAIF standard. Updated entry
/// #3 in $Exceptions section. Moved SPK required reading from
/// $Literature_References to $Required_Reading section.
///
/// - SPICELIB Version 1.0.0, 10-MAR-1995 (KRG)
/// ```
pub fn spkr14(
ctx: &mut SpiceContext,
handle: i32,
descr: &[f64],
et: f64,
record: &mut [f64],
) -> crate::Result<()> {
SPKR14(handle, descr, et, record, ctx.raw_context())?;
ctx.handle_errors()?;
Ok(())
}
//$Procedure SPKR14 ( Read SPK record from segment, type 14 )
pub fn SPKR14(
HANDLE: i32,
DESCR: &[f64],
ET: f64,
RECORD: &mut [f64],
ctx: &mut Context,
) -> f2rust_std::Result<()> {
let DESCR = DummyArray::new(DESCR, 1..);
let mut RECORD = DummyArrayMut::new(RECORD, 1..);
let mut VALUE: f64 = 0.0;
let mut ENDS: i32 = 0;
let mut INDX: i32 = 0;
let mut FOUND: bool = false;
//
// SPICELIB functions
//
//
// Local variables
//
// The number of constant values stored with a type 14 segment
// segment.
//
//
// The beginning location in the output record for the non constant
// segment data.
//
//
// Standard SPICE error handling.
//
if RETURN(ctx) {
return Ok(());
} else {
CHKIN(b"SPKR14", ctx)?;
}
//
// Check the request time against the time bounds in the segment
// descriptor. In order to get the right data back from the generic
// segment calls below, we need to be sure that the desired epoch
// falls within the bounds of the segment, as specified by the
// descriptor. The first two elements of the descriptor are the start
// time for the segment and the stop time for the segment,
// respectively.
//
if ((ET < DESCR[1]) || (ET > DESCR[2])) {
SETMSG(
b"Request time # is outside of descriptor bounds # : #.",
ctx,
);
ERRDP(b"#", ET, ctx);
ERRDP(b"#", DESCR[1], ctx);
ERRDP(b"#", DESCR[2], ctx);
SIGERR(b"SPICE(TIMEOUTOFBOUNDS)", ctx)?;
CHKOUT(b"SPKR14", ctx)?;
return Ok(());
}
//
// Fetch the constants and store them in the first part of
// the output RECORD.
//
SGFCON(
HANDLE,
DESCR.as_slice(),
1,
NCONST,
RECORD.subarray_mut(1),
ctx,
)?;
//
// Locate the time in the file less than or equal to the input ET,
// obtaining its index. This will allow us to retrieve the proper
// record.
//
SGFRVI(
HANDLE,
DESCR.as_slice(),
ET,
&mut VALUE,
&mut INDX,
&mut FOUND,
ctx,
)?;
//
// Fetch the appropriate record from the segment.
//
SGFPKT(
HANDLE,
DESCR.as_slice(),
INDX,
INDX,
RECORD.subarray_mut(BEGEL1),
std::slice::from_mut(&mut ENDS),
ctx,
)?;
CHKOUT(b"SPKR14", ctx)?;
Ok(())
}