Effacer les filtres
Effacer les filtres

how can I read multiple netcdf file using ncread in matlab?

50 vues (au cours des 30 derniers jours)
Deep Shah
Deep Shah le 2 Sep 2017
I have 6000 files netcdf files.In each file there are four variables like latitude, longitude,precipitation and time.
I made a variable
filename=dir(location of folder in which files are)
now I want to read data from each 6000 files for every variable I write
for j=3:length(filename)
a=ncread('\folder',filename(j).name,variable);
end
but files are not being read. can u tell me how can I read these files?
the format of file is 3B42_Daily.20020315.7.nc4
  2 commentaires
KSSV
KSSV le 2 Sep 2017
Are Looking the variable names in all the files same?
Deep Shah
Deep Shah le 2 Sep 2017
yes

Connectez-vous pour commenter.

Réponses (2)

Walter Roberson
Walter Roberson le 2 Sep 2017
ncvars = {'latitude', 'longitude', 'precipitation', 'time'};
projectdir = 'C:\some\location';
dinfo = dir( fullfile(projectdir, '*.nc4') );
num_files = length(dinfo);
filenames = fullfile( projectdir, {dinfo.name} );
lats = cell(num_files, 1);
lons = cell{num_files, 1);
precips = cell(num_files, 1);
times = cell(num_files, 1);
for K = 1 : num_files
this_file = filenames{K};
lats{K} = ncread(this_file, ncvars{1});
lons{K} = ncread(this_file, ncvars{2});
precips{K} = ncread(this_file, ncvars{3});
times{K} = ncread(this_file, ncvars{4});
end
  10 commentaires
Tsehaye Negash
Tsehaye Negash le 14 Juin 2023
This program is not working.
Walter Roberson
Walter Roberson le 14 Juin 2023
What error message are you encountering?

Connectez-vous pour commenter.


Dushantha Sandaruwan WIJENDRA NAIDHELAGE
Modifié(e) : Walter Roberson le 14 Juin 2023
This code can be used to substract the specific region from the set of netcdf files
clear all;clc
av_file=dir('*.nc');
ncdisp(av_file(1).name);
lon1=ncread(av_file(1).name,'lon');
lat1=ncread(av_file(1).name,'lat');
% temp1=zeros(141,121,365);
b=[1982:2021];
count=0;
for i=1:40
a=['sst.day.mean.',num2str(b(i)),'.nc'];
sst0=double(ncread(a,'sst'));
a1=find(lon1>=40&lon1<=110);
b1=find(lat1>=-10&lat1<=30);
lon_ind=lon1(a1);lat_ind=lat1(b1);
arb=sst0(a1,b1,:);
[~,~,nt]=size(lhflx1);
sst_final(:,:,count+1:count+nt)=arb(:,:,:);
count=count+nt;
end
  1 commentaire
DGM
DGM le 21 Fév 2023
How is this an answer to the question?

Connectez-vous pour commenter.

Tags

Community Treasure Hunt

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

Start Hunting!

Translated by