Iphc developer

L615 of tpl was causing the code to break when number of age obs was less than 3. 
Old code:
cout<<setw(4)<<d3_A(nAgears).sub(n_A_nobs(nAgears)-2,n_A_nobs(nAgears))<<endl;

New code:
cout<<setw(4)<<d3_A(nAgears).sub(n_A_nobs(nAgears)-1,n_A_nobs(nAgears))<<endl;