def moist_amplification_factor(Tkelvin, relative_humidity=<NUM_LIT>):
deltaT = <NUM_LIT><EOL>dqsdTs = (qsat(Tkelvin+deltaT/<NUM_LIT:2>, <NUM_LIT>) - qsat(Tkelvin-deltaT/<NUM_LIT:2>, <NUM_LIT>)) / deltaT<EOL>return const.Lhvap / const.cp * relative_humidity * dqsdTs<EOL>
Compute the moisture amplification factor for the moist diffusivity
    given relative humidity and reference temperature profile.