def main():
    ff = evaluate.FitnessFunction()
    base = tools.load_nolan_bedrock()
    surf = tools.load_first_guess_surface()
    for i in range(9):
        mp.plot(range(58), base, "red")
        params = create(58, -500, 500)
        bed = map(operator.add, params, base)
        mp.plot(range(58), bed, "blue")
        run1 = basic_model.isothermalISM(
            58, 1000, 0.0015, 0.0005, 0.00022, bed[:]
        )  # 55 nodes, 1000-meter spacing,  basal slip was .0005

        for i in range(1500):  # 5000 years
            run1.timestep(1)
            if i % 100 == 0:
                print "on timestep", i
        mp.plot(range(58), run1.get_surface_elev()[:58])
        mp.plot(range(57), surf, "red")
        mp.scatter([1, 5, 12, 20, 32], [-70, -200, -500, -600, -250])
        print "surface difference:", tools.calculate_surface_difference(surf, run1.get_surface_elev()[:58])
        print "evaluate surface difference:", ff.evaluate(bed, run1.get_surface_elev())
        mp.show()
import numpy as np

possibilities1 = [ .0008, .00082, .00085, .0009, .00095, .00097, .00099, .001, .0012, .0015]
possibilities2 = [ .0005, .0006, .0007,.00072, .00075, .00077, .0008, .00082, .00085, .0009]
possibilities3 = [.0001, .0002, .00022, .00024, .00026, .00028, .0003, .0004, .0005, .00055]
results = np.zeros([1000, 4])
count = 0
b0 = tools.load_nolan_bedrock()
obs_surface = tools.load_first_guess_surface()
for i in possibilities1:
	for j in possibilities2:
		for k in possibilities3:
			run = basic_model.isothermalISM(58, 1000, i,j,k, b0[:])
			for n in range(1000):
				run.timestep(1)
			difference = tools.calculate_surface_difference(run.get_surface_elev(), obs_surface)
			results[count][0] = i
			results[count][1] = j
			results[count][2] = k
			results[count][3] = difference
			count =+ 1
			print i, j, k, difference

f = file("experiment_results.txt","wb")
np.save(f,results)
f.close()
#best for 2: .00095, .0002

#best for 3: 0.0015 0.0005 0.00022 564.84883375

Beispiel #3
0
	def evaluate(self,bed_elev, surf_elev):
		return self.constrain_to_reality(bed_elev) + self.constrain_to_nolan(bed_elev) + self.smoothness_constraint(bed_elev) + tools.calculate_surface_difference(surf_elev, self.gps_surface)