The use of organic nonlinear optical (ONLO) materials in electro-optic (EO) modulators requires that the active molecular components (chromophores) be acentrically oriented. The fundamental molecular constituents are in a condensed, glassy phase. Molecular orientation in such systems is typically achieved by applying a DC poling field to the glassy material. We are developing efficient coarse-grained classical Monte Carlo (MC) methods to simulate the order of such systems. The most challenging aspects of these simulations are convergence to an experimentally relevant equilibrium ensemble and verification of simulation accuracy. We use a variety of molecular descriptions and a variety of MC methods to achieve proper order in the shortest number of computational cycles possible. Herein, we illustrate a few examples of the types of calculations and compare with experimental results with representative amorphous organic materials, including electro-optic chromophores.