Page MenuHomec4science
Files F124230868

toSBStream
No OneTemporary
Actions

File Metadata

Created
Fri, Aug 1, 03:53

toSBStream
View Options

#!/usr/bin/env python3
import argparse
import os
####################################################################
# option parser
####################################################################
description="""convert a snapshot file to png SB maps
"""
epilog ="""
"""
parser = argparse.ArgumentParser(description=description,epilog=epilog,formatter_class=argparse.RawDescriptionHelpFormatter)
parser.add_argument(action="store",
dest="files",
metavar='FILE',
type=str,
default=None,
nargs='*',
help='a set of files')
parser.add_argument('--mass',
dest="mass",
type=float,
required=True,
help='total mass')
####################################################################
# main
####################################################################
opt = parser.parse_args()
for f in opt.files:
# get the raw name
name = os.path.splitext(f)[0]
name = os.path.basename(name)
# convert to hdf5
cmd = "./erkal2swift DATA_Erkal/{0:s}.txt --mass {1:g} -o hdf5/{0:s}.hdf5".format(name,opt.mass)
print(cmd)
os.system(cmd)
# create fits files
cmd = "./mockimgs_sb_addfields hdf5/{0:s}.hdf5 -o hdf5/{0:s}.hdf5 --minAge 3000 --maxAge 14000 --ref ~/Nextcloud/AMC_Repository/02-Project/Science/Products/WPS6_Cosmological_Simulations/Reprocessed_Simulations_hdf5/Dwarfs_RevazJablonka2018/hdf5_for_sb/h019.hdf5".format(name)
print(cmd)
os.system(cmd)
# create the fits files
cmd = "mockimgs_sb_compute_images hdf5/{0:s}.hdf5 --instrument arrakihs_vis_SDSSg --distance 35 --los 0 0 1 -o fits/sb_{0:s}_SDSSg_CCD273_84_35Mpc.fits".format(name)
print(cmd)
os.system(cmd)
cmd = "mockimgs_sb_compute_images hdf5/{0:s}.hdf5 --instrument arrakihs_vis_SDSSr --distance 35 --los 0 0 1 -o fits/sb_{0:s}_SDSSr_CCD273_84_35Mpc.fits".format(name)
print(cmd)
os.system(cmd)
cmd = "mockimgs_sb_compute_images hdf5/{0:s}.hdf5 --instrument arrakihs_vis_SDSSi --distance 35 --los 0 0 1 -o fits/sb_{0:s}_SDSSi_CCD273_84_35Mpc.fits".format(name)
print(cmd)
os.system(cmd)
# pngs images
cmd = "./mockimgs_sb_display_fits --add_axes --ax_unit kpc --ax_max 150 --sbmin 25 --sbmax 33 --colorbar fits/sb_{0:s}_SDSSg_CCD273_84_35Mpc.fits.gz -o pngs/sb_{0:s}_SDSSg_CCD273_84_35Mpc.png".format(name)
print(cmd)
os.system(cmd)
cmd = "./mockimgs_sb_display_fits --add_axes --ax_unit kpc --ax_max 150 --sbmin 25 --sbmax 33 --colorbar fits/sb_{0:s}_SDSSr_CCD273_84_35Mpc.fits.gz -o pngs/sb_{0:s}_SDSSr_CCD273_84_35Mpc.png".format(name)
print(cmd)
os.system(cmd)
cmd = "./mockimgs_sb_display_fits --add_axes --ax_unit kpc --ax_max 150 --sbmin 25 --sbmax 33 --colorbar fits/sb_{0:s}_SDSSi_CCD273_84_35Mpc.fits.gz -o pngs/sb_{0:s}_SDSSi_CCD273_84_35Mpc.png".format(name)
print(cmd)
os.system(cmd)

Event Timeline

c4science · Help