min_megsb_L1_wave = 33.0d0
max_megsb_L1_wave = 107.0d0
min_megsa1_L1_wave = 3.0d0
max_megsa1_L1_wave = 39.0d0
min_megsa2_L1_wave = 3.0d0
max_megsa2_L1_wave = 39.0d0
megsa_1_cutoff_freq=5.9
megsa_1to2_cutoff_freq = 17.24
megs_AtoB_cutoff_freq = 37.0
num_megsa1_L1_spectral_elements = $
long((max_megsa1_L1_wave - min_megsa1_L1_wave + 1) * 100)
num_megsa2_L1_spectral_elements = $
long((max_megsa2_L1_wave - min_megsa2_L1_wave + 1) * 100)
num_megsa_L1_spectral_elements = $
long((max_megsa2_L1_wave - min_megsa1_L1_wave + 1) * 100)
num_megsb_L1_spectral_elements = $
long((max_megsb_L1_wave - min_megsb_L1_wave) * 50)
megsa_L1_wave = eve_range(min_megsa1_L1_wave,max_megsa2_L1_wave,invdelta=100,num_megsa_L1_spectral_elements)
megsa1_l1_wave=megsa_l1_wave
megsa2_l1_wave=megsa_l1_wave
megsb_L1_wave = eve_range(min_megsb_L1_wave,max_megsb_L1_wave,invdelta=50,num_megsb_L1_spectral_elements)
megs_L2_wave= eve_range(min_megsa1_L1_wave,max_megsb_L1_wave,invdelta=50,num_megs_L2_spectral_elements)
megs_L2_atob_cutoff_bin=max(where(megs_l2_wave lt megs_atob_cutoff_freq))