hdiff output

r29054/rotamer_move.F90 2015-11-17 23:33:46.369126605 +0000 r29053/rotamer_move.F90 2015-11-17 23:33:46.569129286 +0000
  1: MODULE ROTAMER  1: MODULE ROTAMER
  2:    IMPLICIT NONE  2:    IMPLICIT NONE
  3:    LOGICAL                 :: ROTAMER_MOVET  3:    LOGICAL                 :: ROTAMER_MOVET
  4:    CHARACTER(LEN=100)      :: ROTAMER_SCRIPT  4:    CHARACTER(LEN=100)      :: ROTAMER_SCRIPT
  5:   5: 
  6: CONTAINS  6: CONTAINS
  7:   7: 
  8: SUBROUTINE ROTAMER_INITIALISE()  8: SUBROUTINE ROTAMER_INITIALISE()
  9:    USE PORFUNCS 
 10:    IMPLICIT NONE  9:    IMPLICIT NONE
 11:  10: 
 12: ! If the Makefile/CMake have defined _SVN_ROOT_, we can use that to show where the 11: ! If the Makefile/CMake have defined _SVN_ROOT_, we can use that to show where the
 13: ! rotamer python script should be. If not, rely on the user-defined location. 12: ! rotamer python script should be. If not, rely on the user-defined location.
 14: #ifdef _SVN_ROOT_ 13: #ifdef _SVN_ROOT_
 15:    CALL SYSTEM(_SVN_ROOT_ // '/SCRIPTS/AMBER/rotamer/rotamer.py' // ' coords.prmtop') 14:    CALL SYSTEM(_SVN_ROOT_ // '/SCRIPTS/AMBER/rotamer/rotamer.py' // ' coords.prmtop')
 16: #else 15: #else
 17:    CALL SYSTEM(TRIM(ADJUSTL(ROTAMER_SCRIPT)) // ' coords.prmtop') 16:    CALL SYSTEM(TRIM(ADJUSTL(ROTAMER_SCRIPT)) // ' coords.prmtop')
 18: #endif  17: #endif 
 19: END SUBROUTINE ROTAMER_INITIALISE 18: END SUBROUTINE ROTAMER_INITIALISE
 20:  19: 
 21: SUBROUTINE DUMP_COORDS(COORDS) 
 22: ! Dump coordinates in Amber restart format to .coords_before_rotamer.rst for rotamer moves. 
 23:    USE AMBER12_INTERFACE_MOD, ONLY: AMBER12_WRITE_RESTART 
 24:    USE COMMONS, ONLY: AMBERT, AMBER12T, NATOMS 
 25:    IMPLICIT NONE 
 26: ! Arguments 
 27:    DOUBLE PRECISION, INTENT(IN)  :: COORDS(:) 
 28:  
 29:    IF (AMBERT) THEN 
 30:    ! If not using AMBER 12, then use DUMPRST from hbondmatrix.f90 
 31:       CALL DUMPRST(COORDS(:), NATOMS, '.coords_before_rotamer.rst') 
 32:    ELSE IF (AMBER12T) THEN 
 33:       CALL AMBER12_WRITE_RESTART(COORDS(:), & 
 34:                                  '.coords_before_rotamer.rst', & 
 35:                                  LEN('.coords_before_rotamer.rst')) 
 36:    ELSE 
 37:       WRITE(*, *) 'Rotamer moves can only currently be used with the Amber potentials.' 
 38:       STOP 
 39:    END IF 
 40:  
 41: END SUBROUTINE DUMP_COORDS 
 42:  
 43: SUBROUTINE READ_NEW_COORDS(COORDS) 
 44: ! Read the new coordinates (after rotamer moves). 
 45:    IMPLICIT NONE 
 46: ! Arguments 
 47:    DOUBLE PRECISION, INTENT(OUT) :: COORDS(:) 
 48: ! Variables 
 49:    INTEGER                       :: UNIT_NUM 
 50:    CHARACTER(LEN=100)            :: FILE_NAME 
 51:  
 52:    UNIT_NUM = 5827 
 53:    FILE_NAME = '.coords_after_rotamer.rst' 
 54:    OPEN(UNIT=UNIT_NUM, FILE=FILE_NAME, STATUS='UNKNOWN') 
 55:    READ(UNIT_NUM, '(3F20.10)') COORDS(:) 
 56:    CLOSE(UNIT_NUM) 
 57:  
 58: END SUBROUTINE READ_NEW_COORDS 
 59:  
 60: END MODULE ROTAMER 20: END MODULE ROTAMER


legend
Lines Added 
Lines changed
 Lines Removed

hdiff - version: 2.1.0