Documentation of compare_various_grs

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

Help text

fid = fopen('reg_sst.bin', 'r', 'b');

Cross-Reference Information

Listing of script compare_various_grs

get_global
cd /home/disk/tao/data/nmc.reanalysis/monthly
sst = getnc('sst.mon.mean.nc', 'air');
[lat, lon] = getll('sst.mon.mean.nc');
default_global
[sst, clim] = annave(sst);
[ntim, nlat, nlon] = size(sst);
sst = reshape(sst, ntim, nlat*nlon);
cd /home/disk/hayes2/dvimont/data
load gr_ct_pacific_only.mat
grpac = gr;
ttl = ['Smoothed GR time series, calculated from Pacific only (60S to 60N)']
cd /home/disk/hayes2/dvimont/data
load ct_gr.mat
grorig = gr;
gr = rave(rave(gr, 9), 5);
grorig_smooth = (gr - mean(gr)) / std(gr);
cd /home/disk/hayes2/dvimont/data
load gr_ccm.mat
grccm = gr2;
cd /home/disk/tao/dvimont/pub/data
fid = fopen('GR_58.bin', 'r', 'b');
gr_zhang = fread(fid, 'float');
grz = (gr_zhang - mean(gr_zhang)) / std(gr_zhang);

grpat = gr' * sst / ntim;
grpat = reshape(grpat, nlat, nlon);

figure(1)
sd(1)
     gcont(grpat, [-2:.2:2])
     dc
     title(['SST regressed on GR']);
     xlabel(['Contour interval:  0.2 K (std)^-^1'])
sd(2)
     plot(grorig);
     hold on;
     plot(grccm, '--');
     plot(grpac, '-.');
     plot(grlp, '.');
     hold off
     set(gca, 'XTick', [25:60:500], 'XTickLabel', [60:5:100])
     axis([0 481 -2.5 2])
     grid
     xlabel('Year');
     ylabel('STD');
     title(ttl);
     

cd /home/disk/tao/dvimont/matlab/CCM/GR/GR_Plots

subplot(5,1,3)
plot(grpac)
plot(rave(rave(grccm,9),5))
ylabel('gr\_nmc\_ccmtpac')
     set(gca, 'XTick', [1:12:500], 'XTickLabel', [58:100],'YTick',[-5:5])
     axis([0 433 -2.5 2])
     grid

  b = rave(rave(gr1, 9), 5);
  b = grpac;
for i = -15:5;
  a(i+16) = corr(gr_zhang, b(1:395), i);
end

subplot(5,3,13)
bar(-15:5,a);
axis([-15 5 0.5 1])
grid
xlabel('gr\_zhang, grpac')