### import modules ###
#local module for project
from directory_master import Foldermap
folder = Foldermap()
folder.activate_environ()
#GCE calculation and local visualization module
from omega import omega
from visualize import visualize
#get some bestfit-value from 'Eris'
from bestfit_param_omega.find_bestfit_param_v0.set_param_from_eris import *
#get values and filename from document
from read_parameter_space import read_param
#Set cmdline arguments
try:
from sys import argv
n = int(argv[1])
except IndexError:
n = 10
#Run calculations of GCE with 'Omega'
timesteps = n
loa_inflow_vals, save_name = read_param(
"inflow_rate") #get values from 'parameter_space.txt'
print "Using inflow values: ", loa_inflow_vals
bestfit_mgal = 1.0e+10
loa_omega_inst = [
omega(special_timesteps=timesteps,
in_out_control=True,
inflow_rate=inflow,
imf_type=bestfit_imf_type,
"""
Description:
-Make a class that inherits everything from omega or sygma
-change the functions that uses sygma to get yields
-add a funky factor to the yields of a single isotope
"""
######################################
### Various imports and structures ###
######################################
#get folder-handling-script
from directory_master import Foldermap
folder = Foldermap()
folder.activate_environ() #set environment for omega before importing
#NuPyCEE one-zone chemical evolution code
from omega import *
#use pandas for data-storage
import pandas as pd
#import namespace of the default bestfit-parameters
import bestfit_param_omega.bestfit_file as default_bestfit_namespace
#####################################################
### Class inheriting omega for experimental usage ###
#####################################################
class experiment(omega):
def __init__(self, input_isotope='Re-187', input_factor=1.0,
input_timesteps=30, dt=0,
bestfit_namespace=default_bestfit_namespace):
self.experiment_isotope = input_isotope
self.experiment_factor = input_factor
### initialize omega ###
"""
Purpose:
Make a class that that inherits from 'Omega',
takes a 'namespace' as argument and call on 'Omega' with all the
parameters from 'namespace'.
Also add a function for saving relevant data to a 2D numpy-matrix.
Description:
"""
#get folder-handling-script
from directory_master import Foldermap
folder = Foldermap()
folder.activate_environ() #set environment/directory for omega
#NuPyCEE one-zone chemical evolution code
from omega import omega
class child_omega(omega):
def __init__(self, bestfit_namespace):
omega.__init__(
self,
galaxy=bestfit_namespace.bestfit_galaxy,
in_out_control=bestfit_namespace.bestfit_in_out_control,
SF_law=bestfit_namespace.bestfit_SF_law,
DM_evolution=bestfit_namespace.bestfit_DM_evolution,
Z_trans=bestfit_namespace.bestfit_Z_trans,
f_dyn=bestfit_namespace.bestfit_f_dyn,
sfe=bestfit_namespace.bestfit_sfe,
outflow_rate=bestfit_namespace.bestfit_outflow_rate,
inflow_rate=bestfit_namespace.bestfit_inflow_rate,
"""
Plot the current bestfit-parameters, in 'current_bestfit.py'.
Then add the plots to a folder, 'current_bestfit_plots/'.
Input arguments:
n: int
Number of special-timesteps.
A value of zero disables special timesteps,
meaning constant timesteps.
"""
### Get modules ###
import sys, time
from directory_master import Foldermap #structure of directory
folder = Foldermap()
folder.activate_environ() #set NuPyCEE-folder
from omega import omega #GCE
from visualize import visualize #visualization of Omega and Eris
from bestfit_param_omega.current_bestfit import *
### User-made exception ###
class Wtf(Exception):
def __init__(self):
print "What the F**K did you do?!"
sys.exit("Exiting!")
### Set global parameters ###
n = 5
save_dir = "current_bestfit_plots/"
try:
n = int(sys.argv[1])
except:
