Script to calculate bolometric peak