import read_output as rd
from constants import *
from sub_report import *
import cobra_sub as cobra
nargs = len(sys.argv)
if nargs>1:
for i in range(len(sys.argv)):
if sys.argv[i]=='-h' or sys.argv[i]=='h' or sys.argv[i]=='help':
cobra.help_me_emissivities()
sys.exit()
# print out a big cobra sign!
cobra.print_cobra()
# we iterate over blackbody temperatire
t_init = 10000.0 # this is the starting temperture of the blackbody
imax = 10 # maximum steps we iterate over
i_temp_steps = 1000 # step widths in temp space
MACRO = 7
MACRO_ISO=6
# now we want a list of pywind cmds
# we want to make files
# we want electron density
# we want neutral hydrogen density and non neutral
-h print this help message
-c comparison mode, provide more root files to plot (e.g -c 5, provide 5 root files)
'''
# first import modules of use
import csv, sys, os, array, warnings
import matplotlib.pyplot as plt
import numpy as np
import read_output as rd
import cobra_sub as sub
sub.print_cobra()
# filename is provided by command line
print 'Reading your commands from command line'
mode, store = rd.read_args (sys.argv)
print 'Welcome to cobra, the plotting utility for the radiative transfer code, Python.'
rd.setpars() #this just sets some standard parameters, e.g. tex
# if user asked for help we print help message and exit
if mode.help:
sub.help_me_spec()
filename = store.filename
# read the spec file