M4111 Part 1: The Sun

Introduction: L1

Solar structure, solar interior: L2-3     Handwritten stuff1     Handwritten stuff2

MHD derivation, general transport equation: L4     Handwritten stuff

MHD derivation, multi-fluid equations: L5     Handwritten stuff

MHD derivation, entire fluid equations: L6     Handwritten stuff

MHD, applications: L7     Handwritten stuff

MHD, linearization: L8     Handwritten stuff

MHD, waves in plasma: L8_part2_handwritten_stuff

Radiative transport in solar atmosphere: L9





Problems 1

Problems 2

Problems 3

Problems 4




Standard solar Model S (Christensen-Dalsgaard, J. et al, Science, 272, 1286, 1996)



MHD code simulated photospheric data:

Temperature and pressure

Density

Momenta: vxx, vyy, vzz

Magnetic field: bxx, byy, bzz


IDL program to load the MHD data

IDL program to display two-dimensional images

See a note on GDL below.



NOTES:

1) The code works with conservative variables. To get the physical velocity the momenta have to be divided by the density.

2) The code uses normalized units for magnetic field. To get the magnetic field in cgs, it has to be multiplied by sqrt(4*pi)

3) In the code, x and z directions are horizontal. y direction is vertical.

4) The code uses cgs units. Velocity is in cm/s, density is in g/cm^3, temperature is in K, pressure is in dyn/cm^2.



Stokes profiles for the magnetically sensitive absorption line of iron FeI 630.2nm

IDL program to read Stokes profiles




If you have no IDL/Matlab/etc installed, you can try installing and using GDL: http://gnudatalanguage.sourceforge.net/

Note: GDL is slightly :) different from IDL (and much more buggy). I was able to compile and install GDL 0.9.3 using macports and the latest source from the GDL website.

If you use 0.9.3, to display two-dimensional images please download these two source files instead of the one above and put them into the local directory:

GDL program to display two-dimensional images

congrid.pro which is used by the above program

If you have GDL version less than 0.9.3, use the command tvscl (e.g. "tvscl, your_array_name(*,50,*)") to plot 2D picture.