fpath = '/Users/jackie/Desktop/data/access1.3/';
fname = 'CSIRO_access1.3_cli_2000_01.nc';
myFolder = fpath;
if ~isfolder(myFolder)
errorMessage = sprintf('Error: The following folder does not exist:\n%s', myFolder);
uiwait(warndlg(errorMessage));
return;
end
filePattern = fullfile(myFolder,'*.nc');
theFiles = dir(filePattern);
avgClia = zeros(144,91,40);
for i =1:60
baseFileName = theFiles(i).name;
fullFileName = fullfile(myFolder, baseFileName);
ncid=netcdf.open(fullFileName, 'NC_NOWRITE');
vid=netcdf.inqVarID(ncid,'cli');
clia=netcdf.getVar(ncid,vid);
avgClia = avgClia + clia;
netcdf.close(ncid);
end
avgClia = avgClia / 60;
ncid =netcdf.open([fpath,fname],'NC_WRITE');
vid=netcdf.inqVarID(ncid,'lon');
longitude=netcdf.getVar(ncid,vid);
vid=netcdf.inqVarID(ncid,'lat');
latitude=netcdf.getVar(ncid,vid);
vid=netcdf.inqVarID(ncid,'plev');
z1=netcdf.getVar(ncid,vid);
index = find(longitude<0);
longitude(index) = longitude(index)+360;
index1 = find(longitude<0);
longitude(index1) = longitude(index1)+360;
for d = 1:1
lon_b1 = 0;
lon_b2 = 360;
lat_b1 = -90;
lat_b2 = 90;
end
figure(1);
a = find (longitude>=lon_b1 & longitude<=lon_b2);
b = find (latitude>=lat_b1 & latitude<=lat_b2);
z = find (z1>= 897 & z1 <= 912); plev = '900';
lon = longitude(a);
lat1 = latitude(b);
clia = avgClia(a,b,z,d) * 1000000;
[XX,YY]=meshgrid(lon,lat1);
min_axis=floor(min(min(clia))/0.1)*0.1;
max_axis=ceil(max(max(clia))/0.1)*0.1;
contourf(XX,YY,clia',1000,'linestyle','none');
caxis([min_axis,max_axis]);
ScaleTick_step=0.1;
ScaleTick=min_axis:ScaleTick_step:max_axis;
FactorName=cell(1,length(ScaleTick));
for i=1:length(ScaleTick)
FactorName{i}=num2str(ScaleTick(i));
end
hold on
load coast
plot(long+180,lat,'w')
h = colorbar;
set(h, 'ylim', [0 50])
colormap(jet)
