Global Index (short | long) | Local contents | Local Index (short | long)
Round 2
| This script calls | |
|---|---|
clean
cd /home/disk/hayes2/dvimont/ccm/ccm3.6/run/sun/HT/HTRAE/cdtem_e
nc = netcdf('ha0242.nc', 'nowrite');
lat = nc{'lat', 1}(:);
lon = nc{'lon', 1}(:);
z = nc{'Z250', 1}(:,:,:,:);
nc = close(nc);
nc = netcdf('ha0243.nc', 'nowrite');
z = [z; nc{'Z250', 1}(:,:,:,:)];
nc = close(nc);
nc = netcdf('ha0244.nc', 'nowrite');
z = [z; nc{'Z250', 1}(:,:,:,:)];
nc = close(nc);
nc = netcdf('ha0245.nc', 'nowrite');
z = [z; nc{'Z250', 1}(:,:,:,:)];
nc = close(nc);
[z2, clim] = remove_mean(z);
[ntim, nlat, nlon] = size(z);
M = moviein(20);
dg(lat, lon); FRAME = [0 360 20 90];
lat2 = YAX*ones(1, nlon)-mean(diff(YAX))/2;
lon2 = ones(nlat, 1)*XAX'-mean(diff(XAX))/2;
figure(1); fl(1); clf;
subplot(1,1,1);
ind = [ 1:20 ];
for i = 1:length(ind);
tem = squeeze(z(ind(i),:,:));
tem = fliplr(tem);
tem2 = squeeze(z2(ind(i),:,:));
tem2 = fliplr(tem2);
tem2(abs(tem2)<100) = NaN;
[h, c] = mcont3(tem, [9000:200:11000], 'stereo', [90 270]);
hold on;
[hh, cc] = surfacem(lat2, lon2, tem2);
hold off;
zd = get(hh, 'ZData'); set(hh, 'ZData', -0.1*ones(size(zd)));
caxis([-500 500]);
set(c, 'linewidth', 2);
dcmfill(-1, 0.8);
gridm on;
axis_limits(10000);
M(:,i) = getframe;
end
movie(M([1 1 1 1:20 20 20 20]), 4, 2)
cd /home/disk/hayes2/dvimont/ccm/ccm3.6/run/sun/HT/EARTH/cdtem_e
nc = netcdf('ha0170.nc', 'nowrite');
lat = nc{'lat', 1}(:);
lon = nc{'lon', 1}(:);
z = nc{'Z250', 1}(:,:,:,:);
nc = close(nc);
nc = netcdf('ha0171.nc', 'nowrite');
z = [z; nc{'Z250', 1}(:,:,:,:)];
nc = close(nc);
nc = netcdf('ha0172.nc', 'nowrite');
z = [z; nc{'Z250', 1}(:,:,:,:)];
nc = close(nc);
nc = netcdf('ha0173.nc', 'nowrite');
z = [z; nc{'Z250', 1}(:,:,:,:)];
nc = close(nc);
[z2, clim] = remove_mean(z);
[ntim, nlat, nlon] = size(z);
M2 = moviein(20);
dg(lat, lon); FRAME = [0 360 20 90];
lat2 = YAX*ones(1, nlon)-mean(diff(YAX))/2;
lon2 = ones(nlat, 1)*XAX'-mean(diff(XAX))/2;
figure(1); fl(1); clf;
subplot(1,1,1);
ind = [ 1:20 ];
for i = 1:length(ind);
tem = squeeze(z(ind(i),:,:));
tem2 = squeeze(z2(ind(i),:,:));
tem2(abs(tem2)<100) = NaN;
[h, c] = mcont3(tem, [9000:200:11000], 'stereo', [90 270]);
hold on;
[hh, cc] = surfacem(lat2, lon2, tem2);
hold off;
zd = get(hh, 'ZData'); set(hh, 'ZData', -0.1*ones(size(zd)));
caxis([-500 500]);
set(c, 'linewidth', 2);
dcmfill(-1, 0.8);
gridm on;
axis_limits(10000);
M2(:,i) = getframe;
end
movie(M2([1 1 1 1:20 20 20 20]), 4, 4)
movie(M([1 1 1 1:20 20 20 20]), 4, 4)
clean
file_base=[7:48];
fint = [];
for i = 1:1;
fint = [fint 73*(i-1)+file_base];
end
% Start with earth
lims = [0 360 20 90];
tim = 1:5;
cd /home/disk/hayes2/dvimont/ccm/ccm3.6/run/sun/HT/EARTH/cdtem_e
% Get parameters
nc = netcdf('ha0001.nc', 'nowrite');
[lat, lon, yk, xk] = get_nclatlon(lims, nc);
nc = close(nc);
nlat = length(lat); nlon = length(lon);
zdat = repmat(NaN, [5*length(fint) nlat nlon]);;
for i = 1:length(fint);
disp(fint(i))
if fint(i) < 10;
fyr = ['000' num2str(fint(i))];
elseif fint(i) < 100;
fyr = ['00' num2str(fint(i))];
elseif fint(i) < 1000;
fyr = ['0' num2str(fint(i))];
end
fname = ['ha' fyr '.nc'];
nc = netcdf(fname, 'nowrite');
z = nc{'Z500', 1}(:, yk, xk);
nc = close(nc);
zdat(5*(i-1)+[1:5], :, :) = z;
end
% HTRAE run
cd /home/disk/hayes2/dvimont/ccm/ccm3.6/run/sun/HT/HTRAE/cdtem_e
% Get parameters
nc = netcdf('ha0001.nc', 'nowrite');
[lat, lon, yk, xk] = get_nclatlon(lims, nc);
nc = close(nc);
nlat = length(lat); nlon = length(lon);
zdat2 = repmat(NaN, [5*length(fint) nlat nlon]);;
for i = 1:length(fint);
disp(fint(i))
if fint(i) < 10;
fyr = ['000' num2str(fint(i))];
elseif fint(i) < 100;
fyr = ['00' num2str(fint(i))];
elseif fint(i) < 1000;
fyr = ['0' num2str(fint(i))];
end
fname = ['ha' fyr '.nc'];
nc = netcdf(fname, 'nowrite');
z = nc{'Z500', 1}(:, yk, xk);
nc = close(nc);
zdat2(5*(i-1)+[1:5], :, :) = z;
end
zdat(abs(zdat)>1e4)=NaN;
zdat2(abs(zdat2)>1e4)=NaN;
[ntim, nlat, nlon] = size(zdat2);
zdat2 = reshape(shiftdim(zdat2, 2), nlon, ntim*nlat);
zdat2 = flipud(zdat2);
zdat2 = shiftdim(reshape(zdat2, nlon, ntim, nlat), 1);
% Get highpass stats
zlow = rave(zdat, 5);
zhi = zdat - zlow;
zbp = zlow - rave(zdat, 15);
zlow2 = rave(zdat2, 5);
zhi2 = zdat2 - zlow2;
zbp2 = zlow2 - rave(zdat, 15);
zlow2b = rave(zdat, 15);
% Take a look at some movies here
figure_landscape(1);
global_axes(6,6,0,0,1.5);
subplot2(1,1);
global_latlon(lat, lon, lims);
map_axis('stereo', [90 270]);
cd /home/disk/tao/dvimont/matlab/CCM/Htrea/Figs_movie
lag = 60;
M = moviein(90);
for i = 1:90;
[c, h] = map_contour_pn(zlow(i+lag,:,:), 100);
set(h, 'linewidth', 2);
hold on;
map_surface(zhi(i+lag,:,:), -1);
hold off;
caxis([-150 150]);
draw_landmap;
gridm on; framem on; tightmap2;
set(gca, 'visible', 'off');
eval(['print -depsc2 f' num2str(i) '.eps']);
% M(:,i) = getframe;
end
lag = 60;
M2 = moviein(90);
for i = 1:90;
[c, h] = map_contour_pn(zlow2(i+lag,:,:), 100);
set(h, 'linewidth', 2);
hold on;
map_surface(zhi2(i+lag,:,:), -1);
hold off;
caxis([-150 150]);
draw_landmap;
gridm on; framem on; tightmap2;
set(gca, 'visible', 'off');
M2(:,i) = getframe;
end