Prior to this PR, it was not possible to run pyradex on molecules with collision rates for H2 + anything else (electrons, for example) because readdata.f assumes (incorrectly) that >2 collision partners, where one of those partners is H2, always means that the partners are pH2 and oH2.
This PR comments those lines out of readdata.f and leaves the ortho/para handling to pyradex (the python side)
Prior to this PR, it was not possible to run pyradex on molecules with collision rates for H2 + anything else (electrons, for example) because
readdata.fassumes (incorrectly) that >2 collision partners, where one of those partners is H2, always means that the partners are pH2 and oH2.This PR comments those lines out of readdata.f and leaves the ortho/para handling to pyradex (the python side)