code: def instantaneous_tpole(self):
          
          logger.debug("{}.instantaneous_tpole".format(self.component))
  
          if 'tpole' not in self.inst_vals.keys():
              logger.debug("{}.instantaneous_tpole COMPUTING".format(self.component))
  
              if self.mesh is None:
                  raise ValueError("mesh must be computed before determining tpole")
              # Convert from mean to polar by dividing total flux by gravity darkened flux (Ls drop out)
              # see PHOEBE Legacy scientific reference eq 5.20
              self.inst_vals['tpole'] = self.teff*(np.sum(self.mesh.areas) / np.sum(self.mesh.gravs.centers*self.mesh.areas))**(0.25)
  
          return self.inst_vals['tpole']
docs: compute the instantaenous temperature at the pole to achieve the mean
          effective temperature (teff) provided by the user