Documentation of get_monthly_data

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

Help text

cd /home/disk/hayes2/dvimont/ccm/ccm3.6/run/sun/HTRAE/HTRAE/atm/netcdf

Cross-Reference Information

Listing of script get_monthly_data

clean
cd /home/disk/hayes2/dvimont/ccm/ccm3.6/run/sun/HTRAE/EARTH/atm/netcdf

%varn = 'Z3';
%varn = 'Z3SQRD';
varn = 'UU';
varn = 'U';
varn = 'VV';
varn = 'V';
tem = repmat(NaN, [120 18 48 48]);

for i = 9:12;
  ind = i-8;
  if i < 10;
    filin = ['0010-0' num2str(i) '.nc'];
  else
    filin = ['0010-' num2str(i) '.nc'];
  end
  disp(filin)
  nc = netcdf(filin, 'nowrite');
    tem(ind,:,:,:) = nc{varn, 1}(:,:,:,:);
  nc = close(nc);
end

for j = 11:19;
  for i = 1:12; 
    ind = (j-11)*12+4+i;
    if i < 10;
      filin = ['00' num2str(j) '-0' num2str(i) '.nc'];
    else
      filin = ['00' num2str(j) '-' num2str(i) '.nc'];
    end
    disp(filin)
    nc = netcdf(filin, 'nowrite');
      tem(ind,:,:,:) = nc{varn, 1}(:,:,:,:);
    nc = close(nc);
  end
end

for i = 1:8;
  ind = 112+i;
  if i < 10;
    filin = ['0020-0' num2str(i) '.nc'];
  else
    filin = ['0020-' num2str(i) '.nc'];
  end
  disp(filin)
  nc = netcdf(filin, 'nowrite');
    tem(ind,:,:,:) = nc{varn, 1}(:,:,:,:);
  nc = close(nc);
end

nc = netcdf(filin, 'nowrite');
  lat = nc{'lat'}(:);
  lon = nc{'lon'}(:);
  lev = nc{'lev'}(:);
nc = close(nc);

uu = tem;
u = tem;
vv = tem; 
v = tem;

ke = uu+vv;
kep = uu+vv-u.^2-v.^2;
%z3mon = tem;
%z3sqrdmon = tem;

%cd /home/disk/tao/dvimont/matlab/CCM/Htrea/data
%%save z3mon.mat z3mon lat lon lev
%save z3sqrdmon.mat z3sqrdmon lat lon lev


%  Make annual average
cd /home/disk/tao/dvimont/matlab/CCM/Htrea/data
%load z3mon.mat
%load z3sqrdmon.mat

z3psqrd = z3sqrdmon - z3mon.^2;

[z3psqrdmon, z3psqrd] = annave(z3psqrd);

tem = repmat(NaN, size(z3psqrd));
tem(9:12,:,:,:) = z3psqrd(1:4,:,:,:);
tem(1:8,:,:,:) = z3psqrd(5:12,:,:,:);
z3psqrd = tem;

cd /home/disk/tao/dvimont/matlab/CCM/Htrea/data
save z3psqrd_earth.mat z3psqrd z3psqrdmon lat lon lev


cd /home/disk/tao/dvimont/matlab/CCM/Htrea/data

[ke, ke_ann] = annave(ke);
[kep, kep_ann] = annave(kep);

tem = repmat(NaN, size(ke_ann));
tem(9:12,:,:,:) = ke_ann(1:4,:,:,:);
tem(1:8,:,:,:) = ke_ann(5:12,:,:,:);
ke_ann = tem;

tem = repmat(NaN, size(kep_ann));
tem(9:12,:,:,:) = kep_ann(1:4,:,:,:);
tem(1:8,:,:,:) = kep_ann(5:12,:,:,:);
kep_ann = tem;

cd /home/disk/tao/dvimont/matlab/CCM/Htrea/data
%save ke_htrae.mat ke_ann kep_ann lat lon lev;
save ke_earth.mat ke_ann kep_ann lat lon lev;