Here's a Python code of the JRM33 model which could probably be ported to Autoplot's Jython (or used to motivate extensions if necessary): https://github.com/rjwilson-LASP/PSH/blob/main/Jupiter/Python/jovian_jrm33_order13_internal_xyz.py