The model of Beard and Chuang (1987), using the complete form of Laplace’s formula and adjustments to the aerodynamic pressure distribution for the effect of drop distortion, has been extended to raindrop shapes under the influence of vertical electric fields and drop charges. A finite volume method with numerically generated transformation to a boundary‐fitted coordinate system was used to calculate the shape‐dependent electric field. Sufficient constraints (viz, drop volume, overall force balance, and shape‐dependent surface distributions of aerodynamic and electrostatic stresses) allow the calcualtion of a unique shape by integration from the upper to lower pole using a multiple iteration scheme. The model has been verified against solutions for a stationary drop in a uniform electric field (Taylor, 1964; Brazier‐Smith 1971; Zrnicet al. 1984). Numerical shapes of drops falling in electric fields show a pronounced extension of the upper pole. The increased fall speed of electrostatically stretched drops enhances the aerodynamic flattening of the base. The resultant triangular drop profiles are similar to wind tunnel observations (Richards & Dawson 1973; Rasmussenet al. 1985).