search

MingjieJian / pymoog

Python wrapper for MOOG
MIT License
11 stars 4 forks source link

The wavelength of the synthetic spectra is not equidistant after16384.02A #44

Closed MingjieJian closed 1 year ago

MingjieJian commented 2 years ago

The wavelength of the synthetic spectra is not equidistant after16384.02A, as shown in the following figure: image

Such in-equidistant wavelength is not related with linelist, wav_start nor del_wav. I suspect that there is some bug in MOOG.

MingjieJian commented 2 years ago

The wavelength array is calculated in lines 429-439 of Smooth.f and stored in the xsyn. The type of xsyn is defined in Equiv.com as real4. Such precision is not enough for wavelength larger than ~16000A, so the wavelength grid has some problems. Changing the type form real4 to real*8 can solve the problem.

This bug has been solved in version 0.0.28.