读取多个 netcdf 文件 - matlab

Question

我有 40 个名为 'lffd1961_V10M_Iberia.nc'、'lffd1962_V10M_Iberia' ... 'lffd2000_V10M_Iberia' 的 netcdf 文件。

我想在同一个程序中打开所有这些，在所有这些中获取相同的变量，在该变量中应用一些等式并创建 40 个新的 netcdf 文件并修改该变量。

为了打开 40 个文件，我开发了这段代码，它可以工作。

for yr=1961:2000
inFile=strcat('lffd',num2str(yr),'_V10M_Iberia.nc'); disp(inFile)
nc=netcdf.open(inFile,'NC_NOWRITE');
netcdf.close(nc)
end

现在我想在所有变量中获取变量号 4。这个变量是一个 3D 矩阵（70x51x（8760 或 8784）） 像下一行这样的东西，但有 40 年的循环：

ws_1961(1962, etc, etc) = netcdf.getVar(nc,4);

它应该在我的工作区中有 40 个新变量。 在此之后，我想在所有 40 个变量中应用这个方程：

ws_80 = ws.*(log(z/alpha)/log(exp(1))/(log(10/alpha)/log(exp(1))));

最后在每个文件中使用我的新变量创建 40 个新的 netcdf 文件。

类似的东西：

outFile=strcat('year',num2str(yr),'_80m.nc')
ncid = netcdf.create(outFile,'CLOBBER');

dimid0 = netcdf.defDim(ncid,'longitude',nlon); % getvar number 1
dimid1 = netcdf.defDim(ncid,'latitude',nlat); %getvar number 2
dimid2 = netcdf.defDim(ncid,'time',nt); %getvar number 3

varid0 = netcdf.defVar(ncid,'longitude','double', dimid0);
varid1 = netcdf.defVar(ncid,'latitude','double', dimid1);
varid2 = netcdf.defVar(ncid,'time','double', dimid2);
varid3 = netcdf.defVar(ncid,'ws_80','double', [dimid0 dimid1 dimid2]);

netcdf.putAtt(ncid,varid0,'units','degrees_east');
netcdf.putAtt(ncid,varid0,'long_name','longitude');

netcdf.putAtt(ncid,varid1,'units','degrees_north');
netcdf.putAtt(ncid,varid1,'degrees_north','latitude');

netcdf.endDef(ncid);

netcdf.putVar(ncid,varid0,longitude);
netcdf.putVar(ncid,varid1,latitude);
netcdf.putVar(ncid,varid2,time);
netcdf.putVar(ncid,varid3,ws80);

我很难用你们所有人都理解的方式来解释这一点，所以请放心地问任何你想问的问题。

Answer 1

这里有一些“空气代码”可以帮助您入门。

% The variable we want to process
myVarName = 'ws_80';

for yr=1961:2000
  inFile  = strcat('lffd',num2str(yr),'_V10M_Iberia.nc'); 
  disp(inFile);
  outFile = strcat('lffd',num2str(yr),'_V10M_Iberia_processed.nc'); 
  disp(outFile);

  % copy input file to create a template for the output
  copyDone = copyfile(inFile, outFile,'f');
  if(~copyDone)
      error('Could not copy file');
  end

  % Read variable
  inVar = ncread(inFile,myVarName);

  % Replace this line with your equation
  outVar = myProcessFunction(inVar);

  % Write variable
  ncwrite(outFile, myVarName, outVar);

end

您必须对其进行修改以适应您的目标。试一试，然后回复您遇到的问题。