1 Code to create non-Gaussian simulations of the CMB 2 -------------------------------------------------- 3 4Included in this distribution is: 5 6makefile - you will need to edit this with the location of 7your HEALPIX directory, f90 compiler etc 8 9sky_ng_sim.f90 - code to simulate full-sky non-Gaussian maps 10sky_ng_sim_bin.f90 - as above, but using disjoint bins in multipole 11space with a different non-Gaussian being used in each bin. 12sky_sub.f90 - module used by the above programs 13sub_ngpdf_sho.f90 - module containing subroutines to draw 14pixel values from non-Gaussian pdf based on the eigenstates 15of a harmonic oscillator 16sub_ngpdf_powergauss.f90 - module containing subroutines to 17draw pixel values from Gaussian pdf raised to some power 18(has a much greater skewness than SHO) 19 20sky.par - example parameter file for sky_ng_sim 21sky_bin.par - example parameter file for sky_ng_sim_bin 22 23 24Integration in Healpix 25---------------------- 26- in progress, sky_ng_sim and sky_ng_sin_bin will eventually be merged into a single code 27- replacement of idum parameter with iseed, to be consistent with synfast. 28- No pgplot, unless the preprocessing directive PGPLOT is given during compilation 29examples: my_ng_sim.par and my_ng_sim_bin.par