%compare simulated iron with observations 
load ../dat/mod/processed/feMonth_1918.mat x y SFE %unit
load ../dat/obs/TAGLIABUE_FE_DATABASE_OCT2013.asc
oD=TAGLIABUE_FE_DATABASE_OCT2013; clear TAGLIABUE_FE_DATABASE_OCT2013
oD(:,3)=westoff(oD(:,3),360);

MSFE=mean(SFE,3);

for m=1:12
 Io=find(oD(:,1)==m & oD(:,5)<=30 & oD(:,4)<-30);
 D=oD(Io,:);
 M=D(:,1)+NaN;

 for i=1:length(Io)
  Im=find(x>=D(i,3)); 
  Jm=find(y>=D(i,4)); 
  if length(Im)>0 & length(Jm)>0
   M(i)=SFE(Im(1),Jm(1),m);
   M(i)=MSFE(Im(1),Jm(1),1);
  end
 end

 figure
 subplot(2,2,1)
  flat_pcolor(x,y,squee(SFE(:,:,m))'); caxis([0 1]); colorbar
  hold on; plot(D(:,3),D(:,4),'k+');
  set(gca,'ylim',[-70 -30])
 subplot(2,2,2)
  plot(D(:,6),M,'+')
  xlabel('Observed Fe [nmol/l]');
  ylabel('Simulated Fe [nmol/l]');
end
