BC03 and CB13 templates

[steveschulze]

Hi,

I would like to use the BC03 templates and compare the results to the CB13 templates.

The IAP FTP server is not accessible at the moment. The error message is "ftp.iap.fr refused to connect". Is this error due to maintenance or is there something else going on? (I managed to download the BC03 templates from a different location.)

I have a problem with using the BC03 templates in BEAGLE. I commented the TEMPLATES NEBULAR option and set the TEMPLATES = $BEAGLE_TEMPLATES/bc03/bc03_miles_chab_spectra. BEAGLE returned

ERROR while opening file: /BEAGLE_virtual_workdir/templates/bc03/bc03_miles_chab_spectra, error_msg: Cannot open file '/BEAGLE_virtual_workdir/templates/bc03/bc03_miles_chab_spectra': Is a directory ]

Specifying a particular ised file did not help, e.g. TEMPLATES = $BEAGLE_TEMPLATES/bc03/bc03_miles_chab_spectra/bc2003_hr_xmiless_m22_chab_ssp.ised

The new error message is:

At line 392 of file /root/src/github.com/jacopo-chevallard/BEAGLE/build/src/iotemplates.f90
Fortran runtime error: Bad real number in item 1 of list input

Error termination. Backtrace:
#0  0x7f4b0e9f9ef7 in ???
#1  0x7f4b0e9faa95 in ???
#2  0x7f4b0e9fb24c in ???
#3  0x7f4b0eabeba2 in ???
#4  0x7f4b0eac1f66 in ???
#5  0x40f9e7 in __iotemplates_MOD_read_templates_age_z
    at /root/src/github.com/jacopo-chevallard/BEAGLE/build/src/iotemplates.f90:392
#6  0x413dd5 in __iotemplates_MOD_initialize_templates
    at /root/src/github.com/jacopo-chevallard/BEAGLE/build/src/iotemplates.f90:265
#7  0x406c6c in beagle
    at /root/src/github.com/jacopo-chevallard/BEAGLE/src/BEAGLE.f90:321
#8  0x404a44 in main
    at /root/src/github.com/jacopo-chevallard/BEAGLE/src/BEAGLE.f90:6

Could you help me to use to the BC03 templates? Does one always need to choose BC03 template with a specific metallicity? Are the output parameters identical to the CB13 templates?

Is there a paper where the differences between the BC03 and CB13/CB15 templates are discussed? I remember that there was a controversial debate about the contribution of TP-AGB stars.

Cheers,

Steve

[jacopo-chevallard]

Hi Steve

Could you help me to use to the BC03 templates? Does one always need to choose BC03 template with a specific metallicity? Are the output parameters identical to the CB13 templates?

If you look at the Beagle manual (https://www.overleaf.com/read/wkhmzbcwmkqz), Section 3.3, you'll see how to use different *ised templates in Beagle. (If it's still unclear feel free to ask here and/or suggest improvements to the manual)

Is there a paper where the differences between the BC03 and CB13/CB15 templates are discussed?

There isn't yet a full paper discussing BC03 vs newer CBxx templates, but perhaps @scharlot can suggest you some papers in which newer versions of the BC03 models have been tested wrt TP-AGB stars (I remember Zibetti+2013, I don't know if there is anything more recent)

[steveschulze]

Hi Jacopo,

Sorry, I didn't properly read the manual. I still have a problem...

For instance, setting the template to

TEMPLATES = $BEAGLE_TEMPLATES/bc03/bc03_miles_chab_spectra/bc2003_hr_xmiless_m32_chab_ssp.ised,

returns

docker run --rm -it -v $HOME/BEAGLE-general:/BEAGLE_virtual_workdir --env-file $HOME/BEAGLE-general/env.list beagletool/beagle:0.18.9 --parameter-file /BEAGLE_virtual_workdir/params/fit_photometry_ptf_bc03.param --fit

 Reading the file containing the templates (*.ised) file list: /BEAGLE_virtual_workdir/templates/bc03/bc03_miles_chab_spectra/bc2003_hr_xmiless_m32_chab_ssp.ised

At line 392 of file /root/src/github.com/jacopo-chevallard/BEAGLE/build/src/iotemplates.f90
Fortran runtime error: Bad real number in item 1 of list input

Error termination. Backtrace:
#0  0x7f5d689ceef7 in ???
#1  0x7f5d689cfa95 in ???
#2  0x7f5d689d024c in ???
#3  0x7f5d68a93ba2 in ???
#4  0x7f5d68a96f66 in ???
#5  0x40fa37 in __iotemplates_MOD_read_templates_age_z
    at /root/src/github.com/jacopo-chevallard/BEAGLE/build/src/iotemplates.f90:392
#6  0x413e25 in __iotemplates_MOD_initialize_templates
    at /root/src/github.com/jacopo-chevallard/BEAGLE/build/src/iotemplates.f90:265
#7  0x406c6c in beagle
    at /root/src/github.com/jacopo-chevallard/BEAGLE/src/BEAGLE.f90:322
#8  0x404a44 in main
    at /root/src/github.com/jacopo-chevallard/BEAGLE/src/BEAGLE.f90:6

Cheers,

Steve

[steveschulze]

Thanks for the reference to Zibetti. I will look into the paper and write to Stéphane.

Cheers,

Steve

[jacopo-chevallard]

Hi Steve, the TEMPLATES keyword refers to an text (ASCII) file containing a list of the different *ised files, not the ised files themselves. You may just use the example linked in the manual, i.e. ftp://ftp.iap.fr/pub/from_users/chevalla/SSPs/bc03/bc2003_hr_xmiless_chab_ssp

[steveschulze]

Hi Jacopo,

Upps, I overread that. I downloaded the files and stored the file names in a text file. It still refuses to process the files:

Reading the file containing the templates (*.ised) file list: /BEAGLE_virtual_workdir/templates/bc03/bc03_hr_xmiles

At line 392 of file /root/src/github.com/jacopo-chevallard/BEAGLE/build/src/iotemplates.f90
Fortran runtime error: End of file

Error termination. Backtrace:
#0  0x7f7f28f61ef7 in ???
#1  0x7f7f28f62a95 in ???
#2  0x7f7f28f6324c in ???
#3  0x7f7f2902f323 in ???
#4  0x7f7f29029eb9 in ???
#5  0x40f9e7 in __iotemplates_MOD_read_templates_age_z
    at /root/src/github.com/jacopo-chevallard/BEAGLE/build/src/iotemplates.f90:392
#6  0x413dd5 in __iotemplates_MOD_initialize_templates
    at /root/src/github.com/jacopo-chevallard/BEAGLE/build/src/iotemplates.f90:265
#7  0x406c6c in beagle
    at /root/src/github.com/jacopo-chevallard/BEAGLE/src/BEAGLE.f90:321
#8  0x404a44 in main
    at /root/src/github.com/jacopo-chevallard/BEAGLE/src/BEAGLE.f90:6

Cheers,

Steve

[jacopo-chevallard]

Hi Steve, can you let me know which Beagle version are you using? Can you also attach here (in an archive) the parameter file and the file containing the template list?

[steveschulze]

Hi Jacopo,

Thanks a lot. I am using version 0.18.7. The attached archive also contains the filter parameter files and the data file.

Cheers,

Steve

BEAGLE_BC03.zip

[jacopo-chevallard]

Hi Steve, in your file you're missing the column indicating the metallicities, i.e. your file loos like

bc2003_hr_xmiless_m22_chab_ssp.ised
bc2003_hr_xmiless_m32_chab_ssp.ised
bc2003_hr_xmiless_m42_chab_ssp.ised
bc2003_hr_xmiless_m52_chab_ssp.ised
bc2003_hr_xmiless_m62_chab_ssp.ised
bc2003_hr_xmiless_m72_chab_ssp.ised
bc2003_hr_xmiless_m82_chab_ssp.ised

while it should look like

bc2003_hr_xmiless_m22_chab_ssp.ised      0.0001
bc2003_hr_xmiless_m32_chab_ssp.ised      0.0004
bc2003_hr_xmiless_m42_chab_ssp.ised      0.004
bc2003_hr_xmiless_m52_chab_ssp.ised      0.008
bc2003_hr_xmiless_m62_chab_ssp.ised      0.02
bc2003_hr_xmiless_m72_chab_ssp.ised      0.05
bc2003_hr_xmiless_m82_chab_ssp.ised      0.1

[steveschulze]

Blimey. How could I overlook that! I feel so embarrassed. Thanks for solving that.

Cheers,

Steve

[jacopo-chevallard]

No problem Steve ! :) let's see if now everything works fine !

[steveschulze]

One step forward but still a tango. Do I have to set a further option in the parameter file. When I try to fit an SED with the BC03 template, I get a segmentation error. If I comment nebular_logU, BEAGLE does the fit but Multinest does not stop to run.

Cheers,

Steve

[jacopo-chevallard]

Indeed, all the nebular_ parameters require the use of models including nebular emission, this is why you get seg fault (I should insert a more explicit error check, feel free to open a separate issue so that I can keep track of it). With BC03 models you should not use any nebular_ parameter

[steveschulze]

OK. I toggled them off. MultiNest is still not stopping although ln(Z) is barely changing and the acceptance rate is < 7%. With the CB15 templates, the fit would be have been finished.

Is there a possibility to connect the models for nebular emission with the BC03 templates?

[jacopo-chevallard]

Would you mind sending me the log of the fitting?

Is there a possibility to connect the models for nebular emission with the BC03 templates?

Not really, since the nebular emission is tightly coupled to the ionizing emission of the stellar populations used, and therefore to the specific version of the Bruzual & Charlot stellar population synthesis code.

[steveschulze]

The log is in the folder "log". The archive also contains an example in case you want to run it.

BEAGLE_ZTF18abukavn.zip

Thanks for the clarification on the emission-line model.

[steveschulze]

Thanks for running a fitting request on GAZPAR. I used the same parameter file on my computer and MultiNest converged. I can't explain why it didn't work before.

[jacopo-chevallard]

Cool ! Happy that it worked !

[jacopo-chevallard]

(In the README file that I included in the GAZPAR output I reported a few problems that I found in your settings)