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