In [4]:

    

from stack import *

#N=60; L=10e3; dL=0.1*L; n0=1.0; n1=1.6
#my_stack = random_stack(N,n0,n1,dL,L)

my_stack = load('io_files/saved_stack.txt')

wlrange = [1350,1450]
my_stack.scan('transmission',wlrange,400,'')


    

In [2]:

    

threshold = 0.2
new_search = mode_search(my_stack,wlrange,1000,threshold)
new_search.plot_features('mean','sigma')


    

In [3]:

    

new_search.sort('sigma')
new_search.sort('mean',5)
new_search.plot_profiles([1,2,3,4,5])


    

In [4]:

    

new_search.plot_profiles([1,5])


    

In [ ]:

    

save(my_stack,'io_files/saved_stack_new.txt')