IR pumping or not?

Valid from revision 685e6bef.

One of the questions one can pose is whether ro-vibrational levels (CO, water, ...) are populated by collisions and/or are radiatively-pumped. There is an option in ProDiMo to test it. The idea is to scale artificially the collisional rates for a given species by a given factor provided by the user. In Parameter.in,

---------
2            ! NTestColl (Default = 0)
1e-4 CO
1e-5 o-H2O

NTestColl is the number of species for which the collisional rates should be scaled. The next two lines should follow immediately the line with NTestColl 1e-4 is the scaling factor for CO. 1e-5 is the scaling factor for o-H2O. The number of lines should be exactly the same that the value provided (here 2).

The way to test the effect of pumping (mostly IR pumping except for CO do-vibrational levels for which the effects of UV-pumping can be tested) is:

  1. Run your favorite model
  2. Re-run the same model but in Parameter.in:
.false.      ! solve_temp        : solve the heating/cooling balance?
.true.       ! restart
.true.       ! freeze_RT         : freeze radiative transfer results Td/Jv?
.false.      ! freeze_diskstruc  : freeze density/pressure structure?
.true.       ! freeze_Tgas       : freeze gas temperature?
.false.      ! freeze_chemistry  : freeze chemical concentrations?
---------
2            ! NTestColl (Default = 0)
1e-4 CO
1e-5 o-H2O

We do not recompute the gas-temperature via the balance between heating and cooling. You can now compare the line fluxes and line emission area between the two runs.