TOP_OF_MAIN_SECTION // Memory management arrmblsize=500000; // increases memory size gradient_structure::set_GRADSTACK_BUFFER_SIZE(12000000); DATA_SECTION //this section reads in the data init_int Nobs; // number of observations init_int Ntrips;// number of trips init_matrix Data(1,Nobs*100,1,8); //data matrix vector divecell(1,Nobs*100); // creates empty vector vector abundance(1,Nobs*100); vector presence(1,Nobs*100); vector depthmean(1,Nobs*100); vector depthsd(1,Nobs*100); vector tripid(1,Nobs*100); !!tripid = column(Data,3);// extracts first column from data matrix !!presence = column(Data,4); !!abundance = column(Data,5); !!depthmean = column(Data,6); !!depthsd = column(Data,7); !!divecell = column(Data,8); !! cout<<" Nobs "<< Nobs<