updating scripts for cloudy input

This commit is contained in:
caes 2017-06-08 03:59:36 -04:00
parent 36d6294905
commit 18dd4cda3f

View File

@ -1,324 +0,0 @@
////// CLOUDY 13.03 grid run as hypothetical values for AGN observation
//
title ngc 5548 BLR hden_vs_phi .25 dex
//
// single run copied from the standard mehdipour grid
//
////// Output continua
//
save last incident continuum "incident_continuum"
//
////// Formatting
//
//print line column
print line inward
print line iso collapsed
print line faint -2
print last
print line precision 6
//
////// Resolved energy states for atomic models
//
// Defaults: 10, 10, 6
//
atom H-like element hydrogen resolved levels 10
atom H-like element helium resolved levels 10
atom He-like element helium resolved levels 6
//
////// Prevent cloudy from stopping prematurely
//
set nend 5000
//
////// Print diffuse continuum information in cloudy output
// and include additional spectral resolution around Lyman, Briggs, Etc.
//
// diffuse_inward can be calculated with (InwT - InwC = InwDiffuse)
//
set nFnu add 800 angstrom
set nFnu add 840 angstrom
set nFnu add 850 angstrom
set nFnu add 900 angstrom
set nFnu add 950 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 1250 angstrom
set nFnu add 1300 angstrom
set nFnu add 1350 angstrom
set nFnu add 1400 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 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 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 880 angstrom
set nFnu add 910 angstrom
set nFnu add 913 angstrom
set nFnu add 960 angstrom
set nFnu add 1240 angstrom
set nFnu add 1324 angstrom
set nFnu add 1410 angstrom
set nFnu add 3644 angstrom
set nFnu add 3645 angstrom
set nFnu add 3648 angstrom
set nFnu add 8204 angstrom
set nFnu add 8207 angstrom
set nFnu add 8210 angstrom
//
////// Print equivalent widths
//
//normalize to 1215 scale to 1215
//
////// Background radiation
//
background z = 0.0172
//
////// Hydrogen density of cloud is varied
//
hden 11
//hden 11 vary
//grid range from 7 to 14 with .25 dex
//
////// Central AGN source shape generated by external program
//
table SED "reference/magdziarz_sed.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=40