Made the required changes to match local IC check in this case

Examples/ScalarField/params.txt

@@ -7,16 +7,16 @@
 verbosity = 0
 
 # location / naming of output files
-output_path = "/nfs/st01/hpc-gr-epss/eaf49/random-seed-test/config2" # Main path for all files. Must exist!
+output_path = "/nfs/st01/hpc-gr-epss/eaf49/hdf5-super-horizon-test" # Main path for all files. Must exist!
 chk_prefix = ScalarField_
 plot_prefix = ScalarFieldp_
 # restart_file = ScalarField_001000.3d.hdf5
 
 # HDF5files are written every dt = L/N*dt_multiplier*checkpoint_interval
-checkpoint_interval = 300
+checkpoint_interval = 500
 # set to 0 to turn off plot files (except at t=0 and t=stop_time)
 # set to -1 to never ever print plotfiles
-plot_interval = 40
+plot_interval = 80
 num_plot_vars = 8
 plot_vars = h11 h12 h13 h22 h23 h33 Ham Mom
 
@@ -45,7 +45,7 @@ scalar_amplitude = 10.0
 scalar_velocity = -0.00162845913015181
 scalar_mass = 0.01
 N_fine = 128
-tensor_amplitude = 1.e-4
+tensor_amplitude = 1.
 
 #################################################
 # Grid parameters
@@ -57,8 +57,8 @@ tensor_amplitude = 1.e-4
 # NB - if you have a non-cubic grid, you can specify 'N1' or 'N1_full',
 # 'N2' or 'N2_full' and 'N3' or 'N3_full' ( then dx_coarsest = L/N(max) )
 # NB - the N values need to be multiples of the block_factor
-N_full = 64
-L_full = 6.5
+N_full = 128
+L_full = 600.
 
 # Maximum number of times you can regrid above coarsest level
 max_level = 0 # There are (max_level+1) grids, so min is zero (does changing this turn off AMR?)
@@ -126,7 +126,7 @@ extrapolating_vars = phi Pi
 # dt will be dx*dt_multiplier on each grid level
 dt_multiplier = 0.25
 stop_time = 40
-max_steps = 1700
+max_steps = 3500
 
 # Spatial derivative order (only affects CCZ4 RHS)
 max_spatial_derivative_order = 4 # can be 4 or 6

Source/InitialConditions/ScalarFields/RandomField.impl.hpp

@@ -13,8 +13,8 @@
  inline RandomField::RandomField(params_t a_params, InitialScalarData::params_t a_bkgd_params, std::string a_spec_type)
     : m_params(a_params), m_bkgd_params(a_bkgd_params), m_spec_type(a_spec_type)
 {
-    kstar = 16.*(2.*M_PI/m_params.L);
-    epsilon = 0.5;//0.25 * (sqrt(3.)*2.*M_PI/m_params.L); //0.5;
+    kstar = 32.*(2.*M_PI/m_params.L);
+    epsilon = 10.;//0.25 * (sqrt(3.)*2.*M_PI/m_params.L); //0.5;
     H0 = sqrt((8.0 * M_PI/3.0/m_bkgd_params.m_pl/m_bkgd_params.m_pl)
             *(0.5*m_bkgd_params.velocity*m_bkgd_params.velocity + 0.5*pow(m_bkgd_params.m * m_bkgd_params.amplitude, 2.0)));
     norm = pow(m_params.N, 3.);