#!/usr/bin/env python3

import sys
import numpy as np
import phoebe

b = phoebe.default_binary()

b.add_dataset('mesh', compute_times=[0.05], columns=['visibilities', 'loggs'])

b.run_compute()

f = open('twigs.txt', 'w')
for twig in b.twigs:
  f.write("%s\n" % (twig))
f.close()

fig, plt = b.plot(kind='mesh', fc='loggs', show=True)
plt.savefig("test_phoebe0_loggs.png")