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
|
/***********************************************/
/**
* @file graceL1a2Temperature.cpp
*
* @brief Read GRACE L1A data.
*
* @author Beate Klinger
* @date 2017-03-09
*
*/
/***********************************************/
// Latex documentation
#define DOCSTRING docstring
static const char *docstring = R"(
This program converts Level-1A temperature measurments (HRT1B or HRT1A) to the GROOPS instrument file format.
The GRACE Level-1A format is described in GRACE given at \url{http://podaac-tools.jpl.nasa.gov/drive/files/allData/grace/sw/GraceReadSW_L1_2010-03-31.tar.gz}.
Multiple \config{inputfile}s must be given in the correct time order.
The output is one arc of satellite data which can include data gaps.
To split the arc in multiple gap free arcs use \program{InstrumentSynchronize}.
)";
/***********************************************/
#include "programs/program.h"
#include "files/fileInstrument.h"
#include "fileGrace.h"
/***** CLASS ***********************************/
/** @brief Read GRACE L1A data.
* @ingroup programsConversionGroup */
class GraceL1a2Temperature
{
void readFile(const FileName &fileName, MiscValuesArc &arcTemp);
public:
void run(Config &config, Parallel::CommunicatorPtr comm);
};
GROOPS_REGISTER_PROGRAM(GraceL1a2Temperature, SINGLEPROCESS, "read GRACE L1A data (HRT1B or HRT1A)", Conversion, Grace, Instrument)
/***********************************************/
void GraceL1a2Temperature::run(Config &config, Parallel::CommunicatorPtr /*comm*/)
{
try
{
FileName fileNameOut;
std::vector<FileName> fileNameIn;
readConfig(config, "outputfileTemperature", fileNameOut, Config::MUSTSET, "", "MISCVALUES");
readConfig(config, "inputfile", fileNameIn, Config::MUSTSET, "", "HRT1B or HRT1A");
if(isCreateSchema(config)) return;
// =============================================
logStatus<<"read input files"<<Log::endl;
Arc arc;
for(UInt idFile=0; idFile<fileNameIn.size(); idFile++)
{
logStatus<<"read file '"<<fileNameIn.at(idFile)<<"'"<<Log::endl;
UInt numberOfRecords;
FileInGrace file(fileNameIn.at(idFile), numberOfRecords);
for(UInt idEpoch=0; idEpoch<numberOfRecords; idEpoch++)
{
Int32 seconds, microSeconds; // seconds, microseconds part
Char timeRef, GRACE_id; // time reference frame (R = Receiver Time, G = GPS time), GRACE satellite ID
Float TempMEPNegY; // I/F support structure to MEP -y [°C]
Float TempMEPPosY; // I/F support structure to MEP +y [°C]
Float TempMEPm; // I/F support structure to MEP mid [°C]
Float TempICU; // ICU Temperature Reference Point [°C]
Float TempICURed; // ICU Tempeartuer Reference Point (redundant) [°C]
Float TempACCNegZ; // ACC thermal cage -z [°C]
Float TempACCPosZ; // ACC thermal cage +z [°C]
Float TempCFRPPosX, TempCFRPPosXRed; // CFRP Frame at +x I/F to Baseplate (+redundant) [°C]
Float TempCFRPNegX, TempCFRPNegXRed; // CFRP Frame at -x I/F to Baseplate (+redundant) [°C]
Float TempCFRPNegY, TempCFRPNegYRed; // CFRP Frame at -y I/F to Baseplate (+redundant) [°C]
Float TempACCSen; // Harness to ACC sensor [°C]
Float TempICUSpec; // ICU special [°C]
Float TempMWANegY, TempMWANegYOff; // MWA -y Temp. Ref. Point, 2 out 3 (+nominally off) [°C]
Float TempMWAPosY, TempMWAPosYOff; // MWA +y Temp. Ref. Point, 2 out 3 (+nominally off) [°C]
Float TempHornPosX, TempHornPosXRed; // Horn aperture +x (+redundant) [°C]
Float TempHornPl, TempHornPlRed; // Horn / platform I/F (+redundant) [°C]
Float TempHWMANegY; // Harnass to MWA electronics -y [°C]
Float TempHWMAPosY; // Harnass to MWA electronics +y [°C]
Float TempRFSamp; // RF-Sampling unit [°C]
Float TempUSONegY, TempUSONegYRed; // USO Temp. Ref. Point -y (+redundant) [°C]
Float TempUSOPosY, TempUSOPosYRed; // USO Temp. Ref. Point +y (+redundant) [°C]
Byte qualityFlag; // data quality flag
try
{
file>>seconds>>microSeconds>>timeRef>>GRACE_id>>TempMEPNegY>>TempMEPPosY>>TempMEPm>>TempICU>>TempICURed>>TempACCNegZ>>TempACCPosZ;
file>>TempCFRPPosX>>TempCFRPPosXRed>>TempCFRPNegX>>TempCFRPNegXRed>>TempCFRPNegY>>TempCFRPNegYRed>>TempACCSen>>TempICUSpec;
file>>TempMWANegY>>TempMWANegYOff>>TempMWAPosY>>TempMWAPosYOff>>TempHornPosX>>TempHornPosXRed>>TempHornPl>>TempHornPlRed;
file>>TempHWMANegY>>TempHWMAPosY>>TempRFSamp>>TempUSONegY>>TempUSONegYRed>>TempUSOPosY>>TempUSOPosYRed>>FileInGrace::flag(qualityFlag);
}
catch(std::exception &/*e*/)
{
// GRACE-FO number of records issue
logWarning<<arc.back().time.dateTimeStr()<<": file ended at "<<idEpoch<<" of "<<numberOfRecords<<" expected records"<<Log::endl;
break;
}
{
MiscValuesEpoch epoch(21);
epoch.time = mjd2time(51544.5) + seconds2time(seconds) + seconds2time(microSeconds*1e-6);
epoch.values(0) = TempACCNegZ;
epoch.values(1) = TempACCPosZ;
epoch.values(2) = TempACCSen;
epoch.values(3) = TempCFRPPosX;
epoch.values(4) = TempCFRPNegX;
epoch.values(5) = TempCFRPNegY;
epoch.values(6) = TempCFRPNegYRed;
epoch.values(7) = TempMEPNegY;
epoch.values(8) = TempMEPPosY;
epoch.values(9) = TempMEPm;
epoch.values(10) = TempICUSpec;
epoch.values(11) = TempICU;
epoch.values(12) = TempMWANegY;
epoch.values(13) = TempMWAPosY;
epoch.values(14) = TempHWMANegY;
epoch.values(15) = TempHWMAPosY;
epoch.values(16) = TempRFSamp;
epoch.values(17) = TempUSONegY;
epoch.values(18) = TempUSOPosY;
epoch.values(19) = TempHornPosX;
epoch.values(20) = TempHornPl;
arc.push_back(epoch);
}
}// for(idEpoch)
} // for(idFile)
// =============================================
Arc::printStatistics(arc);
if(arc.size() == 0)
return;
if(!fileNameOut.empty())
{
logInfo<<"write temperature data to <"<<fileNameOut<<">"<<Log::endl;
InstrumentFile::write(fileNameOut, arc);
}
}
catch(std::exception &e)
{
GROOPS_RETHROW(e)
}
}
/***********************************************/
|
