tsuite/programs,
get them to work, create dotnet project for them, run them, zip them and have them for download
ETC

hizqso has comment   !Some hydrogenic lines mased, species was H  1, smallest tau was -1.23e+000, transition 10-0

cloudy web site say that chemistry totally updated, in paper say that results of chemistry network are more sensitive than ions 


====================================================================
HIGH PRIORITY

add to AGES - timescale of destruction of grains against sputtering.  "I've become worried 
about sputtering. From Draine & Salpeter (1978), the thermal sputtering 
lifetime is ~2x10^4 yr (1/nH) (a/0.01 microns) for gas with temperatures 
between 10^6 and 10^9 K."

look at top todo in nextdr - will speed up code

change linpack into clapack calls, then able to use std clapack link onto lib for each machine

add estimate of collisional contribution to Ha and Hbeta

Lya pumping of n=4 of HeII - recheck today vs Netzer, Elitzur & Ferland (1985, ApJ, 299, 752), 

remove uneeded var from EmLine - xIntensity could go - val not needed until zone totally converged */
cs2 should go

following not tested in tsuite: fabden, parseagn, dlaw, fileball, ratio, table abund, table table den
cext, cylinder, energy density, force te, hext, illuminate, neutrons

get rid of reference with ITRDIM, malloc when iterations are known 

use new he ct file for helium charge transfer

get rid of existing lowest 16 level feii, use verner instead, pequignot uses 46 level atom - do this

xspec integration

CO turned off and on in grid I ran in india

incorporate tlusty and wm-basic (pouldrach) model atmospheres

dusty torus convergence failures

post peter vh's small sample on .net web site, ask how to speed it up.  run same thing on P3

define differences between cloudy and three reflected spectra of compton thick matter in lex meeting book

CO - csupra now destroys CO but does not create anything in the matrix.  Should it be CO+?
what does high energy light do?  

is rayleigh scattering added as a excitation process for lya - think not but could be important

FeXIII 1216 line under Lya

tarter thesis has good description of diffuse fields in outward only

widla called in pressurechange never does special lya - always uses widlin - could this bew major problem with lya width?

Lya dest prob is nearly always the low-end limit of 1e-8 in kk.in

full test suite with eden solver new, conserv.in failed with excessive eden failures

blister.in fails with impossible bracket with set temp solver new

Fix so that very high temperature exits like very low temperature

get x-ray models strong lines to agree with exact wavelengths/energies

get new temp, eden solvers integrated into code

trace set eeden ot cofirm it still works

newdata has NII and NIII collision strength data - put into code

Stancil et al. O+ ct, A&A sup 140, 225, also Pequignot CT excitation of OI

Tayal 2000, ApJ 544, 575, has collision strengths for Fe10- 11, 13, expand
Fe10 to three level atom 

reinstate conpump routine that Hagai and I played with

Kaastra sent paper with lines from fluorescence, in cloudy/atomic data dir - put them in

way to enter J into continuum array, directly

longer input lines as per Fischer

boron-sequence, put in NeVI - Ca21 

NeVII 465A, Berrington et al, 1985, at dat nuc dat tab 33, 195

Fe15, 283A, better cs in Pradhan ADNDT 40, 335
Mg9 368, Berrington, 1994, ADNDT 57, many articles

hazy_coolcurve.c crashes if init te lowered to 3.4

==================================================================

MEDIUM PRIORITY

stop column density veil does not stop at exact column density 

get rid of NCELL and LIMSPC, they only occur in a few arrays

fileopac tests disabled as is while file opacities, should be reinstated?

get rid of typmatrx.h and dummy calls to non-existent matrix routines 

transfer si 10 606 and S XIII 488.4 then remove expion from source 

hmole, remove H3+ from matrix solver, do it after getting the H2 solution

putline does not need to call linadd

limits to how low an ionization could be done with H level atom were influenced by bug in float/double.  
reconsider all this with harder bounds

level3 does not consider the pumping component before giving up because gas too cool.  
As result it underestimates some lines.

no molecules command - change to set H2 formation rate to be very small

H2 - update photo cs to Yan et al. include H2 photoioniz as dest process 

molecular data base in Millar et al. A&AS 121, 139 - study paper, has lots of details, update cloudy to include all this

look at secondary ionization in Dalgarno et al. 1999, apj sup

put in cobe background tempeature into cloudy/hazy

add electron scattering escape to Lya only, since it is so important and such a problem in grain free case

add option for abundances to be on log scale where H is 12
also to not print abundances as logs

old oi photo into excited states - does it still exist?

new oi ct data

put in mewe line yields for all elements so that all fluor lines are predicted for axaf

keenan atomic data 
!G-bar cooling lines reached 7.29e+001% of the local cooling.  Line=Fe15  282A

go over vanHoof2do.txt and put in his suggested additions

move comments now in linedata.old into data file level1.dat, modify getdata so that it
	counts number of level1 lines and mallocs the right sized array.  skip over comments

make sure that every command in hazy is used at least one in a test case

heii bowen problem should be done correctly

====================================================================

PDR physics

[OI] production in flow from OH, Stoerzer & Hollenbach 1998

new molecular reactions in table 1 of Hollenbach takahashi & tielens apj 377, 192

check CO pops go to lte

TH85 says that charge transfer between C+ and S and Si are important

Maloney & Hollenback 1996, XDR, xrays into pdrs

