Documentation of get_monthly_slp

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

Help text

  Load SLP data from disk

Cross-Reference Information

This script calls
clean name

Listing of script get_monthly_slp


clean
fdir = '/home/cdrom/jisao2';
fname = dir(fdir);
nfil = size(fname, 1);

%  Get parameters
i = 1; flag = 0;
while flag == 0;
  if strcmp(char(fname(i).name(1)), '0');
    nc = netcdf([fdir '/' char(fname(i).name)], 'nowrite');
    lat = nc{'lat',1}(:);
    lon = nc{'lon',1}(:);
    nc = close(nc);
    flag = 1;
  end
  i = i + 1;
end

nlat = length(lat);
nlon = length(lon);

%  Start loading data
slp = repmat(NaN, [12*10, nlat, nlon]);

%  CD1
num = 0;
num2 = 0;
for i = 1:nfil;
  disp([num2str(i) ' of ' num2str(nfil) ';  ' char(fname(i).name)])
  if strcmp(char(fname(i).name(1)), '0');
    nc = netcdf([fdir '/' char(fname(i).name)], 'nowrite');
      tem = nc{'PSL', 1}(:,:,:);
    nc = close(nc);
    slp(i-num+num2,:,:) = tem;
  else
    num = num+1;
  end
end
num2 = num2 + nfil - num;

%  CD2
fname = dir(fdir);
nfil = size(fname, 1);
num = 0;
for i = 1:nfil;
  disp([num2str(i) ' of ' num2str(nfil) ';  ' char(fname(i).name)])
  if strcmp(char(fname(i).name(1)), '0');
    nc = netcdf([fdir '/' char(fname(i).name)], 'nowrite');
      tem = nc{'PSL', 1}(:,:,:);
    nc = close(nc);
    slp(i-num+num2,:,:) = tem;
  else
    num = num+1;
  end
end
num2 = num2 + nfil - num;

%  CD3 (sort of)
fdir = '/home/disk/hayes2/dvimont/ccm/ccm3.6/run/sun/HT/HTRAE/cdtem4';
fname = dir(fdir);
nfil = size(fname, 1);
num = 0;
for i = 1:nfil;
  disp([num2str(i) ' of ' num2str(nfil) ';  ' char(fname(i).name)])
  if strcmp(char(fname(i).name(1)), '0');
    nc = netcdf([fdir '/' char(fname(i).name)], 'nowrite');
      tem = nc{'PSL', 1}(:,:,:);
    nc = close(nc);
    slp(i-num+num2,:,:) = tem;
  else
    num = num+1;
  end
end
num2 = num2 + nfil - num;

%  CD4 (sort of)
fdir = '/home/disk/hayes2/dvimont/ccm/ccm3.6/run/sun/HT/HTRAE/cdtem5';
fname = dir(fdir);
nfil = size(fname, 1);
num = 0;
for i = 1:nfil;
  disp([num2str(i) ' of ' num2str(nfil) ';  ' char(fname(i).name)])
  if strcmp(char(fname(i).name(1)), '0');
    nc = netcdf([fdir '/' char(fname(i).name)], 'nowrite');
      tem = nc{'PSL', 1}(:,:,:);
    nc = close(nc);
    slp(i-num+num2,:,:) = tem;
  else
    num = num+1;
  end
end
num2 = num2 + nfil - num;

%  CD5 (sort of)
fdir = '/home/disk/hayes2/dvimont/ccm/ccm3.6/run/sun/HT/HTRAE/atm/netcdf';
fname = dir(fdir);
nfil = size(fname, 1);
num = 0;
for i = 1:nfil;
  disp([num2str(i) ' of ' num2str(nfil) ';  ' char(fname(i).name)])
  if strcmp(char(fname(i).name(1)), '0');
    nc = netcdf([fdir '/' char(fname(i).name)], 'nowrite');
      tem = nc{'PSL', 1}(:,:,:);
    nc = close(nc);
    slp(i-num+num2,:,:) = tem;
  else
    num = num+1;
  end
end
num2 = num2 + nfil - num;

%  Save SLP data
cd ~/matlab/CCM/Htrea/Data2
save SLP_monthly_htrae.mat slp lat lon





%  Do the same for EARTH
clean

%  Load SLP data from disk
fdir = '/home/cdrom/jisao2';
fname = dir(fdir);
nfil = size(fname, 1);

%  Get parameters
i = 1; flag = 0;
while flag == 0;
  if strcmp(char(fname(i).name(1)), '0');
    nc = netcdf([fdir '/' char(fname(i).name)], 'nowrite');
    lat = nc{'lat',1}(:);
    lon = nc{'lon',1}(:);
    nc = close(nc);
    flag = 1;
  end
  i = i + 1;
end

nlat = length(lat);
nlon = length(lon);

%  Start loading data
slp = repmat(NaN, [12*10, nlat, nlon]);

%  CD1
num = 0;
num2 = 0;
for i = 1:nfil;
  disp([num2str(i) ' of ' num2str(nfil) ';  ' char(fname(i).name)])
  if strcmp(char(fname(i).name(1)), '0');
    nc = netcdf([fdir '/' char(fname(i).name)], 'nowrite');
      tem = nc{'PSL', 1}(:,:,:);
    nc = close(nc);
    slp(i-num+num2,:,:) = tem;
  else
    num = num+1;
  end
end
num2 = num2 + nfil - num;

%  CD2
fname = dir(fdir);
nfil = size(fname, 1);
num = 0;
for i = 1:nfil;
  disp([num2str(i) ' of ' num2str(nfil) ';  ' char(fname(i).name)])
  if strcmp(char(fname(i).name(1)), '0');
    nc = netcdf([fdir '/' char(fname(i).name)], 'nowrite');
      tem = nc{'PSL', 1}(:,:,:);
    nc = close(nc);
    slp(i-num+num2,:,:) = tem;
  else
    num = num+1;
  end
end
num2 = num2 + nfil - num;

%  CD3
fname = dir(fdir);
nfil = size(fname, 1);
num = 0;
for i = 1:nfil;
  disp([num2str(i) ' of ' num2str(nfil) ';  ' char(fname(i).name)])
  if strcmp(char(fname(i).name(1)), '0');
    nc = netcdf([fdir '/' char(fname(i).name)], 'nowrite');
      tem = nc{'PSL', 1}(:,:,:);
    nc = close(nc);
    slp(i-num+num2,:,:) = tem;
  else
    num = num+1;
  end
end
num2 = num2 + nfil - num;

%  CD4
fname = dir(fdir);
nfil = size(fname, 1);
num = 0;
for i = 1:nfil;
  disp([num2str(i) ' of ' num2str(nfil) ';  ' char(fname(i).name)])
  if strcmp(char(fname(i).name(1)), '0');
    nc = netcdf([fdir '/' char(fname(i).name)], 'nowrite');
      tem = nc{'PSL', 1}(:,:,:);
    nc = close(nc);
    slp(i-num+num2,:,:) = tem;
  else
    num = num+1;
  end
end
num2 = num2 + nfil - num;

%  CD5 (sort of)
fdir = '/home/disk/hayes2/dvimont/ccm/ccm3.6/run/sun/HT/EARTH/atm/netcdf';
fname = dir(fdir);
nfil = size(fname, 1);
num = 0;
for i = 1:nfil;
  disp([num2str(i) ' of ' num2str(nfil) ';  ' char(fname(i).name)])
  if strcmp(char(fname(i).name(1)), '0');
    nc = netcdf([fdir '/' char(fname(i).name)], 'nowrite');
      tem = nc{'PSL', 1}(:,:,:);
    nc = close(nc);
    slp(i-num+num2,:,:) = tem;
  else
    num = num+1;
  end
end
num2 = num2 + nfil - num;


%  Save SLP data
cd ~/matlab/CCM/Htrea/Data2
save SLP_monthly_earth.mat slp lat lon