#!/usr/bin/env python3

import sys
import numpy as np
import phoebe
import matplotlib.pyplot as plt

# https://phoebe-project.org/docs/2.4/examples/animation_binary_complete

b = phoebe.default_binary()

times = np.linspace(0,0.25,25)

b.add_dataset('mesh', times=times, columns=['visibilities'], dataset='mesh01')

b.run_compute()

b['mesh@model'].plot(fc='visibilities', ec='None', y='vs')
#b['mesh@model'].plot(fc='visibilities', ec='None', y='ws')

fig = plt.figure(figsize=(11,4))

afig, mplanim = b.savefig('test_phoebe0_mesh.gif', fig=fig, tight_layouot=True, draw_sidebars=False, animate=True, save_kwargs={'writer': 'imagemagick'})