Maybe should also allow for multiplying by another function?  I'm imagining someone modeling dust extinction.  It might be convenient to have dust be a function of wavelength and then subtracting different amounts of extinction from different galaxies:
```
gal_sed -= extict * dust_func * init_gal_sed
```