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]
cfg_ecc_m.imgdir = sscfg.dir + "/viz_ecc____m"
cfg_eng_m.imgdir = sscfg.dir + "/viz_eng____m"
cfg_T___XY_n.imgdir = sscfg.dir + "/viz_T___XY_n"
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) ]
#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"
sscfg.nocpus = 8
sscfg.runargs = "dump 16"
sscfg.keepfiles = "clumps.h5part disk.hdf5 moons.pdf"
eVinK = 11604.505
#axnorm = [-8.0e9 , 8.0e9 , -4.5e9 , 4.5e9 ]
axwide = [-3.2e10, 1.6e10, -1.8e10, 0.9e10]
axmedi = [-1.6e10, 8.0e9 , -9.0e9 , 4.5e9 ]
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_ecc_m.imgdir = sscfg.dir + "/viz_ecc____m"
cfg_eng_m.imgdir = sscfg.dir + "/viz_eng____m"
vizcfg = GIvizConfig()
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)]
#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"
sscfg.nocpus = 8
sscfg.runargs = "dump 16"
sscfg.keepfiles = "clumps.h5part disk.hdf5 moons.pdf"
eVinK = 11604.505
#axnorm = [-8.0e9 , 8.0e9 , -4.5e9 , 4.5e9 ]
axwide = [-3.2e10, 1.6e10, -1.8e10, 0.9e10]
axmedi = [-1.6e10, 8.0e9, -9.0e9, 4.5e9]
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_ecc_m.imgdir = sscfg.dir + "/viz_ecc____m"
cfg_eng_m.imgdir = sscfg.dir + "/viz_eng____m"
vizcfg = GIvizConfig()
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_"
ssconf.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"
ssconf.delcmd = "qdel $JOBID"
ssconf.nocpus = 4
ssconf.runargs = "dump 8"
axnorm = [-1.6e10, 1.6e10, -0.9e10, 0.9e10]
axzoom = [-3.2e9, 3.2e9, -1.8e9, 1.8e9]
eVinK = 11604.505
cfg_mat_XY_n = GIplotConfig()
cfg_mat_XY_n.ax = axnorm
cfg_mat_XY_n.postPlot = plotSimParams
cfg_mat_XY_n.imgdir = ssconf.dir + "/viz_mat_XY_n"
cfg_mat_YZ_n = GIplotConfig()
cfg_mat_YZ_n.ax = axnorm
cfg_mat_YZ_n.X = 1
cfg_mat_YZ_n.Y = 2
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"
cfg_T___XY_n.imgdir = sscfg.dir + "/viz_T___XY_n"
cfg_T___XY_n.scal_min = 150. / eVinK
cfg_T___XY_n.MminPlt = 1.e-5*ME
cfg_T___XY_n.MminPlt = 5.e-4*ME
cfg_T___XY_n.axreluse = True
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"
cfg_T___XY_n.imgdir = sscfg.dir + "/viz_T___XY_n"
cfg_T___XY_n.scal_min = 150. / eVinK
cfg_T___XY_n.MminPlt = 1.e-5 * ME
cfg_T___XY_n.MminPlt = 5.e-4 * ME
cfg_T___XY_n.axreluse = True
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
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]