This code is
dedicated to diagnose turbulence inside interstellar medium. First, the program creates a fractalcube with a turbulent spectrum index &beta from theoretical models. The grid point number can be adjusted to a higher value so that the spatial resolution of calculation is fine enough.
Structure and correlation functions for two points are computed using fast Fourier transform (FFT) and inverse FFT.
This program consists of the following functions and procedures.
The fractal.pro procedure creates a fractal cube of 3 dimensions. The spectrum
index is denoted by &beta. Fluctuation density inside this fractal cube is the
output and stored to calculate correlation and fluctuation functions.
The slicing.pro procedure allows you to handle with the data cube in order to investigate the turtubulence in 3D.
The spektrt.pro returns a power spectrum as a function of wave number.
Procedures cf_fft.pro and sf_fft.pro are used to compute the correlation and
structure function, respectively. Finally, plots of the structure function and
correlation function are set up.
Back to Homepage