code: def wet_bulb_from_db_rh(db_temp, rh, b_press=101325):
      
      es = 6.112 * math.e**((17.67 * db_temp) / (db_temp + 243.5))
      e = (es * rh) / 100
      t_w = 0
      increse = 10.0
      previoussign = 1
      e_d = 1
      while math.fabs(e_d) > 0.005:
          e_wg = 6.112 * (math.e**((17.67 * t_w) / (t_w + 243.5)))
          eg = e_wg - (b_press/100) * (db_temp - t_w) * 0.00066 * (1 + (0.00155 * t_w))
          e_d = e - eg
          if e_d == 0:
              break
          else:
              if e_d < 0:
                  cursign = -1
                  if cursign != previoussign:
                      previoussign = cursign
                      increse = increse / 10
                  else:
                      increse = increse
              else:
                  cursign = 1
                  if cursign != previoussign:
                      previoussign = cursign
                      increse = increse/10
                  else:
                      increse = increse
          t_w = t_w + increse * previoussign
      return t_w
docs: Wet Bulb Temperature (C) at db_temp (C),
      Relative Humidity rh (%), and Pressure b_press (Pa).
  
      Note:
          [1] J. Sullivan and L. D. Sanders. "Method for obtaining wet-bulb temperatures by
          modifying the psychrometric formula." Center for Experiment Design and Data
          Analysis. NOAA - National Oceanic and Atmospheric Administration.
          http://www.srh.noaa.gov/epz/?n=wxcalc_rh