Extract data from NetCDF files to text files (For SWAT Model).

Question

shobhit pipil 2018-12-23

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/437025-extract-data-from-netcdf-files-to-text-files-for-swat-model

编辑： shobhit pipil 2021-3-16

% new code to extract data from netcdf file to SWAT format text file.
data =[];
Files = dir('*.nc');
for k=1:length(Files)
    FilesNames = Files(k).name;
    for t=1:2
        pr = ncread(FilesNames,'rainfall_amount');
        long =  ncread(FilesNames,'x');
        lat = ncread(FilesNames,'y');
    end
    for j = 716:746(lat);
        for i = 328:368(long);
            
            v = pr(i,j,:);
            data2 = (v);
            data2(isnan(data2)) = -99.0;
            
            outfile =sprintf('%d_%d_PCP.txt',lat(j),long(i));
            fopen all
            fid = fopen(outfile,'wt');
            data = {data;data2};
            fprintf(fid,'%5.1f\n',data2);
        end
        fclose all
    end
end
disp([outfile 'created'])
% original code source: 
% http://dficklin.weebly.com/netcdf-to-swat-climate-input-files.html

4 个评论
显示 2更早的评论隐藏 2更早的评论

Image Analyst 2018-12-23

Please read the link I gave you. I know you haven't yet read it because you have not attached the two data files. So please read it.

shobhit pipil 2018-12-23

编辑：shobhit pipil 2018-12-23

problem_new2312.png

outfiles shows 365 days records instead of 730 days, because single netcsd contains 1 year records per lat long.

NetCDF data download link: https://my.pcloud.com/publink/show?code=XZ60yC7ZmQaTSzygQ6FO1Gvlex2yTXYH3piX

one can copy the same NetCDF file in same directory and change year from 1890 to 91 for testing the code for two files. Single file size is 4.5 GB so I have uploaded just one.

Thanks

请先登录，再进行评论。

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

Answer 1

shobhit pipil 2018-12-24

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/437025-extract-data-from-netcdf-files-to-text-files-for-swat-model#answer_353754

编辑：madhan ravi 2018-12-24

在 MATLAB Online 中打开

Now It works for me:-

% new code to extract data from netcdf file to SWAT format text file.
data =[];
Files = dir('*.nc');
for k=1:length(Files)
    FilesNames = Files(k).name;
    for t=1:2
        pr = ncread(FilesNames,'rainfall_amount');
        long =  ncread(FilesNames,'x');
        lat = ncread(FilesNames,'y');
    end
        
    for j = 716:746(lat);
        for i = 328:368(long);
            outfile =sprintf('%d_%d_PCP.txt',lat(j),long(i));
            fopen all;
            
            v = pr(i,j,:);
            data2 = (v);
            data2(isnan(data2)) = -99.0;
            
     %       outfile =sprintf('%d_%d_PCP.txt',lat(j),long(i));
     %       fopen all;
            fid = fopen(outfile,'At');
            data = {data;data2};
            fprintf(fid,'%5.1f\n',data2);
        end
        fclose all;
    end
end
disp([outfile 'created'])
% original code source: 
% http://dficklin.weebly.com/netcdf-to-swat-climate-input-files.html