Adding the script of the biophysical part of Plos CB 2013