code: def show_hydrophobic(self):
          
          hydroph = self.plcomplex.hydrophobic_contacts
          if not len(hydroph.bs_ids) == 0:
              self.select_by_ids('Hydrophobic-P', hydroph.bs_ids, restrict=self.protname)
              self.select_by_ids('Hydrophobic-L', hydroph.lig_ids, restrict=self.ligname)
  
              for i in hydroph.pairs_ids:
                  cmd.select('tmp_bs', 'id %i & %s' % (i[0], self.protname))
                  cmd.select('tmp_lig', 'id %i & %s' % (i[1], self.ligname))
                  cmd.distance('Hydrophobic', 'tmp_bs', 'tmp_lig')
              if self.object_exists('Hydrophobic'):
                  cmd.set('dash_gap', 0.5, 'Hydrophobic')
                  cmd.set('dash_color', 'grey50', 'Hydrophobic')
          else:
              cmd.select('Hydrophobic-P', 'None')
docs: Visualizes hydrophobic contacts.