Documentation of india_sst

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

Cross-Reference Information

Listing of script india_sst

clear
filin = '/home/disk/hayes/dvimont/ccm3.6/data/INDIA/india_sst.nc';
[lat, lon] = getll(filin);
[sst] = getnc(filin, 'SST');
filin = '/home/disk/hayes2/dvimont/ccm/ccm3.6/run/sun/ct/data/wct.nc';
[pcw, plw] = getnc(filin, 'PRECC', 'PRECL'); pw = pcw + plw; clear pcw plw;
filin = '/home/disk/hayes2/dvimont/ccm/ccm3.6/run/sun/ct/data/cct.nc';
[pcw, plw] = getnc(filin, 'PRECC', 'PRECL'); pc = pcw + plw; clear pcw plw;
get_global
define_global
top = squeeze(mean(pw-pc))*1e7;
figure(1)
sp(1)
     gcont(top, [-5:.2:5]);
     dc
limind1 = [40 105 0 10];
FRAME = limind1;
sp(2)
     gcont(top, [-5:.2:5]);
     dc
[xk, yk] = keep_var(limind1, lon, lat);
ssta = top(yk, xk);
ssta = ssta / max(max(ssta));
[ntim, nlat, nlon] = size(sst);
sstnew = sst;
for i = 1:ntim
  sstnew(i,yk,xk) = squeeze(sst(i,yk,xk)) + ssta;
end
FRAME = [0 360 -90 90];
figure(1)
sp(1)
     gcont(squeeze(mean(sst)), [-5:1:40]);
     dc
sp(2)
     gcont(squeeze(mean(sstnew)), [-5:1:40]);
     dc
filin = '/home/disk/hayes/dvimont/ccm3.6/data/INDIA/india_sst.nc';
nc = netcdf(filin, 'write');
  nc{'SST'}(:) = sstnew;
nc = close(nc);