Page Menu
Home
c4science
Search
Configure Global Search
Log In
Files
F90029895
gplot_r-P
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
Mon, Oct 28, 15:38
Size
3 KB
Mime Type
text/x-python
Expires
Wed, Oct 30, 15:38 (1 d, 23 h)
Engine
blob
Format
Raw Data
Handle
21996383
Attached To
rGTOOLS Gtools
gplot_r-P
View Options
#!/home/revaz/local/bin/python
'''
Plot temperature of the model as a function of radius
assuming an ideal gas
Yves Revaz
Thu Feb 23 15:00:11 CET 2006
'''
from
numarray
import
*
from
Nbody
import
*
import
SM
import
string
import
sys
import
os
from
libjeans
import
*
from
Nbody.libutil
import
histogram
from
optparse
import
OptionParser
from
Gtools
import
*
from
Gtools
import
vanderwaals
as
vw
def
parse_options
():
usage
=
"usage: %prog [options] file"
parser
=
OptionParser
(
usage
=
usage
)
parser
=
add_postscript_options
(
parser
)
parser
=
add_color_options
(
parser
)
parser
=
add_limits_options
(
parser
)
parser
=
add_units_options
(
parser
)
parser
=
add_log_options
(
parser
)
parser
=
add_histogram_options
(
parser
)
parser
=
add_reduc_options
(
parser
)
parser
.
add_option
(
"-t"
,
action
=
"store"
,
dest
=
"ftype"
,
type
=
"string"
,
default
=
None
,
help
=
"type of the file"
,
metavar
=
" TYPE"
)
parser
.
add_option
(
"--cgs"
,
action
=
"store_true"
,
dest
=
"cgs"
,
default
=
0
,
help
=
"invert into cgs units"
)
# parser.add_option("--gamma",
# action="store",
# dest="gamma",
# type="float",
# default=5/3.,
# help="adiabatic index")
# parser.add_option("--mu",
# action="store",
# dest="mu",
# type="float",
# default=2,
# help="mean molecular mass")
parser
.
add_option
(
"--av"
,
action
=
"store"
,
dest
=
"av"
,
type
=
"float"
,
default
=
0
,
help
=
"Av consante (in cgs)"
)
parser
.
add_option
(
"--bv"
,
action
=
"store"
,
dest
=
"bv"
,
type
=
"float"
,
default
=
0
,
help
=
"Bv consante (in cgs)"
)
(
options
,
args
)
=
parser
.
parse_args
()
if
options
.
colors
!=
None
:
exec
(
"options.colors = array([
%s
])"
%
(
options
.
colors
))
if
len
(
args
)
==
0
:
print
"you must specify a filename"
sys
.
exit
(
0
)
files
=
args
return
files
,
options
#############################
# graph
#############################
# get options
files
,
options
=
parse_options
()
ps
=
options
.
ps
col
=
options
.
colors
xmin
=
options
.
xmin
xmax
=
options
.
xmax
ymin
=
options
.
ymin
ymax
=
options
.
ymax
log
=
options
.
log
histogram
=
options
.
histogram
reduc
=
options
.
reduc
cgs
=
options
.
cgs
ftype
=
options
.
ftype
#gamma = options.gamma
#mu = options.mu
av
=
options
.
av
bv
=
options
.
bv
localsystem
=
Set_SystemUnits_From_Options
(
options
)
#######################################
# define output system of unit
#######################################
outputunits
=
UnitSystem
(
'mks'
,[
Unit_kpc
,
Unit_Ms
,
Unit_yr
,
Unit_K
,
Unit_mol
,
Unit_C
])
#######################################
# open sm
#######################################
g
=
Graph_Init
(
ps
)
Graph_SetDefaultsGraphSettings
(
g
)
colors
=
Graph_SetColorsForFiles
(
files
,
col
)
#######################################
# LOOP
#######################################
# read files
for
file
in
files
:
nb
=
Nbody
(
file
,
ftype
=
ftype
)
#nb.hdcenter()
if
reduc
!=
None
:
nb
=
nb
.
reduc
(
reduc
)
x
=
nb
.
rxyz
()
.
astype
(
Float
)
y
=
nb
.
u
.
astype
(
Float
)
#
nb
.
localsystem_of_units
=
localsystem
# compute temperature
y
=
nb
.
P
()
# use log
if
log
!=
None
:
x
,
y
=
Graph_UseLog
(
x
,
y
,
log
)
if
file
==
files
[
0
]:
xmin
,
xmax
,
ymin
,
ymax
=
Graph_SetLimits
(
g
,
xmin
,
xmax
,
ymin
,
ymax
,
x
,
y
)
Graph_DrawBox
(
g
,
xmin
,
xmax
,
ymin
,
ymax
,
log
)
g
.
ctype
(
colors
[
file
])
# plot points
if
histogram
:
xb
,
yb
=
Graph_MakeHistrogram
(
g
,
xmin
,
xmax
,
ymin
,
ymax
,
x
,
y
)
g
.
connect
(
xb
,
yb
)
else
:
g
.
points
(
x
,
y
)
g
.
ctype
(
0
)
if
log
==
'xy'
or
log
==
'yx'
:
g
.
xlabel
(
'log R'
)
g
.
ylabel
(
'log T'
)
elif
log
==
'x'
:
g
.
xlabel
(
'log R'
)
g
.
ylabel
(
'T'
)
elif
log
==
'y'
:
g
.
xlabel
(
'R'
)
g
.
ylabel
(
'log T'
)
else
:
g
.
xlabel
(
'R'
)
g
.
ylabel
(
'T'
)
# -- end ---
Graph_End
(
g
,
ps
)
Event Timeline
Log In to Comment