Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

New COMPSP #7

Merged
merged 4 commits into from
May 18, 2016
Merged

New COMPSP #7

merged 4 commits into from
May 18, 2016

Conversation

bd-j
Copy link
Collaborator

@bd-j bd-j commented May 18, 2016

Complete rewrite of COMPSP to do the convolution of the SSPs with SFH using analytic integrals. This means that time_res_inc never needs to be larger than 1 to get accurate spectra. The call signature of the new COMSP is unchanged from the old version. Some variables were added to sps_vars, and a gitignore is added. Behavior has changed slighty in some respects, the most important are

  • nzin > 1 raises an error in compsp. Only single metallicity SFHs are possible.
  • sf_trunc < sf_start causes a warning and sf_trunc to be ignored, instead of giving a zero spectrum.
  • fburst is applied at the same time as const instead of afterwards, so fburst > 0 does not alter the effective const. Related to this, an error is raised if fburst + const > 1.
  • For tage > 0, if tburst > tage then the burst is totally ignored. In old compsp the fluxes and masses were still decreased by a factor (1-fburst) in this case.
  • For tage=0 the returned SFRs integrate to (1.0-fburst) instead of 1.0. That is, the user must add the burst SFR themselves however they want to do that, but the sfrs for the other components are correct for the given stellar masses, etc. M/L and stellar mass history unchanged.
  • Light weighted ages not implemented.
  • Behavior for any age greater than the maximum isochrone age is not defined.

bd-j added 4 commits May 17, 2016 10:06
This commit removes the old compsp and replaces it with a subroutine that performs analytic integrals of the interpolation weights, to increase accuracy for certain SFHs, without increasing time_res_inc.  It also significantly increases the modularity of compsp. intspec and intsfr have been removed, many new functions and subroutines have been added.  However, the API of compsp remains unchanged.  This commit also adds a gitignore, sets time_res_inc=1 and changes some other sps_vars.
This commit implements a tabular SFH option in csp_gen. It moves add_dust from compsp to csp_gen for speed reasons. Treatment of lbol is fixed.  Warnings are added to compsp for certain parameter inputs.  Numerous bugs are fixed in csp_gen and sfhinfo that affected the logic for relative component normalizations.
This commit includes fixes for many bugs in the sfhinfo logic affecting relative normalizations of different components of the SFH.  It also corrects the slope and formed mass calculations used in the tabular SFH portions of csp_gen. sfhinfo returns interpolated sfrs for tabulated SFH.
@cconroy20 cconroy20 merged commit ea39c6e into cconroy20:master May 18, 2016
@benjaminrose
Copy link

Looking at fburst and const, what happens if both are set to 0. I would assume that no SF would happen using the tau-model. But these lines, https://github.com/cconroy20/fsps/blob/master/src/sfhinfo.f90#L109 https://github.com/cconroy20/fsps/blob/master/src/sfhinfo.f90#L119, confuse me.

@cconroy20
Copy link
Owner

cconroy20 commented Dec 12, 2016 via email

@benjaminrose
Copy link

Then I don't understand fburst. I thought that was the fraction of mass formed in the tau-model burst.

@bd-j
Copy link
Collaborator Author

bd-j commented Dec 12, 2016

No, it's the fraction of mass formed in an instantaneous burst at time tburst. The tau-model need not be a burst (since tau can be large)

@cconroy20
Copy link
Owner

cconroy20 commented Dec 12, 2016 via email

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

Successfully merging this pull request may close these issues.

3 participants