[gmx-users] .dat file from gmx densmap does not match prescribed size

[Michele Pellegrino]

Hello,


I am struggling a little bit with gmx densmap. I am trying to obtain a matrix with density values from the .dat file generated by the program; however the number of data in the output file does not match what I am expecting.

I tried to prescribe both the number of bins in each direction (-n1 and -n2 options) and the grid size (-bin option), not at the same, of course:


gmx densmap -f traj.trr -s system.tpr -od densmap.dat -b 800.0 -e 1200.0 -aver y -n1 4772 -n2 3531


OR (equivalently):


gmx densmap -f traj.trr -s system.tpr -od densmap.dat -b 800.0 -e 1200.0 -aver y -bin 0.01


In both cases I get a dimension mismatch (I expect 16849932 values but I get 4359553).

I am using GROMACS 2019.4?.


What am I doing wrong?


Thanks for your help.


Cheers,

Michele