A method for predicting the shape of organic crystals grown from polar solvents is presented. The model is an improvement of the recent method developed by Winn and Doherty (A.I.Ch.E. Journal 44 (1998) 2501) for predicting the shape of organic crystals grown from solvents in which the energy of adhesion at the interfaces is dominated by dispersive forces (e.g., non-polar solvents). The principal characteristic of the new method is that it can account for the role of hydrogen donor and hydrogen acceptor atoms in forming strong bonds at the interface. This technique is a first step towards predicting the shapes of organic crystals grown from polar solvents, and has been applied successfully to predict the shape of adipic acid grown from water, and succinic acid grown from water and from propanol.