Test+Problem+find_cores.pro

=Test Problem: find_cores.pro=

This IDL procedure is completely analogous to the **g2c.x** Fortran utility, but is simpler to play with. To use it, do code IDL> g=read_grid('snapshot.dat') IDL> s=find_cores(g.rho,g.ux,g.uy,g.uz,size=1,aver_n=1e4)                        ; HD IDL> s=find_cores(g.rho,g.ux,g.uy,g.uz,g.bx,g.by,g.bz,size=1,aver_n=1e4)          ; MHD code In the form initially supplied, the procedure makes no use of the velocity and magnetic field information. You are invited to add criteria based on these, to try to improve the success rate of predictions; do potential collapse sites actually collapse, or not?!

The **size** and **aver_n** parameters (size and average number density) are chosen as discussed elsewhere.