Skip to content

Commit

Permalink
Frank's fixes
Browse files Browse the repository at this point in the history
  • Loading branch information
@ivezic
ivezic committed Dec 30, 2017
1 parent a470075 commit 451f84a
Show file tree
Hide file tree
Showing 418 changed files with 58,918 additions and 0 deletions.
64 changes: 64 additions & 0 deletions dusty.mas
View file
@@ -0,0 +1,64 @@
% ------- Master Input File for DUSTY -------------------------------

% Spherical dust distribution with
% constant density profile
% heated by external radiation with a Black Body of 5000K
% for composite dust grain 70% silicates 30% carbon
% for 3 optical depth (30,40,50)

examples/sphere1/sphere1.inp

% Example slab geometry
% directional input of power law spectrum
% left side temperature of 800 K

examples/slab1/slab1.inp

% The Iveciz 97 Benchmark
% Heating source is a 2500K Black Body
% sphercical symmetry
% constant density
% for 4 optical depth (1,10,100,1000)

examples/ivezic97/ivezic97.inp

% dust distriubtuion 1/r^2 profile
% for 4 optical depth (1,10,100,1000)

examples/ivezic97/ivezic97_2.inp

% Spherical dust distribution with
% 1/r density profile
% heated by external raditition with a temperature of 1000K
% no scattering

examples/sphere_BB/sphere_BB.inp


% Slab in a Black Body radiation field
% with temperature of 1000K
% for 4 optical depth (1,10,100,1000)

examples/slab_BB/slab_BB.inp

% Spherical dust distribution with
% constant density profile
% heated by external radiation with a temperature of 5000K
% for composite dust grain 70% silicates 30% carbon
% for 2 optical depth (10,100)

examples/external/external.inp

% Example slab geometry
% isotropic input on both sides
% ISM-OHM dust grain properties

examples/slab2/slab2.inp

% Example slab geometry
% only left side illumination (flux input)
% ISM-OHM dust grain properties

examples/slab3/slab3.inp

% -------------------------------
Binary file added release/dusty.tar
View file
Binary file not shown.
2 changes: 2 additions & 0 deletions release/dusty/Makefile
View file
@@ -0,0 +1,2 @@
all:
gfortran -O3 -lgomp -fopenmp -o dusty dusty.f90
1 change: 1 addition & 0 deletions release/dusty/compile.bat
View file
@@ -0,0 +1 @@
gfortran -O3 -lgomp -fopenmp -o dusty.exe dusty.f90
109 changes: 109 additions & 0 deletions release/dusty/data/AGN-p0.5.dat
View file
@@ -0,0 +1,109 @@
Example of AGN SED, flat in 0.01-0.1mic,
lam*f_lam ~ lam^{-0.5} and RJ-tail beyond 1 micron
lambda lam*flam
>
1.000E-02 2.785E-01
1.500E-02 2.785E-01
2.000E-02 2.785E-01
3.000E-02 2.785E-01
4.000E-02 2.785E-01
5.000E-02 2.785E-01
6.000E-02 2.785E-01
8.000E-02 2.785E-01
1.000E-01 2.785E-01
1.200E-01 2.541E-01
1.500E-01 2.273E-01
2.000E-01 1.969E-01
2.500E-01 1.761E-01
3.000E-01 1.608E-01
3.600E-01 1.467E-01
4.400E-01 1.328E-01
5.500E-01 1.187E-01
7.000E-01 1.052E-01
8.500E-01 9.549E-02
1.000E+00 8.804E-02
1.150E+00 5.789E-02
1.300E+00 4.008E-02
1.700E+00 1.792E-02
2.000E+00 1.101E-02
2.200E+00 8.268E-03
2.700E+00 4.473E-03
3.000E+00 3.261E-03
3.500E+00 2.054E-03
4.000E+00 1.376E-03
4.500E+00 9.662E-04
5.000E+00 7.044E-04
5.500E+00 5.292E-04
6.000E+00 4.076E-04
6.500E+00 3.206E-04
7.000E+00 2.567E-04
7.500E+00 2.087E-04
8.000E+00 1.720E-04
8.500E+00 1.434E-04
9.000E+00 1.208E-04
9.400E+00 1.060E-04
9.550E+00 1.011E-04
9.700E+00 9.647E-05
9.850E+00 9.212E-05
1.000E+01 8.804E-05
1.050E+01 7.606E-05
1.100E+01 6.615E-05
1.130E+01 6.102E-05
1.160E+01 5.641E-05
1.200E+01 5.096E-05
1.250E+01 4.508E-05
1.300E+01 4.008E-05
1.350E+01 3.579E-05
1.400E+01 3.209E-05
1.450E+01 2.888E-05
1.500E+01 2.609E-05
1.600E+01 2.150E-05
1.700E+01 1.792E-05
1.800E+01 1.510E-05
1.900E+01 1.284E-05
2.000E+01 1.101E-05
2.200E+01 8.268E-06
2.300E+01 7.237E-06
2.400E+01 6.369E-06
2.500E+01 5.635E-06
2.600E+01 5.010E-06
2.700E+01 4.473E-06
2.800E+01 4.011E-06
3.000E+01 3.261E-06
3.500E+01 2.054E-06
4.000E+01 1.376E-06
4.500E+01 9.662E-07
5.000E+01 7.044E-07
5.500E+01 5.292E-07
6.000E+01 4.076E-07
6.500E+01 3.206E-07
7.000E+01 2.567E-07
7.500E+01 2.087E-07
8.000E+01 1.720E-07
8.500E+01 1.434E-07
9.000E+01 1.208E-07
9.500E+01 1.027E-07
1.000E+02 8.804E-08
1.050E+02 7.606E-08
1.100E+02 6.615E-08
1.200E+02 5.096E-08
1.300E+02 4.008E-08
1.400E+02 3.209E-08
1.500E+02 2.609E-08
2.000E+02 1.101E-08
2.500E+02 5.635E-09
3.000E+02 3.261E-09
4.000E+02 1.376E-09
5.000E+02 7.044E-10
7.000E+02 2.567E-10
1.000E+03 8.804E-11
1.500E+03 2.609E-11
2.250E+03 7.729E-12
3.350E+03 2.342E-12
5.000E+03 7.044E-13
7.400E+03 2.173E-13
1.100E+04 6.615E-14
1.300E+04 4.008E-14
1.600E+04 2.150E-14
2.400E+04 6.369E-15
3.600E+04 1.887E-15
27 changes: 27 additions & 0 deletions release/dusty/data/Dprofile-1.5.dat
View file
@@ -0,0 +1,27 @@
This is an example for dust density given in a file with
rho = 1/r**1.5. Make sure that this header has three lines.
r(i) rho_dust(i)
>
1.0000 1.0000
1.0200 0.9707
1.0400 0.9429
1.0600 0.9163
1.0800 0.8910
1.1000 0.8668
1.5000 0.5443
2.0000 0.3536
2.5000 0.2530
3.0000 0.1925
3.5000 0.1527
4.0000 0.1250
5.0000 0.0894
7.5000 0.0487
10.0000 0.0316
15.0000 0.0172
20.0000 0.0112
25.0000 8.0000e-3
50.0000 2.8284e-3
75.0000 1.5396e-3
90.0000 1.1712e-3
100.0000 1.0000e-3
101.0000 9.8519e-4
28 changes: 28 additions & 0 deletions release/dusty/data/GrayAbs.dat
View file
@@ -0,0 +1,28 @@
Gray Absorption Cross sections
For tests
lambda <s_abs>/<V> <s_sca>/<V>
>
1.00000E-02 1.00000E+00 0.00000E+00
2.00000E-02 1.00000E+00 0.00000E+00
3.00000E-02 1.00000E+00 0.00000E+00
1.00000E-01 1.00000E+00 0.00000E+00
2.00000E-01 1.00000E+00 0.00000E+00
5.00000E-01 1.00000E+00 0.00000E+00
7.00000E-01 1.00000E+00 0.00000E+00
1.00000E+00 1.00000E+00 0.00000E+00
2.00000E+00 1.00000E+00 0.00000E+00
5.00000E+00 1.00000E+00 0.00000E+00
7.00000E+00 1.00000E+00 0.00000E+00
1.00000E+01 1.00000E+00 0.00000E+00
2.00000E+01 1.00000E+00 0.00000E+00
5.00000E+01 1.00000E+00 0.00000E+00
1.00000E+02 1.00000E+00 0.00000E+00
2.00000E+02 1.00000E+00 0.00000E+00
5.00000E+02 1.00000E+00 0.00000E+00
7.00000E+02 1.00000E+00 0.00000E+00
1.00000E+03 1.00000E+00 0.00000E+00
2.00000E+03 1.00000E+00 0.00000E+00
5.00000E+03 1.00000E+00 0.00000E+00
1.30000E+04 1.00000E+00 0.00000E+00
2.00000E+04 1.00000E+00 0.00000E+00
3.60000E+04 1.00000E+00 0.00000E+00
129 changes: 129 additions & 0 deletions release/dusty/data/ISM-OHM.dat
View file
@@ -0,0 +1,129 @@
Standard MRN mixture (a_min=0.005, amax=0.25)
Silw:0.00 Silc:0.53 SiDL:0.00 amC:0.00 SiC:0.00 Grph:0.47 CrSi:0.00
lambda <s_abs>/<V> <s_sca>/<V>
>
1.00000E-02 1.37692E+01 1.91565E+01
1.50000E-02 1.45889E+01 2.35828E+01
2.00000E-02 1.49011E+01 2.51513E+01
3.00000E-02 1.52609E+01 2.39555E+01
4.00000E-02 1.44957E+01 2.22682E+01
5.00000E-02 1.53590E+01 2.15468E+01
6.00000E-02 1.70872E+01 2.07716E+01
8.00000E-02 2.17842E+01 2.06453E+01
1.00000E-01 1.80401E+01 1.86583E+01
1.20000E-01 1.28086E+01 1.55895E+01
1.50000E-01 9.76698E+00 1.13639E+01
2.00000E-01 1.42963E+01 1.09208E+01
2.50000E-01 1.13794E+01 1.11748E+01
3.00000E-01 9.04540E+00 9.06969E+00
3.60000E-01 7.63721E+00 7.74474E+00
4.40000E-01 6.21257E+00 6.72963E+00
5.50000E-01 4.69712E+00 5.61756E+00
7.00000E-01 3.38964E+00 4.28294E+00
8.50000E-01 2.59384E+00 3.11658E+00
1.00000E+00 2.13997E+00 2.18339E+00
1.15000E+00 1.85926E+00 1.54104E+00
1.30000E+00 1.63576E+00 1.10403E+00
1.70000E+00 1.15204E+00 4.95375E-01
2.00000E+00 9.04046E-01 2.75221E-01
2.20000E+00 7.66565E-01 1.89316E-01
2.70000E+00 5.32395E-01 8.13015E-02
3.00000E+00 4.42378E-01 5.23051E-02
3.50000E+00 3.35433E-01 2.74362E-02
4.00000E+00 2.67999E-01 1.57305E-02
4.50000E+00 2.18639E-01 9.65074E-03
5.00000E+00 1.84746E-01 6.23549E-03
5.50000E+00 1.59698E-01 4.19896E-03
6.00000E+00 1.41654E-01 2.92149E-03
6.50000E+00 1.29064E-01 2.08788E-03
7.00000E+00 1.20243E-01 1.52555E-03
7.50000E+00 1.14356E-01 1.13335E-03
8.00000E+00 1.15293E-01 8.49941E-04
8.20000E+00 1.19055E-01 7.58266E-04
8.40000E+00 1.27507E-01 6.77099E-04
8.60000E+00 1.41320E-01 6.03974E-04
8.80000E+00 1.81709E-01 5.46637E-04
9.00000E+00 2.20006E-01 5.03365E-04
9.20000E+00 2.67792E-01 4.68106E-04
9.30000E+00 2.93935E-01 4.54998E-04
9.40000E+00 3.26493E-01 4.47036E-04
9.45000E+00 3.39166E-01 4.44094E-04
9.50000E+00 3.52246E-01 4.41765E-04
9.55000E+00 3.59885E-01 4.37921E-04
9.60000E+00 3.67566E-01 4.34383E-04
9.65000E+00 3.75499E-01 4.30190E-04
9.70000E+00 3.83548E-01 4.26315E-04
9.75000E+00 3.95550E-01 4.24356E-04
9.80000E+00 4.07848E-01 4.22929E-04
9.85000E+00 4.19681E-01 4.23722E-04
9.90000E+00 4.31504E-01 4.25089E-04
9.95000E+00 4.34045E-01 4.24710E-04
1.00000E+01 4.36102E-01 4.24500E-04
1.02000E+01 4.23210E-01 4.11116E-04
1.05000E+01 3.87015E-01 3.83367E-04
1.10000E+01 3.13667E-01 3.23431E-04
1.13000E+01 2.85787E-01 2.94089E-04
1.16000E+01 2.49705E-01 2.70636E-04
1.20000E+01 1.88450E-01 2.28611E-04
1.25000E+01 1.52398E-01 1.87402E-04
1.30000E+01 1.28633E-01 1.53619E-04
1.35000E+01 1.23224E-01 1.27927E-04
1.40000E+01 1.21880E-01 1.08515E-04
1.45000E+01 1.17825E-01 9.29300E-05
1.50000E+01 1.15639E-01 7.83900E-05
1.55000E+01 1.31753E-01 6.75835E-05
1.60000E+01 1.51647E-01 5.96427E-05
1.65000E+01 1.74045E-01 5.39487E-05
1.70000E+01 1.95426E-01 5.00228E-05
1.75000E+01 2.10369E-01 4.72440E-05
1.80000E+01 2.14152E-01 4.47831E-05
1.85000E+01 2.07661E-01 4.20871E-05
1.90000E+01 1.93983E-01 3.89542E-05
1.95000E+01 1.78618E-01 3.56456E-05
2.00000E+01 1.64394E-01 3.23953E-05
2.10000E+01 1.43682E-01 2.65795E-05
2.20000E+01 1.31951E-01 2.20079E-05
2.30000E+01 1.24228E-01 1.84892E-05
2.40000E+01 1.18967E-01 1.56558E-05
2.50000E+01 1.15890E-01 1.33681E-05
2.60000E+01 1.13977E-01 1.15446E-05
2.70000E+01 1.13469E-01 1.00896E-05
2.80000E+01 1.15037E-01 8.90090E-06
3.00000E+01 1.17369E-01 7.12452E-06
3.50000E+01 9.87381E-02 4.35166E-06
4.00000E+01 6.89579E-02 2.63347E-06
4.50000E+01 4.75022E-02 1.63016E-06
5.00000E+01 3.42660E-02 1.05469E-06
5.50000E+01 2.61953E-02 7.11406E-07
6.00000E+01 2.09203E-02 4.97131E-07
6.50000E+01 1.71519E-02 3.58027E-07
7.00000E+01 1.43341E-02 2.64537E-07
7.50000E+01 1.21488E-02 1.99802E-07
8.00000E+01 1.04382E-02 1.53771E-07
8.50000E+01 9.07464E-03 1.20300E-07
9.00000E+01 7.98339E-03 9.54357E-08
9.50000E+01 7.08557E-03 7.66889E-08
1.00000E+02 6.33709E-03 6.23333E-08
1.05000E+02 5.71597E-03 5.11878E-08
1.10000E+02 5.19395E-03 4.24276E-08
1.20000E+02 4.38606E-03 2.98808E-08
1.30000E+02 3.81568E-03 2.16726E-08
1.40000E+02 3.33885E-03 1.61221E-08
1.50000E+02 2.88017E-03 1.22464E-08
2.00000E+02 1.50716E-03 3.86439E-09
2.50000E+02 9.48225E-04 1.57971E-09
3.00000E+02 6.55054E-04 7.61011E-10
4.00000E+02 3.67628E-04 2.40545E-10
5.00000E+02 2.35461E-04 9.84787E-11
7.00000E+02 1.22200E-04 2.56262E-11
1.00000E+03 5.91449E-05 6.15164E-12
1.50000E+03 2.63333E-05 1.21501E-12
2.25000E+03 1.17136E-05 2.39988E-13
3.35000E+03 5.28074E-06 4.88350E-14
5.00000E+03 2.37255E-06 9.84068E-15
7.40000E+03 1.09182E-06 2.05107E-15
1.10000E+04 5.14910E-07 4.20084E-16
1.30000E+04 4.06692E-07 2.15344E-16
1.60000E+04 3.06126E-07 9.38483E-17
2.40000E+04 1.80672E-07 1.85379E-17
3.60000E+04 1.10043E-07 3.66181E-18

0 comments on commit 451f84a

Please sign in to comment.