web sky issue in tsz simulations

[Debabrataadak]

File "/anaconda3/lib/python3.7/site-packages/pysm3/sky.py", line 179, in get_emission output += comp.getemission(freq, weights=weights) File "/.local/lib/python3.7/site-packages/astropy/units/decorators.py", line 253, in wrapper return = wrapped_function(*func_args, **func_kwargs) File "/anaconda3/lib/python3.7/site-packages/pysm3/models/websky.py", line 264, in get_emission get_sz_emission_numba(freqs, weights, self.m.value, is_thermal) << u.uK_RJ TypeError: expected dtype object, got 'numpy.dtype[float32]'

[zonca]

Please provide minimal runnable code that produces this error.

[Debabrataadak]

import pysm3 import pysm3.units as u sky = pysm3.Sky(nside=128, preset_strings=["cib1", "tsz1"]) map= sky.get_emission(115.3 * u.GHz)

both CIB and SZ has returned same error...

p.s: I am seeing co maps are all zero values as well.

[zonca]

Just to make sure that it is not a problem with Python 3.7 please first try with a newer version of python.

[Debabrataadak]

thanks..yes problem resolved for tsz.

However, none of the co's can generate non-zero maps. Values at all pixels are all zero.

[zonca]

CO are line emissions, you need to specify the exact line frequency up to whatever decimals are in the code.

[Debabrataadak]

Can we simulate multiple realisations of CIB and SZ? I do not find to set seed of the realisation in these classes.

[zonca]

No

[Debabrataadak]

I was trying to simulate band-pass integrated maps of CIB.... sky = pysm3.Sky(nside=512, preset_strings=["cib1"]) map_100GHz_bandpass = sky.get_emission(x[1:]*u.GHz, y[1:])

got following warning with all zero values : /Users/debabrata/anaconda3/envs/py39/lib/python3.9/site-packages/pysm3/models/interpolating.py:103: UserWarning: Frequency not supported, requested 0.505515992641449 Ghz < lower bound 1.0 GHz warnings.warn(

map_100GHz_bandpass <Quantity [[0., 0., 0., ..., 0., 0., 0.], [0., 0., 0., ..., 0., 0., 0.], [0., 0., 0., ..., 0., 0., 0.]] uK_RJ>

[zonca]

Please provide minimal runnable code that produces this error.

[Debabrataadak]

sky = pysm3.Sky(nside=512, preset_strings=["cib1"]) map_100GHz_bandpass = sky.get_emission(x[1:]*u.GHz, y[1:])

x and y contains bandpass frequencies and corresponding bandpass response. map_100GHz_bandpass <Quantity [[0., 0., 0., ..., 0., 0., 0.], [0., 0., 0., ..., 0., 0., 0.], [0., 0., 0., ..., 0., 0., 0.]] uK_RJ>

[zonca]

How do I test without them? Please prepare the simplest x and y that show the problem and share an entirely runnable code I can paste and execute.

[zonca]

2 weeks without follow-up, closing. Please open another issue with a runnable code snippet if still need help