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