py79_pl_loop

This routine carries out a series of thin shell python simulations.

The wind mdot is set to produce a hydrogen density of 1e7 - to change this one had to change the wind_mdot parameter. The loop is carried out over the 2-10kev luminosity of the central source. A luminosity of 1e21 give a mainly neutral plasma (IP=5e-10) whilst 1e37 is almost totally ionized (IP=5e6). The code strips out ion densities and constructs data files with the cloudy type ionization parameter as the first line. If one wants a constant temperature file, then one needs to hack python to stop the temperature changing We currently write out H,He,C,N,O and Fe. We also output heating and cooling mechanisms.

Command line usage (if any):

usage: python_pl_loop optional file suffix

Description:

Primary routines:

Notes:

History:

081214 nsh Commented