Allowed the smoothing scale of the density grid to explicitly be a fr…

…ee parameter in ANAL_PARAMS.H
andreimesinger · Oct 5, 2016 · c43c112 · c43c112
1 parent 156cd3c
commit c43c112
Show file tree

Hide file tree

Showing 2 changed files with 3 additions and 2 deletions.
diff --git a/Parameter_files/ANAL_PARAMS.H b/Parameter_files/ANAL_PARAMS.H
@@ -64,7 +64,7 @@
 /*
   If set to 1, the ZA density field is additionally smoothed (asside from the implicit
   boxcar smoothing performed when re-binning the ICs from DIM to HII_DIM) with a Gaussian
-  filter of width 0.2*BOX_LEN/HII_DIM.  The implicit boxcar smoothing in perturb_field.c
+  filter of width R_smooth_density*BOX_LEN/HII_DIM.  The implicit boxcar smoothing in perturb_field.c
   bins the density field on scale DIM/HII_DIM, similar to what Lagrangian codes do 
   when constructing Eulerian grids. In other words, the density field, \delta,
   is quantized into (DIM/HII_DIM)^3 values. If your usage requires smooth density fields,
@@ -73,6 +73,7 @@
   New in v1.1
  */
 #define SMOOTH_EVOLVED_DENSITY_FIELD (int) (1)
+#define R_smooth_density (float) (0.2)
 
 /*
   Use second-order Lagrangian perturbation theory (2LPT).

diff --git a/Programs/perturb_field.c b/Programs/perturb_field.c
@@ -448,7 +448,7 @@ int main (int argc, char ** argv){
 
     //smooth the field
     if (!EVOLVE_DENSITY_LINEARLY && SMOOTH_EVOLVED_DENSITY_FIELD){
-      HII_filter(updated, 2, 0.2*BOX_LEN/(float)HII_DIM);
+      HII_filter(updated, 2, R_smooth_density*BOX_LEN/(float)HII_DIM);
     }
 
     // save a copy of the k-space density field