hdiff output

r31957/sandbox.f90 2017-02-22 09:30:11.778950481 +0000 r31956/sandbox.f90 2017-02-22 09:30:12.002953475 +0000
907:     else907:     else
908:         sandout_name = 'sandout.xyz'908:         sandout_name = 'sandout.xyz'
909:     end if909:     end if
910:     open(unit=sandout_unit, file=sandout_name, status='unknown')910:     open(unit=sandout_unit, file=sandout_name, status='unknown')
911: 911: 
912:     ! loop over saved minima912:     ! loop over saved minima
913:     do min_num = 1, nsave913:     do min_num = 1, nsave
914:         ! put number of atoms and comment line. for now this we'll do just one atom per molecule914:         ! put number of atoms and comment line. for now this we'll do just one atom per molecule
915:         !write(sandout_unit, '(I0)') num_mols915:         !write(sandout_unit, '(I0)') num_mols
916:         write(sandout_unit, '(I0)') num_atoms916:         write(sandout_unit, '(I0)') num_atoms
917: !        write(sandout_unit, '(A,1x,I0,A,1x,F20.10,1x,A,1x,I0)') 'energy of minimum', min_num, '=', qmin(min_num), 'first found at step', ff(min_num)917:         write(sandout_unit, '(A,1x,I0,A,1x,F20.10,1x,A,1x,I0)') 'energy of minimum', min_num, '=', qmin(min_num), 'first found at step', ff(min_num)
918: 918: 
919:         ! loop over molecules919:         ! loop over molecules
920:         do mol_num = 1, size(molecules)920:         do mol_num = 1, size(molecules)
921:             mol => molecules(mol_num)921:             mol => molecules(mol_num)
922: 922: 
923:             call rmdrvt(qminp(min_num, mol%p_index: mol%p_index + 2), rotmat, dummy, dummy, dummy, .false.)923:             call rmdrvt(qminp(min_num, mol%p_index: mol%p_index + 2), rotmat, dummy, dummy, dummy, .false.)
924: 924: 
925:             call update_sandbox_molecule(.false., mol, qminp(min_num, mol%r_index: mol%r_index + 2),qminp(min_num,mol%p_index: mol%p_index + 2))925:             call update_sandbox_molecule(.false., mol, qminp(min_num, mol%r_index: mol%r_index + 2),qminp(min_num,mol%p_index: mol%p_index + 2))
926: 926: 
927:             !write(sandout_unit, '(A6,3F10.5,A12,1X,3F10.5)') mol%atom_symbol, qminp(min_num, mol%r_index: mol%r_index + 2), &927:             !write(sandout_unit, '(A6,3F10.5,A12,1X,3F10.5)') mol%atom_symbol, qminp(min_num, mol%r_index: mol%r_index + 2), &


legend
Lines Added 
Lines changed
 Lines Removed

hdiff - version: 2.1.0