#!/usr/bin/env python3

import sys
import phoebe

logger = phoebe.logger()

b = phoebe.default_binary()

b.add_dataset('lp', times=[0.25,0.50], wavelengths=phoebe.linspace(549,551,101), flux_densities=phoebe.linspace(0.9,0.9,101), time=0.25)

print("b['lp01@lp@dataset'] = ", b['lp01@lp@dataset'])

b.run_compute(irrad_method='none')

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

print("b['lp01@phoebe01@latest@lp@model'] = ", b['lp01@phoebe01@latest@lp@model'])

b.plot(show=True)