call solenoidn(5,600.936,6.0,1.,'SOL1 ') CALL FIT( 1, 1, 1, 0.0, 1.0, 1) CALL FIT( 1, 3, 3, 0.0, 1.0, 1)