subroutine PrtKin(e, K) c*********************************************************************** c subroutine PrtKin print out the kinematic variables such as c on-shell momentum for every channel and the total energy c and Isheet variable c*********************************************************************** implicit none integer Nifty, Isheet, nch, i double complex e, K(2) c sheet.h: parameters to decide which sheet the energy is on. common /sheet/ Isheet(2) common /control/ Nifty(15) c********************************************************** c nch = Nifty(2) c write out the channel momenta and Isheet write(8, 900) e write(8, 910) (Isheet(i), i = 1, nch) write(8, 920) do 10000 i = 1, nch write(8, 930) i, K(i) 10000 continue 900 format('0', 'Energy = ', 2(f8.3, 1x), 'MeV') 910 format(1x, 'ISHEET: ', 5i2) 920 format(1x, 'Channel momemta (MeV/c):') 930 format(1x, i2, 1x, 2f12.5) return end