Global Index (short | long) | Local contents | Local Index (short | long)
This script calls | |
---|---|
clear lims = [105 285 -6 6]; tim = 1201:1800; cd /home/disk/hayes2/dvimont/csiro/data/Individual_levels nc = netcdf('temp_M_L1_1000_years_new.nc', 'nowrite'); lat = nc{'latitude'}(:); lon = nc{'longitude'}(:); [xk, yk] = keep_var(lims, lon, lat); temp = nc{'temp'}(tim, 1, yk, xk); mv = nc{'temp'}.missing_value(:); nc = close(nc); temp(temp == mv) = NaN; temp = squeeze(temp); temp = squeeze(mean2(shiftdim(temp, 1)))'; lon = lon(xk); [temp, clim] = annave(temp); cd /home/disk/hayes2/dvimont/csiro/data/Individual_levels nc = netcdf('temp_M_L5_1000_years_new.nc', 'nowrite'); temp2 = nc{'temp'}(tim, 1, yk, xk); mv2 = nc{'temp'}.missing_value(:); nc = close(nc); temp2(temp2 == mv2) = NaN; temp2 = squeeze(temp2); temp2 = squeeze(mean2(shiftdim(temp2, 1)))'; [temp2, clim2] = annave(temp2); tim2 = 1201:1440; tim2 = 1201:1261; tim2 = [1201:1440]+240; subplot(1,2,1); contourf(lon, tim2, rave(temp(tim2-1200, :), 5), [-1.5:.15:1.5]); caxis([-1 1]); set(gca, 'YTick', min(tim2):12:max(tim2), ... 'YTickLabel', ((min(tim2)-1)/12):((max(tim2)-1)/12)) shading flat subplot(1,2,2); contourf(lon, tim2, rave(temp2(tim2-1200, :), 5), [-1.5:.15:1.5]); caxis([-1 1]); set(gca, 'YTick', min(tim2):12:max(tim2), ... 'YTickLabel', ((min(tim2)-1)/12):((max(tim2)-1)/12)) shading flat