Extracting Data from Over 600 .nc Files...

2021 4 10
1 个回答
回答已采纳
更新时间：2021 4 13
10 次查看（30 天）

0 个投票

Hi all!
Here's a tough one I'm hoping to receive some guidance on: I'm trying to average a specific variable over six hundred days. That is, I have 600 .nc files from different dates which I have to look into, extract data from, and ultimately work with.
The variable itself is specific humidity (qv), and it is of size 1 x 4, where each column represents: longitude, latitude, pressure level, and time of day (timestep). I have figured out how to read qv for all pressure levels at a specific location and a specific timestep using ONE .nc file, but I would like to automate the process to do this (i.e. read the variable) for all 600 days using all 600 .nc files. Any pointers or suggestions? Thanks in advance!!!
Here's what I have so far:
merra = ('MERRA2_100.inst6_3d_ana_Np.19800115.SUB.nc')
ncdisp(merra)
qv = []
latitude = ncread(merra, 'lat')
longitude = ncread(merra, 'lon')
time = ncread(merra, 'time')
level = ncread(merra, 'lev')
for i = 1:size(level)
i
z = ncread(merra, 'QV', [130, 38, i, 1], [1, 1, 1, 1]) %the 130 is longitude, the 38 is latitude...while 1 is timestep
qv = [qv; z]
end
qv
%the variable 'qv' is an array of size 21 x 1, representing values of
%specific humidity at 21 different pressure levels.

请先登录，再回答此问题。

 采纳的回答

Jan
Jan 2021-4-11

1 个投票

See: FAQ: How to process a sequence of files
Folder = 'C:\Your\Data\Folder';
FileList = dir(fullfile(Folder, '*.nc'));
for iFile = 1:numel(FileList)
merra = fullfile(FileList(iFile).folder, FileList(iFile).name);
... your code comes here
end
Do not reset qv to the empty matrix, but pre-allocate the output and use 2 indices:
Folder = 'C:\Your\Data\Folder';
FileList = dir(fullfile(Folder, '*.nc'));
qv = zeros(21, numel(FileList));
for iFile = 1:numel(FileList)
...
for i = ...
qv(i, iFile) = ncread(merra, 'QV', [130, 38, i, 1], [1, 1, 1, 1]);
end
end
This is tricky:
for i = 1:size(level)
size() replies a vector. Then 1:size(level) uses the 1st element of the size vector and ignores the rest. If this is really wanted, prefer to do this explicitly:
for i = 1:size(level, 1)
If you mean the number of elements, use
for i = 1:numel(level)

3 个评论
显示 1更早的评论 隐藏 1更早的评论

Michelle De Luna
Michelle De Luna 2021-4-11
@Jan: Thank you for taking the time to look into this! Your response was very helpful! It makes sense now, the only issue then is that I receive an error message reading the following: "Unable to perform assignment because the size of the left side is 1x1 and the size of the right side is 1x21." I know this is a follow-up question and understand if you don't have time to respond, but any suggestions would be greatly appreciated. (Please note, I changed the range of my latitude to be [130, 38, i, 1] [1, 21, 1, 1]...) Thanks!
folderData = 'C:\Your\Data\Folder'
fileList = dir(fullfile(folderData, '*.nc'))
qv = zeros(21, numel(fileList))
for iFile = 1:numel(fileList)
merra = fullfile(fileList(iFile).folder, fileList(iFile).name)
ncdisp(merra);
lon = ncread(merra, 'lon');
lat = ncread(merra, 'lat');
time = ncread(merra, 'time');
level = ncread(merra, 'lev')
for i = 1:numel(level)
i
qv(i, iFile) = ncread(merra, 'QV', [130, 38, i, 1], [1, 21, 1, 1])
end
end
Jan
Jan 2021-4-11
Then expand qv:
qv = zeros(21, numel(fileList), 21);
...
qv(i, iFile, :) = ncread(merra, 'QV', [130, 38, i, 1], [1, 21, 1, 1])
Michelle De Luna
Michelle De Luna 2021-4-13
@Jan: Thank you! You've been very helpful. Wishing you a great week!

请先登录，再进行评论。

更多回答（0 个）

类别

在 帮助中心 和 File Exchange 中查找有关 Loops and Conditional Statements 的更多信息

标签

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by