I was able to read an fchk with IOdata, use gbasis to compute integrals, then write to an FCIDUMP with IOdata.
Please add a warning to the gbasis documentation examples: IOdata assumes integrals are stored internally in physicists' notation, hence transposes to chemists' notation upon dumping an FCIDUMP. Hence, one should not compute the 2 electron integrals with notation='chemist'.
Everything makes sense, and I am sorry for bringing this up. It would help.