A previous initial value coupled coherent state (CCS) representation is applied to Gaussian wave packet propagation on multidimensional Henon Heiles potentials. Solutions of the time-dependent integro-differential Schrodinger equation are obtained in a basis of trajectory guided Frozen Gaussian Coherent States, with Monte Carlo sampling to ensure a unique capability for propagating multidimensional wave functions. Results, which are obtained for up to 14 D, are compared with those derived by the Herrnan-Kluk semiclassical initial value representation (IVR) wave packet method.