Global Index (short | long) | Local contents | Local Index (short | long)
filin = '36yearly.cdf';
This script calls | |
---|---|
clear if 0; cd /home/disk/hayes2/werner/matlab filin = 'histape.cdf.36'; tim = 1:240; nc = netcdf(filin, 'nowrite'); ts = nc{'TS'}(tim,:,:); oro = nc{'ORO'}(tim,:,:); precc = nc{'PRECC'}(tim,:,:); precl = nc{'PRECL'}(tim,:,:); nc = close(nc); [lat, lon] = getll('histape.cdf.36'); else cd /home/disk/hayes2/dvimont/ccm/ccm3.6/run/sun/clim_5093/data filin = 'PRECC_5093.nc'; [precc, precl] = getnc(filin, 'PRECC', 'PRECL'); end [lat, lon] = getll(filin); default_global prec_tot = (precc+precl)*(3600*24*1000); for i = 1:27; prec(i,:,:) = mean(prec_tot(((3*(i-1))+[3:5]),:,:)); end colormap('default'); a = colormap; b = interp2(1:3, 1:size(a,1), a, 1:3, [1:(size(a,1)/9):size(a,1)]'); for i = 1:10; b = [b; [1 0 0]]; end colormap(b); prec_tot(:,1,1) = zeros(size(prec_tot(:,1,1))); prec_tot(:,48,48) = 36*ones(size(prec_tot(:,1,1))); num=1 figure(1) sp(1) gshade(squeeze(prec(num,:,:)),[0:2:50]) sp(2) gshade(squeeze(prec(num+1,:,:)),[0:2:50]) sp(1) n=48; M = moviein(n); for j = 1:n; gshade(squeeze(prec_tot(j,:,:)),[0:2:50]); dc M(:,j) = getframe; end movie(M,5,6)