function oaciplot
%casco
x =1.0e+05 *[0.7649  1.6633];
y =1.0e+05 *[1.5102    0.6899];

%penobscot
%x =1.0e+05 *[ 1.7369 2.3182];
%y =1.0e+05 *[ 2.2926  1.2886];

p=[x(1),y(1),x(2),y(2)];

run /home/ruoying/Software/opnml_matlab5/startup
g=loadgrid('../mesh/stw2');
e=g.e';

  dpost=read_s3r('../output/Tanalysis.s3r');
  dprior=read_s3r('../obs/ICS.NLLC4.MJ.T.s3r');
  sdpost=read_s3r('../output/Sanalysis.s3r');
  sdprior=read_s3r('../obs/ICS.NLLC4.MJ.S.s3r');
  sca=[2 12]; anm=[0 2];
  nn=length(g.x);
  nnv=21;
  z=reshape(dpost(:,2),nnv,nn)';
  tpost=reshape(dpost(:,3),nnv,nn)';
  tprior=reshape(dprior(:,3),nnv,nn)';
  spost=reshape(sdpost(:,3),nnv,nn)';
  sprior=reshape(sdprior(:,3),nnv,nn)';
  whos
  [zi,di,tpriori]=zslice1(g,z,tprior,p,100);
  [zi,di,tposti]=zslice1(g,z,tpost,p,100);
  [zi,di,spriori]=zslice1(g,z,sprior,p,100);
  [zi,di,sposti]=zslice1(g,z,spost,p,100);
  clf;
  
  tax=[-2,12];
  sax=[28,34];
  subplot(2,2,1);
  pcolor(di,zi,tpriori);caxis(tax);shading interp;colorbar 'h';
  subplot(2,2,2);
  pcolor(di,zi,spriori);caxis(sax);shading interp;colorbar 'h';
  subplot(2,2,3);
  pcolor(di,zi,tposti);caxis(tax);shading interp;colorbar 'h';
  subplot(2,2,4);
  pcolor(di,zi,sposti);caxis(sax);shading interp;colorbar 'h';

  