search

cconroy20 / fsps

Flexible Stellar Population Synthesis
MIT License
84 stars 48 forks source link

Why is nzinit for BPASS set to 1 and not 12 in sps_setup.py ? #47

Closed prerakgarg07 closed 3 years ago

prerakgarg07 commented 3 years ago

I was wondering why is the nzinit for BPASS set to 1 by default (line 232 in sps_setup.py ).

Is that intended? Should it not be set to the length of zlegend.dat which in this case 12.

bd-j commented 3 years ago

I think you mean sps_vars.f90? Anyway the metallicity grid is set by the isochrones. For other isochrone libraries, the spectra are constructed from independent stellar libraries that must be interpolated to the isochrone metallicities, and nzinit describes these stellar libraries. As described in the manual (S1.16) the BPASS models are included as complete SSPs, so the the only relevant variable, used in sps_setup.f90 is the number of isochrone metallicities, which is set at line 224 of sps_vars.f90.

prerakgarg07 commented 3 years ago

Thanks a lot for explaining this. I understand the reasoning now and yes I meant sps_vars.f90.