code: def _fill_loggs(self, mesh=None, ignore_effects=False):
          
          logger.debug("{}._fill_loggs".format(self.component))
  
          if mesh is None:
              mesh = self.mesh
  
          loggs = np.log10(mesh.normgrads.for_computations * g_rel_to_abs(self.masses[self.ind_self], self.sma))
  
          if not ignore_effects:
              for feature in self.features:
                  if feature.proto_coords:
                      loggs = feature.process_loggs(loggs, mesh.roche_coords_for_computations, s=self.polar_direction_xyz, t=self.time)
                  else:
                      loggs = feature.process_loggs(loggs, mesh.coords_for_computations, s=self.polar_direction_xyz, t=self.time)
  
          mesh.update_columns(loggs=loggs)
  
          if not self.needs_recompute_instantaneous:
              logger.debug("{}._fill_loggs: copying loggs to standard mesh".format(self.component))
              theta = 0.0
              self._standard_meshes[theta].update_columns(loggs=loggs)
docs: TODO: add documentation
  
          Calculate local surface gravity
  
          GMSunNom = 1.3271244e20 m**3 s**-2
          RSunNom = 6.597e8 m