示例#1
0
def grid_radex(radexfilename, outfile, flow, fhigh):
    """
    Use the less-hackey version of read_radex to read a radex.out file and turn it into
    a machine-readable .dat table
    """
    grid = open(outfile, 'w')
    fmt = '%10.3e %10.3e %10.3e %10.3e %10.3e %10.3e %10.3e %10.3e %10.3e %10.3e %10.3e \n'
    grid.write(
        fmt.replace('.3e', 's') %
        ("Temperature", "log10(dens)", "log10(col)", "Tex_low", "Tex_hi",
         "TauLow", "TauUpp", "TrotLow", "TrotUpp", "FluxLow", "FluxUpp"))

    radexfile = open(radexfilename)

    rmin = 100
    rmax = 0.1

    while (1):
        radex_out = read_radex.read_radex(radexfile, flow, fhigh)
        if radex_out == 0:
            break
        else:
            temp, dens, col, tlow, tupp, taulow, tauupp, trotlow, trotupp, fluxlow, fluxupp = radex_out

        grid.write(fmt % (temp, math.log10(dens), math.log10(col), tlow, tupp,
                          taulow, tauupp, trotlow, trotupp, fluxlow, fluxupp))

    grid.close()
    radexfile.close()
示例#2
0
def grid_radex(radexfilename,outfile,flow,fhigh):
    """
    Use the less-hackey version of read_radex to read a radex.out file and turn it into
    a machine-readable .dat table
    """
    grid = open(outfile,'w')
    fmt  = '%10.3e %10.3e %10.3e %10.3e %10.3e %10.3e %10.3e %10.3e %10.3e %10.3e %10.3e \n'
    grid.write(fmt.replace('.3e','s') % ("Temperature","log10(dens)",
        "log10(col)","Tex_low","Tex_hi","TauLow","TauUpp","TrotLow","TrotUpp","FluxLow","FluxUpp"))

    radexfile  = open(radexfilename)

    rmin = 100
    rmax = 0.1

    while(1):
        radex_out = read_radex.read_radex(radexfile,flow,fhigh)
        if radex_out == 0:
            break
        else:
            temp,dens,col,tlow,tupp,taulow,tauupp,trotlow,trotupp,fluxlow,fluxupp = radex_out

        grid.write(fmt %(temp, math.log10(dens), math.log10(col),
            tlow, tupp, taulow, tauupp, trotlow,trotupp,fluxlow,fluxupp))

    grid.close()
    radexfile.close()
示例#3
0
文件: radex_grid.py 项目: Fade89/agpy
    if verbose > 0: print "Processor %i: Beginning output parsing." % mpirank
    grid = open(gfil,'w')
    fmt  = '%10.3e %10.3e %10.3e %10.3e %10.3e %10.3e %10.3e %10.3e %10.3e %10.3e %10.3e \n'
    grid.write(fmt.replace('.3e','s') % ("Temperature","log10(dens)",
        "log10(col)","Tex_low","Tex_hi","TauLow","TauUpp","TrotLow","TrotUpp","FluxLow","FluxUpp"))

    outfile  = open('radex.out')

    rmin = 100
    rmax = 0.1

    for temp in temperatures:
        for dens in densities:
            for col in columns:

                temp,dens,col,tlow,tupp,taulow,tauupp,trotlow,trotupp,fluxlow,fluxupp = read_radex(outfile,lowfreq,uppfreq)

                grid.write(fmt %(temp, math.log10(dens), math.log10(col),
                    tlow, tupp, taulow, tauupp, trotlow,trotupp,fluxlow,fluxupp))

    grid.close()
    outfile.close()

stop = time.time()
dure = stop - start
if verbose > 0: print "Processor %i Run time = %f seconds" % (mpirank,dure)
if mpisize > 1:
    os.chdir(pwd)

MPI.COMM_WORLD.Barrier()
if mpisize > 1 and mpirank == 0:
示例#4
0
    fmt = '%10.3e %10.3e %10.3e %10.3e %10.3e %10.3e %10.3e %10.3e %10.3e %10.3e %10.3e \n'
    grid.write(
        fmt.replace('.3e', 's') %
        ("Temperature", "log10(dens)", "log10(col)", "Tex_low", "Tex_hi",
         "TauLow", "TauUpp", "TrotLow", "TrotUpp", "FluxLow", "FluxUpp"))

    outfile = open('radex.out')

    rmin = 100
    rmax = 0.1

    for temp in temperatures:
        for dens in densities:
            for col in columns:

                temp, dens, col, tlow, tupp, taulow, tauupp, trotlow, trotupp, fluxlow, fluxupp = read_radex(
                    outfile, lowfreq, uppfreq)

                grid.write(
                    fmt % (temp, math.log10(dens), math.log10(col), tlow, tupp,
                           taulow, tauupp, trotlow, trotupp, fluxlow, fluxupp))

    grid.close()
    outfile.close()

stop = time.time()
dure = stop - start
if verbose > 0: print "Processor %i Run time = %f seconds" % (mpirank, dure)
if mpisize > 1:
    os.chdir(pwd)

MPI.COMM_WORLD.Barrier()