; ; This IDL procedure outputs model solution data cube with the correct format for the SHINE validation study ; ; Written : Peter MacNeice March 2012 ;-------------------------------------------------------------------------- ; ; Instruction from FORWARD/MODELS/DATACUBE/make_my_cube.pro ; ; ; ; The input should be a save set named 'simname'.dat which should include ; at a minimum (for 3D simulations) these ; six quantities with these same names: ; ; r = fltarr(NR) (r=0 is center of sun) ; th = fltarr(NTH) ; ph = fltarr(NPH) ; ; (NOTE ALL GRIDS MUST BE REGULAR AND ORDERED THEY DO NOT HAVE TO BE ; UNIFORM) Theta and phi should be in RADIANS, where (0 LE theta GE ; pi, zero is the north pole ) and (-pi GE phi LE pi). ; The central longitude of your simulation ; should be zero. ; ; Br = fltarr(NR, NTH, NPH) ; Bth = fltarr(NR, NTH, NPH) ; Bph = fltarr(NR, NTH, NPH) ; ; (note, float used as example; double is fine too): ; ; The FORWARD codes also require specification of electron number ; density, temperature, and pressure. There are several ways to ; proceed. ; ; A) - If you want to use the plasma data from your cube, include the ; following three arrays ; ; dens = fltarr(NR, NTH, NPH), electron # density ; pres = fltarr(NR, NTH, NPH), dyne/cm^2 ; temp = fltarr(NR, NTH, NPH), kelvin ; ; If you do not have all three, you can use the ideal gas law, ; P=2NKT where Boltzman's constant is 1.387e-16 [erg/K] to ; calculate the others. ; ; B) - If the simulation does NOT have a density distribution, or if ; one wants to ignore that distribution and just consider the ; magnetic field in a spherically symmetric hydrostatic ; atmosphere (reasonable for low beta plasmas), ; ; keyword HYDRO should be set to: ; ; HYDRO = 1; exponential isothermal atmosphere ; HYDRO = 2; power law density/pressure in hydrostatic balance ; (note HYDRO=0 by default, and it will look for ; your dens,pres,temp arrays) ;-------------------------------------------------------------------------- CODE TO BE ADDED