Documentation of topography

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

Help text

tem2 = mean(tem2);

Cross-Reference Information

Listing of script topography

clear
cd /home/disk/hayes/dvimont/ccm3.6/data 
nc = netcdf('t21_CSIRO_resolution_IC.nc', 'nowrite');
  lat = nc{'lat'}(:);
  lon = nc{'lon'}(:);
  oro = nc{'PHIS'}(:);
nc = close(nc);
cd /home/disk/hayes/dvimont/ccm3.6/tools/definesurf
nc = netcdf('topo.nc', 'nowrite');
  oro1 = nc{'htopo'}(:);
nc = close(nc);
[nlat, nlon] = size(oro1);
lat1 = 180*[0:(nlat-1)]'/nlat-90;
lon1 = 360*[0:(nlon-1)]'/nlon;
lims = [60 150  -10 50];
[xk, yk] = keep_var(lims, lon, lat);
lat2 = lat(yk); lon2 = lon(xk)'; tem2 = oro(yk,xk);
[lon2, lat2] = meshgrid(lon2, lat2);
tem2 = tem2/9.8;
[xk, yk] = keep_var(lims, lon1, lat1);
lat3 = lat1(yk); lon3 = lon1(xk)'; tem3 = oro1(yk,xk);
[lon3, lat3] = meshgrid(lon3, lat3);
tem4 = zeros(size(tem2));
tem5 = zeros(size(tem3));
clev = [2 100 250 500 750 1000:500:(max(max(tem3))+500)];
for i = 2:length(clev)
     tem4(tem2 >= clev(i-1) & tem2 < clev(i)) = i-1;
end
for i = 2:length(clev)
     tem5(tem3 >= clev(i-1) & tem3 < clev(i)) = i-1;
end
sptalk(2,1,1);
  surf(lon2, lat2, tem2/9.8, tem4);
  axis([lims -5000 5000]);
  shading interp
  caxis([0 length(clev)])
  view(-40, 60);
  colorbar
sptalk(2,1,2);
  surf(lon3, lat3, tem3, tem5);
  axis([lims -5000 5000]);
  shading interp
  caxis([0 length(clev)])
  view(-40, 60);
  colorbar
tem3 = mean(tem3);

figure(1); fo; clf;
sptalk(2,1,1);
plot(lon2, tem2/9.8, '-k', lon3, tem3, '-k');