[molpro-user] Density in the cube file

[Rene Spada]

Good day,

                I 'm a graduate student and I'm using molpro to write cube
files for atoms, and then I've wrote a fortran code to read it, because i
want to recover the radial solution.  But when I plot the cube diagonal it
shows something like a gaussian function with the maximum near the nucleus
position. I was expecting to see something like the Hydrogen radial
solution with zero value in the nucleus position. Is there any
transformation that I should carry out to recover the hydrogen radial
solution or am I missing something? The fortran code is attached, as the
input file for molpro.



PS. I've compiled the file using gfortran and to use it is just
"./executable input" where input is the cube file.

-- 
*Rene Felipe Keidel Spada
Físico por Formação
Humano por Natureza
Instituto Tecnológico de Aeronáutica*
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://www.molpro.net/pipermail/molpro-user/attachments/20120306/2f93e272/attachment.html>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: Densidade.f
Type: text/x-fortran
Size: 3480 bytes
Desc: not available
URL: <http://www.molpro.net/pipermail/molpro-user/attachments/20120306/2f93e272/attachment.f>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: h-hf-pVQZ
Type: application/octet-stream
Size: 103 bytes
Desc: not available
URL: <http://www.molpro.net/pipermail/molpro-user/attachments/20120306/2f93e272/attachment.obj>