diff --git a/atomic/do_mesh.f90 b/atomic/do_mesh.f90
index 1a15c7bbb..d9414e6b8 100644
--- a/atomic/do_mesh.f90
+++ b/atomic/do_mesh.f90
@@ -15,7 +15,7 @@
       implicit none
       integer::  ibound, ndm, mesh, i
       real(DP) :: zmesh, xmin, xmax, x, dx, rmax, &
-             r(ndm), r2(ndm), rab(mesh), sqr(ndm)
+             r(ndm), r2(ndm), rab(ndm), sqr(ndm)
 !      
       xmax=log(rmax*zmesh)
       mesh=(xmax-xmin)/dx+1