def correlation(location, independentName, dependentName, locationName, sampleNumbers, timeSteps):
location=boolean(location)
name=independentName + '_' + dependentName + '_' + locationName
tssFileIntercept = file("%s%s.tss" % (name,'_int'), "w")
tssFileSlope = file("%s%s.tss" % (name,'_slope'), "w")
tssFileRSquared = file("%s%s.tss" % (name,'_rSq'), "w")
for step in timeSteps:
values=[]
for sample in sampleNumbers:
smallValue=0.0000000000000000001
fileNameOne=generateNameST(dependentName,sample,step)
valueOne=generalfunctions.getCellValueAtBooleanLocation(location,scalar(fileNameOne))
pairList=[valueOne + smallValue]
fileNameTwo=generateNameST(independentName,sample,step)
valueTwo=generalfunctions.getCellValueAtBooleanLocation(location,scalar(fileNameTwo))
pairList.append(valueTwo + smallValue)
values.append(pairList)
#print valueOne + smallValue, valueTwo + smallValue
reg=regression.linearRegression(values, 1)
rSq=regression.linearRSquared(values,reg)
#print step, reg, rSq
tssFileIntercept.write("%d %g\n" % (step, reg[0]))
tssFileSlope.write("%d %g\n" % (step, reg[1]))
tssFileRSquared.write("%d %g\n" % (step, rSq))
tssFileIntercept.close()
tssFileSlope.close()
tssFileRSquared.close()
def correlation(location, independentName, dependentName, locationName, sampleNumbers, timeSteps):
location=boolean(location)
name=independentName + '_' + dependentName + '_' + locationName
tssFileIntercept = file("%s%s.tss" % (name,'_int'), "w")
tssFileSlope = file("%s%s.tss" % (name,'_slope'), "w")
tssFileRSquared = file("%s%s.tss" % (name,'_rSq'), "w")
for step in timeSteps:
values=[]
for sample in sampleNumbers:
smallValue=0.0000000000000000001
fileNameOne=generateNameST(dependentName,sample,step)
valueOne=generalfunctions.getCellValueAtBooleanLocation(location,scalar(fileNameOne))
pairList=[valueOne + smallValue]
fileNameTwo=generateNameST(independentName,sample,step)
valueTwo=generalfunctions.getCellValueAtBooleanLocation(location,scalar(fileNameTwo))
pairList.append(valueTwo + smallValue)
values.append(pairList)
#print valueOne + smallValue, valueTwo + smallValue
reg=regression.linearRegression(values, 1)
rSq=regression.linearRSquared(values,reg)
#print step, reg, rSq
tssFileIntercept.write("%d %g\n" % (step, reg[0]))
tssFileSlope.write("%d %g\n" % (step, reg[1]))
tssFileRSquared.write("%d %g\n" % (step, rSq))
tssFileIntercept.close()
tssFileSlope.close()
tssFileRSquared.close()
def createtimeseries(names, nameExtension, locations,sampleNumbers,timeSteps):
if deterministicInput(sampleNumbers):
for name in names:
tssFile = file(name + nameExtension + '.tss', "w")
tssFile.write("timeseries scalar\n")
tssFile.write("2\n")
tssFile.write("timestep\n")
tssFile.write("%s\n" % (name))
for step in timeSteps:
timeseriesValue=mapmaximum(ifthen(locations,generateNameT(name,step)))
value, valid = cellvalue(timeseriesValue, 1, 1); assert valid
tssFile.write("%d %g\n" % (step, value))
tssFile.close()
else:
for name in names:
for sample in sampleNumbers:
tssFile = file(generateNameS("%s%s.tss" % (name,nameExtension), sample), "w")
tssFile.write("timeseries scalar\n")
tssFile.write("2\n")
tssFile.write("timestep\n")
tssFile.write("%s\n" % (name))
for step in timeSteps:
timeseriesValue=mapmaximum(ifthen(locations,generateNameST(name,sample,step)))
value, valid = cellvalue(timeseriesValue, 1, 1); assert valid
tssFile.write("%d %g\n" % (step, value))
tssFile.close()
def createtimeseries(names, nameExtension, locations,sampleNumbers,timeSteps):
if deterministicInput(sampleNumbers):
for name in names:
tssFile = file(name + nameExtension + '.tss', "w")
tssFile.write("timeseries scalar\n")
tssFile.write("2\n")
tssFile.write("timestep\n")
tssFile.write("%s\n" % (name))
for step in timeSteps:
timeseriesValue=mapmaximum(ifthen(locations,generateNameT(name,step)))
value, valid = cellvalue(timeseriesValue, 1, 1); assert valid
tssFile.write("%d %g\n" % (step, value))
tssFile.close()
else:
for name in names:
for sample in sampleNumbers:
tssFile = file(generateNameS("%s%s.tss" % (name,nameExtension), sample), "w")
tssFile.write("timeseries scalar\n")
tssFile.write("2\n")
tssFile.write("timestep\n")
tssFile.write("%s\n" % (name))
for step in timeSteps:
timeseriesValue=mapmaximum(ifthen(locations,generateNameST(name,sample,step)))
value, valid = cellvalue(timeseriesValue, 1, 1); assert valid
tssFile.write("%d %g\n" % (step, value))
tssFile.close()
def createAllGstatRealizations(nameCommandFile, nameOutMapList,
nrRealPerGstatCall, sampleNumbers, timeSteps):
# create list names with filenames required
names = []
names.append([])
i = 0
j = 0
for sample in sampleNumbers:
for step in timeSteps:
if j == nrRealPerGstatCall:
j = 1
i = i + 1
names.append([])
else:
j = j + 1
namesOneSampleOneTimeStep = []
for nameOutMap in nameOutMapList:
if staticInput(timeSteps):
fileName = generateNameS(nameOutMap, sample) + '.map'
else:
fileName = generateNameST(nameOutMap, sample, step)
namesOneSampleOneTimeStep.append(fileName)
names[i].append(namesOneSampleOneTimeStep)
for setOfRealizations in names:
createGstatRealizations(setOfRealizations, nameCommandFile,
nameOutMapList)
def createAllGstatRealizations(nameCommandFile,nameOutMapList,nrRealPerGstatCall,sampleNumbers,timeSteps):
# create list names with filenames required
names=[]
names.append([])
i = 0
j = 0
for sample in sampleNumbers:
for step in timeSteps:
if j == nrRealPerGstatCall:
j = 1
i = i+1
names.append([])
else:
j = j+1
namesOneSampleOneTimeStep=[]
for nameOutMap in nameOutMapList:
if staticInput(timeSteps):
fileName=generateNameS(nameOutMap,sample) + '.map'
else:
fileName=generateNameST(nameOutMap,sample,step)
namesOneSampleOneTimeStep.append(fileName)
names[i].append(namesOneSampleOneTimeStep)
for setOfRealizations in names:
createGstatRealizations(setOfRealizations,nameCommandFile, nameOutMapList)