def METHOD_NAME( contraDC, simulation, plot = False):
    thruAmplitude = 10*np.log10(np.abs(contraDC.E_Thru[0,:])**2)
    dropAmplitude = 10*np.log10(np.abs(contraDC.E_Drop[0,:])**2)
    if plot == True:
        import matplotlib.pyplot as plt
        plt.figure()
        plt.plot(contraDC.wavelength[:]*1e9, thruAmplitude, label='Through Port')
        plt.plot(contraDC.wavelength[:]*1e9, dropAmplitude, label='Drop Port')
        plt.legend()
        plt.ylabel('Response (dB)')
        plt.xlabel('Wavelength (nm)')
    return [thruAmplitude, dropAmplitude]