This example is ready to go, if only we had the raw data to compare to.  Note that sqrt(iarea/pi) is r(t) to compare to data in Figure 4(a) in Sayag & Worster (2012).