Documentation of velocity_plots

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

Help text

%%%%%%%%%%%

Cross-Reference Information

Listing of script velocity_plots

clear
cd /home/disk/hayes2/dvimont/csiro/data
tim = 101:450; lev = 5;
nc = netcdf('ul_L1-10.nc', 'nowrite');
  lat = nc{'latitude'}(:);
  lon = nc{'longitude'}(:);
  u = nc{'u'}(tim, lev, :, :);
  mv = nc{'u'}.missing_value(:);
nc = close(nc);
nc = netcdf('vl_L1-10.nc', 'nowrite');
  v = nc{'v'}(tim, lev, :, :);
nc = close(nc);
nc = netcdf('temp_A_L1-10.nc', 'nowrite');
  lat1 = nc{'latitude'}(:);
  lon1 = nc{'longitude'}(:);
  temp = nc{'temp'}(tim, lev, :, :);
nc = close(nc);
u(u == mv) = NaN;
v(v == mv) = NaN;
ubar = squeeze(mean(u));
vbar = squeeze(mean(v));
tbar = squeeze(mean(temp));
default_global;
figure(1); figure_orient;
sp(2)
     [h, a, c] = gquiv(ubar, vbar, 1.5);

clear
cd /home/disk/hayes2/dvimont/csiro/data

tim = 101:550; lev = 1:10;
nc = netcdf('temp_A_L1-10.nc', 'nowrite');
  lat1 = nc{'latitude'}(:);
  lon1 = nc{'longitude'}(:);
  depth = nc{'depth'}(:);
  temp = nc{'temp'}(tim, lev, :, :);
  mv = nc{'temp'}.missing_value(:);
nc = close(nc);
temp = squeeze(mean(temp));
temp(temp == mv) = NaN;

nc = netcdf('ul_L1-10.nc', 'nowrite');
  lat = nc{'latitude'}(:);
  lon = nc{'longitude'}(:);
  u = nc{'u'}(tim, lev, :, :);
  mv = nc{'u'}.missing_value(:);
nc = close(nc);
u = squeeze(mean(u));
u(u == mv) = NaN;

nc = netcdf('vl_L1-10.nc', 'nowrite');
  v = nc{'v'}(tim, lev, :, :);
nc = close(nc);
v = squeeze(mean(v));
v(v == mv) = NaN;

[nlev, nlat, nlon] = size(temp);
for i = 1:nlev;
  temp(i,:,:) = interp2(lon1, lat1, squeeze(temp(i,:,:)), lon, lat');
end

temp(isnan(u)) = NaN;
u(isnan(temp)) = NaN;
v(isnan(temp)) = NaN;

temp = reshape(temp, nlev, nlat*nlon);
u = reshape(u, nlev, nlat*nlon);
v = reshape(v, nlev, nlat*nlon);

u20 = NaN * ones(1, nlat*nlon);
v20 = NaN * ones(1, nlat*nlon);
d20 = NaN * ones(1, nlat*nlon);
for i = 1:nlat*nlon;
  if isempty(find(diff(temp(:,i)) > 0)) & isempty(find(isnan(temp(:,i))));
    u20(i) = interp1(temp(:,i), u(:,i), 20);
    v20(i) = interp1(temp(:,i), v(:,i), 20);
    d20(i) = interp1(temp(:,i), depth, 20);
  end
end

u20 = reshape(u20, nlat, nlon);
v20 = reshape(v20, nlat, nlon);
d20 = reshape(d20, nlat, nlon);

default_global

figure(1); figure_orient;
sp(1)
     gcont(d20/100, [0:25:600]);
     dc
sp(2)
     gquiv(u20, v20, 2.5, 1);
     dc