docs: def InitDetectorNetwork(numDets, detNetwork, LIGO3FLAG=0):
      """ InitDetectorNetwork - function to initialise desired detector network
          specified as input.
  
  	numDets - Number of IFOs in the detector network to be considered, can
          take any value in the range [1,5].
          detNetwork - 5 character string to signify which detector network should
          be considered.  1st, 2nd, 3rd, 4th and 5th characters represent H1, L1,
          V1, I1 and K1 detectors respectively.
  	LIGO3FLAG - Set to 1 to use LIGO3 PSD instead of aLIGO PSD for H1, L1
          and I1 detectors.
  
  	Returns Noise - (N_fd,numDets) array containing simulated Gaussian detector
  			noise for each detector in the network.
  		PSD - (N_fd,numDets) array containing the noise PSD for each
  		      detector in the network.
  		detNames - (numDets) list containing names of detectors in the
  			   network.
  		detPosVector - (3,numDets) array containing position vector of
  			       detector vertex for each detector in the network
  			       in geocentric coordinates.
  		detRespTensor - (3,3,numDets) array containing detector response
  				tensor for each detector in the network. 
  
          Sarah Gossan 2012. Last updated 02/18/14. """
  
      # Create arrays for output
      Noise = np.zeros((N_fd,numDets),complex)
      PSD = np.zeros((N_fd,numDets),float)
      detPosVector = np.zeros((3,numDets),float)
      detRespTensor = np.zeros((3,3,numDets),float)
  
      # Create detector names array
      detNames = []
      if detNetwork[0] == 'H':
  	detNames.append('H1')
      if detNetwork[1] == 'L':
  	detNames.append('L1')
      if detNetwork[2] == 'V':
  	detNames.append('V1')
      if detNetwork[3] == 'I':
  	detNames.append('I1')
      if detNetwork[4] == 'K':
  	detNames.append('K1')
  
      # Get detector PSD, Noise, position vector, response tensor and time delay
      # from geocenter dependent on working domain
      for detIdx in range(0,numDets):
  	# PSDs
  	PSD[:,detIdx] = GetDetectorPSD(detNames[detIdx])
          Noise[:,detIdx] = GenerateGaussianNoise(PSD[:,detIdx])
         	# Position vector and response tensor
  	[detRespTensor[:,:,detIdx],detPosVector[:,detIdx]] = \
  	GetDetectorData(detNames[detIdx])
  
      return Noise,PSD,detNames,detPosVector,detRespTensor