hdiff output

r29254/chirality.F90 2015-11-17 23:32:47.880342206 +0000 r29253/chirality.F90 2015-11-17 23:32:48.084344942 +0000
 78: ! doi = {10.1002/(SICI)1096-987X(19960715)17:9<1132::AID-JCC5>3.0.CO;2-T}, 78: ! doi = {10.1002/(SICI)1096-987X(19960715)17:9<1132::AID-JCC5>3.0.CO;2-T},
 79: ! pages = {1132--1141}, 79: ! pages = {1132--1141},
 80: ! year = {1996}, 80: ! year = {1996},
 81: ! } 81: ! }
 82:    angle = atan2(ddot(3, b1xb2_x_b2xb3, 1, b2_norm, 1), ddot(3, b1xb2, 1, b2xb3, 1)) 82:    angle = atan2(ddot(3, b1xb2_x_b2xb3, 1, b2_norm, 1), ddot(3, b1xb2, 1, b2xb3, 1))
 83:  83: 
 84: end function dihedral 84: end function dihedral
 85:  85: 
 86: function cis_trans(coords, tolerance_in) result(state) 86: function cis_trans(coords, tolerance_in) result(state)
 87: ! Works out whether a molecule is cis, trans or neither, to within a certain tolerance 87: ! Works out whether a molecule is cis, trans or neither, to within a certain tolerance
 88: ! (in degrees). If tolerance isn't specified, uses the default of 30.0 degrees (cis_trans_tol). 88: ! (in degrees). If tolerance isn't specified, uses the default of 5.0 degrees.
 89:    double precision, dimension(12), intent(in)  :: coords 89:    double precision, dimension(12), intent(in)  :: coords
 90:    double precision, intent(in), optional       :: tolerance_in 90:    double precision, intent(in), optional       :: tolerance_in
 91:    double precision                             :: tolerance 91:    double precision                             :: tolerance
 92:    character                                    :: state 92:    character                                    :: state
 93:    double precision                             :: angle_degrees 93:    double precision                             :: angle_degrees
 94:    double precision                             :: pi 94:    double precision                             :: pi
 95:  95: 
 96: ! Set tolerance appropriately to 30 degrees if not specified 96: ! Set tolerance appropriately to 5 degrees if not specified
 97:    if (present(tolerance_in)) then 97:    if (present(tolerance_in)) then
 98:       tolerance = tolerance_in 98:       tolerance = tolerance_in
 99:    else 99:    else
100:       tolerance = cis_trans_tol100:       tolerance = cis_trans_tol
101:    end if101:    end if
102: 102: 
103: 103: 
104: ! Calculate the dihedral angle in degrees104: ! Calculate the dihedral angle in degrees
105:    pi = 4.0d0 * atan(1.0d0) 105:    pi = 4.0d0 * atan(1.0d0) 
106:    angle_degrees = 180.0d0 * (dihedral(coords) / pi)106:    angle_degrees = 180.0d0 * (dihedral(coords) / pi)


legend
Lines Added 
Lines changed
 Lines Removed

hdiff - version: 2.1.0