2
c $Id: print_output.f,v 1.1 2008-04-18 17:48:13 marat Exp $
5
SUBROUTINE print_output(istep,nprnt,ekin,ecoul,eshrt,ebond,eshel)
15
real*8 ekin,ecoul,eshrt,ebond,eshel
16
real*8 stptmp,stpke,stppe,stpte
24
stptmp=2.0*ekin/boltzmann/degfree
26
stppe=ecoul+eshrt+ebond+eshel
29
if(lnew)write(output,'(/,a12,8a12)')' STEP ',
30
$ ' TOTAL E. ',' KINETIC E. ','POTENTIAL E.',' COUL. E. ',
31
$ ' VDW E. ',' BOND E. ',' SHELL E. ',' TEMPERATURE'
34
if(mod(istep,nprnt).eq.0)write(output,'(i12,8f12.4)')
35
$ istep,stpte,stpke,stppe,ecoul,eshrt,ebond,eshel,stptmp