0

我有一个初始化 netcdf 文件的函数,然后是一个写入字段的函数。他们是:

void initNc(string name){
  NcFile myfile(name.c_str(), NcFile::Replace);
  NcDim *xdim = myfile.add_dim("xdim", NX);
  NcDim *tdim = myfile.add_dim("tdim");
  myfile.add_var("x", ncDouble, tdim, xdim);
  myfile.add_var("vx", ncDouble, tdim, xdim);
}

void writeToNcFile(double *x, string name, string varname){
  NcFile myfile(name.c_str(), NcFile::Write);
  NcDim *tdim = myfile.get_dim("tdim");
  NcVar *xvar = myfile.get_var(varname.c_str());
  long sz = tdim->size();
  xvar->set_rec(sz);
  xvar->put_rec(x);
  myfile.close();
}

对于一个变量(在我的 init 函数示例中我有两个),并多次调用 writeToNcFile,这工作正常,在记录维度之后附加结果。但是,一旦我有两个变量(如在本例中,称为 x 和 vx),并且我为每个变量调用一个函数(因此,两次调用),它就会为每个变量 PER 函数调用记录!以下是文件的内容,例如 NX=10:

dimensions:
xdim = 10 ;
tdim = UNLIMITED ; // (2 currently)
variables:
    double x(tdim, xdim) ;
    double vx(tdim, xdim) ;
data:
 x =
  0, 1, 2, 3, 4, 5, 6, 7, 8, 9,
  _, _, _, _, _, _, _, _, _, _ ;
 vx =
  _, _, _, _, _, _, _, _, _, _,
  0, 2, 4, 6, 8, 10, 12, 14, 16, 18 ;}

我该如何解决这个问题?

4

1 回答 1

1

错误来自这样一个事实,即为每个变量调用一次写入函数,在每次调用之间关闭文件。这样每次调用函数时都会为 nc 文件中的所有变量放置一条记录。一种解决方案是在一个函数中写入一条记录的所有字段。

于 2013-07-16T14:28:34.257 回答