Modifications to enable using the roman prism bandpass in the roman m…

…odule

Co-authored-by: Michael Troxel <[email protected]>

galsim/roman/__init__.py

@@ -65,7 +65,7 @@
 # F184, K213
 longwave_bands = ['F184', 'K213']
 # R062, Z087, Y106, J129, H158, W146
-shortwave_bands = ['R062', 'Z087', 'Y106', 'J129', 'H158', 'W146']
+shortwave_bands = ['R062', 'Z087', 'Y106', 'J129', 'H158', 'W146','SNPrism']
 
 stray_light_fraction = 0.1
 

galsim/roman/roman_bandpass.py

@@ -135,7 +135,7 @@ def getBandpasses(AB_zeropoint=True, default_thin_trunc=True, **kwargs):
     # Loop over the bands.
     for index, bp_name in enumerate(data.dtype.names[1:]):
         # Need to skip the prism and grism (not used for weak lensing imaging).
-        if bp_name=='SNPrism' or bp_name=='Grism_1stOrder' or bp_name=='Grism_0thOrder':
+        if bp_name=='Grism_1stOrder' or bp_name=='Grism_0thOrder':
             continue
 
         # Initialize the bandpass object.
@@ -154,9 +154,10 @@ def getBandpasses(AB_zeropoint=True, default_thin_trunc=True, **kwargs):
             bp = bp.withZeropoint('AB')
 
         # Store the sky level information as an attribute.
-        bp._ecliptic_lat = ecliptic_lat
-        bp._ecliptic_lon = ecliptic_lon
-        bp._sky_level = sky_data[2+index, :]
+        if bp_name!='SNPrism':
+            bp._ecliptic_lat = ecliptic_lat
+            bp._ecliptic_lon = ecliptic_lon
+            bp._sky_level = sky_data[2+index, :]
 
         # Add it to the dictionary.
         bp.name = bp_name if bp_name != 'W149' else 'W146'