listData={'seals*','XBT*','argo_in*','argo_at*','argo_pa*'}; listVar={'T','S'};

for ii=2:3;
    if ii==1; dirData='firstguess/profiles/input/';
    elseif ii==2; dirData='firstguess/profiles/output/';
    elseif ii==3; dirData='solution/profiles/output/';
    end;
    
    MITprof_stats_cost(dirData,listData,listVar);    
    MITprof_stats_misfit(dirData,listData,listVar,2004,2008,[10 25 45],[[3 3 1];[0.5 0.5 0.2]]);
    drawnow;    
end;