tmp update
This commit is contained in:
parent
541ef075f2
commit
367277189b
@ -310,21 +310,22 @@ namespace gctl
|
||||
* @param[in] xname 文件中 x 坐标轴的变量名，默认为 x
|
||||
* @param[in] yname 文件中 y 坐标轴的变量名，默认为 y
|
||||
* @param[in] dataname 文件中网格数据的变量名，默认为 z
|
||||
* @param[in] cor 文件中网格数据排序方式（默认是从左下到右上）
|
||||
*/
|
||||
template <typename T, typename P>
|
||||
void save_netcdf_grid(std::string filename, const array<T> &out_data, int xnum, int ynum,
|
||||
P xmin, P dx, P ymin, P dy, std::string xname = "x", std::string yname = "y",
|
||||
std::string dataname = "z", T zmin = NAN, T zmax = NAN)
|
||||
P xmin, P dx, P ymin, P dy, matrix_corner_e cor = BtmLeft, std::string xname = "x",
|
||||
std::string yname = "y", std::string dataname = "z", T zmin = NAN, T zmax = NAN)
|
||||
{
|
||||
// 准备写出数据
|
||||
if (out_data.empty())
|
||||
if (out_data.empty() || xnum*ynum != out_data.size())
|
||||
{
|
||||
throw runtime_error("The output data is empty. From save_netcdf_grid(...)");
|
||||
throw std::runtime_error("[gctl::save_netcdf_grid] Invalid data dimensions.");
|
||||
}
|
||||
|
||||
if (filename == "")
|
||||
{
|
||||
throw domain_error("The input filename is empty. From save_netcdf_grid(...)");
|
||||
throw std::domain_error("[gctl::save_netcdf_grid] Empty filename empty.");
|
||||
}
|
||||
|
||||
// 判断文件名结尾是否为.nc如果不是则添加.nc结尾
|
||||
@ -338,7 +339,7 @@ namespace gctl
|
||||
NcFile nc_outfile(full_name.c_str(), NcFile::Replace);
|
||||
if (!nc_outfile.is_valid())
|
||||
{
|
||||
throw runtime_error("Can't not open file: " + full_name + ". From save_netcdf_grid(...)");
|
||||
throw runtime_error("[gctl::save_netcdf_grid] Can't not open file: " + full_name + ".");
|
||||
}
|
||||
|
||||
if (std::isnan(zmin))
|
||||
@ -418,7 +419,21 @@ namespace gctl
|
||||
dVar->add_att("valid_range", 2, &data_range[0]);
|
||||
dVar->add_att("missing_value", GCTL_BDL_MAX);
|
||||
|
||||
array<T> tmp_data(out_data);
|
||||
array<T> tmp_data;
|
||||
if (cor == BtmLeft) tmp_data = out_data;
|
||||
else if (cor == TopLeft)
|
||||
{
|
||||
tmp_data.resize(xnum*ynum);
|
||||
for (int i = 0; i < ynum; i++)
|
||||
{
|
||||
for (int j = 0; j < xnum; j++)
|
||||
{
|
||||
tmp_data[j + i*xnum] = out_data[j + (ynum - 1 - i)*xnum];
|
||||
}
|
||||
}
|
||||
}
|
||||
else throw std::runtime_error("[gctl::save_netcdf_grid] Invalid data ordering type.");
|
||||
|
||||
for (int i = 0; i < tmp_data.size(); i++)
|
||||
{
|
||||
if (tmp_data[i] < zmin || tmp_data[i] > zmax)
|
||||
@ -458,7 +473,7 @@ namespace gctl
|
||||
out_data.reform(tmpdata);
|
||||
|
||||
save_netcdf_grid(filename, tmpdata, out_data.col_size(), out_data.row_size(),
|
||||
xmin, dx, ymin, dy, xname, yname, dataname, zmin, zmax);
|
||||
xmin, dx, ymin, dy, TopLeft, xname, yname, dataname, zmin, zmax);
|
||||
return;
|
||||
}
|
||||
|
||||
@ -474,21 +489,21 @@ namespace gctl
|
||||
* @param[in] xname 文件中 x 坐标轴的变量名，默认为 x
|
||||
* @param[in] yname 文件中 y 坐标轴的变量名，默认为 y
|
||||
* @param[in] dataname 文件中网格数据的变量名，默认为 z
|
||||
* @param[in] cor 文件中网格数据排序方式（默认是从左下到右上）
|
||||
*/
|
||||
template <typename T, typename P>
|
||||
void save_netcdf_grid(std::string filename, const array<T> &out_data, const array<P> &out_x,
|
||||
const array<P> &out_y, std::string xname = "x", std::string yname = "y", std::string dataname = "z",
|
||||
T zmin = NAN, T zmax = NAN)
|
||||
const array<P> &out_y, matrix_corner_e cor = BtmLeft, std::string xname = "x",
|
||||
std::string yname = "y", std::string dataname = "z", T zmin = NAN, T zmax = NAN)
|
||||
{
|
||||
// 准备写出数据
|
||||
if (out_data.empty())
|
||||
if (out_data.empty() || out_x.size()*out_y.size() != out_data.size())
|
||||
{
|
||||
throw runtime_error("The output data is empty. From save_netcdf_grid(...)");
|
||||
throw std::runtime_error("[gctl::save_netcdf_grid] Invalid data dimensions.");
|
||||
}
|
||||
|
||||
if (filename == "")
|
||||
{
|
||||
throw domain_error("The input filename is empty. From save_netcdf_grid(...)");
|
||||
throw std::domain_error("[gctl::save_netcdf_grid] Empty filename empty.");
|
||||
}
|
||||
|
||||
// 判断文件名结尾是否为.nc如果不是则添加.nc结尾
|
||||
@ -502,7 +517,7 @@ namespace gctl
|
||||
NcFile nc_outfile(full_name.c_str(), NcFile::Replace);
|
||||
if (!nc_outfile.is_valid())
|
||||
{
|
||||
throw runtime_error("Can't not open file: " + full_name + ". From save_netcdf_grid(...)");
|
||||
throw runtime_error("[gctl::save_netcdf_grid] Can't not open file: " + full_name + ".");
|
||||
}
|
||||
|
||||
if (std::isnan(zmin))
|
||||
@ -571,7 +586,21 @@ namespace gctl
|
||||
dVar->add_att("valid_range", 2, &data_range[0]);
|
||||
dVar->add_att("missing_value", GCTL_BDL_MAX);
|
||||
|
||||
array<T> tmp_data(out_data.get(), out_data.size());
|
||||
array<T> tmp_data;
|
||||
if (cor == BtmLeft) tmp_data = out_data;
|
||||
else if (cor == TopLeft)
|
||||
{
|
||||
tmp_data.resize(xnum*ynum);
|
||||
for (int i = 0; i < ynum; i++)
|
||||
{
|
||||
for (int j = 0; j < xnum; j++)
|
||||
{
|
||||
tmp_data[j + i*xnum] = out_data[j + (ynum - 1 - i)*xnum];
|
||||
}
|
||||
}
|
||||
}
|
||||
else throw std::runtime_error("[gctl::save_netcdf_grid] Invalid data ordering type.");
|
||||
|
||||
for (int i = 0; i < tmp_data.size(); i++)
|
||||
{
|
||||
if (tmp_data[i] < zmin || tmp_data[i] > zmax)
|
||||
@ -608,7 +637,7 @@ namespace gctl
|
||||
gctl::array<T> tmpdata;
|
||||
out_data.reform(tmpdata);
|
||||
|
||||
save_netcdf_grid(filename, tmpdata, out_x, out_y, xname, yname, dataname, zmin, zmax);
|
||||
save_netcdf_grid(filename, tmpdata, out_x, out_y, TopLeft, xname, yname, dataname, zmin, zmax);
|
||||
return;
|
||||
}
|
||||
|
||||
@ -624,10 +653,12 @@ namespace gctl
|
||||
* @param[in] dataname 保存的数据名称，注意不要与已有数据重名
|
||||
* @param[in] zmin 网格数据的最小值，默认为 NAN 。 此时函数会自动计算数组的最小值
|
||||
* @param[in] zmax 网格数据的最大值，默认为 NAN 。 此时函数会自动计算数组的最大值
|
||||
* @param[in] cor 文件中网格数据排序方式（默认是从左下到右上）
|
||||
*/
|
||||
template <typename T>
|
||||
void append_netcdf_grid(std::string filename, const array<T> &out_data, std::string xname,
|
||||
std::string yname, std::string dataname, T zmin = NAN, T zmax = NAN)
|
||||
std::string yname, std::string dataname, matrix_corner_e cor = BtmLeft,
|
||||
T zmin = NAN, T zmax = NAN)
|
||||
{
|
||||
// 准备写出数据
|
||||
if (out_data.empty())
|
||||
@ -730,14 +761,27 @@ namespace gctl
|
||||
dVar->add_att("valid_range", 2, &data_range[0]);
|
||||
dVar->add_att("missing_value", GCTL_BDL_MAX);
|
||||
|
||||
array<T> tmp_data(out_data.get(), out_data.size());
|
||||
array<T> tmp_data;
|
||||
if (cor == BtmLeft) tmp_data = out_data;
|
||||
else if (cor == TopLeft)
|
||||
{
|
||||
tmp_data.resize(xnum*ynum);
|
||||
for (int i = 0; i < ynum; i++)
|
||||
{
|
||||
for (int j = 0; j < xnum; j++)
|
||||
{
|
||||
tmp_data[j + i*xnum] = out_data[j + (ynum - 1 - i)*xnum];
|
||||
}
|
||||
}
|
||||
}
|
||||
else throw std::runtime_error("[gctl::append_netcdf_grid] Invalid data ordering type.");
|
||||
|
||||
for (int i = 0; i < tmp_data.size(); i++)
|
||||
{
|
||||
if (tmp_data[i] < zmin || tmp_data[i] > zmax)
|
||||
tmp_data[i] = GCTL_BDL_MAX;
|
||||
}
|
||||
dVar->put(tmp_data.get(), ynum, xnum);
|
||||
|
||||
return;
|
||||
}
|
||||
|
||||
|
||||
Loading…
Reference in New Issue
Block a user