Nt = generate_data([800],[0.8 0.8],1e-1,0.15,'../cod_admb/cod',1,15,2,1*14*2,1980);
% (...run ADMB on the data set somehow, either from matlab/octave or in a separate terminal window...)
datpar2plot('../cod_admb/cod'); 
add_N_to_plot(Nt,1980,gcf)