Saturday, 28 February 2015

1D Solar Atmosphere Models in IDL: Gingerich et al (HSRA)

The Harvard-Smithsonian Reference Atmosphere (HSRA) by Gingerich et al (1971, 1971SoPh...18..347G) is another widely used semi-empirical model of the idealized plane-parallel solar atmosphere in hydrostatic equilibrium. This model also includes the chromosphere where the hydrogen ionization is solved using the statistical equilibrium equations. The model is still often used as the reference solar atmosphere in 1D and its number of citations steadily grows (> 820 so far). The paper is very clearly written and it's a very recommendable reading especially if you use this model. Beside the model itself, there is several tables with computed intensities, brightness temperature, optical depths at different wavelengths, etc.


The depth-dependent variables are listed in two tables. The first one (Table IV, p.359) contains the variables:
  • optical depth,
  • geometrical depth,
  • temperature,
  • pressure,
  • electron pressure,
  • opacity,
  • density, 
  • hydrogen ionization ratio,
  • distribution of the main electron donors (Mg, Si, Fe, C),
  • radiative and adiabatic gradients.
All these variables are also in my files below. The other table (Table V, p.360) contains the depth distribution of the following variables:
  • Rayleigh scattering contribution in the total absorption,
  • electron scattering in the total absorption,
  • opacity per negative hydrogen ion and per neutral hydrogen,
  • number density of CO,
  • departure coefficient for hydrogen ground state, $b_1$,
  • ionization fraction of the metals that are the main electron donors. 
The data from this table I did not digitized.

The data in my files is OCR-ed and manually double checked. There are two versions of the model: .dat (ascii) and .sav (IDL). In the IDL files the model is stored in the structure variables hsra that has the following tags:
hsra = {h:h, tau5:tau5, t:t, pgas:pgas, pel:pel, kappa:kappa, rho:rho, hion:hion, 
  el_metals:el_metals, el_si:el_si, el_mg:el_mg, el_fe:el_fe, el_c:el_c, 
  rgrad:rgrad, agrad:agrad, variables:variables, reference:reference}
where variables has two self-explanatory tags, varnames and varunits.


The files are all in a tarball:

hsra.tar


Disclaimer: This data is publicly available in Gingerich et al (1971, 1971SoPh...18..347G). If you use the data in a publication please cite it properly. I double-checked that the digitized version here matches the original, still typos are possible no matter how unlikely. If you have any doubts, compare the numbers to those in the paper.

2 comments:

  1. Hi Nikola,
    The tags listed here are not what I get when I download the tarball. Where can I get the .sav file with the original tags? Those are much more useful to me...

    ReplyDelete
    Replies
    1. Hi Erica,

      Sorry for a late reply. I was on vacation.

      I've just checked the file. It does contain the actual data from the Table IV of Gingerich et al, but you are right - the list of tags in the post was wrong (copied from a different model - MACKKL). I corrected the tags listed in the post so that they correspond to the data in the files. I hope this resolves the confusion. Thanks for this correction. N.

      Delete