Merge netCDF files to have sequence of daily precipitation of several weather stations

Please, I posted this question long time ago and I didn't get any help. Please, I appreciate any help and suggestions! The question posted here: https://www.mathworks.com/matlabcentral/answers/343582-merge-netcdf-files-to-have-sequence-of-daily-precipitation-of-several-weather-stations

Thanks, Majid

[KSSV]

You can read the respective nc files....and write them into new nc file....what difficulty you facing?

[Majid Mohamod]

When I run the code which you have written, it give me the following error:

   Error in netcdf.getAtt (line 43)
          xtype = netcdflib('inqAtt',ncid,varid,attname);

Error in internal.matlab.imagesci.nc/readAttribute (line 174)
            attValue    = netcdf.getAtt(gid, varid, attName);

Error in ncreadatt (line 27)
attValue = ncObj.readAttribute(location, attName);

Error in Untitled2 (line 9)
    date{i,1} = ncreadatt(files(i).name,'/','HDF5_GLOBAL.BeginDate')  ;

I'm good if I get the script work as the following: 1- extract data from specific latitude and longitude 2- export them to excel as separate excel file (each nc file export to excel file).

Note, I want to run this script automatically on thousand of nc file and export them automatically too!

Thanks KSSV for all your help

[KSSV]

That case....attach the nc file which shows error....and the code....

[Majid Mohamod]

Sure thing, I herewith attached everything what you want. Please, feel free to ask you for more data!

I appreciate your help!

[KSSV]

Dear friend....you need to have a look into nc files....that code is trying to read date, time, precipitation from the nc file. Your nc files doesn't have that information....

[Majid Mohamod]

I tried to enlarge the study area and now I am getting results " the script works".

You had helped me someday to export the array as excel worksheets as shown below:

files = dir('*.nc') ;

nfiles = length(files) ;
P = cell(nfiles,1) ; % precipitation of all files
date = cell(nfiles,2) ;
time = cell(nfiles,2) ;
for i = 1:nfiles
    %%Read dat
    date{i,1} = ncreadatt(files(i).name,'/','HDF5_GLOBAL.BeginDate')  ;
    date{i,2} = ncreadatt(files(i).name,'/','HDF5_GLOBAL.EndDate')  ;
    %%Read time
    time{i,1} = ncreadatt(files(i).name,'/','HDF5_GLOBAL.BeginTime')  ;
    time{i,2} = ncreadatt(files(i).name,'/','HDF5_GLOBAL.EndTime')  ;
    %%REad precipitation
    P{i} = ncread(files(i).name,'precipitation') ;
end

myfile = 'myfile.xlsx' ;
for i = 1:nfiles
    xlswrite(myfile,P{i},i)
end

Could you change the part of exporting from worksheet to excel file individually. Which mean each nc file will be excel file instead excel sheet?

I appreciate it man!

Majid