The modelling of the habit modifying effect of water on crystalline substances is still problematic, since the approaches for morphology prediction are investigated for organic molecules rather than for inorganic molecules. Therefore, selected well known 3-, 4- and 5-site water models have been implemented as the solvent in a molecular modelling method. By means of the surface docking method, the investigations on crystalline benzoic acid aimed to model the effect of water as a solvent. The modelling comprised additional geometry constraints of water as a flexible and a rigid molecule. By proof of experimental results, one 3-site water model was favoured for the application in morphology predictions. Future modelling work will be simplified, since it can focus on one single model of water. Furthermore, corresponding to the desired level of supersaturation, the best model can be taken from a modelling toolbox. [ABSTRACT FROM AUTHOR]