check_H2Odens.m

Check the water density
Assumes that the water resname is 'SOL'
atom is the atom struct
Box_dim is the box dimension vector [x y z]

Contents

Version
Contact
Examples

Version

2.11

Contact

Please report problems/bugs to michael.holmboe@umu.se

Examples

H2Odens = check_H2Odens(atom,Box_dim)

function H2Odens = check_H2Odens(atom,Box_dim)
nSOL=sum(strncmpi([atom.type],'Ow',2))
nSOLatoms=sum(ismember([atom.resname],{'Water' 'water' 'WAT' 'wat' 'SOL' 'sol'}))
nAtomsperSOL=nSOLatoms/nSOL;
V = Box_dim(1)*Box_dim(2)*Box_dim(3);
H2Odens=nSOL/nAtomsperSOL*18.09/6.022E23*1E24/V;