Page Menu
Home
c4science
Search
Configure Global Search
Log In
Files
F92693181
ghelp
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, Nov 22, 20:19
Size
6 KB
Mime Type
text/x-python
Expires
Sun, Nov 24, 20:19 (1 d, 19 h)
Engine
blob
Format
Raw Data
Handle
22488673
Attached To
rGTOOLS Gtools
ghelp
View Options
#!/usr/bin/env python
print
"""
################################################
gscripts/
################################################
################################################
# run preparation
galloc : determine the memory allocation for gadget
gmktimesteps : set output timesteps for cosmological simulations
gmktimesteps -p 10 > outputtimesteps.txt
gread_header : return some info on the header (incomplete, actually return the atime)
gread_header snap.dat
gheader : return some info on the header
gheader snap.dat
addgas : add gas to a cosmological box
################################################
# global statitics on the run
gcpu : plot cpu info
gcpu cpu.txt
gcpu -o budget --legend cpu.txt
gcpu -o gravity --legend cpu.txt
gcpu -o hydro --legend cpu.txt
genergy : plot info relative to the energy
genergy --relative energy.txt
genergy -o EnergyKineticFeedback energy.txt
gplot_energy : plot info relative to the energy budget
gplot_energy --legend energy.txt
gsteps : plot info relative to the steps
gsteps info.txt
################################################
pscripts/
################################################
################################################
# global output analysis
# star formation rate and stellar mass (sfr.txt)
pSfrvsTime --param params -o SFR sfr.txt
pSfrvsTime --param params -o MSTARS sfr.txt
pSfrvsTime --param params -o MNEWSTARS sfr.txt
pSfrvsTime --param params -o sfr sfr.txt
pSfrvsTime --param params -o mnewstars --integrate --derive --rf 100 sfr.txt
pSfrvsTime --param params -o mstars sfr.txt
# supernovae (chimie.txt)
gSNvsTime -o NSNII,NSNIa --rf 1000 chimie.txt
gSNvsTime -o EgySNThermal,EgySNKinetic --rf 1000 chimie.txt
gSNvsTime -o EgySNThermal,EgySNKinetic --integrate chimie.txt
gSNvsTime -o Nflag --rf 1000 chimie.txt
gSNvsTime -o Nflag --integrate chimie.txt
# accretion (accretion.txt)
gAccretionvsTime -o MT accretion.txt
gAccretionvsTime -o MT_th accretion.txt
gAccretionvsTime -o DMDT accretion.txt
gAccretionvsTime -o DMDT_th accretion.txt
################################################
# snapshot analysis
# logRho - logT
pplot --param params --select=gas --x logrho --y logT --xmin -3 --xmax 0 --ymin 1 --ymax 5 snap.dat
pplot --param params --select=gas --x logrho --y logT --xmin -3 --xmax 0 --ymin 1 --ymax 5 --map snap.dat
pplot --param params --select=gas --colorbar --log z --x logrho --y logT --z Tcool --xmin -3 --xmax 0 --ymin 1 --ymax 5 --zmin 1 --zmax 1e16 snap.dat
# rho - T
pplot --param params --select=gas --x rho --y T --log xy --xmin 1e-3 --xmax 1e0 --ymin 1e1 --ymax 1e5 snap.dat
# R - T
pplot --param params --select=gas --x R --y logT snap.dat
# R - Fe
pplot --param params --select=gas --x R --y Fe snap.dat
# MgFe - Fe
pplot --param params --select=stars1 --x Fe --y MgFe --xmin=-4 --xmax=0.5 --ymin=-1 --ymax=1.5 snap.dat
pplot --param params --select=stars1 --x Fe --y MgFe --xmin=-4 --xmax=0.5 --ymin=-1 --ymax=1.5 --map snap.dat
# Fe - Age
pplot --param params --select=stars1 --x Age --y Fe --xmin=0 --xmax=14 --ymin=-4 --ymax=0.5 snap.dat
# Tcool
pplot --param params --select=gas --x logrho --y Tcool --xmin=-5 --xmax=2 --ymin 0 --ymax 1e16 --log y snap.dat
################################
# metallicite
################################
# MgFevsFe
pMgFevsFe --select=stars1 --xmin=-4 --xmax=0.5 --ymin=-1 --ymax=1.5 --map snap.dat
# BaFevsFe
pBaFevsFe --select=stars1 --xmin=-4 --xmax=0.5 --ymin=-1 --ymax=1.5 --map snap.dat
# [E1/E2] vs [E3/E4]
pplot --select=stars1 --xmin=-4 --xmax=0.5 --ymin=-1 --ymax=1.5 --map --x Fe/H --y Ba/Fe snap.dat
# NvsAge
pNvsX --param params --select stars1 --x Age --nx 14 --xmin=0 --xmax=14 snap.dat
# NvsFe
pNvsX --param params --select stars1 --x Fe --nx 40 --xmin=-4 --xmax=0.5 --ymin=0 --ymax=0.15 snap.dat
################################
# profiles
################################
# grille spherique
# density
pSphericalProfile -o density --xmax 10 --rmax 10 --nr 64 --center histocenter --log xy --param params snap.dat
# masse
pSphericalProfile -o mass --xmax 10 --rmax 10 --nr 64 --center histocenter --log xy --param params snap.dat
# integrated mass
pSphericalProfile -o imass --xmax 10 --rmax 10 --nr 64 --center histocenter --log xy --param params snap.dat
# velocity disp. in z
pSphericalProfile -o sigmaz --xmax 10 --rmax 10 --nr 64 --center histocenter --param params snap.dat
# circular velocity
pSphericalProfile -o vcirc --xmax 10 --rmax 10 --nr 64 --center histocenter --param params --eps 0.1 snap.dat
# grille cylindrique
# density
pCylindricalProfile -o sdens --xmax 10 --rmax 10 --nr 64 --center histocenter --log xy --param params snap.dat
# masse
pCylindricalProfile -o mass --xmax 10 --rmax 10 --nr 64 --center histocenter --log xy --param params snap.dat
# integrated mass
pCylindricalProfile -o imass --xmax 10 --rmax 10 --nr 64 --center histocenter --log xy --param params snap.dat
# velocity disp. in z
pCylindricalProfile -o sigmaz --xmax 10 --rmax 10 --nr 64 --center histocenter --param params snap.dat
# circular velocity
pCylindricalProfile -o vcirc --xmax 10 --rmax 10 --nr 64 --center histocenter --param params --eps 0.1 snap.dat
# luminosity profile
pCylindricalProfile -o Lum --xmax 10 --rmax 10 --nr 64 --center histocenter --select stars1 --param params snap.dat
############################################################
# velocities and velocity dispertions in cylindrical coords
############################################################
# extract all
gExtractCylindrical --Rmax 10 --components="('gas','stars1','halo1')" --disk="('gas','stars1')" --nR 32 --eps=0.01 --nmin 3 --param params snap.dat
# number of points per bins
pplotdmp --mode=Numbers gas.dmp
# circular velocity per component
pplotdmp --mode=Velocities --legend gas.dmp stars1.dmp halo1.dmp
# velocities and velocity disp
pplotdmp --mode=Dispersions --legend gas.dmp
# stability Toomre paramter Q
pplotdmp --mode=Q gas.dmp
"""
Event Timeline
Log In to Comment