-
Notifications
You must be signed in to change notification settings - Fork 0
/
Copy pathcalculate_products.h
106 lines (89 loc) · 5.29 KB
/
calculate_products.h
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
/*
SLIDEM Processor: calculate_products.h
Copyright (C) 2024 Johnathan K Burchill
This program is free software: you can redistribute it and/or modify
it under the terms of the GNU General Public License as published by
the Free Software Foundation, either version 3 of the License, or
(at your option) any later version.
This program is distributed in the hope that it will be useful,
but WITHOUT ANY WARRANTY; without even the implied warranty of
MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
GNU General Public License for more details.
You should have received a copy of the GNU General Public License
along with this program. If not, see <https://www.gnu.org/licenses/>.
*/
#ifndef _CALCULATE_PRODUCTS_H_
#define _CALCULATE_PRODUCTS_H_
#include <stdint.h>
#include <stdbool.h>
#include "modified_oml.h"
void calculateProducts(const char satellite, uint8_t **hmDataBuffers, double *fpCurrent, double *vn, double *ve, double *vc, double *dipLatitude, double *faceplateVoltage, double f107Adj, int dayOfYear, double *ionEffectiveMass, double *ionDensity, double *ionDriftRaw, double *ionDrift, double *IonEffectiveMassError, double *ionDensityError, double *ionDriftError, double *fpAreaOML, double *rProbeOML, double *electronTemperature, double *spacecraftPotential, uint32_t *electronTemperatureSource, uint32_t *spacecraftPotentialSource, double *ionEffectiveMassTTS, uint32_t *mieffFlags, uint32_t *viFlags, uint32_t *niFlags, uint16_t *iterationCount, long nHmRecs, probeParams sphericalProbeParams, long *numberOfSlidemEstimates);
void getTeVs(const char satellite, uint8_t **hmDataBuffers, long hmTimeIndex, double *te, uint32_t *teSource, double *vs, uint32_t *vsSource);
int iterateEquations(double *niIO, double nil1b, double *vionsIO, double *mieffIO, uint32_t *viFlagIO, uint32_t *mieffFlagIO, uint32_t *niFlagIO, double *fpAreaIO, double *rProbeIO, double te, double vs, double faceplateVoltage, probeParams sphericalProbeParams, double ifp, double di, double vionsram, double mieffmodel, double qdlat, bool postProcessing, const char satellite);
void updateFlags(int iterations, double *mieffIO, double *mieffErrorIO, double *viIO, double *viErrorIO, double *niIO, double *niErrorIO, double *fpAreaIO, double *rProbeIO, double te, double vs, uint32_t teSource, uint32_t vsSource, double vionsram, double dipLat, double *vn, double *ve, double *vc, uint32_t *mieffFlagIO, uint32_t *viFlagIO, uint32_t *niFlagIO, long *slidemEstimatesIO, uint8_t **hmDataBuffers, long hmTimeIndex);
enum LP_FLAGS {
LP_HGN_OVERFLOW_LINEAR_BIAS = 1 << 2,
LP_LGN_OVERFLOW_LINEAR_BIAS = 1 << 3,
LP_HGN_OVERFLOW_RETARDED_BIAS = 1 << 4,
LP_LGN_OVERFLOW_RETARDED_BIAS = 1 << 5,
LP_HGN_ZERO_TRACKING_FAILED = 1 << 6,
LP_LGN_ZERO_TRACKING_FAILED = 1 << 7,
LP_HGN_LINEAR_BIAS_LESS_THAN_RETARDED_BIAS = 1 << 9,
LP_LGN_LINEAR_BIAS_LESS_THAN_RETARDED_BIAS = 1 << 10,
LP_HGN_LINEAR_BIAS_GREATER_THAN_5V_16BIT_OVERFLOW = 1 << 11,
LP_LGN_LINEAR_BIAS_GREATER_THAN_5V_16BIT_OVERFLOW = 1 << 12,
LP_HGN_ION_ADMITTANCE_GREATER_THAN_RETARDED_ADMITTANCE = 1 << 13,
LP_LGN_ION_ADMITTANCE_GREATER_THAN_RETARDED_ADMITTANCE = 1 << 14,
LP_HGN_ION_CURRENT_GREATER_THAN_RETARDED_CURRENT = 1 << 15,
LP_LGN_ION_CURRENT_GREATER_THAN_RETARDED_CURRENT = 1 << 16,
LP_HGN_RETARDED_ADMITTANCE_GREATER_THAN_LINEAR_ADMITTANCE = 1 << 17,
LP_LGN_RETARDED_ADMITTANCE_GREATER_THAN_LINEAR_ADMITTANCE = 1 << 18,
LP_HGN_RETARDED_CURRENT_GREATER_THAN_LINEAR_CURRENT = 1 << 19,
LP_LGN_RETARDED_CURRENT_GREATER_THAN_LINEAR_CURRENT = 1 << 20,
LP_NE_FROM_LGN_PROBE = 1 << 21
};
enum LP_SOURCE_PROBE {
LP_HGN_PROBE = 0b01,
LP_LGN_PROBE = 0b10,
LP_NO_PROBE = 0b00,
LP_BLENDED_PROBE = 0b11
};
enum LP_VALIDITY_MASKS {
LP_TE_HGN_MASK = (
LP_HGN_OVERFLOW_LINEAR_BIAS
| LP_HGN_OVERFLOW_RETARDED_BIAS
| LP_HGN_ZERO_TRACKING_FAILED
| LP_HGN_ION_ADMITTANCE_GREATER_THAN_RETARDED_ADMITTANCE
| LP_HGN_ION_CURRENT_GREATER_THAN_RETARDED_CURRENT
),
LP_TE_LGN_MASK = (
LP_LGN_OVERFLOW_LINEAR_BIAS
| LP_LGN_OVERFLOW_RETARDED_BIAS
| LP_LGN_ZERO_TRACKING_FAILED
| LP_LGN_ION_ADMITTANCE_GREATER_THAN_RETARDED_ADMITTANCE
| LP_LGN_ION_CURRENT_GREATER_THAN_RETARDED_CURRENT
),
LP_VS_HGN_MASK = (
LP_HGN_OVERFLOW_LINEAR_BIAS
| LP_HGN_OVERFLOW_RETARDED_BIAS
| LP_HGN_ZERO_TRACKING_FAILED
| LP_HGN_LINEAR_BIAS_LESS_THAN_RETARDED_BIAS
| LP_HGN_LINEAR_BIAS_GREATER_THAN_5V_16BIT_OVERFLOW
| LP_HGN_ION_ADMITTANCE_GREATER_THAN_RETARDED_ADMITTANCE
| LP_HGN_ION_CURRENT_GREATER_THAN_RETARDED_CURRENT
| LP_HGN_RETARDED_ADMITTANCE_GREATER_THAN_LINEAR_ADMITTANCE
| LP_HGN_RETARDED_CURRENT_GREATER_THAN_LINEAR_CURRENT
),
LP_VS_LGN_MASK = (
LP_LGN_OVERFLOW_LINEAR_BIAS
| LP_LGN_OVERFLOW_RETARDED_BIAS
| LP_LGN_ZERO_TRACKING_FAILED
| LP_LGN_LINEAR_BIAS_LESS_THAN_RETARDED_BIAS
| LP_LGN_LINEAR_BIAS_GREATER_THAN_5V_16BIT_OVERFLOW
| LP_LGN_ION_ADMITTANCE_GREATER_THAN_RETARDED_ADMITTANCE
| LP_LGN_ION_CURRENT_GREATER_THAN_RETARDED_CURRENT
| LP_LGN_RETARDED_ADMITTANCE_GREATER_THAN_LINEAR_ADMITTANCE
| LP_LGN_RETARDED_CURRENT_GREATER_THAN_LINEAR_CURRENT
)
};
#endif // _CALCULATE_PRODUCTS_H_