Fully working model; added separate code for calculation of theta parameter for HN given distance at which detection is 50% of max; yet to incorporate in model parameterisation.