Comments (9)
I was able to fix this by ingesting the model like
database.add_model(model=model_choice, teff_range=(700., 1300.), wavel_sampling=2000., wavel_range=(0.9,5))
from species.
Hm I am not sure what caused the warning to show up. I am not able reproduce it.
It is indeed anyway a good idea to use wavel_sampling
and wavel_range
(if the data allows) with a high-resolution grid, to limit the disk storage and speed up the interpolation at the start of FitModel
.
from species.
Hello Tomas and William! I get the same as William and thought it was due to the spacing of my data. (By the way, I ignored the warnings, especially since the deviation is small, and it seems ok.)
For what it is worth, I downloaded exo-rem-highres.tgz
a few days ago and it has (~20 s to compute, instead of re-downloading the file):
~$ md5sum exo-rem-highres.tgz
0149900c72bf44fb3342b0892e94d38d exo-rem-highres.tgz
@wbalmer, maybe you can compare against this and/or the timestamp on Tomas' website (https://home.strw.leidenuniv.nl/~stolker/species/):
[ ] exo-rem-highres.tgz 2023-07-23 15:14 16G
from species.
After some testing, I found that the issue was with the grid. I had used a smaller number of decimals to store the data compared to other grids.
The wavelength precision had therefore changed between the resampling of the grid (on my end) and storing the data in the text files that are used by species
.
I have created and new grid for 'exo-rem'
and 'exo-rem-highres'
. It is uploading but will take another hour since it is 22 GB...
from species.
Thanks for noticing that! Reopening this issue because the warning should not have occurred because all grid have been resampled to a constant lambda/Delta_lambda.
from species.
I think it could be for two reasons, either because my SPHERE spectrum has some nans (e.g. some upper limit values) or because it is extracted on an uneven wavelength grid that is more true to the instrument's wavelength calibration (the extraction comes from our SPHERE-DC friends). Probably the latter? I can send you the file if you'd like to use it to bug test.
from species.
Any NaN fluxes are removed from the spectra when using add_object
I think, so that should hopefully not cause a problem. In case you want to treat these wavelengths as upper limits in the fit, then you could also set the fluxes to zero and the uncertainty to the upper limit.
For the smoothing, the wavelength sampling of the data will not matter, only of the model spectra. The model spectra are first smoothed to the instrument resolution and then resampled to the data wavelengths. The latter is just a linear interpolation.
Perhaps you could try redownloading the grid in case you downloaded it a while ago and I may have fixed an issue that I forgot about. But the grid is very large so might not be worth it if it is running fine now.
from species.
Ah! right, I should have noticed (current version):
# Wavelength (um) - Flux (W m-2 um-1)
6.700000e-01 5.828044e+02
6.700166e-01 5.827946e+02
6.700332e-01 5.825590e+02
6.700498e-01 5.821541e+02
Two more decimal places for lambda should do it. The flux however has enough precision (1e-7 in relative terms).
A bigger issue is the non-monotonicity over Teff=1100–1200 K [spectral features change smoothly up to and down to those temperatures but not smoothly across that range], however, which is either physical and not caught by the Teff spacing, or numerical (which would be due to non-converged structures and/or spectra).
from species.
Ah, good to note! Thanks for updating this Tomas.
from species.
Related Issues (20)
- Upper limits on photometric measurements? HOT 1
- Change fontsize of axis labels with plot_spectrum? HOT 3
- Problem in add_custom_model with data_path HOT 1
- Tutorial "Fitting data with a grid of model spectra" (and a few other pages): small things HOT 5
- database overview: list_models() and verbosity control? HOT 4
- Database problem: OSError: Unable to open file (file is already open for read-only) HOT 4
- Making unpacking/storing of atmospheric models more memory-efficient HOT 10
- Adding dynesty support (especially to retrievals) HOT 8
- Various small things HOT 7
- Fitting with nested sampling (UltraNest or MultiNest) HOT 4
- Enhancing plot_spectrum a bit HOT 21
- Running multinest or ultranest in parallel? HOT 11
- Higher-order interpolation? HOT 3
- Retrieval with radial velocity / rotational broadening vsini HOT 2
- Installation error: No matching distribution found for matplotlib~=3.8.0 HOT 3
- Using wavel_range with database.add_model() HOT 5
- pypi package version HOT 1
- Different fsed for cloud species in retrieval HOT 4
- sonora elf-owl as successor to bobcat and cholla HOT 18
Recommend Projects
-
React
A declarative, efficient, and flexible JavaScript library for building user interfaces.
-
Vue.js
🖖 Vue.js is a progressive, incrementally-adoptable JavaScript framework for building UI on the web.
-
Typescript
TypeScript is a superset of JavaScript that compiles to clean JavaScript output.
-
TensorFlow
An Open Source Machine Learning Framework for Everyone
-
Django
The Web framework for perfectionists with deadlines.
-
Laravel
A PHP framework for web artisans
-
D3
Bring data to life with SVG, Canvas and HTML. 📊📈🎉
-
Recommend Topics
-
javascript
JavaScript (JS) is a lightweight interpreted programming language with first-class functions.
-
web
Some thing interesting about web. New door for the world.
-
server
A server is a program made to process requests and deliver data to clients.
-
Machine learning
Machine learning is a way of modeling and interpreting data that allows a piece of software to respond intelligently.
-
Visualization
Some thing interesting about visualization, use data art
-
Game
Some thing interesting about game, make everyone happy.
Recommend Org
-
Facebook
We are working to build community through open source technology. NB: members must have two-factor auth.
-
Microsoft
Open source projects and samples from Microsoft.
-
Google
Google ❤️ Open Source for everyone.
-
Alibaba
Alibaba Open Source for everyone
-
D3
Data-Driven Documents codes.
-
Tencent
China tencent open source team.
from species.