Currently the gravitational field of a celestial body can be customised only in latitude (by defining the coefficients of a zonal expansion). I'm still planning (at some point) to implement a full spherical harmonics expansion in both angular variables (which will allow the definition of arbitrary local variations, but may require a very high order expansion). This hasn't happened yet. See Doc\technotes\gravity.pdf for the current implementation.