print V, "(%d)" % len(ts)
print " ... model run",
pm = ParcelModel(initial_aerosols, V, T0, S0, P0, console=False)
parcel, aerosols = pm.run(z_top, ts=ts)
parcel = parcel.ix[parcel.index % 1. == 0]
aero_subset = {}
for key in aerosols:
aerosol = aerosols[key]
subset = aerosol.ix[aerosol.index % 1. == 0]
aero_subset[key] = subset
aerosols = aero_subset
pm.write_csv(parcel, aerosols, "./temp_data")
print "done"
print " ... activation"
params = OrderedDict()
with open("./temp_data/parcel.csv", "r") as param_file:
for i, line in enumerate(param_file.readlines()):
if i == 0: continue
key, P, S, T, wc, wv = line.strip().split(',')
params[key] = {'T': float(T), 'S': float(S)}
## Compute Activation stats
from razzak_disco import ActivatedFractionJob
for species in aer_species:
if species == "Mode 2": continue
