Test+Problem+Grid+Data

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


 * NOTE:** Consider first if accessing the precomputed initial value data files and reading them with IDL is sufficient for your purposes!

//**THE WEEK OF DECEMBER 3-7 THE SRB SYSTEM (AND HENCE ALSO /gpfs-wan) IS OFF-LINE FOR UPGRADES!**//

To reproduce the generation of one of the data files (and to test that you have the correct work area and UNIX path setup), do code cd ~/KITPcompare07 g2g.x ICS/hd256.g2g code This creates, in the **data/** subdirectory, a raw binary file **hd256.dat**, and a small text file **hd256.dim** describing the data. The binary file contains the density and the three velocity components. To access the data from you own code you probably need to convert the data to a different format.

Other aspects of the data, such as numerical resolution and centering of the velocities may be controlled with parameters in the input file. To change these, copy the file to your own directory, edit it, and run the program again: code cp ICS/hd256.g2g mydata.g2g $EDITOR mydata.g2g g2g.x mydata.g2g code Here's the content of the **hd256.g2g** file: code &in file='cm1000b/d1000t07000b0.dat' offset=-0.5 do_exp=t / &out file='output/hd256.dat' m=256 mv=4 offset=.0 do_smooth=f / code and here's the content of the m**hd256.g2g** file: code &in file='cm1000b/d1000t06600.dat' offset=-0.5 do_exp=t / &out file='output/mhd256.dat' m=256 mv=7 offset=.0 do_smooth=f / code The **&in** and **&out** namelists describe the input and output data, as follows: code &in file='cm1000b/d1000t07000b0.dat'      ! input file (1000x1000x1000 HD) offset=-0.5                           ! velocities centered on left cube faces do_exp=t                              ! input data has ln(rho), so do rho=exp(ln(rho)) after reading / &out file='data/hd256.dat'                 ! output file m=256                                 ! dimension of output data mv=4                                  ! 4 variables (rho,ux,uy,uz), mv=7 for MHD (rho,ux,uy,uz,bx,by,bz) offset=0.0                            ! velocities centered at cube centers do_smooth=t                           ! apply a 3-point (1/4/1) smoothing to the output (should be =f for 1024 case!) / code In the MHD case the pre-computed initial value data files all have offset=-0.5 and div(B) defined as in the Stagger Code (6th order derivatives).