Homec4science
Diffusion pNbody (master) Browse

pNbody/src/PyGadget/test_Parameters.py
021ca11fcae0master

Last Change
Raw File

test_Parameters.py
View Options

#!/usr/bin/env python
from mpi4py import MPI
from pNbody import ic
from numpy import *
import gadget
import Ptools as pt
import sys
import time
random.seed(MPI.COMM_WORLD.Get_rank()) # 2 for two points
n = 100
nb = ic.plummer(n,1,1,1,eps=0.1,rmax=1.)
gadget.InitMPI() # init MPI
gadget.InitDefaultParameters() # init default parameters
gadget.Info()
params = gadget.GetParameters()
keys = ['TimeLimitCPU',
'ResubmitOn',
'ICFormat',
'SnapFormat',
'ComovingIntegrationOn',
'TypeOfTimestepCriterion',
'OutputListOn',
'PeriodicBoundariesOn',
'TimeBegin',
'TimeMax',
'Omega0',
'OmegaLambda',
'OmegaBaryon',
'HubbleParam',
'BoxSize',
'TimeBetSnapshot',
'TimeOfFirstSnapshot',
'CpuTimeBetRestartFile',
'TimeBetStatistics',
'NumFilesPerSnapshot',
'NumFilesWrittenInParallel',
'ErrTolIntAccuracy',
'MaxRMSDisplacementFac',
'CourantFac',
'MaxSizeTimestep',
'MinSizeTimestep',
'ErrTolTheta',
'TypeOfOpeningCriterion',
'ErrTolForceAcc',
'TreeDomainUpdateFrequency',
'DesNumNgb',
'MaxNumNgbDeviation',
'ArtBulkViscConst',
'InitGasTemp',
'MinGasTemp',
'PartAllocFactor',
'TreeAllocFactor',
'BufferSize',
'UnitLength_in_cm',
'UnitMass_in_g',
'UnitVelocity_in_cm_per_s',
'GravityConstantInternal',
'MinGasHsmlFractional',
'SofteningGas',
'SofteningHalo',
'SofteningDisk',
'SofteningBulge',
'SofteningStars',
'SofteningBndry',
'SofteningGasMaxPhys',
'SofteningHaloMaxPhys',
'SofteningDiskMaxPhys',
'SofteningBulgeMaxPhys',
'SofteningStarsMaxPhys',
'SofteningBndryMaxPhys']
for key in keys:
print key,params[key]
i=0
for key in keys:
params[key] = i
i = i + 1
print
for key in keys:
print key,params[key]
gadget.SetParameters(params)
params = gadget.GetParameters()
#gadget.LoadParticles(array(nb.npart),nb.pos,nb.vel,nb.mass,nb.num,nb.tpe)
c4science · Help