Beispiel #1
0
from sim_admin import SimAdmin, SimParam, SimSetConfig
from simulation import resolvePath
from gi_plot import *
from gi_viz  import GIviz, GIvizConfig
from clumps import ClumpsPlotConfig
import os
import stat
import numpy as np
import commands

nan = float('nan')
lem = 3.40000e41

sscfg = SimSetConfig()
sscfg.name = "r1"
sscfg.T = 3.50e11

sscfg.auxfiles = "./aux/ANEOS.QUARTZ ./aux/aneos_tables.hdf5"
sscfg.srcdir   = "./repos/sphlatch/apps/simple_sph"

sscfg.maketarg = "simple_sph_GHUmD_"
sscfg.binary   = "simple_sph_GHUmD_"
sscfg.attr     = [ ("umin", 1.e9), ("rhominclump", 1.5000), ("diskrmax",5.e10) ]

sscfg.nocpus  = 4
sscfg.runargs = "dump 6"

simadm = SimAdmin(sscfg)
sims = simadm._sims
Beispiel #2
0
from sim_admin import SimAdmin, SimParam, SimSetConfig
from simulation import resolvePath
from gi_plot import *
from gi_viz import GIviz, GIvizConfig
from clumps import ClumpsPlotConfig
import os
import stat
import numpy as np
import commands

nan = float('nan')
lem = 3.40000e41

sscfg = SimSetConfig()
sscfg.name = "r2"
sscfg.T = 3.50e11

sscfg.auxfiles = "./aux/ANEOS.QUARTZ ./aux/aneos_tables.hdf5"
sscfg.srcdir = "./repos/sphlatch/apps/simple_sph"

sscfg.maketarg = "simple_sph_GCSHUmD_"
sscfg.binary = "simple_sph_GCSHUmD_"
sscfg.attr = [("umin", 1.e9), ("rhominclump", 1.5000), ("diskrmax", 1.5e10)]

sscfg.nocpus = 4
sscfg.runargs = "dump 6"

simadm = SimAdmin(sscfg)
sims = simadm._sims
Beispiel #3
0
import matplotlib.pyplot as plt

from sim_admin import SimAdmin, SimParam, SimSetConfig
from gi_plot import *
from gi_viz import GIviz, GIvizConfig
from clumps import ClumpsPlotConfig
import os
import stat
import numpy as np
import commands

nan = float('nan')
lem = 3.40000e41

sscfg = SimSetConfig()
sscfg.name = "m0"
sscfg.bodtol = 0.2
sscfg.T = 3.50e11
sscfg.tcolanim = 0.1

sscfg.attr.append(("umin", 1.e9))
sscfg.auxfiles = "./moon_colls/aux/aneos.input ./moon_colls/aux/aneos_tables.hdf5"
sscfg.srcdir = "/home0/areufer/repos/sphlatch/apps/simple_sph"

sscfg.maketarg = "simple_sph_GHUaE_"
sscfg.binary = "simple_sph_GHUaE_"
sscfg.attr = [("rmaxunbound", 1.e10), ("rmaxbound", 2.e10), ("umin", 1.e9),
              ("rhominclump", 1.9875)]

#sscfg.subcmd  = "qsub -M [email protected] -b y -cwd -l h_cpu=240:00:00 -l h_vmem=512M -pe smp $NOCPUS -R y -N $SIMNAME ./$BINARY initial.h5part $SAVETIME $STOPTIME $RUNARGS"
sscfg.subcmd = "qsubPSMP $NOCPUS $SIMNAME ./$BINARY initial.h5part $SAVETIME $STOPTIME $RUNARGS"
Beispiel #4
0
from gi_plot import *
from gi_viz import GIviz, GIvizConfig
from clumps import ClumpsPlotConfig
import os
import stat
import time
import numpy as np
import commands
import socket
import datetime

nan = float('nan')
eVinK = 11604.505

sscfg = SimSetConfig()
sscfg.name = "v0"

sscfg.T = 3.50e11
sscfg.tcolanim = 0.1

sscfg.attr.append(("umin", 1.e9))

sscfg.maketarg = "simple_sph_GCSHUmD_"
sscfg.binary = "simple_sph_GCSHUmD_"
sscfg.attr = [("umin", 1.e9), ("rhominclump", 0.7000), ("diskrmax", 1.5e9),
              ("courant", 0.2)]

#sscfg.subcmd  = "qsub -M [email protected] -b y -cwd -l h_cpu=240:00:00 -l h_vmem=512M -pe smp $NOCPUS -R y -N $SIMNAME ./$BINARY initial.h5part $SAVETIME $STOPTIME $RUNARGS"
#sscfg.subcmd  = "qsubPSMP $NOCPUS $SIMNAME ./$BINARY initial.h5part $SAVETIME $STOPTIME $RUNARGS"
sscfg.nocpus = 4
sscfg.runargs = "dump 6"
Beispiel #5
0
from gi_plot import *
from gi_viz  import GIviz, GIvizConfig
from clumps import ClumpsPlotConfig
import os
import stat
import time
import numpy as np
import commands
import socket
import datetime

nan = float('nan')
eVinK = 11604.505

sscfg = SimSetConfig()
sscfg.name = "f3"

sscfg.T = 3.50e11
sscfg.tcolanim = 0.1

sscfg.attr.append( ("umin",1.e9) )

sscfg.maketarg = "simple_sph_GCSHUmD_"
sscfg.binary   = "simple_sph_GCSHUmD_"
sscfg.attr     = [ ("umin", 1.e9), ("rhominclump", 0.7000), ("diskrmax",1.5e9), ("courant",0.2) ]

#sscfg.subcmd  = "qsub -M [email protected] -b y -cwd -l h_cpu=240:00:00 -l h_vmem=512M -pe smp $NOCPUS -R y -N $SIMNAME ./$BINARY initial.h5part $SAVETIME $STOPTIME $RUNARGS"
#sscfg.subcmd  = "qsubPSMP $NOCPUS $SIMNAME ./$BINARY initial.h5part $SAVETIME $STOPTIME $RUNARGS"
sscfg.nocpus  = 4
sscfg.runargs = "dump 6"
sscfg.keepfiles = "clumps.h5part disk.hdf5 moons.pdf"
Beispiel #6
0
from sim_admin import SimAdmin, SimParam, SimSetConfig
from simulation import resolvePath
from gi_plot import *
from gi_viz import GIviz, GIvizConfig
from clumps import ClumpsPlotConfig
import os
import stat
import numpy as np
import commands

nan = float('nan')
lem = 3.40000e41

sscfg = SimSetConfig()
sscfg.name = "f0"
sscfg.T = 3.50e11

sscfg.auxfiles = "./aux/ANEOS.QUARTZ ./aux/aneos_tables.hdf5"
sscfg.srcdir = "./repos/sphlatch/apps/simple_sph"

sscfg.maketarg = "simple_sph_GHUmD_"
sscfg.binary = "simple_sph_GHUmD_"
sscfg.attr = [("umin", 1.e9), ("rhominclump", 1.5000), ("diskrmax", 5.e10)]

sscfg.nocpus = 8
sscfg.runargs = "dump 12"

simadm = SimAdmin(sscfg)
sims = simadm._sims
Beispiel #7
0
import matplotlib.pyplot as plt

from sim_admin import SimAdmin, SimParam, SimSetConfig
from gi_plot import *
from gi_viz  import GIviz, GIvizConfig
from clumps import ClumpsPlotConfig
import os
import stat
import numpy as np
import commands

nan = float('nan')
lem = 3.40000e41

sscfg = SimSetConfig()
sscfg.name = "m0"
sscfg.bodtol = 0.2
sscfg.T = 3.50e11
sscfg.tcolanim = 0.1

sscfg.attr.append( ("umin",1.e9) )
sscfg.auxfiles = "./moon_colls/aux/aneos.input ./moon_colls/aux/aneos_tables.hdf5"
sscfg.srcdir   = "/home0/areufer/repos/sphlatch/apps/simple_sph"

sscfg.maketarg = "simple_sph_GHUaE_"
sscfg.binary   = "simple_sph_GHUaE_"
sscfg.attr     = [ ("rmaxunbound",1.e10), ("rmaxbound",2.e10), ("umin", 1.e9), ("rhominclump", 1.9875) ]

#sscfg.subcmd  = "qsub -M [email protected] -b y -cwd -l h_cpu=240:00:00 -l h_vmem=512M -pe smp $NOCPUS -R y -N $SIMNAME ./$BINARY initial.h5part $SAVETIME $STOPTIME $RUNARGS"
sscfg.subcmd  = "qsubPSMP $NOCPUS $SIMNAME ./$BINARY initial.h5part $SAVETIME $STOPTIME $RUNARGS"
sscfg.delcmd  = "qdel $JOBID"
Beispiel #8
0
import pylab as pl

from plot_func import *

#rcParams['text.usetex']=True
mp.rc('text', usetex=True)
mp.rc('text.latex', preamble = '\usepackage{amssymb}, \usepackage{wasysym}')

nan = float('nan')

ssname = "c1"
if len(sys.argv) > 1:
  ssname = sys.argv[1]

sscfg = SimSetConfig()
sscfg.name = ssname
sscfg.dir = "/Users/areufer/Documents/UniTAPS/gi_sims"

simadm = SimAdmin(sscfg)

fig = plt.figure( figsize=(8, 6) )

x0 = 0.12
y0 = 0.10

dx = 0.21
dy = 0.27

params = []
if ssname == "c1":
  params = [ \
Beispiel #9
0
from sim_admin import SimAdmin, SimParam, SimSetConfig
from simulation import resolvePath
from gi_plot import *
from gi_viz import GIviz, GIvizConfig
from clumps import ClumpsPlotConfig
import os
import stat
import numpy as np
import commands
import time

nan = float('nan')
lem = 3.40000e41

sscfg = SimSetConfig()
sscfg.name = "m1H"
#sscfg.attr     = [ ("umin", 1.e9), ("rhominclump", 2.0000), ("diskrmax",5.e10) ]
sscfg.attr = [("umin", 1.e9), ("rhominclump", 2.0000), ("diskrmax", 1.5e10),
              ("rmaxunbound", 1.0e10), ("rmaxbound", 1.5e10)]
sscfg.nocpus = 12
sscfg.runargs = "dump 32"
sscfg.desnop = 1000000
sscfg.bodtol = 0.22
simadm = SimAdmin(sscfg)

sims = simadm._sims

eVinK = 11604.505

cfg_orb_XY_m.imgdir = sscfg.dir + "/viz_orb_XY_m"
#cfg_orb_XY_m.ax        = [-3.2e10, 1.6e10, -1.8e10, 0.9e10]
Beispiel #10
0
from gi_viz import GIviz, GIvizConfig
from clumps import ClumpsPlotConfig
import os
import stat
import time
import numpy as np
import commands
import socket
import datetime

nan = float('nan')
lem = 3.40000e41
eVinK = 11604.505

sscfg = SimSetConfig()
sscfg.name = "c2_P03h"

sscfg.maketarg = "simple_sph_GCSHUmD_"
sscfg.binary = "simple_sph_GCSHUmD_"
sscfg.attr = [("umin", 1.e9), ("rhominclump", 1.5000), ("diskrmax", 1.5e10),
              ("courant", 0.2)]

sscfg.nocpus = 4
sscfg.runargs = "dump 6"

simadm = SimAdmin(sscfg)
sims = simadm._sims

cfg_orb_XY_m.imgdir = sscfg.dir + "/viz_orb_XY_m"

cfg_ecc_m.imgdir = sscfg.dir + "/viz_ecc____m"
Beispiel #11
0
import matplotlib.pyplot as plt

from sim_admin import SimAdmin, SimParam, SimSetConfig
from gi_plot import *
from gi_viz  import GIviz, GIvizConfig
from clumps import ClumpsPlotConfig
import os
import stat
import numpy as np
import commands

nan = float('nan')
lem = 3.40000e41

sscfg = SimSetConfig()
sscfg.name = "i0"
sscfg.dir = "./" + sscfg.name
sscfg.bodiesdb = "./icy_colls/bodies"
sscfg.bodtol = 0.5
sscfg.T = 3.50e11
sscfg.tcolmax = 100.

sscfg.attr.append( ("umin",1.e9) )
sscfg.auxfiles = "./icy_colls/aux/aneos.input ./icy_colls/aux/aneos_tables.hdf5"
sscfg.srcdir   = "/home0/areufer/repos/sphlatch/apps/simple_sph"

sscfg.maketarg = "simple_sph_GHUaE_"
sscfg.binary   = "simple_sph_GHUaE_"
#sscfg.maketarg = "simple_sph_GHUAC_"
#sscfg.binary   = "simple_sph_GHUAC_"
Beispiel #12
0
from gi_viz import GIviz, GIvizConfig
from clumps import ClumpsPlotConfig
import os
import stat
import time
import numpy as np
import commands
import socket
import datetime

nan = float('nan')
lem = 3.40000e41
eVinK = 11604.505

sscfg = SimSetConfig()
sscfg.name = "c2_P04h"

sscfg.maketarg = "simple_sph_GCSHUmD_"
sscfg.binary = "simple_sph_GCSHUmD_"
sscfg.attr = [("umin", 1.e9), ("rhominclump", 1.5000), ("diskrmax", 1.5e10),
              ("courant", 0.2)]

sscfg.nocpus = 4
sscfg.runargs = "dump 6"

simadm = SimAdmin(sscfg)
sims = simadm._sims

cfg_orb_XY_m.imgdir = sscfg.dir + "/viz_orb_XY_m"

cfg_ecc_m.imgdir = sscfg.dir + "/viz_ecc____m"
Beispiel #13
0
from sim_admin import SimAdmin, SimParam, SimSetConfig
from simulation import resolvePath
from gi_plot import *
from gi_viz import GIviz, GIvizConfig
from clumps import ClumpsPlotConfig
import os
import stat
import numpy as np
import commands
import time

nan = float('nan')
lem = 3.40000e41

sscfg = SimSetConfig()
sscfg.name = "m2H"
#sscfg.attr     = [ ("umin", 1.e9), ("rhominclump", 2.0000), ("diskrmax",5.e10) ]
sscfg.attr = [("umin", 1.e9), ("rhominclump", 2.0000), ("diskrmax", 1.5e10),
              ("rmaxunbound", 1.5e10), ("rmaxbound", 1.5e10)]
sscfg.nocpus = 12
sscfg.runargs = "dump 32"
sscfg.desnop = 1000000
sscfg.bodtol = 0.22
simadm = SimAdmin(sscfg)

sims = simadm._sims

eVinK = 11604.505

cfg_orb_XY_m.imgdir = sscfg.dir + "/viz_orb_XY_m"
#cfg_orb_XY_m.ax        = [-3.2e10, 1.6e10, -1.8e10, 0.9e10]
Beispiel #14
0
from gi_viz import GIviz, GIvizConfig
from clumps import ClumpsPlotConfig
import os
import stat
import time
import numpy as np
import commands
import socket
import datetime

nan = float('nan')
lem = 3.40000e41
eVinK = 11604.505

sscfg = SimSetConfig()
sscfg.name = "c2_P06h"

sscfg.maketarg = "simple_sph_GCSHUmD_"
sscfg.binary = "simple_sph_GCSHUmD_"
sscfg.attr = [("umin", 1.e9), ("rhominclump", 1.5000), ("diskrmax", 1.5e10),
              ("courant", 0.2)]

sscfg.nocpus = 4
sscfg.runargs = "dump 6"

simadm = SimAdmin(sscfg)
sims = simadm._sims

cfg_orb_XY_m.imgdir = sscfg.dir + "/viz_orb_XY_m"

cfg_ecc_m.imgdir = sscfg.dir + "/viz_ecc____m"
Beispiel #15
0
from gi_viz  import GIviz, GIvizConfig
from clumps import ClumpsPlotConfig
import os
import stat
import time
import numpy as np
import commands
import socket
import datetime

nan = float('nan')
lem = 3.40000e41
eVinK = 11604.505

sscfg = SimSetConfig()
sscfg.name = "c2_P03h"

sscfg.maketarg = "simple_sph_GCSHUmD_"
sscfg.binary   = "simple_sph_GCSHUmD_"
sscfg.attr     = [ ("umin", 1.e9), ("rhominclump", 1.5000), ("diskrmax",1.5e10), ("courant",0.2) ]

sscfg.nocpus  = 4
sscfg.runargs = "dump 6"

simadm = SimAdmin(sscfg)
sims = simadm._sims

cfg_orb_XY_m.imgdir = sscfg.dir + "/viz_orb_XY_m"

cfg_ecc_m.imgdir = sscfg.dir + "/viz_ecc____m"
cfg_eng_m.imgdir = sscfg.dir + "/viz_eng____m"
Beispiel #16
0
import matplotlib as mp
mp.use('Agg')
import matplotlib.pyplot as plt

from sim_admin import SimAdmin, SimParam, SimSetConfig
from gi_plot import GIplot, GIplotConfig, colorScalarLinear, colorScalarLog10, colorClumps, colorBodyAndPhaseANEOS, plotGrid, plotSimParams
from gi_viz import GIviz, GIvizConfig
from clumps import ClumpsPlotConfig
import os
import numpy as np

nan = float('nan')

ssconf = SimSetConfig()
ssconf.name = "r0"
ssconf.dir = "./" + ssconf.name
ssconf.bodiesdb = "./bodies"
ssconf.T = 1500.

ssconf.mimp = np.array([0.1, 1.0])
ssconf.mtar = np.array([0.1, 1.0])
ssconf.vimp = np.array([1.0, 1.1, 1.15, 1.2, 1.25, 1.3, 1.5])
ssconf.impa = np.array([0, 15, 30, 45, 60, 75, 80, 85])

ssconf.attr.append(("UMIN", 1.e9))
ssconf.auxfiles = "./aux/aneos.input"
ssconf.srcdir = "~/repos/sphlatch/apps/simple_sph"
ssconf.maketarg = "simple_sph_GHUaC_"
ssconf.binary = "simple_sph_GHUaC_"
Beispiel #17
0
from gi_plot import *
from gi_viz import GIviz, GIvizConfig
from clumps import ClumpsPlotConfig
import os
import stat
import time
import numpy as np
import commands
import socket
import datetime

nan = float('nan')
eVinK = 11604.505

sscfg = SimSetConfig()
sscfg.name = "f3"

sscfg.T = 3.50e11
sscfg.tcolanim = 0.1

sscfg.attr.append(("umin", 1.e9))

sscfg.maketarg = "simple_sph_GCSHUmD_"
sscfg.binary = "simple_sph_GCSHUmD_"
sscfg.attr = [("umin", 1.e9), ("rhominclump", 0.7000), ("diskrmax", 1.5e9),
              ("courant", 0.2)]

#sscfg.subcmd  = "qsub -M [email protected] -b y -cwd -l h_cpu=240:00:00 -l h_vmem=512M -pe smp $NOCPUS -R y -N $SIMNAME ./$BINARY initial.h5part $SAVETIME $STOPTIME $RUNARGS"
#sscfg.subcmd  = "qsubPSMP $NOCPUS $SIMNAME ./$BINARY initial.h5part $SAVETIME $STOPTIME $RUNARGS"
sscfg.nocpus = 4
sscfg.runargs = "dump 6"
Beispiel #18
0
from sim_admin import SimAdmin, SimParam, SimSetConfig
from simulation import resolvePath
from gi_plot import *
from gi_viz  import GIviz, GIvizConfig
from clumps import ClumpsPlotConfig
import os
import stat
import numpy as np
import commands
import time

nan = float('nan')
lem = 3.40000e41

sscfg = SimSetConfig()
sscfg.name = "m1H"
#sscfg.attr     = [ ("umin", 1.e9), ("rhominclump", 2.0000), ("diskrmax",5.e10) ]
sscfg.attr     = [ ("umin", 1.e9), ("rhominclump", 2.0000), ("diskrmax",1.5e10), ("rmaxunbound", 1.0e10), ("rmaxbound", 1.5e10) ]
sscfg.nocpus  = 12
sscfg.runargs = "dump 32"
sscfg.desnop  = 1000000
sscfg.bodtol  = 0.22
simadm = SimAdmin(sscfg)

sims = simadm._sims

eVinK = 11604.505

cfg_orb_XY_m.imgdir = sscfg.dir + "/viz_orb_XY_m"
#cfg_orb_XY_m.ax        = [-3.2e10, 1.6e10, -1.8e10, 0.9e10]
cfg_orb_XY_m.ax        = [-1.2e10, 0.4e10, -0.6e10, 0.3e10]
Beispiel #19
0
import matplotlib.pyplot as plt

from sim_admin import SimAdmin, SimParam, SimSetConfig
from gi_plot import *
from gi_viz import GIviz, GIvizConfig
from clumps import ClumpsPlotConfig
import os
import stat
import numpy as np
import commands

nan = float('nan')
lem = 3.40000e41

sscfg = SimSetConfig()
sscfg.name = "f2"
sscfg.dir = "./" + sscfg.name
sscfg.bodiesdb = "fluffy_colls/bodies"
sscfg.resultsdb = sscfg.dir + "/results"
sscfg.bodtol = 0.2
sscfg.T = 3.50e11
sscfg.tcolanim = 0.1
sscfg.relsep = 10.

sscfg.auxfiles = "./aux/aneos.input ./aux/aneos_tables.hdf5"
sscfg.srcdir = "/home0/areufer/repos/sphlatch/apps/simple_sph"

sscfg.maketarg = "simple_sph_GHUaE_"
sscfg.binary = "simple_sph_GHUaE_"
sscfg.attr = [("rmaxunbound", 5.e10), ("rmaxbound", 5.e10), ("umin", 1.e9),
              ("rhominclump", 1.9875), ("diskrmax", 5.e10)]
Beispiel #20
0
from gi_plot import *
from gi_viz  import GIviz, GIvizConfig
from clumps import ClumpsPlotConfig
import os
import stat
import time
import numpy as np
import commands
import socket
import datetime

nan = float('nan')
lem = 3.40000e41

sscfg = SimSetConfig()
sscfg.name = "i1"
sscfg.T = 3.50e11
sscfg.desnop = 100000

sscfg.auxfiles = "./aux/ANEOS.QUARTZ ./aux/aneos_tables.hdf5"
sscfg.srcdir   = "./repos/sphlatch/apps/simple_sph"

sscfg.maketarg = "simple_sph_GCSHUmD_"
sscfg.binary   = "simple_sph_GCSHUmD_"
sscfg.attr     = [ ("umin", 1.e9), ("rhominclump", 1.5000), ("diskrmax",1.5e10), ("courant",0.2) ]

sscfg.nocpus  = 4
sscfg.runargs = "dump 6"

simadm = SimAdmin(sscfg)
sims = simadm._sims
Beispiel #21
0
from sim_admin import SimAdmin, SimParam, SimSetConfig
from simulation import resolvePath
from gi_plot import *
from gi_viz  import GIviz, GIvizConfig
from clumps import ClumpsPlotConfig
import os
import stat
import numpy as np
import commands
import time

nan = float('nan')
lem = 3.40000e41

sscfg = SimSetConfig()
sscfg.name = "m1G"
#sscfg.attr     = [ ("umin", 1.e9), ("rhominclump", 2.0000), ("diskrmax",5.e10) ]
sscfg.attr     = [ ("umin", 1.e9), ("rhominclump", 2.0000), ("diskrmax",1.5e10), ("rmaxunbound", 1.0e10), ("rmaxbound", 1.5e10) ]
sscfg.nocpus  = 16
sscfg.runargs = "dump 32"
sscfg.desnop  = 2000000
sscfg.bodtol  = 0.22
simadm = SimAdmin(sscfg)

sims = simadm._sims

eVinK = 11604.505

cfg_orb_XY_m.imgdir = sscfg.dir + "/viz_orb_XY_m"
#cfg_orb_XY_m.ax        = [-3.2e10, 1.6e10, -1.8e10, 0.9e10]
cfg_orb_XY_m.ax        = [-1.2e10, 0.4e10, -0.6e10, 0.3e10]
Beispiel #22
0
from sim_admin import SimAdmin, SimParam, SimSetConfig
from simulation import resolvePath
from gi_plot import *
from gi_viz  import GIviz, GIvizConfig
from clumps import ClumpsPlotConfig
import os
import stat
import numpy as np
import commands
import time

nan = float('nan')
lem = 3.40000e41

sscfg = SimSetConfig()
sscfg.name = "m3H"
#sscfg.attr     = [ ("umin", 1.e9), ("rhominclump", 2.0000), ("diskrmax",5.e10) ]
sscfg.attr     = [ ("umin", 1.e9), ("rhominclump", 2.0000), ("diskrmax",1.5e10), ("rmaxunbound", 1.5e10), ("rmaxbound", 1.5e10) ]
sscfg.nocpus  = 12
sscfg.runargs = "dump 32"
sscfg.desnop  = 1000000
sscfg.bodtol  = 0.22
sscfg.maketarg = "simple_sph_GMSHUmD_"
sscfg.binary   = "simple_sph_GMSHUmD_"

simadm = SimAdmin(sscfg)

sims = simadm._sims

eVinK = 11604.505
Beispiel #23
0
from sim_admin import SimAdmin, SimParam, SimSetConfig
from simulation import resolvePath
from gi_plot import *
from gi_viz import GIviz, GIvizConfig
from clumps import ClumpsPlotConfig
import os
import stat
import numpy as np
import commands
import time

nan = float('nan')
lem = 3.40000e41

sscfg = SimSetConfig()
sscfg.name = "m3H"
#sscfg.attr     = [ ("umin", 1.e9), ("rhominclump", 2.0000), ("diskrmax",5.e10) ]
sscfg.attr = [("umin", 1.e9), ("rhominclump", 2.0000), ("diskrmax", 1.5e10),
              ("rmaxunbound", 1.5e10), ("rmaxbound", 1.5e10)]
sscfg.nocpus = 12
sscfg.runargs = "dump 32"
sscfg.desnop = 1000000
sscfg.bodtol = 0.22
sscfg.maketarg = "simple_sph_GMSHUmD_"
sscfg.binary = "simple_sph_GMSHUmD_"

simadm = SimAdmin(sscfg)

sims = simadm._sims

eVinK = 11604.505
Beispiel #24
0
from gi_plot import *
from gi_viz  import GIviz, GIvizConfig
from clumps import ClumpsPlotConfig
import os
import stat
import time
import numpy as np
import commands
import socket
import datetime

nan = float('nan')
eVinK = 11604.505

sscfg = SimSetConfig()
sscfg.name = "f4"

sscfg.T = 3.50e11
sscfg.tcolanim = 0.1

sscfg.attr.append( ("umin",1.e9) )

sscfg.maketarg = "simple_sph_GCSHUmD_"
sscfg.binary   = "simple_sph_GCSHUmD_"
sscfg.attr     = [ ("umin", 1.e9), ("rhominclump", 0.7000), ("diskrmax",1.5e9), ("courant",0.2) ]

sscfg.nocpus  = 4
sscfg.runargs = "dump 6"
sscfg.keepfiles = "clumps.h5part disk.hdf5 moons.pdf"

Beispiel #25
0
from sim_admin import SimAdmin, SimParam, SimSetConfig
from simulation import resolvePath
from gi_plot import *
from gi_viz  import GIviz, GIvizConfig
from clumps import ClumpsPlotConfig
import os
import stat
import numpy as np
import commands
import time

nan = float('nan')
lem = 3.40000e41

sscfg = SimSetConfig()
sscfg.name = "m2H"
#sscfg.attr     = [ ("umin", 1.e9), ("rhominclump", 2.0000), ("diskrmax",5.e10) ]
sscfg.attr     = [ ("umin", 1.e9), ("rhominclump", 2.0000), ("diskrmax",1.5e10), ("rmaxunbound", 1.5e10), ("rmaxbound", 1.5e10) ]
sscfg.nocpus  = 12
sscfg.runargs = "dump 32"
sscfg.desnop  = 1000000
sscfg.bodtol  = 0.22
simadm = SimAdmin(sscfg)

sims = simadm._sims

eVinK = 11604.505

cfg_orb_XY_m.imgdir = sscfg.dir + "/viz_orb_XY_m"
#cfg_orb_XY_m.ax        = [-3.2e10, 1.6e10, -1.8e10, 0.9e10]
cfg_orb_XY_m.ax        = [-1.2e10, 0.4e10, -0.6e10, 0.3e10]
Beispiel #26
0
from sim_admin import SimAdmin, SimParam, SimSetConfig
from simulation import resolvePath
from gi_plot import *
from gi_viz import GIviz, GIvizConfig
from clumps import ClumpsPlotConfig
import os
import stat
import numpy as np
import commands

nan = float('nan')
lem = 3.40000e41

sscfg = SimSetConfig()
sscfg.name = "r1"
sscfg.T = 3.50e11

sscfg.auxfiles = "./aux/ANEOS.QUARTZ ./aux/aneos_tables.hdf5"
sscfg.srcdir = "./repos/sphlatch/apps/simple_sph"

sscfg.maketarg = "simple_sph_GHUmD_"
sscfg.binary = "simple_sph_GHUmD_"
sscfg.attr = [("umin", 1.e9), ("rhominclump", 1.5000), ("diskrmax", 5.e10)]

sscfg.nocpus = 4
sscfg.runargs = "dump 6"

simadm = SimAdmin(sscfg)
sims = simadm._sims
Beispiel #27
0
import pylab as pl

from plot_func import *

#rcParams['text.usetex']=True
mp.rc('text', usetex=True)
mp.rc('text.latex', preamble='\usepackage{amssymb}, \usepackage{wasysym}')

nan = float('nan')

ssname = "c1"
if len(sys.argv) > 1:
    ssname = sys.argv[1]

sscfg = SimSetConfig()
sscfg.name = ssname
sscfg.dir = "/Users/areufer/Documents/UniTAPS/gi_sims"

simadm = SimAdmin(sscfg)

fig = plt.figure(figsize=(8, 6))

x0 = 0.12
y0 = 0.10

dx = 0.21
dy = 0.27

params = []
if ssname == "c1":
    params = [ \
Beispiel #28
0
from gi_viz  import GIviz, GIvizConfig
from clumps import ClumpsPlotConfig
import os
import stat
import time
import numpy as np
import commands
import socket
import datetime

nan = float('nan')
lem = 3.40000e41
eVinK = 11604.505

sscfg = SimSetConfig()
sscfg.name = "c1"

sscfg.maketarg = "simple_sph_GCSHUmD_"
sscfg.binary   = "simple_sph_GCSHUmD_"
sscfg.attr     = [ ("umin", 1.e9), ("rhominclump", 1.5000), ("diskrmax",1.5e10), ("courant",0.2) ]

sscfg.nocpus  = 4
sscfg.runargs = "dump 6"

simadm = SimAdmin(sscfg)
sims = simadm._sims

cfg_orb_XY_m.imgdir = sscfg.dir + "/viz_orb_XY_m"

cfg_ecc_m.imgdir = sscfg.dir + "/viz_ecc____m"
cfg_eng_m.imgdir = sscfg.dir + "/viz_eng____m"
Beispiel #29
0
from gi_viz  import GIviz, GIvizConfig
from clumps import ClumpsPlotConfig
import os
import stat
import time
import numpy as np
import commands
import socket
import datetime

nan = float('nan')
lem = 3.40000e41
eVinK = 11604.505

sscfg = SimSetConfig()
sscfg.name = "c2_P04h"

sscfg.maketarg = "simple_sph_GCSHUmD_"
sscfg.binary   = "simple_sph_GCSHUmD_"
sscfg.attr     = [ ("umin", 1.e9), ("rhominclump", 1.5000), ("diskrmax",1.5e10), ("courant",0.2) ]

sscfg.nocpus  = 4
sscfg.runargs = "dump 6"

simadm = SimAdmin(sscfg)
sims = simadm._sims

cfg_orb_XY_m.imgdir = sscfg.dir + "/viz_orb_XY_m"

cfg_ecc_m.imgdir = sscfg.dir + "/viz_ecc____m"
cfg_eng_m.imgdir = sscfg.dir + "/viz_eng____m"
Beispiel #30
0
from gi_viz  import GIviz, GIvizConfig
from clumps import ClumpsPlotConfig
from sim_machine import *
import os
import stat
import time
import numpy as np
import commands
import socket
import datetime

nan = float('nan')
lem = 3.40000e41

sscfg = SimSetConfig()
sscfg.name = "i2"
sscfg.T = 3.50e11
sscfg.desnop = 200000
sscfg.bodtol = 0.1

#sscfg.auxfiles = "./aux/ANEOS.QUARTZ ./aux/aneos_tables.hdf5"
#sscfg.srcdir   = "./repos/sphlatch/apps/simple_sph"
sscfg.srcdir   = srcdir


sscfg.maketarg = "simple_sph_GCSHUmD_"
sscfg.binary   = "simple_sph_GCSHUmD_"
sscfg.attr     = [ ("umin", 1.e9), ("rhominclump", 0.7000), ("diskrmax",1.5e9), ("courant",0.2) ]

sscfg.nocpus  = 6
sscfg.runargs = "dump 8"
Beispiel #31
0
from gi_viz  import GIviz, GIvizConfig
from clumps import ClumpsPlotConfig
from sim_machine. import *
import os
import stat
import time
import numpy as np
import commands
import socket
import datetime

nan = float('nan')
lem = 3.40000e41

sscfg = SimSetConfig()
sscfg.name = "r3"
sscfg.T = 3.50e11
sscfg.desnop = 200000

sscfg.auxfiles = "./aux/ANEOS.QUARTZ ./aux/aneos_tables.hdf5"
sscfg.srcdir   = srcdir

sscfg.maketarg = "simple_sph_GCSHUmD_"
sscfg.binary   = "simple_sph_GCSHUmD_"
sscfg.attr     = [ ("umin", 1.e9), ("rhominclump", 1.5000), ("diskrmax",1.5e10), ("courant",0.2) ]

sscfg.nocpus  = 4
sscfg.runargs = "dump 6"

simadm = SimAdmin(sscfg)
sims = simadm._sims
Beispiel #32
0
from gi_viz  import GIviz, GIvizConfig
from clumps import ClumpsPlotConfig
import os
import stat
import time
import numpy as np
import commands
import socket
import datetime

nan = float('nan')
lem = 3.40000e41
eVinK = 11604.505

sscfg = SimSetConfig()
sscfg.name = "c2_P06h"

sscfg.maketarg = "simple_sph_GCSHUmD_"
sscfg.binary   = "simple_sph_GCSHUmD_"
sscfg.attr     = [ ("umin", 1.e9), ("rhominclump", 1.5000), ("diskrmax",1.5e10), ("courant",0.2) ]

sscfg.nocpus  = 4
sscfg.runargs = "dump 6"

simadm = SimAdmin(sscfg)
sims = simadm._sims

cfg_orb_XY_m.imgdir = sscfg.dir + "/viz_orb_XY_m"

cfg_ecc_m.imgdir = sscfg.dir + "/viz_ecc____m"
cfg_eng_m.imgdir = sscfg.dir + "/viz_eng____m"
Beispiel #33
0
import matplotlib.pyplot as plt

from sim_admin import SimAdmin, SimParam, SimSetConfig
from gi_plot import *
from gi_viz  import GIviz, GIvizConfig
from clumps import ClumpsPlotConfig
import os
import stat
import numpy as np
import commands

nan = float('nan')
lem = 3.40000e41

sscfg = SimSetConfig()
sscfg.name = "f2"
sscfg.dir = "./" + sscfg.name
sscfg.bodiesdb = "fluffy_colls/bodies"
sscfg.resultsdb = sscfg.dir + "/results"
sscfg.bodtol = 0.2
sscfg.T = 3.50e11
sscfg.tcolanim = 0.1
sscfg.relsep = 10.

sscfg.auxfiles = "./aux/aneos.input ./aux/aneos_tables.hdf5"
sscfg.srcdir   = "/home0/areufer/repos/sphlatch/apps/simple_sph"

sscfg.maketarg = "simple_sph_GHUaE_"
sscfg.binary   = "simple_sph_GHUaE_"
sscfg.attr     = [ ("rmaxunbound",5.e10), ("rmaxbound",5.e10), ("umin", 1.e9), ("rhominclump", 1.9875), ("diskrmax",5.e10) ]
Beispiel #34
0
from sim_admin import SimAdmin, SimParam, SimSetConfig
from simulation import resolvePath
from gi_plot import *
from gi_viz  import GIviz, GIvizConfig
from clumps import ClumpsPlotConfig
import os
import stat
import numpy as np
import commands

nan = float('nan')
lem = 3.40000e41

sscfg = SimSetConfig()
sscfg.name = "f0"
sscfg.T = 3.50e11

sscfg.auxfiles = "./aux/ANEOS.QUARTZ ./aux/aneos_tables.hdf5"
sscfg.srcdir   = "./repos/sphlatch/apps/simple_sph"

sscfg.maketarg = "simple_sph_GHUmD_"
sscfg.binary   = "simple_sph_GHUmD_"
sscfg.attr     = [ ("umin", 1.e9), ("rhominclump", 1.5000), ("diskrmax",5.e10) ]

sscfg.nocpus  = 8
sscfg.runargs = "dump 12"

simadm = SimAdmin(sscfg)
sims = simadm._sims
Beispiel #35
0
from gi_plot import *
from gi_viz  import GIviz, GIvizConfig
from clumps import ClumpsPlotConfig
import os
import stat
import time
import numpy as np
import commands
import socket
import datetime

nan = float('nan')
eVinK = 11604.505

sscfg = SimSetConfig()
sscfg.name = "v0"

sscfg.T = 3.50e11
sscfg.tcolanim = 0.1

sscfg.attr.append( ("umin",1.e9) )

sscfg.maketarg = "simple_sph_GCSHUmD_"
sscfg.binary   = "simple_sph_GCSHUmD_"
sscfg.attr     = [ ("umin", 1.e9), ("rhominclump", 0.7000), ("diskrmax",1.5e9), ("courant",0.2) ]

#sscfg.subcmd  = "qsub -M [email protected] -b y -cwd -l h_cpu=240:00:00 -l h_vmem=512M -pe smp $NOCPUS -R y -N $SIMNAME ./$BINARY initial.h5part $SAVETIME $STOPTIME $RUNARGS"
#sscfg.subcmd  = "qsubPSMP $NOCPUS $SIMNAME ./$BINARY initial.h5part $SAVETIME $STOPTIME $RUNARGS"
sscfg.nocpus  = 4
sscfg.runargs = "dump 6"
sscfg.keepfiles = "clumps.h5part disk.hdf5 moons.pdf"
Beispiel #36
0
import matplotlib.pyplot as plt

from sim_admin import SimAdmin, SimParam, SimSetConfig
from gi_plot import *
from gi_viz import GIviz, GIvizConfig
from clumps import ClumpsPlotConfig
import os
import stat
import numpy as np
import commands

nan = float('nan')
lem = 3.40000e41

sscfg = SimSetConfig()
sscfg.name = "f1"
sscfg.dir = "./" + sscfg.name
sscfg.bodiesdb = "fluffy_colls/bodies"
sscfg.resultsdb = sscfg.dir + "/results"
sscfg.bodtol = 0.2
sscfg.T = 3.50e11
sscfg.tcolanim = 0.1
sscfg.relsep = 10.

sscfg.auxfiles = "./aux/aneos.input ./aux/aneos_tables.hdf5"
sscfg.srcdir = "/home0/areufer/repos/sphlatch/apps/simple_sph"

sscfg.maketarg = "simple_sph_GHUaE_"
sscfg.binary = "simple_sph_GHUaE_"
sscfg.attr = [("rmaxunbound", 5.e10), ("rmaxbound", 5.e10), ("umin", 1.e9),
              ("rhominclump", 1.9875), ("diskrmax", 5.e10)]
Beispiel #37
0
from gi_plot import *
from gi_viz import GIviz, GIvizConfig
from clumps import ClumpsPlotConfig
import os
import stat
import time
import numpy as np
import commands
import socket
import datetime

nan = float('nan')
eVinK = 11604.505

sscfg = SimSetConfig()
sscfg.name = "f4"

sscfg.T = 3.50e11
sscfg.tcolanim = 0.1

sscfg.attr.append(("umin", 1.e9))

sscfg.maketarg = "simple_sph_GCSHUmD_"
sscfg.binary = "simple_sph_GCSHUmD_"
sscfg.attr = [("umin", 1.e9), ("rhominclump", 0.7000), ("diskrmax", 1.5e9),
              ("courant", 0.2)]

sscfg.nocpus = 4
sscfg.runargs = "dump 6"
sscfg.keepfiles = "clumps.h5part disk.hdf5 moons.pdf"
Beispiel #38
0
import matplotlib.pyplot as plt

from sim_admin import SimAdmin, SimParam, SimSetConfig
from gi_plot import *
from gi_viz  import GIviz, GIvizConfig
from clumps import ClumpsPlotConfig
import os
import stat
import numpy as np
import commands

nan = float('nan')
lem = 3.40000e41

sscfg = SimSetConfig()
sscfg.name = "f1"
sscfg.dir = "./" + sscfg.name
sscfg.bodiesdb = "fluffy_colls/bodies"
sscfg.resultsdb = sscfg.dir + "/results"
sscfg.bodtol = 0.2
sscfg.T = 3.50e11
sscfg.tcolanim = 0.1
sscfg.relsep = 10.

sscfg.auxfiles = "./aux/aneos.input ./aux/aneos_tables.hdf5"
sscfg.srcdir   = "/home0/areufer/repos/sphlatch/apps/simple_sph"

sscfg.maketarg = "simple_sph_GHUaE_"
sscfg.binary   = "simple_sph_GHUaE_"
sscfg.attr     = [ ("rmaxunbound",5.e10), ("rmaxbound",5.e10), ("umin", 1.e9), ("rhominclump", 1.9875), ("diskrmax",5.e10) ]
Beispiel #39
0
from sim_admin import SimAdmin, SimParam, SimSetConfig
from simulation import resolvePath
from gi_plot import *
from gi_viz  import GIviz, GIvizConfig
from clumps import ClumpsPlotConfig
import os
import stat
import numpy as np
import commands

nan = float('nan')
lem = 3.40000e41

sscfg = SimSetConfig()
sscfg.name = "r2"
sscfg.T = 3.50e11

sscfg.auxfiles = "./aux/ANEOS.QUARTZ ./aux/aneos_tables.hdf5"
sscfg.srcdir   = "./repos/sphlatch/apps/simple_sph"

sscfg.maketarg = "simple_sph_GCSHUmD_"
sscfg.binary   = "simple_sph_GCSHUmD_"
sscfg.attr     = [ ("umin", 1.e9), ("rhominclump", 1.5000), ("diskrmax",1.5e10) ]

sscfg.nocpus  = 4
sscfg.runargs = "dump 6"

simadm = SimAdmin(sscfg)
sims = simadm._sims