code: def _D2O_Sublimation_Pressure(T):
      
      if 210 <= T <= 276.969:
          Tita = T/276.969
          suma = 0
          ai = [-0.1314226e2, 0.3212969e2]
          ti = [-1.73, -1.42]
          for a, t in zip(ai, ti):
              suma += a*(1-Tita**t)
          return exp(suma)*0.00066159
      else:
          raise NotImplementedError("Incoming out of bound")
docs: Sublimation Pressure correlation for heavy water
  
      Parameters
      ----------
      T : float
          Temperature, [K]
  
      Returns
      -------
      P : float
          Pressure at sublimation line, [MPa]
  
      Notes
      ------
      Raise :class:`NotImplementedError` if input isn't in limit:
  
          * 210 ≤ T ≤ 276.969
  
      Examples
      --------
      >>> _Sublimation_Pressure(245)
      3.27390934e-5
  
      References
      ----------
      IAPWS, Revised Release on the IAPWS Formulation 2017 for the Thermodynamic
      Properties of Heavy Water, http://www.iapws.org/relguide/Heavy.html.