Test+Problem+Fortran+analysis

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

A (growing) number of analysis procedures is available as part of the Test Problem Utilities. You can either choose to use the precompiled binaries; this requires converting your data back to the same (raw binary) format as the initial state, or else you can modify and recompile the utilities, to make them able to read your private data format. The Fortran analysis procedures are complemented with IDL procedures for reading and plotting the output from the analysis, so normally you only need to modify the I/O read routines of the Fortran procedures, and then possibly details of the plotting aspect of the IDL procedures.

To analyze SPH results, use either the test problem utility **s2g.x** or your own corresponding procedure for converting SPH data to grid data. In the latter case you may want do call the grid output routines in the Test Problem Utility file **io_grid.f90** into your own code, to produce output files ready to be processed by the available analysis procedures.

Fortran procedures (with supporting IDL procedures)

 * **g2p.x (grid2power.f90):** Computes a number of different power spectra (including those of velocity and magnetic field) from a single snapshot
 * **g2sf.x (grid2struturefn.f90):** Computes structure functions (yet to be released)
 * **g2c.x (grid2cores.f90):** Find the potential collapse sites, by using various criteia

IDL procedures
For reference, here is a list of IDL analysis procedures. For more information, see the IDL analysis page.
 * **plot_power,'data/file.pow'**: plots, on the current IDL graphics device, the power spectra computed by the **g2p.x** Fortran procedure
 * **ps_power,'data/file.pow':** plots, as a full page PostScript file, the power spectra computed by the **g2p.x** Fortran procedure
 * **power3d.pro,ux,uy,uz**: Computes and plots 3-D power spectra, optionally compensated (**comp=2**)
 * **plhist,rho,/pdf,/xlog,/ylog**: Plot the density PDF on log-log scales
 * **sf=structurefn(ux,uy,uz)**:Compute the structure functions (latitudinal and longitudinal) of the velocity field
 * **find_cores(rho,ux,uy,uz,bx,by,bz):** Algorithm for finding potential collapse sites