
 Will generate initial conditions on grid of size          256         256         256

 Select type of initial power spectrum (matter transfer function):
     1 for T(k) from linger.dat
     2 for T(k) from approx. analytical fit (BBKS)
     3 for T(k)=1 (scale-free)
     4 for T(k) from approx. analytical fit (E&Hu)
 Enter case (1,2,3 or 4) now:
           4
 Enter Omega_m, Omega_v, H0 (km/s/Mpc, set H0=1 for scale-free case)
  0.300000012      0.699999988       70.0000000    
 Enter Omega_b
   5.00000007E-02
  Initializing power spectrum for Omega_m, Omega_v, H0=  0.300000012      0.699999988       70.0000000    
 Enter long-wave spectral index n (scale-invariant is n=1)
   1.00000000    
 Enter desired normalization at a=1:
   Qrms-ps/micro-K (if > 0) or -sigma_8 (if < 0) = ?
 Computing Delta_2(k) using Sachs-Wolfe approxmation
 This may take several minutes.
 Qrms-ps/micro-K=   20.5477428    

 I can write the power spectrum of delta_rho/rho to disk (power.dat).
 If you would like this, please enter kmin and kmax (1/Mpc)
 Enter 0,0 if you don't want the power spectrum written to a file
   9.99999978E-03   100.000000    
 Enter dx (initial particle spacing in Mpc, not Mpc/h)
    or enter -boxlength in Mpc/h
    i.e. dx in Mpc if > 0, or -boxlength in Mpc/h if < 0
  -50.0000000    
   dx=  0.279017866      Mpc
 Enter 1 if these initial conditions will not be refined, otherwise
   enter the ultimate refinement factor for the smallest grid spacing
   (This is used only to set astart.)
           1
 Scaling initial conditions to starting a=   3.00068725E-02
   when sigma(delta) at ultimate refinement scale=  0.200000003    
 vfact=   221.333496    
 Enter 0 for final output or 1 for further refinement
           0
 Setting output grid to          256         256         256
 Enter <RETURN> to skip output grid size
 Enter <RETURN> to skip output grid offset
 Enter <RETURN> to skip final grid size
 Enter <RETURN> to skip final grid offset
 Will produce output files ic_deltab, ic_velb[xyz], ic_velc[xyz]

 Subgrid white noise:
   Choose irand (1 or 2) from the following list:
     irand=1 to generate new noise and save to file
     irand=2 to read noise from existing file
 Enter irand
           1
   Enter random number seed (9-digit integer, ignored if irand=2)
        1234
   Enter filename of white noise file (or <RETURN> if irand=0)
 white-outfile.dat                                                                                                               
  Enter 1 if padding of large scales is required, 0 otherwise
           0
  Enter padding file name (white noise format, ignored if ipad=0)
 noisefile_dummy.dat                                                                                                             
 Warning: Generating new random numbers in mpic4!
 myid,seeds:            0   822569775
 Normalizing rms, subtracting average
 Mean value of the white noise box =   -2.3314084748130102E-004
 Will be removed
 Standard deviation of the white noise box =   0.99985838683469952     
 Will be factored out so that input stdev=1
 ic4 white noise: chisq, dof, nu=   16777216.0                 16777215   2.44140610E-04
 Statistics of ic4 for idim, itide=           0           0
    Mean sigma, sampled sigma, maximum=  0.243928477      0.243429005       1.40875602    
 Reading random numbers used in ic4 from white-outfile.dat
   Random numbers generated with iseed=        1234
 Normalizing rms, subtracting average
 Mean value of the white noise box =   -2.3314085176806781E-004
 Will be removed
 Standard deviation of the white noise box =   0.99985838683251660     
 Will be factored out so that input stdev=1
 ic4 white noise: chisq, dof, nu=   16777216.0                 16777215   2.44141935E-04
 Statistics of ic4 for idim, itide=           1           0
    Mean sigma, sampled sigma, maximum=  0.219509587      0.218785495      0.914245963    
 Reading random numbers used in ic4 from white-outfile.dat
   Random numbers generated with iseed=        1234
 Normalizing rms, subtracting average
 Mean value of the white noise box =   -2.3314085176806781E-004
 Will be removed
 Standard deviation of the white noise box =   0.99985838683251660     
 Will be factored out so that input stdev=1
 ic4 white noise: chisq, dof, nu=   16777216.0                 16777215   2.44141935E-04
 Statistics of ic4 for idim, itide=           2           0
    Mean sigma, sampled sigma, maximum=  0.219509587      0.204993233      0.792325854    
 Reading random numbers used in ic4 from white-outfile.dat
   Random numbers generated with iseed=        1234
 Normalizing rms, subtracting average
 Mean value of the white noise box =   -2.3314085176806781E-004
 Will be removed
 Standard deviation of the white noise box =   0.99985838683251660     
 Will be factored out so that input stdev=1
 ic4 white noise: chisq, dof, nu=   16777216.0                 16777215   2.44141935E-04
 Statistics of ic4 for idim, itide=           3           0
    Mean sigma, sampled sigma, maximum=  0.219509587      0.175664037      0.743286073    
 Reading random numbers used in ic4 from white-outfile.dat
   Random numbers generated with iseed=        1234
 Normalizing rms, subtracting average
 Mean value of the white noise box =   -2.3314085176806781E-004
 Will be removed
 Standard deviation of the white noise box =   0.99985838683251660     
 Will be factored out so that input stdev=1
 ic4 white noise: chisq, dof, nu=   16777216.0                 16777215   2.44141935E-04
 Statistics of ic4 for idim, itide=           1           0
    Mean sigma, sampled sigma, maximum=  0.219509587      0.218785346      0.922476888    
 Reading random numbers used in ic4 from white-outfile.dat
   Random numbers generated with iseed=        1234
 Normalizing rms, subtracting average
 Mean value of the white noise box =   -2.3314085176806781E-004
 Will be removed
 Standard deviation of the white noise box =   0.99985838683251660     
 Will be factored out so that input stdev=1
 ic4 white noise: chisq, dof, nu=   16777216.0                 16777215   2.44141935E-04
 Statistics of ic4 for idim, itide=           2           0
    Mean sigma, sampled sigma, maximum=  0.219509587      0.204993084      0.794799328    
 Reading random numbers used in ic4 from white-outfile.dat
   Random numbers generated with iseed=        1234
 Normalizing rms, subtracting average
 Mean value of the white noise box =   -2.3314085176806781E-004
 Will be removed
 Standard deviation of the white noise box =   0.99985838683251660     
 Will be factored out so that input stdev=1
 ic4 white noise: chisq, dof, nu=   16777216.0                 16777215   2.44141935E-04
 Statistics of ic4 for idim, itide=           3           0
    Mean sigma, sampled sigma, maximum=  0.219509587      0.175663874      0.748545229    
