Ejemplo n.º 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
Ejemplo n.º 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
Ejemplo n.º 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"
Ejemplo n.º 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"
Ejemplo n.º 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"
Ejemplo n.º 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
Ejemplo n.º 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"
Ejemplo n.º 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 = [ \
Ejemplo n.º 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]
Ejemplo n.º 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"
Ejemplo n.º 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_"
Ejemplo n.º 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"
Ejemplo n.º 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]
Ejemplo n.º 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"
Ejemplo n.º 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"
Ejemplo n.º 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_"
Ejemplo n.º 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"
Ejemplo n.º 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]
Ejemplo n.º 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)]
Ejemplo n.º 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
Ejemplo n.º 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]
Ejemplo n.º 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
Ejemplo n.º 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
Ejemplo n.º 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"

Ejemplo n.º 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]
Ejemplo n.º 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
Ejemplo n.º 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 = [ \
Ejemplo n.º 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"
Ejemplo n.º 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"
Ejemplo n.º 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"
Ejemplo n.º 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
Ejemplo n.º 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"
Ejemplo n.º 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) ]
Ejemplo n.º 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
Ejemplo n.º 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"
Ejemplo n.º 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)]
Ejemplo n.º 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"
Ejemplo n.º 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) ]
Ejemplo n.º 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