thacker : add export_value_at_point

benchmarks/wetting-drying/thacker/thacker.py

@@ -62,11 +62,14 @@ slimPre.write_file('g.nc', region=None, time=None, data=[('g', g)])
 ###########
 ### run ###
 ###########
+x = numpy.arange(0,R,5000)
+y = 0*x
 
 domain = slim.Domain(meshfile,('bath.nc','bath'), g=("g.nc","g"))
 eq = slim.ShallowWater2dWD(domain, "implicit", final_time=tend, hlimt=hlimt, hlimd=hlimd, linDrag=linDrag, manning=manning)
 eq.set_initial_condition(('init.nc','eta'),('init.nc','u'),('init.nc','v'))   
 eq.set_boundary_coast('Wall')
+eq.export_value_at_point("sol.txt", x, y)
 loop = slim.Loop(maximum_time_step=dt, export_time=600)
 loop.add_equation(eq)
 loop.run()