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
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.