PHASE 2: FINAL — CF4 DENSITY GRID # Using official CosmicFlows-4 density field
# ============================================================ # PHASE 2: FINAL — CF4 DENSITY GRID # Using official CosmicFlows-4 density field # ============================================================ import numpy as np import pandas as pd from astropy.io import fits from scipy.interpolate import RegularGridInterpolator from scipy.stats import ks_2samp, ttest_ind import matplotlib.pyplot as plt from astropy.coordinates import SkyCoord import astropy.units as u print("=" * 60) print("PHASE 2: FINAL — CF4 DENSITY GRID") print("=" * 60) # ------------------------------------------------------------ # STEP 1: Load CF4 density grid # ------------------------------------------------------------ print("\n1. Loading CF4 density grid...") hdul = fits.open('CF4_new_64-z008_delta.fits') density_grid = hdul[0].data header = hdul[0].header print(f" Grid shape: {density_grid.shape}") print(f" Density range: {density_grid.mi...