clean lra test and remove py script

tests/ncdm_linearresponse.lua

@@ -28,26 +28,27 @@ T_cmb   = 2.7255
 N_eff   = 3.046
 N_nu    = 3                 -- number of neutrinos species (including massless species)
 m_ncdm  = {0.12, 0.06, 0.02}
-ncdm_linearresponse = true
-ncdm_transfer_redshift = z_i
 n_shell = 0
 every_ncdm = 0              -- this defines the ratio of the cdm grid number to the ncdm grid number
-ncdm_freestreaming = false  -- choose whether to treat ncdm as free-streaming for the growth ODE and Poisson source terms
+ncdm_freestreaming = true   -- choose whether to treat ncdm as free-streaming for the growth ODE and Poisson source terms
 ncdm_matterlike = false     -- choose whether to approximate ncdm as matter-like in the background
-
+--lra params
+ncdm_linearresponse = true
+ncdm_transfer_redshift = z_i
 ncdm_transfer_nu_file = "lra_trans.txt"
+
 -------- Perturbations --------
 -- Input powerspectrum and growth rate
 -- We advise using REPS: https://github.com/matteozennaro/reps
 -- Power spectrum units: k P(k) in Mpc/h units
 -- Must be compatible with the cosmology parameters above
 read_powerspectrum = "Pcb.txt"
-read_powerspectrum_ncdm = "Pncdm.txt"
-read_linear_growth_rate = "fcb.txt"
-read_linear_growth_rate_ncdm = "fncdm.txt"
+--read_powerspectrum_ncdm = "Pncdm.txt"
+--read_linear_growth_rate = "fcb.txt"
+--read_linear_growth_rate_ncdm = "fncdm.txt"
 
 linear_density_redshift = z_i      -- the redshift of the input cdm files
-linear_density_redshift_ncdm = z_i -- the redshift of the input ncdm files
+--linear_density_redshift_ncdm = z_i -- the redshift of the input ncdm files
 
 random_seed = 100
 particle_fraction = 1.0
@@ -64,6 +65,6 @@ np_alloc_factor= 4.0      -- Amount of memory allocated for particle
 
 -------- Output ---------------
 -- destination of particle outputs (all particles)
-write_snapshot= "ncdm/fastpm"
+write_snapshot= "ncdm_lra/fastpm"
 -- destination of 1d power spectrum (raw)
-write_powerspectrum = "ncdm/powerspec"
+write_powerspectrum = "ncdm_lra/powerspec"