Documentation of som_prec

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

Help text

  filin = '36yearly.cdf';

Cross-Reference Information

Listing of script som_prec

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)