Skip to content

This repository contains some supplementary codes for data analysis and visualization of Jet ( A Relativistic Moving-Mesh Hydro Code) originally written by Paul Duffell.

Notifications You must be signed in to change notification settings

ctiede/Jet_Supplement

Repository files navigation

Jet_Supplement

PyJet.py provides routines to read checkpoint file generated by Jet with (HDF5 output). And output variable values column wise. pyjet_prim output primitive variable value ( 0 for Angle, 1 radius, 2 density, 3 pressure, 4 radial velocity, 5 azimuthal velocity, 6 (Dim=3) rotational velocity and passive scalar. pyjet_cons output conserved variable value ( 0 for angle, 1 radius, 2 mass, 3 energy, 4 radial momentum, 5 azimuthal angular momentum, 6 rotational angular momentum (Dim=3). write.py provides routines to rewrite checkpoint file utilizing h5py.

import PyJet
import sys
filename = sys.argv[1]
Prim = PyJet.pyjet_prim(filename)
Cons = PyJet.pyjet_cons(filename, Dim=3)
Den = np.array(Prim[:,2])

video_parallel.cpp will loop through a given directory to find datafiles matching Regular expression checkpoint.\.h5* and use a plotting script to generate correspoinding image in parallel. It uses C++ boost library. An example PBS script is included.

About

This repository contains some supplementary codes for data analysis and visualization of Jet ( A Relativistic Moving-Mesh Hydro Code) originally written by Paul Duffell.

Resources

Stars

Watchers

Forks

Releases

No releases published

Packages

No packages published

Languages