code: def trim_occluded_throats(network, mask='all'):
      r
      occluded_ts = network['throat.area'] == 0
      if sp.sum(occluded_ts) > 0:
          occluded_ts *= network["throat."+mask]
          trim(network=network, throats=occluded_ts)
docs: r"""
      Remove throats with zero area from the network and also remove
      pores that are isolated (as a result or otherwise)
  
      Parameters
      ----------
      network : OpenPNM Network Object
  
      mask : string
          Applies routine only to pores and throats with this label