I need to calculate the density maps (or preferably density data) on xy-plane for each "slice"
along z, but g_densmap and g_density seem to calculate only the" average" density over z direction. 
I searched and found the following discussion about a similar topic last September. According to
Berk in the discussion, g_densmap in gromacs-3.3.2 should now be able to calculate the density map
per slice. Does anyone know if that is indded so and, if so, how do I use that option? I searched
the manual on g_densmap in gromacs-3.3.2 but couldn't find the option. 
BTW my gromacs version is 3.3.1, I assume I cannot use that new option. Is there any simple way to
do the job on the old version 3.3.1? I am a gromacs beginner and have little programming skill. Any 
help and/or advice would be greatly appreciated.

Sung


---------------------------------------------------------------------------------------------------------------------------
    From: "Berk Hess" <[EMAIL PROTECTED]>
    Reply-To: Discussion list for GROMACS users <gmx-users@gromacs.org>
    To: gmx-users@gromacs.org
    Subject: RE: [gmx-users] g_densmap per slice
    Date: Fri, 21 Sep 2007 11:16:21 +0200

        From: gil claudio <[EMAIL PROTECTED]>
        Reply-To: Discussion list for GROMACS users <gmx-users@gromacs.org>
        To: gmx-users@gromacs.org
        Subject: [gmx-users] g_densmap per slice
        Date: Fri, 21 Sep 2007 01:37:45 -0700 (PDT)

        Is there a way to use g_densmap for a slice of the
        box?  For example, g_densmap can average the density
        of a box on the x-y plane, and thus average through
        the z-axis.  If, for example, the z-axis is 10 nm
        deep, can I use g_densmap to get a slice along the
        z-axis, for example, along the slice between 1nm to
        2nm?

    Not without editing the code.
    But it should not be do difficult to edit the code to make
    g_densmap do what you want.

    Berk.

This seems a useful option, so I have implemented it in CVS,
also for the 3.3.2 release, which should hopefully come out today.

Berk.