Test+Problem+Sanity+Checks

=Test Problem Sanity Checks= Pages: Workshop, Overview, CADAC, Data Access, Work Area Setup, Path Setup, Utilities, Grid Data, SPH Data, Running, Reporting, LOG

As a way of checking that you data import / conversion of the initial data was successful, here are some statistics to compare with (these lines are just copied from the conversion utility **g2g.x** output -- version compiled november 21, 2007:

HD case
code FILE=cm1000b/d1000t07000b0.dat 1 INPUT: average,rms = -1.319       1.568 2 INPUT: average,rms = 0.4601       4.967 3 INPUT: average,rms = -0.6169      5.561 4 INPUT: average,rms =  1.252       5.395 code The four lines show the statistics of **ln(rho)**, **Ux**, **Uy,** and **Uz**, respectively, in the input file **cm1000b/d1000t07000b0.dat**.

In the interpolated HD output file **hd256.dat** the corresponding values are (with the .g2g input file shown): code &in file='cm1000b/d1000t07000b0.dat' vx=2,5 offset=-.5 do_exp=t / &out file='hd256.dat' m=256 mv=4 offset=.0 do_smooth=t / 1 OUTPUT: average,rms =  1.000       3.060 2 OUTPUT: average,rms = 0.4600       4.924 3 OUTPUT: average,rms = -0.6169      5.525 4 OUTPUT: average,rms =  1.252       5.360 code The same, for the **hd512.dat** case: code &in file='cm1000b/d1000t07000b0.dat' vx=2,5 offset=-.5 do_exp=t / &out file='hd512.dat' m=512 mv=4 offset=.0 do_smooth=t / 1 OUTPUT: average,rms =  1.000       3.620 2 OUTPUT: average,rms = 0.4600       4.951 3 OUTPUT: average,rms = -0.6169      5.547 4 OUTPUT: average,rms =  1.252       5.382 code and for the **hd1024.dat** case: code &in file='cm1000b/d1000t07000b0.dat' vx=2,5 offset=-.5 do_exp=t / &out file='hd1024.dat' m=1024 mv=4 offset=.0 do_smooth=f / 1 OUTPUT: average,rms =  1.000       4.242 2 OUTPUT: average,rms = 0.4600       4.965 3 OUTPUT: average,rms = -0.6169      5.559 4 OUTPUT: average,rms =  1.252       5.393 code
 * NOTE:** In the 1024^3 smoothing should //**not**// be on; smoothing is on in the other cases to compensate for the lack of dissipation range when downsampling high resolution data.

MHD case
code &in file='cm1000b/d1000t06600.dat' offset=-.5 do_exp=t / 1 INPUT: average,rms = -0.9667      1.431 2 INPUT: average,rms = 0.4204       4.843 3 INPUT: average,rms = -0.3820      5.442 4 INPUT: average,rms = 0.8081       5.187 5 INPUT: average,rms = 0.4274E-10   1.611 6 INPUT: average,rms = 0.2414E-11   1.707 7 INPUT: average,rms = 0.3000       1.667 code The 7 lines show the statistics of **ln(rho)**, **Ux**, **Uy,** **Uz**, **Bx**, **By,** and **Bz**, respectively, in the input file **cm1000b/d1000t06600.dat**.

In the interpolated MHD output file **mhd256.dat** the corresponding values are (with the .g2g input file shown): code &in file='cm1000b/d1000t06600.dat' offset=-.5 do_exp=t / &out file='data/mhd256.dat' m=256 mv=7 offset=-.5 do_smooth=t / 1 OUTPUT: average,rms =  1.000       2.142 2 OUTPUT: average,rms = 0.4205       4.797 3 OUTPUT: average,rms = -0.3820      5.401 4 OUTPUT: average,rms = 0.8081       5.146 5 OUTPUT: average,rms = -0.2397E-04  1.470 6 OUTPUT: average,rms = 0.2498E-04   1.572 7 OUTPUT: average,rms = 0.3000       1.532 code The same, for the **mhd512.dat** case: code &in file='cm1000b/d1000t06600.dat' offset=-.5 do_exp=t / &out file='data/mhd512.dat' m=512 mv=7 offset=-.5 do_smooth=t / 1 OUTPUT: average,rms =  1.000       2.392 2 OUTPUT: average,rms = 0.4204       4.825 3 OUTPUT: average,rms = -0.3820      5.426 4 OUTPUT: average,rms = 0.8081       5.171 5 OUTPUT: average,rms = 0.4207E-05   1.553 6 OUTPUT: average,rms = 0.3006E-05   1.651 7 OUTPUT: average,rms = 0.3000       1.612 code and for the **mhd1024.dat** case: code &in file='cm1000b/d1000t06600.dat' offset=-.5 do_exp=t / &out file='data/mhd1024.dat' m=1024 mv=7 offset=-.5 do_smooth=f / TBD code