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
|
/***********************************************/
/**
* @file viennaMappingFunctionGrid2File.cpp
*
* @brief Converts VMF grid time series to one file.
*
* @author Torsten Mayer-Guerr
* @date 2010-04-03
*
*/
/***********************************************/
// Latex documentation
#define DOCSTRING docstring
static const char *docstring = R"(
This program converts the gridded time series of the Vienna Mapping Functions (VMF) into
the \file{GROOPS file format}{griddedDataTimeSeries}.
Gridded VMF data is available at: \url{https://vmf.geo.tuwien.ac.at/trop_products/GRID/}
)";
/***********************************************/
#include "programs/program.h"
#include "base/griddedData.h"
#include "inputOutput/file.h"
#include "files/fileGriddedDataTimeSeries.h"
#include "classes/timeSeries/timeSeries.h"
/***** CLASS ***********************************/
/** @brief Converts VMF grid time series to one file.
* @ingroup programsConversionGroup */
class ViennaMappingFunctionGrid2File
{
public:
void run(Config &config, Parallel::CommunicatorPtr comm);
};
GROOPS_REGISTER_PROGRAM(ViennaMappingFunctionGrid2File, SINGLEPROCESS, "converts VMF grid time series to one file.", Conversion)
/***********************************************/
void ViennaMappingFunctionGrid2File::run(Config &config, Parallel::CommunicatorPtr /*comm*/)
{
try
{
FileName fileNameCoefficients;
std::vector<FileName> fileNamesIn;
TimeSeriesPtr timeSeries;
Angle deltaL, deltaB;
Bool isCellRegistered;
readConfig(config, "outputfileVmfCoefficients", fileNameCoefficients, Config::OPTIONAL, "", "");
readConfig(config, "inputfile", fileNamesIn, Config::MUSTSET, "", "files must be given for each point in time");
readConfig(config, "timeSeries", timeSeries, Config::MUSTSET, "", "times of input files");
readConfig(config, "deltaLambda", deltaL, Config::DEFAULT, "1", "[deg] sampling in longitude");
readConfig(config, "deltaPhi", deltaB, Config::DEFAULT, "1", "[deg] sampling in latitude");
readConfig(config, "isCellRegistered", isCellRegistered, Config::DEFAULT, "1", "grid points represent cells (VMF3), not grid corners (VMF1)");
if(isCreateSchema(config)) return;
// ======================================================
const std::vector<Time> times = timeSeries->times();
if(times.size()!=fileNamesIn.size())
throw(Exception("fileCount("+fileNamesIn.size()%"%i) != timeCount("s+times.size()%"%i)"s));
std::vector<Angle> longitudes(static_cast<UInt>(std::round(2*PI/deltaL))); //!< Longitude (columns)
for(UInt i=0; i<longitudes.size(); i++)
longitudes.at(i) = Angle((i+(isCellRegistered ? 0.5 : 0))*deltaL);
std::vector<Angle> latitudes(static_cast<UInt>(std::round(PI/deltaB))+(isCellRegistered ? 0 : 1)); //!< Latitude (rows)
for(UInt i=0; i<latitudes.size(); i++)
latitudes.at(i) = Angle(PI/2-(i+(isCellRegistered ? 0.5 : 0))*deltaB);
GriddedDataRectangular gridRectangular;
gridRectangular.longitudes = longitudes;
gridRectangular.latitudes = latitudes;
gridRectangular.heights.resize(latitudes.size(), 0.);
GriddedData grid(gridRectangular);
// ======================================================
logStatus<<"read coefficients from files"<<Log::endl;
std::vector<Matrix> data(times.size(), Matrix(grid.points.size(), 8));
Single::forEach(fileNamesIn.size(), [&](UInt i)
{
InFile file(fileNamesIn.at(i));
std::string line;
Double lon, lat;
for(UInt k=0; k<grid.points.size(); k++)
{
while(std::getline(file, line) && line.size() && (line.at(0) == '!' || line.at(0) == '#')); // skip header/comment lines
std::stringstream ss(line);
ss>>lat>>lon;
ss>>data.at(i)(k, 0)>>data.at(i)(k, 1)>>data.at(i)(k, 2)>>data.at(i)(k, 3); // ah, aw, zhd, zwd
ss>>data.at(i)(k, 4)>>data.at(i)(k, 5)>>data.at(i)(k, 6)>>data.at(i)(k, 7); // gradients
}
data.at(i).column(4, 4) *= 0.001; // gradients mm --> m
});
// ======================================================
if(!fileNameCoefficients.empty())
{
logStatus<<"write VMF coefficients to <"<<fileNameCoefficients<<">"<<Log::endl;
writeFileGriddedDataTimeSeries(fileNameCoefficients, 1, times, grid, data);
}
}
catch(std::exception &e)
{
GROOPS_RETHROW(e)
}
}
/***********************************************/
/***********************************************/
|
