Home

Here is a simple example of how to use the code:

import pymsz
import os
import numpy as np
from readsnapsgl import readsnapsgl

cc, rr = np.array([500000,500000,500000]), 3000  # cluster centre and radius
simd = pymsz.load_data("snapshot file name", snapshot=True, center=cc, radius=rr)

pj = pymsz.TT_model(simd, npixel=500, axis="z", redshift=0.177105, AR=5.3887782, zthick=2269)

pj.write_fits_image("output file name .fits", clobber=True)