hdiff output

r22146/mdread.f 2016-07-06 15:30:06.186769603 +0100 r22145/mdread.f 2016-07-06 15:30:06.598775178 +0100
889:    mdin_pb=.false.889:    mdin_pb=.false.
890:    mdin_lmod=.false.890:    mdin_lmod=.false.
891:    mdin_amoeba=.false.891:    mdin_amoeba=.false.
892:    iamoeba = 0892:    iamoeba = 0
893: 893: 
894: ! sf344> additions for a smooth non-bonded cutoff894: ! sf344> additions for a smooth non-bonded cutoff
895:    ifswitch=0895:    ifswitch=0
896:    fswitchbeta=10896:    fswitchbeta=10
897:    nrespa2=1897:    nrespa2=1
898: 898: 
899:    rewind 550 
900:  
901:    call nmlsrc('cntrl',550,ifind)899:    call nmlsrc('cntrl',550,ifind)
902:    if (ifind /= 0) mdin_cntrl=.true.900:    if (ifind /= 0) mdin_cntrl=.true.
903:     
904:    if ( mdin_cntrl ) then 
905:       read(550,nml=cntrl) 
906:    else 
907:       write(MYUNITNEW, '(1x,a,/)') 'Could not find cntrl namelist' 
908:       call mexit(MYUNITNEW,1) 
909:    end if 
910: 901: 
911:    call nmlsrc('ewald',550,ifind)902:    call nmlsrc('ewald',550,ifind)
912:    if (ifind /= 0) mdin_ewald=.true.903:    if (ifind /= 0) mdin_ewald=.true.
913: 904: 
914:    call nmlsrc('pb',550,ifind)905:    call nmlsrc('pb',550,ifind)
915:    if (ifind /= 0) mdin_pb=.true.906:    if (ifind /= 0) mdin_pb=.true.
916: 907: 
917:    call nmlsrc('lmod',550,ifind)908:    call nmlsrc('lmod',550,ifind)
918:    if (ifind /= 0) mdin_lmod=.true.909:    if (ifind /= 0) mdin_lmod=.true.
919: 910: 
920:    call nmlsrc('amoeba',550,ifind)911:    call nmlsrc('amoeba',550,ifind)
921:    if (ifind /= 0) mdin_amoeba=.true.912:    if (ifind /= 0) mdin_amoeba=.true.
922: 913: 
 914:    rewind 550
 915:    if ( mdin_cntrl ) then
 916:       read(550,nml=cntrl)
 917:    else
 918:       write(MYUNITNEW, '(1x,a,/)') 'Could not find cntrl namelist'
 919:       call mexit(MYUNITNEW,1)
 920:    end if
923:    if (ifqnt>0) then921:    if (ifqnt>0) then
924:      qmmm_nml%ifqnt = .true.922:      qmmm_nml%ifqnt = .true.
925:      if (saltcon /= 0.0d0) then923:      if (saltcon /= 0.0d0) then
926:        qm_gb%saltcon_on = .true.924:        qm_gb%saltcon_on = .true.
927:      else925:      else
928:        qm_gb%saltcon_on = .false.926:        qm_gb%saltcon_on = .false.
929:      end if927:      end if
930:      if (alpb == 1) then928:      if (alpb == 1) then
931:        qm_gb%alpb_on = .true.929:        qm_gb%alpb_on = .true.
932:      else930:      else
1094: end subroutine mdread1 1092: end subroutine mdread1 
1095: 1093: 
1096: !+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++1094: !+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
1097: !+ Read constant pH file and initialize.1095: !+ Read constant pH file and initialize.
1098: subroutine cnstphread(stateinf,resstate,protcnt,trescnt,&1096: subroutine cnstphread(stateinf,resstate,protcnt,trescnt,&
1099:       statene,chrgdat,charge)1097:       statene,chrgdat,charge)
1100:    use modamber9, only : myunitnew, mdcrd_unit, mdinfo_unit1098:    use modamber9, only : myunitnew, mdcrd_unit, mdinfo_unit
1101:    use constants, only : AMBER_ELECTROSTATIC1099:    use constants, only : AMBER_ELECTROSTATIC
1102:    implicit none1100:    implicit none
1103: 1101: 
 1102: 
 1103: 
 1104: 
 1105: 
 1106: 
1104: type :: const_ph_info1107: type :: const_ph_info
1105:    sequence1108:    sequence
1106:    integer :: num_states, first_atom, num_atoms, first_state, first_charge1109:    integer :: num_states, first_atom, num_atoms, first_state, first_charge
1107: end type const_ph_info1110: end type const_ph_info
1108: 1111: 
1109:    type (const_ph_info) :: stateinf(0:50-1)1112:    type (const_ph_info) :: stateinf(0:50-1)
1110: !+ Specification and control of Amber's Input/Output1113: !+ Specification and control of Amber's Input/Output
1111: 1114: 
1112: ! File names1115: ! File names
1113: character(len=4096) groupbuffer1116: character(len=4096) groupbuffer


legend
Lines Added 
Lines changed
 Lines Removed

hdiff - version: 2.1.0