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 ()