TREX/fparser
1
0
Fork 0
mirror of https://github.com/TREX-CoE/fparser.git synced 2024-08-24 22:21:45 +02:00
Code Issues Projects Releases Wiki Activity

xyz file comment line ignored

Browse Source
This commit is contained in:
Ravindra Shinde 2021-02-18 23:21:04 +01:00
parent f393df0756
commit 1f6b634d06
2 changed files with 13 additions and 18 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

21
src/iochamp.f90
View File

@ -165,12 +165,7 @@ PROGRAM iochamp
write(*,*) "number of atoms", number_of_atoms
endif
na = number_of_atoms
if (pline%ntokens == 1) then
! molecule_name = fdf_bnames(pline, 1)
molecule_name = fdf_get('molecule', 'Unknown molecule')
endif
if (pline%ntokens == 4) then
symbol(ia) = fdf_bnames(pline, 1)
do i= 1, 3
@ -224,18 +219,14 @@ PROGRAM iochamp
write(6,*) 'Reading an inline_xyz block '
ia = 1
do while((fdf_bline(bfdf, pline)) .and. (ia .le. na))
do while((fdf_bline(bfdf, pline)))
if (pline%ntokens == 1) then
if (pline%ntokens == 1) then
number_of_atoms = fdf_bintegers(pline, 1)
write(*,*) "Number of atoms", number_of_atoms
endif
na = number_of_atoms
if (pline%ntokens == 1) then
molecule_name = fdf_string('', 'Unknown molecule')
endif
if (pline%ntokens == 4) then
symbol(ia) = fdf_bnames(pline, 1)
do i= 1, 3
@ -264,11 +255,7 @@ PROGRAM iochamp
write(*,*) "Number of atoms", number_of_atoms
endif
na = number_of_atoms
if (pline%ntokens == 1) then
molecule_name = fdf_string('', 'Unknown molecule')
endif
if (pline%ntokens == 4) then
symbol(ia) = fdf_bnames(pline, 1)
do i= 1, 3

10
src/test-champ.inp
View File

@ -9,6 +9,8 @@ sr_eps 0.001
sr_adiag 0.05
energy_tol 0.004
%block optimization_flags
optimize_wavefunction no
optimize_ci No
@ -16,6 +18,8 @@ energy_tol 0.004
optimize_orbitals False
%endblock
# The sub blocks currently does not parse.
%block general
@ -54,7 +58,7 @@ energy_tol 0.004
%block inline_xyz
3
molecule name unknown
hypothetical molecule
Nh 1.50 1.108000 -3
Ox 0.333000 0 45
Pd -1.431000 -1.108000 20
@ -76,4 +80,8 @@ H 0 -0.75545 -0.47116
%block molecule < "caffeine.xyz"
### Two blocks with the same name not supported. Ignored without warning
%block molecule < "benzene.xyz"