Documentation of plot_monthly_ct_time

Global Index (short | long) | Local contents | Local Index (short | long)

Cross-Reference Information

Listing of script plot_monthly_ct_time

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