f1001 and f1010 in twiss table

[ldeniau]

These RDTs are computed there: https://github.com/MethodicalAcceleratorDesign/MAD-X/blob/88211ac4a0022d29b24a64374592fb72cadc72cf/src/twiss.f90#L2524 And it would be useful to have them in the twiss table for coupling calculation...

Alternatively, it would be more MAD-X oriented to save the coupling matrix Cbar (aka Ctmp in the code), since f1001 and f1010 can be easily recalculated from it.