"KeyError: 'calculateRamp'" in hst.py

[pmccreery]

Hello!

When running jdi.run_pandexo() with all instruments, or specifically with WFC3 G141, a KeyError returns.

Appears that in line 652 of hst.py, calRamp = pandeia_input['strategy']['calculateRamp'], the key 'calculateRamp' doesn't exist. Running the run_test.py script with the 'RUN ALL' option reproduces the error. Not sure if I should be adjusting something for HST/WFC use, but I couldn't find anything.

Thank you!

[natashabatalha]

Apologies for the delay. Can you post the two code snippets you are using?

[lheinke]

I encountered the same issue. It can be replicated even when running the example from the docs:

import pandexo.engine.justdoit as jdi

exo_dict = jdi.load_exo_dict()

#WASP-43
exo_dict['star']['mag']      = 9.397                # H magnitude of the system
#WASP-43b
exo_dict['planet']['type']    = 'user'               # user specified inputs
exo_dict['planet']['exopath'] = 'WASP43b-Eclipse_Spectrum.txt' # filename for model spectrum
exo_dict['planet']['w_unit']  = 'um'                 # wavelength unit
exo_dict['planet']['f_unit']  = 'fp/f*'              # flux ratio unit (can also put "rp^2/r*^2")
exo_dict['planet']['depth']   = 4.0e-3               # flux ratio
exo_dict['planet']['i']       = 82.6                 # Orbital inclination in degrees
exo_dict['planet']['ars']     = 5.13                 # Semi-major axis / stellar radius
exo_dict['planet']['period']  = 0.8135               # Orbital period in days
exo_dict['planet']['transit_duration']= 4170.0/60/60/24#(optional if given above info) transit duration in days
exo_dict['planet']['w']       = 90                   #(optional) longitude of periastron. Default is 90
exo_dict['planet']['ecc']     = 0                    #(optional) eccentricity. Default is 0

inst_dict = jdi.load_mode_dict('WFC3 G141')

exo_dict['observation']['noccultations']               = 5            # Number of transits/eclipses
inst_dict['configuration']['detector']['subarray']     = 'GRISM256'   # GRISM256 or GRISM512
inst_dict['configuration']['detector']['nsamp']        = 10           # WFC3 N_SAMP, 1..15
inst_dict['configuration']['detector']['samp_seq']     = 'SPARS5'     # WFC3 SAMP_SEQ, SPARS5 or SPARS10
inst_dict['strategy']['norbits']                       = 4            # Number of HST orbits
inst_dict['strategy']['nchan']                       = 15           # Number of spectrophotometric channels
inst_dict['strategy']['scanDirection']               = 'Forward'    # Spatial scan direction, Forward or Round Trip
inst_dict['strategy']['schedulability']              = 30           # 30 for small/medium program, 100 for large program
inst_dict['strategy']['windowSize']                  = 20           # (optional) Observation start window size in minutes. Default is 20 minutes.

foo = jdi.run_pandexo(exo_dict, inst_dict, output_file='wasp43b.p')

The error occurs on the last line. Here is the full traceback:

Running Single Case w/ User Instrument Dict
---------------------------------------------------------------------------
KeyError                                  Traceback (most recent call last)
Cell In[6], [line 1](vscode-notebook-cell:?execution_count=6&line=1)
----> [1](vscode-notebook-cell:?execution_count=6&line=1) foo = jdi.run_pandexo(exo_dict, inst_dict, output_file='wasp43b.p')

File [~/.local/mambaforge/envs/pandexo/lib/python3.11/site-packages/pandexo/engine/justdoit.py:379](https://file+.vscode-resource.vscode-cdn.net/home/lheinke/Downloads/~/.local/mambaforge/envs/pandexo/lib/python3.11/site-packages/pandexo/engine/justdoit.py:379), in run_pandexo(exo, inst, param_space, param_range, save_file, output_path, output_file, num_cores, verbose)
    [377](https://file+.vscode-resource.vscode-cdn.net/home/lheinke/Downloads/~/.local/mambaforge/envs/pandexo/lib/python3.11/site-packages/pandexo/engine/justdoit.py:377) if type(inst) == dict:
    [378](https://file+.vscode-resource.vscode-cdn.net/home/lheinke/Downloads/~/.local/mambaforge/envs/pandexo/lib/python3.11/site-packages/pandexo/engine/justdoit.py:378)     if verbose: print("Running Single Case w/ User Instrument Dict")
--> [379](https://file+.vscode-resource.vscode-cdn.net/home/lheinke/Downloads/~/.local/mambaforge/envs/pandexo/lib/python3.11/site-packages/pandexo/engine/justdoit.py:379)     results =wrapper({"pandeia_input": inst , "pandexo_input":exo}, verbose=verbose)
    [380](https://file+.vscode-resource.vscode-cdn.net/home/lheinke/Downloads/~/.local/mambaforge/envs/pandexo/lib/python3.11/site-packages/pandexo/engine/justdoit.py:380)     if output_file == '':
    [381](https://file+.vscode-resource.vscode-cdn.net/home/lheinke/Downloads/~/.local/mambaforge/envs/pandexo/lib/python3.11/site-packages/pandexo/engine/justdoit.py:381)         output_file = 'singlerun.p'

File [~/.local/mambaforge/envs/pandexo/lib/python3.11/site-packages/pandexo/engine/pandexo.py:39](https://file+.vscode-resource.vscode-cdn.net/home/lheinke/Downloads/~/.local/mambaforge/envs/pandexo/lib/python3.11/site-packages/pandexo/engine/pandexo.py:39), in wrapper(dictinput, verbose)
     [37](https://file+.vscode-resource.vscode-cdn.net/home/lheinke/Downloads/~/.local/mambaforge/envs/pandexo/lib/python3.11/site-packages/pandexo/engine/pandexo.py:37) elif telescope=='hst':
     [38](https://file+.vscode-resource.vscode-cdn.net/home/lheinke/Downloads/~/.local/mambaforge/envs/pandexo/lib/python3.11/site-packages/pandexo/engine/pandexo.py:38)     from .hst import compute_sim_hst
---> [39](https://file+.vscode-resource.vscode-cdn.net/home/lheinke/Downloads/~/.local/mambaforge/envs/pandexo/lib/python3.11/site-packages/pandexo/engine/pandexo.py:39)     return compute_sim_hst(dictinput, verbose=verbose)
     [40](https://file+.vscode-resource.vscode-cdn.net/home/lheinke/Downloads/~/.local/mambaforge/envs/pandexo/lib/python3.11/site-packages/pandexo/engine/pandexo.py:40) elif telescope=='wfirst':
     [41](https://file+.vscode-resource.vscode-cdn.net/home/lheinke/Downloads/~/.local/mambaforge/envs/pandexo/lib/python3.11/site-packages/pandexo/engine/pandexo.py:41)     return

File [~/.local/mambaforge/envs/pandexo/lib/python3.11/site-packages/pandexo/engine/hst.py:652](https://file+.vscode-resource.vscode-cdn.net/home/lheinke/Downloads/~/.local/mambaforge/envs/pandexo/lib/python3.11/site-packages/pandexo/engine/hst.py:652), in compute_sim_hst(dictinput, verbose)
    [649](https://file+.vscode-resource.vscode-cdn.net/home/lheinke/Downloads/~/.local/mambaforge/envs/pandexo/lib/python3.11/site-packages/pandexo/engine/hst.py:649) disperser = pandeia_input['configuration']['instrument']['disperser'].lower()
    [651](https://file+.vscode-resource.vscode-cdn.net/home/lheinke/Downloads/~/.local/mambaforge/envs/pandexo/lib/python3.11/site-packages/pandexo/engine/hst.py:651) # add a switch for ramp calculation
--> [652](https://file+.vscode-resource.vscode-cdn.net/home/lheinke/Downloads/~/.local/mambaforge/envs/pandexo/lib/python3.11/site-packages/pandexo/engine/hst.py:652) calRamp = pandeia_input['strategy']['calculateRamp']
    [654](https://file+.vscode-resource.vscode-cdn.net/home/lheinke/Downloads/~/.local/mambaforge/envs/pandexo/lib/python3.11/site-packages/pandexo/engine/hst.py:654) if not pandeia_input['strategy']['useFirstOrbit']:
    [655](https://file+.vscode-resource.vscode-cdn.net/home/lheinke/Downloads/~/.local/mambaforge/envs/pandexo/lib/python3.11/site-packages/pandexo/engine/hst.py:655)     if verbose:print("Dropping first orbit designed by observation strategy")

KeyError: 'calculateRamp'