Documentation of interp_hgt_domsom

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

Help text

a = 'hi'
function make_time_plots

Cross-Reference Information

Listing of script interp_hgt_domsom

clear

for i = 2:3;
i
cd /home/disk/hayes2/dvimont/ccm/ccm3.6/data

if i == 1;
  in1 = 'Z3'
elseif i == 2;
  in1 = 'U'
elseif i == 3;
  in1 = 'V'
end

varid1 = in1; varid2 = 'PS';
[dw1, dw2] = getnc('time_wgr_5093.nc', varid1, varid2);
[dc1, dc2] = getnc('time_cgr_5093.nc', varid1, varid2);
[sw1, sw2] = getnc('time_wgr_hard.nc', varid1, varid2);
[sc1, sc2] = getnc('time_cgr_hard.nc', varid1, varid2);

[lat, lon, lev, hyam, hybm, P0] = getll('wgr_SOM_hard.nc');
[ntim, nlev, nlat, nlon] = size(dw1);

ind = 5:18;
tim = sort([1:120]);  % 2:12:120 3:12:120 4:12:120 11:12:120 12:12:120]);
hyam = hyam(ind); hybm = hybm(ind);
dw1 = dw1(tim,ind,:,:); dw2 = dw2(tim,:,:);
dc1 = dc1(tim,ind,:,:); dc2 = dc2(tim,:,:);
sw1 = sw1(tim,ind,:,:); sw2 = sw2(tim,:,:);
sc1 = sc1(tim,ind,:,:); sc2 = sc2(tim,:,:);
lev = lev(ind);

lev2 = [150 250 500 850 1000];
dwz = atlev(dw1, lev2, dw2, hyam, hybm, P0);
dcz = atlev(dc1, lev2, dc2, hyam, hybm, P0);
swz = atlev(sw1, lev2, sw2, hyam, hybm, P0);
scz = atlev(sc1, lev2, sc2, hyam, hybm, P0);

eval([in1 'dw = dwz;']);
eval([in1 'dc = dcz;']);
eval([in1 'sw = swz;']);
eval([in1 'sc = scz;']);

savestr = [in1 'dw ' in1 'dc ' in1 'sw ' in1 'sc lev2 lat lon'];
savestr = ['save ' in1 '_som_dom.mat ' savestr];

cd /home/disk/hayes2/dvimont/ccm/ccm3.6/data
%save hgt_som_dom.mat dwz dcz swz scz lev2 lat lon
eval(savestr);

end;

exit