function oaciplot

load /home/ruoying/MERHAB/CTD/2003oc391.mat
x=-1*cruise.lon;
y=cruise.lat;
[x,y]=ll2xy(y,x);


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

  dpost=read_s3r('../output/Tresidual.s3r');
  dprior=read_s3r('../obs/ICS.NLLC4.MJ.T.s3r');
  sca=[2 12]; anm=[0 2];

  tpost=dpost(21:21:end,3);
  tprior=dprior(21:21:end,3);
 
  subplot(2,2,1);
  patch(g.x(e),g.y(e),tprior(e)); 
  set(gca,'xtick',[],'ytick',[]);
  shading interp;caxis(sca); colorbar('horiz')
  title('Prior','fontsize',12,'fontweight','bold')
  ylabel('Surface temperature','fontsize',12,'fontweight','bold');
  hold on; plot(x,y,'*')
  axis equal

  subplot(2,2,2);
  patch(g.x(e),g.y(e),tpost(e));
  set(gca,'xtick',[],'ytick',[]);
  shading interp;caxis([-3 3]); colorbar('horiz')
  title('Posterior','fontsize',12,'fontweight','bold')
  hold on; plot(x,y,'*')
  axis equal

  dpost=read_s3r('../output/Sresidual.s3r');
  dprior=read_s3r('../obs/ICS.NLLC4.MJ.S.s3r');
  sca=[30 33]; anm=[0 5];
 
  tpost=dpost(21:21:end,3);
  tprior=dprior(21:21:end,3);
 
  subplot(2,2,3);
  patch(g.x(e),g.y(e),tprior(e));
  set(gca,'xtick',[],'ytick',[]);
  shading interp;caxis(sca); colorbar('horiz')
  title('Prior','fontsize',12,'fontweight','bold')
  ylabel('Surface salinity','fontsize',12,'fontweight','bold');
  hold on; plot(x,y,'*')
  axis equal

  subplot(2,2,4);
  patch(g.x(e),g.y(e),tpost(e));
  set(gca,'xtick',[],'ytick',[]);
  shading interp;caxis([-2 2]); colorbar('horiz')
  title('Posterior','fontsize',12,'fontweight','bold')
  hold on; plot(x,y,'*')
  axis equal
