Documentation of compare_prec

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

Cross-Reference Information

Listing of script compare_prec

cd /home/disk/hayes2/dvimont/ccm/ccm3.6/data
[a, b] = getnc('wgr_SOM_hard.nc', 'PRECC', 'PRECL');
pw = (24*3600*1000)*(a+b);
[a, b] = getnc('cgr_SOM_hard.nc', 'PRECC', 'PRECL');
pc = (24*3600*1000)*(a+b);
[lat, lon] = getll('clim_ice.nc');
get_global
default_global
FRAME = [0 360 -40 40];
top = pw-pc;
oro = getnc('wgr_SOM_relax.nc', 'ORO');
cd /home/disk/hayes2/dvimont/ccm/ccm3.6/data
var1 = 'TS';
a = getnc('wgr_SOM_relax.nc', var1);
b = getnc('som_clim5093.nc', var1);
top = a-b;
top(find(oro ~= 0)) = NaN*ones(size(find(oro ~= 0)));
figure(1)
sp(1)
     gcont(squeeze(mean(top([12 1 2],:,:))), [-10:1:10]);
     dc
clf
cd /home/disk/hayes2/dvimont/ccm/ccm3.6/data
num = 11;
var1 = getnc('wgr_SOM_hard.nc', num, 'Z3');
var2 = getnc('cgr_SOM_hard.nc', num, 'Z3');
var3 = getnc('wgr_ice.nc, num, 'Z3');
var4 = getnc('cgr_ice.nc, num, 'Z3');
var1 = getnc('wgr_SOM_hard.nc', 'PSL');
var2 = getnc('cgr_SOM_hard.nc', 'PSL');
top = squeeze(var1 - var2)/200;
figure(1)
sp(1)
     clma
     mcont(squeeze(mean(top([12 1 2],:,:))), [-20:1:20])
sp(2)
     clma
     mcont(squeeze(mean(top([12 1 2],:,:,:))), [-20:1:20])