code: def dodirot_V(di_block, Dbar, Ibar):
      
      N = di_block.shape[0]
      DipDir, Dip = np.ones(N, dtype=np.float).transpose(
      )*(Dbar-180.), np.ones(N, dtype=np.float).transpose()*(90.-Ibar)
      di_block = di_block.transpose()
      data = np.array([di_block[0], di_block[1], DipDir, Dip]).transpose()
      drot, irot = dotilt_V(data)
      drot = (drot-180.) % 360.  #
      return np.column_stack((drot, irot))
docs: Rotate an array of dec/inc pairs to coordinate system with Dec,Inc as 0,90
  
      Parameters
      ___________________
      di_block : array of [[Dec1,Inc1],[Dec2,Inc2],....]
      Dbar : declination of desired center
      Ibar : inclination of desired center
  
      Returns
      __________
      array of rotated decs and incs: [[rot_Dec1,rot_Inc1],[rot_Dec2,rot_Inc2],....]