from astropy.io import ascii
from astropy import units as u
from poliastro.atmosphere.models import COESA62
from prettytable import PrettyTable
import numpy as np

coesa62 = COESA62()

t = PrettyTable()
t.align='l'
t.border=False

t.add_column("Z [km]", coesa62._Zb_table.to(u.km).round(0).value)
t.add_column("H [km]", coesa62._Hb_table.to(u.km).round(0).value)
t.add_column("Tm [K]", coesa62._Tb_table.to(u.K).round(2).value)
t.add_column("p[mbar]", coesa62._pb_table.to(u.mbar).value)

atmosphere_file = open('coesa62.dat', 'w')
atmosphere_file.write(str(t))