subroutine grid_get_size (iatom,nsize)
implicit double precision (a-h,o-z)
Return in nsize the size of the grid. If iatom is zero,
then this is for the whole grid, otherwise for the specified atom.
Will force evaluation of
the grid if necessary.