Documentation of calc_clim

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

Help text

  Check the curl(tau) method

Cross-Reference Information

Listing of script calc_clim

clear
lev = 1:10;
tim = 101:550;
lims = [100 300 8 70];
[v, lat, lon, depth, middepth] = getheat(lims, lev, tim, 'v');
vclim = squeeze(mean(v/100));
[u, lat, lon, depth, middepth] = getheat(lims, lev, tim, 'u');
uclim = squeeze(mean(u/100));
[strf, lat_out, lon_out] = vel_to_strfnc(vclim, lat, lon);
bc = strf;
default_global; XAX = lon_out; YAX = lat_out; FRAME = lims;
figure(1); figure_orient;
subplot(3,1,1);
  gcont(1e-5*strf, [-5:.2:5]);
  dc2(strf)
XAX = lon; YAX = lat; FRAME = lims;
subplot(3,1,2);
  gcont(vclim, [-.1:.01:.1]);
  dc2(vclim);
subplot(3,1,3);
  gcont(uclim, [-.1:.01:.1]);
  dc2(uclim);

[taux, tauy] = getnc('taux', 'tauy', lims, 1, tim);
[lat, lon] = getll('taux', lims);

xclim = squeeze(mean(taux/10));
yclim = squeeze(mean(tauy/10));

H = middepth(max(lev)+1) - middepth(min(lev));
[strf2, lat_out2, lon_out2] = tau_to_strf(xclim, yclim, lat, lon, bc, H);

default_global; XAX = lon_out2; YAX = lat_out2; FRAME = lims;

figure(2); figure_orient;
subplot(3,1,1);
  gcont(1e-5*strf2, [-5:.2:5]);
  dc2(strf2)

XAX = lon; YAX = lat; FRAME = lims;
subplot(3,1,2);
  gcont(vclim, [-.1:.01:.1]);
  dc2(vclim);

subplot(3,1,3);
  gcont(uclim, [-.1:.01:.1]);
  dc2(uclim);