def test(
        array):
    numpy2pcr(pcraster.Scalar, array, mv=6)
  index0= min(mapStack.shape[0]-1,int(percentile*mapStack.shape[0]))
  index1= min(mapStack.shape[0]-1,int(percentile*mapStack.shape[0])+1)
  x0= float(index0)/mapStack.shape[0]
  x1= float(index1)/mapStack.shape[0]
  if x0 <> x1:
    y= mapStack[index0,:]+(percentile-x0)*\
       (mapStack[index1,:]-mapStack[index0,:])/(x1-x0)
  else:
    y= mapStack[index0,:]
  #-convert a slice of the stack into an array
  mapArray= np.ones((mapArray.shape))*MV
  for iCnt in xrange(coordinates.shape[0]):
    row= coordinates[iCnt,0]-1
    col= coordinates[iCnt,1]-1
    mapArray[row,col]= y[iCnt]
  pcr.report(numpy2pcr(pcr.Scalar,mapArray,MV),'q%03d_cumsec.map' % percent)
print
#-retrieve average recharge
R3AVG= pcr.scalar(0.)
for year in yearList:
  fileName= rootR3AVGFileName % ('%04d' % year)
  tempFile= open(fileName,'wb')
  tempFile.write(zipArchive.read(fileName))
  tempFile.close()  
  R3AVG+= pcr.readmap(fileName)
  os.remove(fileName)
R3AVG/= len(yearList)
pcr.report(R3AVG,rootR3AVGFileName % '')
print 'retrieving average runoff components'
for rootQSpecFileName in rootQSpecFileNames:
  specificRunoffQ= pcr.scalar(0.)