It appears that the NBTH and NBPH parameters are based on the Penepma 2012 "bilayer.in" file which is the default input file for thin film calculations. Apparently the boys in Barcelona must have changed this parameter in Penepma 2016.
So the best solution is to simply edit the bilayer.in file for the new parameter. I've attached a new version below that should work with Penepma 2016 and will upload a new Penepma16.zip file with this modification.
All I did was remove these two lines:
NBTH 45 [No. of bins for the polar angle THETA]
NBPH 30 [No. of bins for the azimuthal angle PHI]
And replace them with this line:
NBANGL 45 30 [Nos. of bins for the angles THETA and PHI]
That should work.
To post in-line images, login and click on the Gallery link at the top