docs: def fishrot(k=20, n=100, dec=0, inc=90, di_block=True):
      """
      Generates Fisher distributed unit vectors from a specified distribution
      using the pmag.py fshdev and dodirot functions.
  
      Parameters
      ----------
      k : kappa precision parameter (default is 20)
      n : number of vectors to determine (default is 100)
      dec : mean declination of distribution (default is 0)
      inc : mean inclination of distribution (default is 90)
      di_block : this function returns a nested list of [dec,inc,1.0] as the default
      if di_block = False it will return a list of dec and a list of inc
  
      Returns
      ---------
      di_block : a nested list of [dec,inc,1.0] (default)
      dec, inc : a list of dec and a list of inc (if di_block = False)
  
      Examples
      --------
      >>> ipmag.fishrot(k=20, n=5, dec=40, inc=60)
      [[44.766285502555775, 37.440866867657235, 1.0],
       [33.866315796883725, 64.732532250463436, 1.0],
       [47.002912770597163, 54.317853800896977, 1.0],
       [36.762165614432547, 56.857240672884252, 1.0],
       [71.43950604474395, 59.825830945715431, 1.0]]
      """
      directions = []
      declinations = []
      inclinations = []
      if di_block == True:
          for data in range(n):
              d, i = pmag.fshdev(k)
              drot, irot = pmag.dodirot(d, i, dec, inc)
              directions.append([drot, irot, 1.])
          return directions
      else:
          for data in range(n):
              d, i = pmag.fshdev(k)
              drot, irot = pmag.dodirot(d, i, dec, inc)
              declinations.append(drot)
              inclinations.append(irot)
          return declinations, inclinations