hello Augusto
this code will loop over all specified monthes (it determines automatically how many days are in the current month)
remains now what you want to do will the daily data ...
% Download the CPC data used in the script below
which_year = 2000; % year to download
httpsUrl = "https://www.ncei.noaa.gov/data/cmorph-high-resolution-global-precipitation-estimates/access/daily/0.25deg/";
% filename_tmp="temp.nc"; % temporary filename (for loop); data will be always stored under this temp name
%% main loop
mm = (1:12); % all monthes in year
% mm = [1 6]; % process only some specific monthes in year
for ci = 1:numel(mm) % loop over monthes
days_per_month = datenum(which_year, mm(ci)+1, 1) - datenum(which_year, mm(ci), 1); % compute how many days this month/year
for ck = 1:days_per_month % loop over days
ymd_str = strcat(sprintf('%04d', which_year), sprintf('%02d', mm(ci)), sprintf('%02d', ck)); % year month day string
Url = strcat(httpsUrl, sprintf('%04d', which_year), "/", sprintf('%02d', mm(ci)),"/");
filename = strcat("CMORPH_V1.0_ADJ_0.25deg-DLY_00Z_", ymd_str ,".nc"); % create filename according to list above
disp(strcat(filename, " is processed"));
dataUrl = strcat(Url, filename);
data = webread(dataUrl);
filename_tmp=strcat("data_", ymd_str ,".nc"); % temporary filename (for loop); data will be always stored under this temp name
filename_out = websave(filename_tmp,dataUrl); % save downloaded data
% ncdisp(filename_tmp,'/','min'); % commented, just to avoid filling the command window
precip=ncread(filename_tmp,'cmorph'); % Dimensions: lon,lat,time
long=ncread(filename_tmp,'lon');
lat=ncread(filename_tmp,'lat');
time=ncread(filename_tmp,'time');
end
end
