《热带太平洋温度.doc》由会员分享,可在线阅读,更多相关《热带太平洋温度.doc(6页珍藏版)》请在金锄头文库上搜索。
1、ncload t2007.1.cdf;temp(temp-10)=nan; %陆地设为非数for i=1:1:20 %一次绘出、保存20层图像l=temp(i,103:198,197:564); %纬度从23.75S到23.75N,经度从98.75E到77.75Wlev=squeeze(l); lon1=lon(197:564);lat1=lat(103:198);c,h=contourf(lon1,lat1,lev,8,10,15,20,22,24,26,28,29,30,31,32);colormap(jet)caxis(8,32); %所有图的colorbar统一colorbar(hor
2、iz);clabel(c,h,10,20,24,28,29,30,31);xlabel( longitude );ylabel( latitude );str1 = strcat( Temp of Tropical Pacific of Tropical Pacific 2007 Jan lev, num2str(i);title(str1); %用字符串str1作为标题str2 = strcat(f:2007t2007_1_lev, num2str(i), .jpg);saveas(gcf, str2, jpg); %用字符串str2命名并保存图象到该路径end; clf;clear;ncl
3、oad t2007.2.cdf;temp(temp-10)=nan; %陆地设为非数for i=1:1:20 %一次绘出、保存20层图像l=temp(i,103:198,197:564); %纬度从23.75S到23.75N,经度从98.75E到77.75Wlev=squeeze(l); lon1=lon(197:564);lat1=lat(103:198);c,h=contourf(lon1,lat1,lev,8,10,15,20,22,24,26,28,29,30,31,32);colormap(jet)caxis(8,32); %所有图的colorbar统一colorbar(horiz)
4、;clabel(c,h,10,20,24,28,29,30,31);xlabel( Longitude );ylabel( Latitude );str1 = strcat( Temp of Tropical Pacific 2007 Feb lev, num2str(i);title(str1); %用字符串str1作为标题str2 = strcat(f:2007t2007_2_lev, num2str(i), .jpg);saveas(gcf, str2, jpg); %用字符串str2命名并保存图象到该路径end; clf;clear;ncload t2007.3.cdf;temp(te
5、mp-10)=nan; %陆地设为非数for i=1:1:20 %一次绘出、保存20层图像l=temp(i,103:198,197:564); %纬度从23.75S到23.75N,经度从98.75E到77.75Wlev=squeeze(l); lon1=lon(197:564);lat1=lat(103:198);c,h=contourf(lon1,lat1,lev,8,10,15,20,22,24,26,28,29,30,31,32);colormap(jet)caxis(8,32); %所有图的colorbar统一colorbar(horiz);clabel(c,h,10,20,24,28
6、,29,30,31);xlabel( Longitude );ylabel( Latitude );str1 = strcat( Temp of Tropical Pacific 2007 Mar lev, num2str(i);title(str1); %用字符串str1作为标题str2 = strcat(f:2007t2007_3_lev, num2str(i), .jpg);saveas(gcf, str2, jpg); %用字符串str2命名并保存图象到该路径end; clf;clear;ncload t2007.4.cdf;temp(temp-10)=nan; %陆地设为非数for
7、i=1:1:20 %一次绘出、保存20层图像l=temp(i,103:198,197:564); %纬度从23.75S到23.75N,经度从98.75E到77.75Wlev=squeeze(l); lon1=lon(197:564);lat1=lat(103:198);c,h=contourf(lon1,lat1,lev,8,10,15,20,22,24,26,28,29,30,31,32);colormap(jet)caxis(8,32); %所有图的colorbar统一colorbar(horiz);clabel(c,h,10,20,24,28,29,30,31);xlabel( Long
8、itude );ylabel( Latitude );str1 = strcat( Temp of Tropical Pacific 2007 Apr lev, num2str(i);title(str1); %用字符串str1作为标题str2 = strcat(f:2007t2007_4_lev, num2str(i), .jpg);saveas(gcf, str2, jpg); %用字符串str2命名并保存图象到该路径end; clf;clear;ncload t2007.5.cdf;temp(temp-10)=nan; %陆地设为非数for i=1:1:20 %一次绘出、保存20层图像l
9、=temp(i,103:198,197:564); %纬度从23.75S到23.75N,经度从98.75E到77.75Wlev=squeeze(l); lon1=lon(197:564);lat1=lat(103:198);c,h=contourf(lon1,lat1,lev,8,10,15,20,22,24,26,28,29,30,31,32);colormap(jet)caxis(8,32); %所有图的colorbar统一colorbar(horiz);clabel(c,h,10,20,24,28,29,30,31);xlabel( Longitude );ylabel( Latitud
10、e );str1 = strcat( Temp of Tropical Pacific 2007 May lev, num2str(i);title(str1); %用字符串str1作为标题str2 = strcat(f:2007t2007_5_lev, num2str(i), .jpg);saveas(gcf, str2, jpg); %用字符串str2命名并保存图象到该路径end; clf;clear;ncload t2007.6.cdf;temp(temp-10)=nan; %陆地设为非数for i=1:1:20 %一次绘出、保存20层图像l=temp(i,103:198,197:564
11、); %纬度从23.75S到23.75N,经度从98.75E到77.75Wlev=squeeze(l); lon1=lon(197:564);lat1=lat(103:198);c,h=contourf(lon1,lat1,lev,8,10,15,20,22,24,26,28,29,30,31,32);colormap(jet)caxis(8,32); %所有图的colorbar统一colorbar(horiz);clabel(c,h,10,20,24,28,29,30,31);xlabel( Longitude );ylabel( Latitude );str1 = strcat( Temp
12、 of Tropical Pacific 2007 Jun lev, num2str(i);title(str1); %用字符串str1作为标题str2 = strcat(f:2007t2007_6_lev, num2str(i), .jpg);saveas(gcf, str2, jpg); %用字符串str2命名并保存图象到该路径end; clf;clear;ncload t2007.7.cdf;temp(temp-10)=nan; %陆地设为非数for i=1:1:20 %一次绘出、保存20层图像l=temp(i,103:198,197:564); %纬度从23.75S到23.75N,经度
13、从98.75E到77.75Wlev=squeeze(l); lon1=lon(197:564);lat1=lat(103:198);c,h=contourf(lon1,lat1,lev,8,10,15,20,22,24,26,28,29,30,31,32);colormap(jet)caxis(8,32); %所有图的colorbar统一colorbar(horiz);clabel(c,h,10,20,24,28,29,30,31);xlabel( Longitude );ylabel( Latitude );str1 = strcat( Temp of Tropical Pacific 20
14、07 Jul lev, num2str(i);title(str1); %用字符串str1作为标题str2 = strcat(f:2007t2007_7_lev, num2str(i), .jpg);saveas(gcf, str2, jpg); %用字符串str2命名并保存图象到该路径end; clf;clear;ncload t2007.8.cdf;temp(temp-10)=nan; %陆地设为非数for i=1:1:20 %一次绘出、保存20层图像l=temp(i,103:198,197:564); %纬度从23.75S到23.75N,经度从98.75E到77.75Wlev=squeeze(l); lon1=lon(197:564);lat1=lat(103:198);c,h=contourf(lon1,lat1,lev,8,10,15,20,22,24,26,28,29,30,31,32);colormap(jet)caxis(8,32); %所有图的colorbar统一colo