D.7 Example: A model with hydrogen and helium

T = 100000 K, log g = 7 and number fractions H = 0.9 and He = 0.1
Basic H and He model atoms are provided by TMAD ( http://astro.uni-tuebingen.de/~TMAD). These are combined then in ${AA}/H+He.

D.7.1 ATOMS2

• Copy the job file in /home/<user>/jobs/atoms2 to H+He.job. It needs a file H+He with the atomic data in /home/<user>/adaten. Edit both files the test_atmos.job and the test_atoms_lte.job (with emacs). Make sure you use the atomic data file H+He.
• After editing, execute the job file and create an output file for both jobs and check the output files for error messages.
  

D.7.2 SETF2

• Edit both files H+He_lte_setf2 and H+He_setf2 in /home/<user>/jobs/setf2. Set the temperature on 100000K and make sure that it refers to H+He_lte.atoms2.

  ---

  Table 6: Example for a setf2.job in the case of LTE and NLTE model atmosphere calculations LTE model NLTE model #!/bin/sh #!/bin/sh set +x;. ${HOME}/.jobstart set +x;. ${HOME}/.jobstart # # TT=0100000 TT=0100000 # # type='_lte' type=' ' # # cp ${AA}/H+He${type}.atoms2 ATOMS cp ${AA}/H+He${type}.atoms2 ATOMS # # ${BI}/setf2${sys} << % ${BI}/setf2${sys} << % $TT $TT PRINT CHECK PRINT CHECK 1.30E+16 1.30E+16 -1 -1 10 10 10 10 % % # # if test -s FGRID if test -s FGRID then then cp FGRID ${FF}/H+He_${TT}${type}.setf2 cp FGRID ${FF}/H+He_${TT}${type}.setf2 else else echo 'no FGRID created' echo 'no FGRID created' fi fi # # set +x; ${HOME}/.jobend ${TMPDIR} set +x; ${HOME}/.jobend ${TMPDIR}

  ---

• Run both jobs and create output files for them. The files H+He_100000(_lte).setf2 and are written to /home/<user>/fgrids.

D.7.3 Parameter files

• Extract the parameters from the H+He_atoms.out file with
  grep para H+He_atoms.out and from the output file H+He_setf2.out of the
  frequency grid If you want to save the parameters you can write them in a file with
  grep para H+He_atoms.out > object.out
• the parameter files for PRO2 can be found in /home/<user>/parameter/pro2
• copy the existing parameter files to PARA_H+He.INC, PARA1_H+He.INC, PARA3_H+He.INC and adjust them by using the parameters from H+He_atoms.out and H+He_setf2.out.
  NOTICE: the values in the parameter files should be never 0!
• compile with ssh ait320 load pro2 H+He [P4/x64] (use P4 for 32bit machines, x64 for 64bit machines). The program /home/<user>/bimod/pro2_H+He.Linux_[P4/x64] is created
  

D.7.4 LTE2

• Edit H+He_lte.job with emacs. Set T=100000K, logg=7 and H=0.9 and He=0.1. The input files are H+He_lte.atoms2 and H+He_100000_lte.setf2.
• Start the job (an example is given below) and write its output into a file.
  nice +19 H+He_lte.job > H+He_lte.out The lte model is written to /home/<user>/model/H+He/lte. The file is called
  0100000_7.00_0.900_0.100.
  

D.7.5 PRO2

• Edit the H+He_nlte.job in emacs. Adjust the temperature, log gand the abundances of the used elements as before and check if the used files have the right name. The Input files are H+He.atoms2, H+He_100000.setf2 and the model from LTE2 is the input model.
• Start the job (an example is given in table below) with
  nice +19 H+He_nlte.job > H+He_nlte.out
• The output model is written to /home/$<$user$>$/model/H+He/.