hdiff output

r26555/mdread.f 2017-03-30 13:31:00.476917331 +0100 r26554/mdread.f 2017-03-30 13:31:01.052924987 +0100
702: 702: 
703:    ! Define default water residue name and the names of water oxygen & hydrogens703:    ! Define default water residue name and the names of water oxygen & hydrogens
704:    704:    
705:    data watdef/'WAT ','O   ','H1  ','H2  '/705:    data watdef/'WAT ','O   ','H1  ','H2  '/
706:    706:    
707:    !     ----- READ THE CONTROL DATA AND OPEN DIFFERENT FILES -----707:    !     ----- READ THE CONTROL DATA AND OPEN DIFFERENT FILES -----
708: 708: 
709: ! sf344> commented out, although not necessary later on709: ! sf344> commented out, although not necessary later on
710: !   if (mdout /= "stdout" ) &710: !   if (mdout /= "stdout" ) &
711: !         call amopen(6,mdout,owrite,'F','W')711: !         call amopen(6,mdout,owrite,'F','W')
712: ! sf344> changed 5 to 550712:    call amopen(550,mdin,'O','F','R')    ! sf344> changed 5 to 550
713:    call amopen(550,mdin,'O','F','R')     
714:    if(master) write(MYUNITNEW,9308)713:    if(master) write(MYUNITNEW,9308)
715:    call date_and_time( DATE=date, TIME=time )714:    call date_and_time( DATE=date, TIME=time )
716:    if (master) write(MYUNITNEW,'(12(a))') '| Run on ', date(5:6), '/', date(7:8), '/',  &715:    if (master) write(MYUNITNEW,'(12(a))') '| Run on ', date(5:6), '/', date(7:8), '/',  &
717:         date(1:4), ' at ', time(1:2), ':', time(3:4), ':', time(5:6)716:         date(1:4), ' at ', time(1:2), ':', time(3:4), ':', time(5:6)
718:    if (master.and.owrite /= 'N') write(MYUNITNEW, '(2x,a)') '[-O]verwriting output'717:    if (master.and.owrite /= 'N') write(MYUNITNEW, '(2x,a)') '[-O]verwriting output'
719:    718:    
720:    ! Echo the file assignments to the user:719:    ! Echo the file assignments to the user:
721:    720:    
722:    if (master) write(MYUNITNEW,9700) 'MDIN'   ,mdin(1:70)  , 'MDOUT' ,mdout(1:70) , &721:    if (master) write(MYUNITNEW,9700) 'MDIN'   ,mdin(1:70)  , 'MDOUT' ,mdout(1:70) , &
723:          'INPCRD' ,inpcrd(1:70), 'PARM'  ,parm(1:70)  , &722:          'INPCRD' ,inpcrd(1:70), 'PARM'  ,parm(1:70)  , &


r26555/runmd.f 2017-03-30 13:31:00.764921159 +0100 r26554/runmd.f 2017-03-30 13:31:01.768934501 +0100
1221:       if (nstep <= 0 .and. master .and. facc /= 'A') then1221:       if (nstep <= 0 .and. master .and. facc /= 'A') then
1222:          if(tsgld)then1222:          if(tsgld)then
1223:            ener(48)=sgft1223:            ener(48)=sgft
1224:            ener(49)=tempsg1224:            ener(49)=tempsg
1225:          endif1225:          endif
1226:          rewind(7)1226:          rewind(7)
1227:          call prntmd(nstep,nitp,nits,t,ener,onefac,7,.false.)1227:          call prntmd(nstep,nitp,nits,t,ener,onefac,7,.false.)
1228: 1228: 
1229: !--- BEGIN DIPOLE PRINTING CODE ---1229: !--- BEGIN DIPOLE PRINTING CODE ---
1230:          ! See code further on for comments-explanations1230:          ! See code further on for comments-explanations
1231:          call nmlsrc('dipoles',550,prndipfind)1231:          call nmlsrc('dipoles',5,prndipfind)
1232:          if(prndipfind /= 0 ) then1232:          if(prndipfind /= 0 ) then
1233:             write(MYUNITNEW,*) '------------------------------- DIPOLE INFO ----------------------------------'1233:             write(MYUNITNEW,*) '------------------------------- DIPOLE INFO ----------------------------------'
1234:             write(MYUNITNEW,9018) nstep,t1234:             write(MYUNITNEW,9018) nstep,t
1235:             read (550,'(a)') prndiptest1235:             read (5,'(a)') prndiptest
1236:             call rgroup(natom,natc,nres,prndipngrp,ix(i02),ih(m02), &1236:             call rgroup(natom,natc,nres,prndipngrp,ix(i02),ih(m02), &
1237:                  ih(m04),ih(m06),ih(m08),ix(icnstrgp), &1237:                  ih(m04),ih(m06),ih(m08),ix(icnstrgp), &
1238:                  jgroup,index,npdec, &1238:                  jgroup,index,npdec, &
1239:                  xx(l60),xx(lcrdr),0,0,0,idecomp,5,.false.)1239:                  xx(l60),xx(lcrdr),0,0,0,idecomp,5,.false.)
1240:             rewind(550)1240:             rewind(550)
1241:             if(prndipngrp > 0) then1241:             if(prndipngrp > 0) then
1242:                call printdip(prndipngrp,ix(icnstrgp),xx(lcrd), &1242:                call printdip(prndipngrp,ix(icnstrgp),xx(lcrd), &
1243:                     xx(l15),xx(linddip),xx(Lmass), natom)1243:                     xx(l15),xx(linddip),xx(Lmass), natom)
1244:             end if1244:             end if
1245:             write(6,*) '----------------------------- END DIPOLE INFO --------------------------------'1245:             write(6,*) '----------------------------- END DIPOLE INFO --------------------------------'
2036:          2036:          
2037: !--- BEGIN DIPOLE PRINTING CODE ---2037: !--- BEGIN DIPOLE PRINTING CODE ---
2038: 2038: 
2039:          ! RCW 2nd Dec 2003 - also output dipole information if2039:          ! RCW 2nd Dec 2003 - also output dipole information if
2040:          ! the dipoles namelist has been specified and corresponding2040:          ! the dipoles namelist has been specified and corresponding
2041:          ! groups defined.2041:          ! groups defined.
2042: 2042: 
2043:          ! Check input unit 5 for namelist dipoles2043:          ! Check input unit 5 for namelist dipoles
2044:          ! We expect to find &dipoles followed by a group2044:          ! We expect to find &dipoles followed by a group
2045:          ! specification of the dipoles to output.2045:          ! specification of the dipoles to output.
2046:          call nmlsrc('dipoles',550,prndipfind)2046:          call nmlsrc('dipoles',5,prndipfind)
2047: 2047: 
2048:         if(prndipfind /= 0 ) then2048:         if(prndipfind /= 0 ) then
2049:            !We calculate the dipoles2049:            !We calculate the dipoles
2050:            write(MYUNITNEW,*) '------------------------------- DIPOLE INFO ----------------------------------'2050:            write(MYUNITNEW,*) '------------------------------- DIPOLE INFO ----------------------------------'
2051:            write(MYUNITNEW,9018) nstep,t2051:            write(MYUNITNEW,9018) nstep,t
2052:            9018 format(/1x, 'NSTEP =',i7,1x,'TIME(PS) =',f10.3)2052:            9018 format(/1x, 'NSTEP =',i7,1x,'TIME(PS) =',f10.3)
2053: 2053: 
2054:            !Get the groups for the dipoles - Ideally we only really want2054:            !Get the groups for the dipoles - Ideally we only really want
2055:            !to call this the once but for the time being I will call it2055:            !to call this the once but for the time being I will call it
2056:            !every time2056:            !every time
2057:          2057:          
2058:            read (550,'(a)') prndiptest2058:            read (5,'(a)') prndiptest
2059: 2059: 
2060:            call rgroup(natom,natc,nres,prndipngrp,ix(i02),ih(m02), &2060:            call rgroup(natom,natc,nres,prndipngrp,ix(i02),ih(m02), &
2061:                 ih(m04),ih(m06),ih(m08),ix(icnstrgp), &2061:                 ih(m04),ih(m06),ih(m08),ix(icnstrgp), &
2062:                 jgroup,index,npdec, &2062:                 jgroup,index,npdec, &
2063:                 xx(l60),xx(lcrdr),0,0,0,idecomp,5,.false.)2063:                 xx(l60),xx(lcrdr),0,0,0,idecomp,5,.false.)
2064: 2064: 
2065:           ! Need to rewind input file after rgroup so it is available2065:           ! Need to rewind input file after rgroup so it is available
2066:           ! when we next loop through2066:           ! when we next loop through
2067: 2067: 
2068:            rewind(550)2068:            rewind(550)


legend
Lines Added 
Lines changed
 Lines Removed

hdiff - version: 2.1.0