This IDL code creates a thick magneto-static structure with parameters of a typical sunspot in a solar like photosphere - chromosphere. The variable parameters are field strength on the axis, radius, and Wilson depression (displacement of the atmosphere on the axis with respect to the field-free atmosphere). Output are magnetic field vector, pressure and density distributions with radius and height. The structure has azimuthal symmetry. The codes are relatively self explanatory and the download packages contain README files.
Top: Contours of constant density for the MHS solution. The density deficit associated with the Wilson depression is evident.
Bottom: Contours of constant B. The labels indicate the magnetic field strength in units of kG. The thick line denotes the surface at which . The dotted lines indicate magnetic field lines. The red box corresponds to the domain taken for the dynamical simulations and the red arrow indicates the location of the driver. Note that, for the sake of clarity, the vertical axis has been expanded, giving an appearance of reduced inclination of the magnetic field lines.
If you use this code and/or its results, please contact schunker [at] mps.mpg.de, for the appropriate acknowledgement.
This website is best viewed using Firefox, Safari or Internet Explorer browsers.