;Name: ; NUMCUBE ; ;Purpose: ; Main wrapper program to call the data cube and get back out the ; parameters of the cube: pressure, density, temp, br, bth, bphi ; ;Keyword inputs: ; ; CUBENAME: Filename (inluding extention and path) of the numerical datacube ; you want to look at. File must be in the correct format, ; use make_my_cube.pro in this directory to convert your ; cube. ; ; BOUT: how to handle the magnetic field ouside the ; cube. The options are: ; ; BOUT=1; no magnetic field outside the cube (DEFAULT) ; BOUT=2; radial open field. default strength. ; Changes sign at equator. ; note that ; there will probably be a discontionuity at the ; edge of the cube. ; BOUT=3; dipole field. default strength. note that ; there will probably be a discontionuity at the ; edge of the cube. ; ; BO: magnitude of external field for BOUT = 2 or 3 ; DEFAULT 10 Gauss ; ; HYDRO: how to handle the plasma properties outside the cube ; volume. The options are: ; ; HYDRO=0; no plasma outside the cube (DEFAULT) ; HYDRO=1; exponential isothermal hydrostatic equilibrium ; HYDRO=2; radial hydrostatic equilibrium ; HYDRO=3; Vasquez et al 2003 hydrostatic atmosphere ; DEFAULT 0 ; ; DENSPROF (OUTSIDE cube volume) ; For HYDRO=1, DENSPROF = density at coronal base in CGS units of 5d8 ; For HYDRO=2, DENSPROF can be array [A,B,C,D,E,F] ; or multiplier of array [densprof*A,B,densprof*C,D,densprof*E,F] ; (note for widget it has to be multiplier, so, scalar) ; dens = A*r^-B + C*r^-D + E*r^-F ; and inside FOR_INTERP_CUBE program A,B,C,D,E,F defaults to WSM streamer values ; For HYDRO=3, DENSPROF can be array [A1,A2,A3,A4,A5,aa,bb,alpha,beta] ; or multiplier of array as above ; dens=A1*exp(A2/r)*r^-2*(1+A3/r+A4/r^2+A5/r^3) ; temp=TE*(aa+1)/(aa+bb*r^alpha + (1-bb)*r^-beta) ; and inside FOR_DIPOLE program A1,A2,A3,A4,A5,aa,bb,alpha,beta ; defaults to Vasquez streamer values ; Note- default=1 for all ; ; TE ; isothermal temperature value ; DEFAULT 1.5D6 ; ; ODENSPROF, OT0, CDENSPROF, CT0 (INSIDE cube volume) ; - only used if TOPOLOGY tag in datacube. If set: ; odensprof, oT0 sets plasma for open field ; cdensprof, cT0 sets plasma for closed field ; follows same form as DENSPROF exactly for CDENSPROF ; but HYDRO=1 ODENSPROF is CGS units of 1d8 ; and inside FOR_INTERP_CUBE A,B,C,D,E,F defaults to WSM coronal hole values ; or A1,A2,A3,A4,A5,aa,bb,alpha,beta to Vasquez 2003 values ; **NOTE for HYDRO = 0, TOPOLOGY profiles will apply HYDRO=3** ; ; VELIMPOSE - impose a velocity of constant magnitude VELIMPOSE directed along the field ; overwrites any velocity field already loaded in if nonzero ; DEFAULT 0.d0 ; ; ; DTHRES - DOUBLE : Identify POP2 coronal points that are dense and can be treated ; as different population ; If unset or set to 0, will be ignored; if set to non-zero ; will be interpreted as density above which points are ; designated population 2. (If set to 1 will default to 5d8 cgs) ; This will result in creation of POP2DENS zeroed everywhere except ; D>THRESD - this will be set up in NUMCUBE.PRO ; This will overwrite any POP2DENS tag already in CUBE. ; Note if COLDPTS also set, DTHRES will be ignored. ; DEFAULT 0 ; COLDPTS - DOUBLE : Identify points that are cold as different population ; If unset or set to 0, will be ignored; if set to non-zero ; will be interpreted as temperature below which points are ; designated population 2. (If set to 1 will default to 3d4 K) ; This will result in creation of POP2DENS zeroed everywhere except ; T