import numpy as np import os vc_A = 2.99792458e+18 # speed of light: A/s vc_m = 2.99792458e+8 # speed of light: m/s h_Plank = 6.626196e-27 # Plank constant: erg s def photonEnergy(lambd): """ The energy of photon at a given wavelength Parameter: lambd: the wavelength in unit of Angstrom Return: eph: energy of photon in unit of erg """ nu = vc_A / lambd eph = h_Plank * nu return eph