5 ! Contributed by Elliott Sales de Andrade
9 dimension xpos
(5000), ypos
(5000)
10 namelist
/geometry
/ xpos
, ypos
15 open
(unit
=4,file
='geometry.in')
16 write(4,'(a)') '$geometry'
17 write(4,'(a)') ' xpos(1)= 0.00, 0.10, 0.20, 0.30, 0.40,'
18 write(4,'(a)') ' ypos(1)= 0.50, 0.60, 0.70, 0.80, 0.90,'
23 open
(unit
=4,file
='geometry.in',status
='old',form
='formatted')
25 close
(4, status
='delete')
27 !print
*, 'xpos', xpos
(1:10), 'ypos', ypos
(1:10)
29 if (any
(xpos
(1:5) /= [0.00, 0.10, 0.20, 0.30, 0.40]))call abort
()
30 if (any
(ypos
(1:5) /= [0.50, 0.60, 0.70, 0.80, 0.90]))call abort
()
31 if (any
(xpos
(6:) /= -huge
(xpos
))) call abort
()
32 if (any
(ypos
(6:) /= -huge
(ypos
))) call abort
()