Exemplo n.º 1
0
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
Exemplo n.º 2
0
    -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