In [1]:

    

using spec_io
using HDF5
using InterpSpec


    

    




loaded grid
These parameters are stored in this grid
{5600,5700,5800,5900,6000,6100,6200,6300,6400,6500,6600}
{2.5,3.0,3.5,4.0,4.5,5.0,5.5,6.0}
{-1.5,-1.0,-0.5,0.0,0.5,1.0}
Size and strides
(11,8,6,216631)
(1,11,88,528)
(11,8,6)(1,11,88)
216631

Will interpolate over  5600 1.0 5700 2.0 5800 3.0 5900 4.0 6000 5.0 6100 6.0 6200 7.0 6300 7.999999999999999 6400 9.0 6500 10.0 6600 11.0
Will interpolate over  2.5 1.0 3.0 2.0 3.5 3.0 4.0 4.0 4.5 5.0 5.0 6.0 5.5 7.0 6.0 8.0
Will interpolate over  -1.5 1.0 -1.0 2.0 -0.5 3.0000000000000004 0.0 3.9999999999999996 0.5 5.000000000000001 1.0 6.0completed interp_invert!
InterpGridCoefs done

In [3]:

    

fid = h5open("libraries/PHOENIX_F_julia_hires.hdf5", "r")

(temps, loggs, Zs), grid = get_grid(fid)
println("loaded grid")


    

    




loaded grid

In [2]:

    

using PyPlot


    

    




INFO: Loading help data...

In [3]:

    

for temp in 5600:10:6600
    plot(temp, interp_spec(temp, 4.0, -0.5)[4000], "b.")
end

temps = [temp for temp in 5600:100:6600]
for i=1:11
    plot(temps[i], rawgrid[i, 4, 3, 4000], "ko")
end

title("pixel 4000, Quadratic")
xlabel("temperature")
ylabel("flux")

savefig("temp_4000_quad.png")


    

    




Setting position to[1.0,4.0,3.0000000000000004]
Setting position to[1.1,4.0,3.0000000000000004]
Setting position to[1.2,4.0,3.0000000000000004]
Setting position to[1.3,4.0,3.0000000000000004]
Setting position to[1.4,4.0,3.0000000000000004]
Setting position to[1.5,4.0,3.0000000000000004]
Setting position to[1.6,4.0,3.0000000000000004]
Setting position to[1.7,4.0,3.0000000000000004]
Setting position to[1.8,4.0,3.0000000000000004]
Setting position to[1.9,4.0,3.0000000000000004]
Setting position to[2.0,4.0,3.0000000000000004]
Setting position to[2.1,4.0,3.0000000000000004]
Setting position to[2.2,4.0,3.0000000000000004]
Setting position to[2.3,4.0,3.0000000000000004]
Setting position to[2.4,4.0,3.0000000000000004]
Setting position to[2.5,4.0,3.0000000000000004]
Setting position to[2.6,4.0,3.0000000000000004]
Setting position to[2.7,4.0,3.0000000000000004]
Setting position to[2.8,4.0,3.0000000000000004]
Setting position to[2.9,4.0,3.0000000000000004]
Setting position to[3.0,4.0,3.0000000000000004]
Setting position to[3.1,4.0,3.0000000000000004]
Setting position to[3.2,4.0,3.0000000000000004]
Setting position to[3.3,4.0,3.0000000000000004]
Setting position to[3.4,4.0,3.0000000000000004]
Setting position to[3.5,4.0,3.0000000000000004]
Setting position to[3.6,4.0,3.0000000000000004]
Setting position to[3.7,4.0,3.0000000000000004]
Setting position to[3.8,4.0,3.0000000000000004]
Setting position to[3.9,4.0,3.0000000000000004]
Setting position to[4.0,4.0,3.0000000000000004]
Setting position to[4.1,4.0,3.0000000000000004]
Setting position to[4.2,4.0,3.0000000000000004]
Setting position to[4.3,4.0,3.0000000000000004]
Setting position to[4.4,4.0,3.0000000000000004]
Setting position to[4.5,4.0,3.0000000000000004]
Setting position to[4.6,4.0,3.0000000000000004]
Setting position to[4.7,4.0,3.0000000000000004]
Setting position to[4.8,4.0,3.0000000000000004]
Setting position to[4.9,4.0,3.0000000000000004]
Setting position to[5.0,4.0,3.0000000000000004]
Setting position to[5.1,4.0,3.0000000000000004]
Setting position to[5.2,4.0,3.0000000000000004]
Setting position to[5.3,4.0,3.0000000000000004]
Setting position to[5.4,4.0,3.0000000000000004]
Setting position to[5.5,4.0,3.0000000000000004]
Setting position to[5.6,4.0,3.0000000000000004]
Setting position to[5.7,4.0,3.0000000000000004]
Setting position to[5.8,4.0,3.0000000000000004]
Setting position to[5.9,4.0,3.0000000000000004]
Setting position to[6.0,4.0,3.0000000000000004]
Setting position to[6.1,4.0,3.0000000000000004]
Setting position to[6.2,4.0,3.0000000000000004]
Setting position to[6.3,4.0,3.0000000000000004]
Setting position to[6.4,4.0,3.0000000000000004]
Setting position to[6.5,4.0,3.0000000000000004]
Setting position to[6.6,4.0,3.0000000000000004]
Setting position to[6.7,4.0,3.0000000000000004]
Setting position to[6.8,4.0,3.0000000000000004]
Setting position to[6.9,4.0,3.0000000000000004]
Setting position to[7.0,4.0,3.0000000000000004]
Setting position to[7.1,4.0,3.0000000000000004]
Setting position to[7.2,4.0,3.0000000000000004]
Setting position to[7.299999999999999,4.0,3.0000000000000004]
Setting position to[7.4,4.0,3.0000000000000004]
Setting position to[7.5,4.0,3.0000000000000004]
Setting position to[7.599999999999999,4.0,3.0000000000000004]
Setting position to[7.699999999999999,4.0,3.0000000000000004]
Setting position to[7.799999999999999,4.0,3.0000000000000004]
Setting position to[7.8999999999999995,4.0,3.0000000000000004]
Setting position to[7.999999999999999,4.0,3.0000000000000004]
Setting position to[8.099999999999998,4.0,3.0000000000000004]
Setting position to[8.2,4.0,3.0000000000000004]
Setting position to[8.3,4.0,3.0000000000000004]
Setting position to[8.399999999999999,4.0,3.0000000000000004]
Setting position to[8.5,4.0,3.0000000000000004]
Setting position to[8.6,4.0,3.0000000000000004]
Setting position to[8.7,4.0,3.0000000000000004]
Setting position to[8.8,4.0,3.0000000000000004]
Setting position to[8.9,4.0,3.0000000000000004]
Setting position to[9.0,4.0,3.0000000000000004]
Setting position to[9.1,4.0,3.0000000000000004]
Setting position to[9.2,4.0,3.0000000000000004]
Setting position to[9.3,4.0,3.0000000000000004]
Setting position to[9.4,4.0,3.0000000000000004]
Setting position to[9.5,4.0,3.0000000000000004]
Setting position to[9.6,4.0,3.0000000000000004]
Setting position to[9.7,4.0,3.0000000000000004]
Setting position to[9.8,4.0,3.0000000000000004]
Setting position to[9.9,4.0,3.0000000000000004]
Setting position to[10.0,4.0,3.0000000000000004]
Setting position to[10.1,4.0,3.0000000000000004]
Setting position to[10.2,4.0,3.0000000000000004]
Setting position to[10.3,4.0,3.0000000000000004]
Setting position to[10.4,4.0,3.0000000000000004]
Setting position to[10.5,4.0,3.0000000000000004]
Setting position to[10.6,4.0,3.0000000000000004]
Setting position to[10.7,4.0,3.0000000000000004]
Setting position to[10.8,4.0,3.0000000000000004]
Setting position to[10.9,4.0,3.0000000000000004]
Setting position to[11.0,4.0,3.0000000000000004]






    

In [5]:

    

for z in -1.:0.1:0.5  
    plot(z, interp_spec(6000, 4.0, z)[4000],  "b.")
end

zs = [z for z in -1.5:0.5:1.0]
for i=1:6
    plot(zs[i], rawgrid[5, 4, i, 4000], "ko")
end

title("pixel 1")
xlabel("Z")
ylabel("flux")

savefig("Z4000.png")


    

    




Setting position to[5.0,4.0,2.0]
Setting position to[5.0,4.0,2.2]
Setting position to[5.0,4.0,2.4000000000000004]
Setting position to[5.0,4.0,2.6000000000000005]
Setting position to[5.0,4.0,2.8000000000000003]
Setting position to[5.0,4.0,3.0000000000000004]
Setting position to[5.0,4.0,3.2000000000000006]
Setting position to[5.0,4.0,3.4]
Setting position to[5.0,4.0,3.5999999999999996]
Setting position to[5.0,4.0,3.7999999999999994]
Setting position to[5.0,4.0,3.9999999999999996]
Setting position to[5.0,4.0,4.2]
Setting position to[5.0,4.0,4.4]
Setting position to[5.0,4.0,4.6]
Setting position to[5.0,4.0,4.800000000000001]
Setting position to[5.0,4.0,5.000000000000001]






    

In [7]:

    

plt.plot(vec(grid[1, 4, 3, 1:50]), "bo")
plt.plot(vec(grid[2, 4, 3, 1:50]), "go")
plt.plot(vec(grid[3, 4, 3, 1:50]), "ro")
plt.show()


    

In [7]:

    

plot(interp_spec(6100, 4.0, -0.5)[3950:4050])
plot(interp_spec(6150, 4.0, -0.5)[3950:4050])
plot(interp_spec(6200, 4.0, -0.5)[3950:4050])


    

    




Setting position to[6.0,4.0,3.0000000000000004]
Setting position to[6.5,4.0,3.0000000000000004]
Setting position to[7.0,4.0,3.0000000000000004]






    













    
Out[7]:





1-element Array{Any,1}:
 PyObject <matplotlib.lines.Line2D object at 0x7f9fd8db15c0>