import sys, os, fnmatch
import numpy as np
from galacticus.utils.progress import Progress
from galacticus.io import GalacticusHDF5
from galacticus.dust.Ferrara2000 import dustAtlas
from galacticus.dust.CharlotFall2000 import CharlotFall2000
from galacticus.dust.screens import dustScreen
hdf5File = sys.argv[1]
ATLAS = dustAtlas()
CHARLOTFALL = CharlotFall2000()
SLAB = dustScreen()
# Open HDF5 file for updating
GALACTICUS = GalacticusHDF5(hdf5File, 'a')
# Loop over outputs
redshifts = GALACTICUS.outputs.z
for i, z in enumerate(redshifts):
print "-----> Processing z = " + str(z) + " (iz = " + str(
i + 1) + "/" + str(len(redshifts)) + ")"
# Get list of available properties at this redshift
allDatasets = GALACTICUS.availableDatasets(z)
if len(allDatasets) == 0:
continue
visibleDatasets = list(np.unique(fnmatch.filter(allDatasets,"*LineLuminosity*"))) + \
list(np.unique(fnmatch.filter(allDatasets,"*LuminositiesStellar*")))