Page Menu
Home
c4science
Search
Configure Global Search
Log In
Files
F85249118
plot.py
No One
Temporary
Actions
Download File
Edit File
Delete File
View Transforms
Subscribe
Mute Notifications
Award Token
Subscribers
None
File Metadata
Details
File Info
Storage
Attached
Created
Fri, Sep 27, 19:19
Size
1 KB
Mime Type
text/x-python
Expires
Sun, Sep 29, 19:19 (1 d, 21 h)
Engine
blob
Format
Raw Data
Handle
21146575
Attached To
rGEAR Gear
plot.py
View Options
#!/usr/bin/env python
from
numpy
import
*
import
sys
import
coolinglib
##############################################################################
# init object
##############################################################################
ct
=
coolinglib
.
CTable
(
"HM01/IonRates_Hydrogen.HM01.hdf5"
)
ion
=
0
# recombinaison rate
Ts
=
10
**
arange
(
1
,
9
,
0.1
)
As
=
zeros
(
len
(
Ts
))
Bs
=
zeros
(
len
(
Ts
))
Cools
=
zeros
(
len
(
Ts
))
for
i
,
T
in
enumerate
(
Ts
):
As
[
i
]
=
ct
.
RecombinasionRateCoefficient
(
ion
,
T
)
for
i
,
T
in
enumerate
(
Ts
):
Bs
[
i
]
=
ct
.
CollisionalIonisationRateCoefficient
(
ion
,
T
)
for
i
,
T
in
enumerate
(
Ts
):
Cools
[
i
]
=
ct
.
CoolingEfficiency
(
ion
,
T
)
# photo heating efficiency
zs
=
arange
(
0
,
20
,
0.01
)
Hs
=
zeros
(
len
(
zs
))
Cs
=
zeros
(
len
(
zs
))
for
i
,
z
in
enumerate
(
zs
):
Hs
[
i
]
=
ct
.
PhotoHeatingEfficiency
(
ion
,
z
)
for
i
,
z
in
enumerate
(
zs
):
Cs
[
i
]
=
ct
.
PhotoIonisationRate
(
ion
,
z
)
##############################################################################
# some plots
##############################################################################
import
Ptools
as
pt
#pt.scatter(ct.Temperatures,ct.Alpha[ion])
#pt.plot(Ts,As)
#pt.scatter(ct.Temperatures,ct.BetaColl[ion])
#pt.plot(Ts,Bs)
#pt.scatter(ct.Temperatures,ct.Cool[ion])
#pt.plot(Ts,Cools)
#pt.scatter(ct.Redshifts,ct.Heat[ion])
#pt.plot(zs,Hs)
pt
.
scatter
(
ct
.
Redshifts
,
ct
.
Heat
[
ion
])
pt
.
plot
(
zs
,
Hs
)
pt
.
show
()
Event Timeline
Log In to Comment