Documentation of TAU_curl

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

Help text

 print -dps2 TAU_prec.ps

Cross-Reference Information

Listing of script TAU_curl

clear
tim = 101:1000;
lims = [-0.1 360 -90 90];
[taux, tauy] = getnc('taux', 'tauy', lims, 1, tim);
[lat, lon] = getll('taux', lims);
taux = squeeze(mean(taux));
tauy = squeeze(mean(tauy));
[lat1, lon1, depth, lm] = getll('temp', lims);
[curlt, lat2, lon2] = sph_curl1(taux, tauy, lat, lon, 1);
default_global; FRAME = [0 360 -90 90];
figure(1); fo(1); clf;
spthes(1);
  gquiv(taux, tauy, -2, 1, 'dynes cm^-^2');
  XAX = lon1; YAX = lat1;
  dc2(lm, .55, 1);
  title('Wind Stress');
prec = getflx('rnd', lims, tim);
prec = squeeze(mean(prec));
spthes(2);
  XAX = lon1; YAX = lat1;
  gcont(prec, 1.5);
  dc2(lm, .55, -1);
  title('Precipitation');
  xlabel('Contour Interval:  1.5 mm day^-^1');
cd ~/Thesis/Clim

clear

tim = 101:1000;
lims = [-0.1 360 -90 90];

psl = getflx('psl', lims, tim);
z250 = getnc('z250', lims, 1, tim);

psl = squeeze(mean(psl));
z250 = squeeze(mean(z250));

[lat, lon] = getll('z250', lims);
[lat2, lon2, depth, lm] = getll('temp', lims);

figure(2); fo(1); clf;
spthes(1);
  XAX = lon2; YAX = lat2;
  gcont(psl, 4);
  dc2(lm, .55, -1);
  title('SLP');
  xlabel('Contour Interval:  4 hPa');

tem = psl - (mean(psl, 2)*ones(1, 64));
spthes(2);
  gcont(tem, 2);
  dc2(lm, .55, -1);
  title('SLP - Zonal mean of SLP');
  xlabel('Contour Interval:  2 hPa (the 0 contour is omitted)');

cd ~/Thesis/Clim
% print -dps2 SLP.ps


default_global; FRAME = [0 360 -90 90];

figure(1); clf;

sptalk(3,2,1); FRAME = [0 360 0 90];
  mcont3(1e-3*z250, .15, 'stereo', [90 270]);
  XAX = lon2; YAX = lat2;
  drawmap2(lm, 0.55, -1);
  XAX = lon; YAX = lat;
  axis_limits; gridm on;
sptalk(3,2,2); FRAME = [0 360 -90 0];
  mcont3(1e-3*z250, .15, 'stereo', [-90 270]);
  XAX = lon2; YAX = lat2;
  drawmap2(lm, 0.55, -1);
  XAX = lon; YAX = lat;
  axis_limits; gridm on;

tem = z250 - mean(z250, 2)*ones(1, length(lon));
sptalk(3,2,3); FRAME = [0 360 0 90];
  mcont3(tem, 30, 'stereo', [90 270]);
  XAX = lon2; YAX = lat2;
  drawmap2(lm, 0.55, -1);
  XAX = lon; YAX = lat;
  axis_limits; gridm on;
sptalk(3,2,4); FRAME = [0 360 -90 0];
  mcont3(tem, 30, 'stereo', [-90 270]);
  XAX = lon2; YAX = lat2;
  drawmap2(lm, 0.55, -1);
  XAX = lon; YAX = lat;
  axis_limits; gridm on;

sptalk(3,2,1);
  title('NH:  Z250')
  xlabel('Contour Interval:  0.15 km');
sptalk(3,2,2);
  title('SH:  Z250')
  xlabel('Contour Interval:  0.15 km');
sptalk(3,2,3);
  title('NH:  Z250 - zonal mean Z250')
  xlabel('Contour Interval:  30 m');
sptalk(3,2,4);
  title('SH:  Z250 - zonal mean Z250')
  xlabel('Contour Interval:  30 m');


cd ~/Thesis/Clim
% print -dps2 Z250.ps