write (iout,*) i,ib,beta_h(ib,iparm),&
1.0d0/(1.987d-3*beta_h(ib,iparm)),iperm(i),&
potE(iperm(i),iparm),&
-entfac(iperm(i)),fdimless(i),sumprob
write (iout,*) i,ib,beta_h(ib,iparm),&
1.0d0/(1.987d-3*beta_h(ib,iparm)),iperm(i),&
potE(iperm(i),iparm),&
-entfac(iperm(i)),fdimless(i),sumprob