mirror of
https://asciireactor.com/otho/cloudy-agn.git
synced 2024-12-05 01:35:08 +00:00
374 lines
9.9 KiB
Plaintext
374 lines
9.9 KiB
Plaintext
////// Cloudy 17.00 4thdex grid using mehdipour ngc5548 SED
|
|
//
|
|
title mehdipour ngc5548 4thdex grid
|
|
//
|
|
////// Output continua and included lines
|
|
//
|
|
//save last incident continuum "incident_continuum"
|
|
//save line labels "line_labels"
|
|
//save continuum bins "continuum_bins"
|
|
//save diffuse continuum "diffuse_continuum"
|
|
save emitted continuum "emitted_continuum"
|
|
//
|
|
////// Report almost everything
|
|
//
|
|
print line column
|
|
print line inward
|
|
print line iso collapsed on
|
|
print line faint -3
|
|
print last
|
|
print line precision 6
|
|
Database print
|
|
//
|
|
////// Resolved energy states for atomic models
|
|
//
|
|
// Defaults: 10, 10, 6
|
|
//
|
|
//database H-like element hydrogen resolved levels 18
|
|
//database H-like element helium resolved levels 15
|
|
//database He-like element helium resolved levels 15
|
|
database H-like element hydrogen resolved levels 18 collapsed levels 30
|
|
database H-like element helium resolved levels 15 collapsed levels 30
|
|
database He-like element helium resolved levels 15 collapsed levels 40
|
|
//
|
|
////// low-metal abundance for AGN testing ~2017
|
|
// metals 0.5
|
|
// element scale factor carbon 2
|
|
// element scale factor nitrogen 2
|
|
//
|
|
////// Prevent cloudy from stopping prematurely
|
|
//
|
|
set nend 5000
|
|
//
|
|
////// Print diffuse continuum information in cloudy output
|
|
// and include additional spectral resolution around Lyman, Briggs, Etc.
|
|
//
|
|
// Continuum labels by default refer to the following fields:
|
|
// nFnu = Total Diffuse (in & out) + Reflected Inward
|
|
// nInu = not of interest, described in hazy 1
|
|
// InwT = Inward Diffuse + Reflected Inward
|
|
// InwC = Reflected Inward
|
|
//
|
|
// possible derived quantities added in Otho's program ~2017
|
|
// InwD = InwT - InwC = Inward Diffuse
|
|
// nFnD = nFnu - InwC = Total Diffuse
|
|
//
|
|
// "print continuum" redundant with set nFnu in c17
|
|
//
|
|
print continuum
|
|
//
|
|
// Jumps Added from Inspection, with mean wavelength as identifier
|
|
// Hydrogen Balmer, Paschen, Brackett
|
|
set nFnu add 3641 angstrom
|
|
set nFnu add 3653 angstrom
|
|
set nFnu add 8184 angstrom
|
|
set nFnu add 8211 angstrom
|
|
set nFnu add 8239 angstrom
|
|
set nFnu add 1456 angstrom
|
|
set nFnu add 1461 angstrom
|
|
// Possible Helium lines
|
|
set nFnu add 2046 angstrom
|
|
set nFnu add 2059 angstrom
|
|
set nFnu add 2592 angstrom
|
|
set nFnu add 2600 angstrom
|
|
set nFnu add 2609 angstrom
|
|
set nFnu add 1408 angstrom
|
|
set nFnu add 1413 angstrom
|
|
set nFnu add 1521 angstrom
|
|
set nFnu add 1826 angstrom
|
|
//
|
|
// Spread of continuum points
|
|
//
|
|
set nFnu add 450 angstrom
|
|
set nFnu add 650 angstrom
|
|
set nFnu add 875 angstrom
|
|
set nFnu add 800 angstrom
|
|
set nFnu add 840 angstrom
|
|
set nFnu add 850 angstrom
|
|
set nFnu add 880 angstrom
|
|
set nFnu add 900 angstrom
|
|
set nFnu add 910 angstrom
|
|
set nFnu add 913 angstrom
|
|
set nFnu add 950 angstrom
|
|
set nFnu add 960 angstrom
|
|
set nFnu add 1000 angstrom
|
|
set nFnu add 1050 angstrom
|
|
set nFnu add 1100 angstrom
|
|
set nFnu add 1150 angstrom
|
|
set nFnu add 1200 angstrom
|
|
set nFnu add 1240 angstrom
|
|
set nFnu add 1250 angstrom
|
|
set nFnu add 1300 angstrom
|
|
set nFnu add 1324 angstrom
|
|
set nFnu add 1350 angstrom
|
|
set nFnu add 1400 angstrom
|
|
set nFnu add 1410 angstrom
|
|
set nFnu add 1450 angstrom
|
|
set nFnu add 1500 angstrom
|
|
set nFnu add 1550 angstrom
|
|
set nFnu add 1600 angstrom
|
|
set nFnu add 1650 angstrom
|
|
set nFnu add 1700 angstrom
|
|
set nFnu add 1750 angstrom
|
|
set nFnu add 1800 angstrom
|
|
set nFnu add 1850 angstrom
|
|
set nFnu add 1900 angstrom
|
|
set nFnu add 1950 angstrom
|
|
set nFnu add 2000 angstrom
|
|
set nFnu add 2050 angstrom
|
|
set nFnu add 2100 angstrom
|
|
set nFnu add 2150 angstrom
|
|
set nFnu add 2200 angstrom
|
|
set nFnu add 2250 angstrom
|
|
set nFnu add 2300 angstrom
|
|
set nFnu add 2350 angstrom
|
|
set nFnu add 2400 angstrom
|
|
set nFnu add 2450 angstrom
|
|
set nFnu add 2500 angstrom
|
|
set nFnu add 2550 angstrom
|
|
set nFnu add 2600 angstrom
|
|
set nFnu add 2650 angstrom
|
|
set nFnu add 2700 angstrom
|
|
set nFnu add 2750 angstrom
|
|
set nFnu add 2800 angstrom
|
|
set nFnu add 2850 angstrom
|
|
set nFnu add 2900 angstrom
|
|
set nFnu add 2950 angstrom
|
|
set nFnu add 3000 angstrom
|
|
set nFnu add 3050 angstrom
|
|
set nFnu add 3100 angstrom
|
|
set nFnu add 3150 angstrom
|
|
set nFnu add 3200 angstrom
|
|
set nFnu add 3250 angstrom
|
|
set nFnu add 3300 angstrom
|
|
set nFnu add 3350 angstrom
|
|
set nFnu add 3400 angstrom
|
|
set nFnu add 3450 angstrom
|
|
set nFnu add 3500 angstrom
|
|
set nFnu add 3550 angstrom
|
|
set nFnu add 3600 angstrom
|
|
set nFnu add 3644 angstrom
|
|
set nFnu add 3645 angstrom
|
|
set nFnu add 3648 angstrom
|
|
set nFnu add 3650 angstrom
|
|
set nFnu add 3700 angstrom
|
|
set nFnu add 3750 angstrom
|
|
set nFnu add 3800 angstrom
|
|
set nFnu add 3850 angstrom
|
|
set nFnu add 3900 angstrom
|
|
set nFnu add 3950 angstrom
|
|
set nFnu add 4000 angstrom
|
|
set nFnu add 4050 angstrom
|
|
set nFnu add 4100 angstrom
|
|
set nFnu add 4150 angstrom
|
|
set nFnu add 4200 angstrom
|
|
set nFnu add 4250 angstrom
|
|
set nFnu add 4300 angstrom
|
|
set nFnu add 4350 angstrom
|
|
set nFnu add 4400 angstrom
|
|
set nFnu add 4450 angstrom
|
|
set nFnu add 4500 angstrom
|
|
set nFnu add 4550 angstrom
|
|
set nFnu add 4600 angstrom
|
|
set nFnu add 4650 angstrom
|
|
set nFnu add 4700 angstrom
|
|
set nFnu add 4750 angstrom
|
|
set nFnu add 4800 angstrom
|
|
set nFnu add 4850 angstrom
|
|
set nFnu add 4900 angstrom
|
|
set nFnu add 4950 angstrom
|
|
set nFnu add 5000 angstrom
|
|
set nFnu add 5050 angstrom
|
|
set nFnu add 5100 angstrom
|
|
set nFnu add 5150 angstrom
|
|
set nFnu add 5200 angstrom
|
|
set nFnu add 5250 angstrom
|
|
set nFnu add 5300 angstrom
|
|
set nFnu add 5350 angstrom
|
|
set nFnu add 5400 angstrom
|
|
set nFnu add 5450 angstrom
|
|
set nFnu add 5500 angstrom
|
|
set nFnu add 5550 angstrom
|
|
set nFnu add 5600 angstrom
|
|
set nFnu add 5650 angstrom
|
|
set nFnu add 5700 angstrom
|
|
set nFnu add 5750 angstrom
|
|
set nFnu add 5800 angstrom
|
|
set nFnu add 5850 angstrom
|
|
set nFnu add 5900 angstrom
|
|
set nFnu add 5950 angstrom
|
|
set nFnu add 6000 angstrom
|
|
set nFnu add 6050 angstrom
|
|
set nFnu add 6100 angstrom
|
|
set nFnu add 6150 angstrom
|
|
set nFnu add 6200 angstrom
|
|
set nFnu add 6250 angstrom
|
|
set nFnu add 6300 angstrom
|
|
set nFnu add 6350 angstrom
|
|
set nFnu add 6400 angstrom
|
|
set nFnu add 6450 angstrom
|
|
set nFnu add 6500 angstrom
|
|
set nFnu add 6550 angstrom
|
|
set nFnu add 6600 angstrom
|
|
set nFnu add 6650 angstrom
|
|
set nFnu add 6700 angstrom
|
|
set nFnu add 6750 angstrom
|
|
set nFnu add 6800 angstrom
|
|
set nFnu add 6850 angstrom
|
|
set nFnu add 6900 angstrom
|
|
set nFnu add 6950 angstrom
|
|
set nFnu add 7000 angstrom
|
|
set nFnu add 7050 angstrom
|
|
set nFnu add 7100 angstrom
|
|
set nFnu add 7150 angstrom
|
|
set nFnu add 7200 angstrom
|
|
set nFnu add 7250 angstrom
|
|
set nFnu add 7300 angstrom
|
|
set nFnu add 7350 angstrom
|
|
set nFnu add 7400 angstrom
|
|
set nFnu add 7450 angstrom
|
|
set nFnu add 7500 angstrom
|
|
set nFnu add 7550 angstrom
|
|
set nFnu add 7600 angstrom
|
|
set nFnu add 7650 angstrom
|
|
set nFnu add 7700 angstrom
|
|
set nFnu add 7750 angstrom
|
|
set nFnu add 7800 angstrom
|
|
set nFnu add 7850 angstrom
|
|
set nFnu add 7900 angstrom
|
|
set nFnu add 7950 angstrom
|
|
set nFnu add 8000 angstrom
|
|
set nFnu add 8050 angstrom
|
|
set nFnu add 8100 angstrom
|
|
set nFnu add 8150 angstrom
|
|
set nFnu add 8200 angstrom
|
|
set nFnu add 8204 angstrom
|
|
set nFnu add 8205 angstrom
|
|
set nFnu add 8207 angstrom
|
|
set nFnu add 8210 angstrom
|
|
set nFnu add 8250 angstrom
|
|
set nFnu add 8300 angstrom
|
|
set nFnu add 8350 angstrom
|
|
set nFnu add 8400 angstrom
|
|
set nFnu add 8450 angstrom
|
|
set nFnu add 8500 angstrom
|
|
set nFnu add 8550 angstrom
|
|
set nFnu add 8600 angstrom
|
|
set nFnu add 8650 angstrom
|
|
set nFnu add 8700 angstrom
|
|
set nFnu add 8750 angstrom
|
|
set nFnu add 8800 angstrom
|
|
set nFnu add 8850 angstrom
|
|
set nFnu add 8900 angstrom
|
|
set nFnu add 8950 angstrom
|
|
set nFnu add 9000 angstrom
|
|
set nFnu add 9050 angstrom
|
|
set nFnu add 9100 angstrom
|
|
set nFnu add 9150 angstrom
|
|
set nFnu add 9200 angstrom
|
|
set nFnu add 9250 angstrom
|
|
set nFnu add 9300 angstrom
|
|
set nFnu add 9350 angstrom
|
|
set nFnu add 9400 angstrom
|
|
set nFnu add 9450 angstrom
|
|
set nFnu add 9500 angstrom
|
|
set nFnu add 9550 angstrom
|
|
set nFnu add 9600 angstrom
|
|
set nFnu add 9650 angstrom
|
|
set nFnu add 9700 angstrom
|
|
set nFnu add 9750 angstrom
|
|
set nFnu add 9800 angstrom
|
|
set nFnu add 9850 angstrom
|
|
set nFnu add 9900 angstrom
|
|
set nFnu add 9950 angstrom
|
|
set nFnu add 10000 angstrom
|
|
set nFnu add 10050 angstrom
|
|
set nFnu add 10100 angstrom
|
|
set nFnu add 10150 angstrom
|
|
set nFnu add 10200 angstrom
|
|
set nFnu add 10250 angstrom
|
|
set nFnu add 10300 angstrom
|
|
set nFnu add 10350 angstrom
|
|
set nFnu add 10400 angstrom
|
|
set nFnu add 10450 angstrom
|
|
set nFnu add 10500 angstrom
|
|
set nFnu add 10550 angstrom
|
|
set nFnu add 10600 angstrom
|
|
set nFnu add 10650 angstrom
|
|
set nFnu add 10700 angstrom
|
|
set nFnu add 10750 angstrom
|
|
set nFnu add 10800 angstrom
|
|
set nFnu add 10850 angstrom
|
|
set nFnu add 10900 angstrom
|
|
set nFnu add 10950 angstrom
|
|
set nFnu add 11000 angstrom
|
|
set nFnu add 11050 angstrom
|
|
set nFnu add 11100 angstrom
|
|
set nFnu add 11150 angstrom
|
|
set nFnu add 11200 angstrom
|
|
set nFnu add 11250 angstrom
|
|
set nFnu add 11300 angstrom
|
|
set nFnu add 11350 angstrom
|
|
set nFnu add 11400 angstrom
|
|
set nFnu add 11450 angstrom
|
|
set nFnu add 11500 angstrom
|
|
set nFnu add 11550 angstrom
|
|
set nFnu add 11600 angstrom
|
|
set nFnu add 11650 angstrom
|
|
set nFnu add 11700 angstrom
|
|
set nFnu add 11750 angstrom
|
|
set nFnu add 11800 angstrom
|
|
set nFnu add 11850 angstrom
|
|
set nFnu add 11900 angstrom
|
|
set nFnu add 11950 angstrom
|
|
set nFnu add 12000 angstrom
|
|
set nFnu add 12050 angstrom
|
|
set nFnu add 12100 angstrom
|
|
set nFnu add 12150 angstrom
|
|
set nFnu add 12200 angstrom
|
|
set nFnu add 12250 angstrom
|
|
set nFnu add 12300 angstrom
|
|
set nFnu add 12350 angstrom
|
|
set nFnu add 12400 angstrom
|
|
set nFnu add 12450 angstrom
|
|
set nFnu add 12500 angstrom
|
|
set nFnu add 12550 angstrom
|
|
set nFnu add 12600 angstrom
|
|
set nFnu add 12650 angstrom
|
|
set nFnu add 12700 angstrom
|
|
set nFnu add 12750 angstrom
|
|
set nFnu add 12800 angstrom
|
|
set nFnu add 12850 angstrom
|
|
set nFnu add 12900 angstrom
|
|
set nFnu add 12950 angstrom
|
|
set nFnu add 13000 angstrom
|
|
set nFnu add 800000 angstrom
|
|
set nFnu add 1000000 angstrom
|
|
//
|
|
////// Print equivalent widths
|
|
//
|
|
normalize to "Inci" 1215 scale to 1215
|
|
//
|
|
////// Background radiation
|
|
//
|
|
background z = 0.0172
|
|
//
|
|
////// Hydrogen density of cloud is varied
|
|
//
|
|
hden 11 vary
|
|
grid range from 7 to 14 with .25 dex
|
|
//
|
|
////// mehdipour et al. 2013 NGC 5548 SED
|
|
//
|
|
table SED "mehdipour2013.tab"
|
|
//
|
|
////// Flux is varied
|
|
//
|
|
//phi(h) 23
|
|
phi(h) 23 vary
|
|
grid range from 17 to 24 with .25 dex
|
|
//
|
|
////// Stopping condition and run
|
|
//
|
|
stop total column density 23
|
|
iterate to convergence max=15 |