'''
Author: Alice Radcliffe
'''
import numpy as np
import h5py

teff = 500
log = 4.0
co = 0.55
met = 1.00
fsed = 5

filename = f"/travail/aradcliffe/exo_k/HR_grid/all_chunks/sorted/500-600K/spect_{teff:.0f}K_logg{log}_met{met:.2f}_CO{co:.2f}_fsed{fsed}.h5" 
filename_wavenumber = "path/to/wavenumber/file/wavenumber.h5"                               #only one wavenumber file for the whole grid 

with h5py.File(filename, "r") as f:
    fk = f["flux"][:]                       #flux as a function of wavenumber [W/m2/cm]

with h5py.File(filename_wavenumber, "r") as k:
    wn = k["wavenumber"][:]                  # wavenumber, [cm-1]