nt = 12;
t = 1:nt;

nx = 1440;
ny = 320;
x0 = 0.125;
y0 = -39.875;
dx = 0.25;
dy = 0.25;
lon = x0 + dx*((1:nx)-1);
lat = y0 + dy*((1:ny)-1);

nvar = 6;
A = NaN+zeros(nx,ny,nvar,nt);
% TMI monthly climatology
A = load('/home/honuea/hafner/xie/CLIM/trmmC.dat');
