astrofrog / astrofrog/sedfitter
Error
- Dominant language
- Python
- Stars
- 21
- Forks
- 23
- PR merge metrics
- No merged PRs in 30d
Description
i tried to run the first example from documantation and i got an Error:
```
------------------------------------------------------------
=> Model parameters
------------------------------------------------------------
Models : YSO models
Log[d] stepping : 0.01
Number of distances : 61
------------------------------------------------------------
=> Reading in convolved fluxes
------------------------------------------------------------
Reading models\models_r06/convolved/2J.fits
WARNING: VerifyWarning: Error validating header for HDU #2 (note: Astropy uses zero-based indexing).
'ascii' codec can't decode byte 0xf8 in position 2: ordinal not in range(128)
There may be extra bytes after the last HDU or the file is corrupted. [astropy.io.fits.hdu.hdulist]
---------------------------------------------------------------------------
ValueError Traceback (most recent call last)
in ()
21 'output.fitinfo',
22 extinction_law=extinction,
---> 23 distance_range=[1., 4.] * u.kpc)
C:\Users\Karolina\Desktop\sedfitter\sedfitter\fit.py in fit(data, filter_names, apertures, model_dir, output, n_data_min, extinction_law, av_range, distance_range, output_format, output_convolved, remove_resolved)
177 extinction_law=extinction_law, av_range=av_range,
178 distance_range=distance_range,
--> 179 remove_resolved=remove_resolved)
180
181 print(" ------------------------------------------------------------")
C:\Users\Karolina\Desktop\sedfitter\sedfitter\fit.py in __init__(self, filter_names, apertures, model_dir, extinction_law, av_range, distance_range, remove_resolved)
83
84 # Read in models
---> 85 self.models = Models.read(model_dir, self.filters, distance_range=distance_range, remove_resolved=remove_resolved)
86
87 # Add wavelength to filters
C:\Users\Karolina\Desktop\sedfitter\sedfitter\models.py in read(cls, directory, filters, distance_range, remove_resolved)
140 return cls._read_version_1(directory, filters,
141 distance_range=distance_range,
--> 142 remove_resolved=remove_resolved)
143 else:
144 return cls._read_version_2(directory, filters,
C:\Users\Karolina\Desktop\sedfitter\sedfitter\models.py in _read_version_1(cls, directory, filters, distance_range, remove_resolved)
196 print(" Reading " + filename)
197
--> 198 conv = ConvolvedFluxes.read(filename)
199
200 if ifilt == 0:
C:\Users\Karolina\Desktop\sedfitter\sedfitter\convolved_fluxes\convolved_fluxes.py in read(cls, filename)
228 conv.flux = tc['TOTAL_FLUX'].data.reshape(tc['TOTAL_FLUX'].shape[0], 1) * tc['TOTAL_FLUX'].unit
229 else:
--> 230 conv.flux = tc['TOTAL_FLUX'].data * tc['TOTAL_FLUX'].unit
231
232 if tc['TOTAL_FLUX_ERR'].ndim == 1 and conv.n_ap == 1:
C:\Users\Karolina\Desktop\sedfitter\sedfitter\convolved_fluxes\convolved_fluxes.py in flux(self, value)
114 self._flux = validate_array('flux', value, ndim=2,
115 shape=(self.n_models, self.n_ap),
--> 116 physical_type=('power', 'flux', 'spectral flux density'))
117
118 @property
C:\Users\Karolina\Desktop\sedfitter\sedfitter\utils\validator.py in validate_array(name, value, domain, ndim, shape, physical_type)
66 raise ValueError("{0} has incorrect length (expected {1} but found {2})".format(name, shape[0], value.shape[0]))
67 else:
---> 68 raise ValueError("{0} has incorrect shape (expected {1} but found {2})".format(name, shape, value.shape))
69
70 return value
ValueError: flux has incorrect shape (expected (200710, 1) but found (200710, 50))
```
Contributor guide
No contributing guide indexed for this repository
Research direction
Start by rerunning the documentation example and inspect the traceback locations in sedfitter/fit.py, sedfitter/models.py, sedfitter/convolved_fluxes/convolved_fluxes.py, and sedfitter/utils/validator.py. Investigate the warning for models/models_r06/convolved/2J.fits and the expected-versus-found flux shapes. Done means the documented example reads the model data without the reported exception, with a regression test if the relevant behavior can be isolated.
Written by the indexing model from the issue text.
Assessment
- Tech stack
- python
- Domain
- data
- Issue type
- Bug
- Difficulty
- 4/5
- Estimated time
- 3-5 days
- Activity status
- Stale
- Clarity
- Needs clarification
- Newbie friendliness
- 35/100