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

.ncnetcdfswat

% 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

Best Answer

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
Related Question