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