crash in prt_comment.cpp
Nick Abel found that the attached sim crashes on an FP exception with the trunk@11275. Closer inspection revealed that the culprit is this line in prt_comment.cpp:
differ = fabs(1.-iso_sp[ipH_LIKE][nelem].trans(ipH3p,ipH2s).Emis().TauIn()* rt.DoubleTau/iso_sp[ipH_LIKE][nelem].trans(ipH3p,ipH2s).Emis().TauTot())*100.;
since iso_sp[ipH_LIKE][nelem].trans(ipH3p,ipH2s).Emis().TauTot() is zero.
Replacing the "case B" command with "database H-like Lyman pumping off" seems to make the bug go away, so this command maybe the culprit. But then, this bug seems to be quite sensitive to details in the input script, so that may be a coincidence too.
The sim crashes at the end of the first grid point. It takes a few minutes to get there.