docs: def moist_amplification_factor(Tkelvin, relative_humidity=0.8):
      '''Compute the moisture amplification factor for the moist diffusivity
      given relative humidity and reference temperature profile.'''
      deltaT = 0.01
      #  slope of saturation specific humidity at 1000 hPa
      dqsdTs = (qsat(Tkelvin+deltaT/2, 1000.) - qsat(Tkelvin-deltaT/2, 1000.)) / deltaT
      return const.Lhvap / const.cp * relative_humidity * dqsdTs