title He atom in LTE at high densities 
c
c commands controlling continuum =========
blackbody 50000 STE
c
c commands for density & abundances =========
hden 18
* set eden 18
c a high he abundance so that it dominates 
element helium abundance 1
init file "hheonly.ini"
c
c commands controlling geometry  =========
stop zone 1
c
c other commands for details     =========
atom he-like resolved levels 7
atom he-like collapsed levels 30
c atom h2 matrix all
iterate
c
c commands controlling output    =========
print he-like departure helium 
print he-like population helium 
print h-like departure helium 
print h-like population helium 
print h-like departure hydrogen 
normalise to "He 1" 4471
print line faint -2
save overview "limit_lte_he1_ste.ovr"
save performance "limit_lte_he1_ste.per"
save heating "limit_lte_he1_ste.het" 
save cooling "limit_lte_he1_ste.col" 
save dr "limit_lte_he1_ste.dr" 
save dominant rates "limit_lte_he1_ste.h2" "H2"
c
c commands giving the monitors    =========
monitor he-like excited helium departures 1 
monitor temperature helium 1 50000
c limit_lte_he1_ste.in 
c class limit 
c ========================================
c 

test whether a gas dominated by He goes to LTE in high-density limit.
The level populations should be in LTE, the departure coefficients
should be unity, and the temperature equal to the BB temp.