コード例 #1
0
ファイル: mysql.py プロジェクト: PascalSteger/darcoda
def snap_exists(n):
    exists=my.file_exists(d(n)+"/info_"+str(n).zfill(5)+".txt");
    return exists
コード例 #2
0
ファイル: finish_all.py プロジェクト: hasanmoudud/darcoda
if i != 3:
    print "Usage: finish_all.py snapstart snapstop"
    exit(1)

os.nice(1)

nstart = int(sys.argv[1])
nstop = int(sys.argv[2])

show = True
run = True
radius = 500  #[kpc/h]
simdir = "/scratch/psteger/sim/"
for ncounter in range(nstop - nstart + 1):
    nc = ncounter + nstart
    num5 = str(nc).zfill(5)
    d = simdir + "/output_" + num5 + "/"

    if (not my.file_exists(d + "info_" + num5 + ".txt")):
        continue
    os.system("head " + d + "/info_" + num5 +
              ".txt | tail -n1 | cut -d'=' -f2>" + d + "aexp")
    af = open(d + "aexp", "r")
    a = af.readline()
    print float(a)

    cmd = "finish_gas_box.py " + simdir + "mov/gas/gas_boxall_" + num5 + ".png " + simdir + "mov/gas/fgas_boxall_" + num5 + ".png " + str(
        nc) + " " + str(float(a)) + " " + str(radius)
    if (show): print cmd
    if (run): my.thread(cmd)
コード例 #3
0
ファイル: finish_all.py プロジェクト: PascalSteger/darcoda
# check syntax
i=len(sys.argv)
if i!=3:
    print "Usage: finish_all.py snapstart snapstop"
    exit(1)

os.nice(1)

nstart = int(sys.argv[1])
nstop = int(sys.argv[2])

show=True
run=True
radius=500 #[kpc/h]
simdir = "/scratch/psteger/sim/"
for ncounter in range(nstop-nstart+1):
    nc = ncounter + nstart
    num5 = str(nc).zfill(5)
    d  = simdir + "/output_"+num5+"/"

    if(not my.file_exists(d+"info_"+num5+".txt")):
        continue
    os.system("head "+d+"/info_"+num5+".txt | tail -n1 | cut -d'=' -f2>"+d+"aexp")
    af = open(d+"aexp","r")
    a  = af.readline()
    print float(a)

    cmd = "finish_gas_box.py "+simdir+"mov/gas/gas_boxall_"+num5+".png "+simdir+"mov/gas/fgas_boxall_"+num5+".png "+str(nc)+" "+str(float(a))+" "+str(radius)
    if(show): print cmd
    if(run):  my.thread(cmd)