function []=insitu_diags(dirModel,dirMat,doComp,dirTex);
%object:     driver for insitu_misfit and insitu_cost
%inputs:     dirModel is the model run directory     
%            dirMat is the directory where diagnozed .mat files will be saved
%                     -> set it to '' to use the default [dirModel 'mat/']
%            doComp states whether to compute (1) or display (0)
%            dirTex is the directory where tex and fig files will be created

if dirMat(1)~='/'; dirMat=[pwd '/' dirMat]; end; %make full path
subdirMat={'subdirMat',[dirMat '/']};

subdirData={'subdirData',[dirMat '/profiles/output/']};

if isempty(who('dirTex'));
  addToTex={'addToTex',0}; subdirTex={'subdirTex',''};
else;
  if ~ischar(dirTex); error('mis-specified dirTex'); end;
  if dirTex(1)~='/'; dirTex=[pwd '/' dirTex]; end; %make full path
  addToTex={'addToTex',1}; subdirTex={'subdirTex',[dirTex '/']};
end;

if doComp;
insitu_misfit('',1,subdirData,subdirMat);
insitu_cost('',1,subdirData,subdirMat);
else;
insitu_misfit('',0,subdirData,subdirMat,addToTex,subdirTex);
insitu_cost('',0,subdirData,subdirMat,addToTex,subdirTex);
end;

return
%old example that may or may not work ...

dirModel='./apr1alpha_it0/';
listData={'listData','WOD09_XBT*'};
listVar={'listVar','T'};
suffMat={'suffMat','XBT'};
addToTex={'addToTex',1};

insitu_misfit(dirModel,1,listData,listVar,suffMat);
insitu_cost(dirModel,1,listData,listVar,suffMat);

mkdir([dirModel 'tex']);
write2tex([dirModel 'tex/myPlots.tex'],0);
insitu_misfit(dirModel,0,suffMat,addToTex);
insitu_cost(dirModel,0,suffMat,addToTex);
write2tex([dirModel 'tex/myPlots.tex'],4);