示例#1
0
文件: f1.py 项目: optimux/sphlatch
mp.use('Agg')
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),
示例#2
0
import matplotlib.pyplot as plt

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
示例#3
0
import numpy as np
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":
示例#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')
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"
示例#5
0
文件: r0.py 项目: optimux/sphlatch
#!/usr/bin/env ipython

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_"
示例#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
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]
示例#7
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
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"
示例#8
0
from simulation import resolvePath
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
示例#9
0
import matplotlib.pyplot as plt

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
示例#10
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
示例#11
0
from simulation import resolvePath
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"
示例#12
0
文件: r3.py 项目: optimux/sphlatch
from gi_plot import *
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)
示例#13
0
import numpy as np
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":
示例#14
0
文件: m3H.py 项目: optimux/sphlatch
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
示例#15
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]
示例#16
0
from gi_plot import *
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
示例#17
0
mp.use('Agg')
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) ]
示例#18
0
文件: i0.py 项目: optimux/sphlatch
mp.use('Agg')
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_"