| 1 | title check that grains equilibriate at correct temp in ste limit |
|---|
| 2 | c this also tests that the code works when H He ionization is ZERO |
|---|
| 3 | c |
|---|
| 4 | c commands controlling continuum ========= |
|---|
| 5 | trace hmole 0 2 |
|---|
| 6 | atom h2 trace full |
|---|
| 7 | c print dominant rates "h2" |
|---|
| 8 | c black 3000K |
|---|
| 9 | c ionization parameter -15 |
|---|
| 10 | black 3 lte |
|---|
| 11 | c |
|---|
| 12 | c commands for density & abundances ========= |
|---|
| 13 | hden 5 |
|---|
| 14 | eden 5 |
|---|
| 15 | init file "hheonly.ini" |
|---|
| 16 | c need to normalize to some line, oxygen will provide some |
|---|
| 17 | c element oxygen on |
|---|
| 18 | c ism grains with high dust to gas ratio |
|---|
| 19 | grain abund 1 |
|---|
| 20 | c |
|---|
| 21 | c commands controlling geometry ========= |
|---|
| 22 | stop zone 1 |
|---|
| 23 | set dr 0 |
|---|
| 24 | c |
|---|
| 25 | c other commands for details ========= |
|---|
| 26 | iterate |
|---|
| 27 | atom h2 |
|---|
| 28 | c the equilibrium temp would not be correct since small H_2 molecule |
|---|
| 29 | c does not go to ste, being an equivalent two-level system |
|---|
| 30 | c constant temperature 3 |
|---|
| 31 | c |
|---|
| 32 | c commands controlling output ========= |
|---|
| 33 | normalize to "O 1" 63.17m |
|---|
| 34 | punch asserts last "grains_lte.asr" |
|---|
| 35 | c |
|---|
| 36 | c commands giving the asserts ========= |
|---|
| 37 | assert temperature H2 3 error 0.01 |
|---|
| 38 | assert grain temperature 1 3 error 0.005 |
|---|
| 39 | assert grain temperature 2 3 error 0.005 |
|---|
| 40 | assert grain temperature 3 3 error 0.005 |
|---|
| 41 | assert grain temperature 4 3 error 0.005 |
|---|
| 42 | assert grain temperature 5 3 error 0.005 |
|---|
| 43 | assert grain temperature 6 3 error 0.005 |
|---|
| 44 | assert grain temperature 7 3 error 0.005 |
|---|
| 45 | assert grain temperature 8 3 error 0.005 |
|---|
| 46 | assert grain temperature 9 3 error 0.005 |
|---|
| 47 | assert grain temperature 10 3 error 0.005 |
|---|
| 48 | assert grain temperature 11 3 error 0.005 |
|---|
| 49 | assert grain temperature 12 3 error 0.005 |
|---|
| 50 | assert grain temperature 13 3 error 0.005 |
|---|
| 51 | assert grain temperature 14 3 error 0.005 |
|---|
| 52 | assert grain temperature 15 3 error 0.005 |
|---|
| 53 | assert grain temperature 16 3 error 0.005 |
|---|
| 54 | assert grain temperature 17 3 error 0.005 |
|---|
| 55 | assert grain temperature 18 3 error 0.005 |
|---|
| 56 | assert grain temperature 19 3 error 0.005 |
|---|
| 57 | assert grain temperature 20 3 error 0.005 |
|---|
| 58 | c grains_lte.in |
|---|
| 59 | c class limit |
|---|
| 60 | c ======================================== |
|---|
| 61 | c |
|---|
| 62 | |
|---|
| 63 | This test irradiates a set of grains with a true blackbody in |
|---|
| 64 | strict thermodynamic equilibrium. We expect the grains |
|---|
| 65 | (and everything else) to equilibriate at the blackbody temperature. |
|---|
| 66 | The gas temperature is forced to the radiation temperature because |
|---|
| 67 | the current molecule network (based on ISM approximations) does not |
|---|
| 68 | go to LTE in the high-radiation density limit. The calculation |
|---|
| 69 | asserts that all grain temperatures are very close to the radiation |
|---|
| 70 | temperature. |
|---|
| 71 | |
|---|
| 72 | |
|---|
| 73 | |
|---|
