# make the directory if needed :
out_dir = './results/'
set_log_active(True)
#set_log_level(PROGRESS)
thklim = 10.0
# collect the raw data :
searise = DataFactory.get_searise(thklim = thklim)
measure = DataFactory.get_gre_measures()
#meas_shf = DataFactory.get_shift_gre_measures()
bamber = DataFactory.get_bamber(thklim = thklim)
fm_qgeo = DataFactory.get_gre_qgeo_fox_maule()
#sec_qgeo = DataFactory.get_gre_qgeo_secret()
merged = DataFactory.get_gre_merged()
# define the meshes :
mesh = Mesh('../../../meshes/greenland/greenland_coarse_mesh.xml')
flat_mesh = Mesh('../../../meshes/greenland/greenland_coarse_mesh.xml')
mesh.coordinates()[:,2] /= 1000.0
flat_mesh.coordinates()[:,2] /= 1000.0
# create data objects to use with varglas :
dsr = DataInput(None, searise, mesh=mesh)
dbm = DataInput(None, bamber, mesh=mesh)
#dms = DataInput(None, measure, mesh=mesh)
#dmss = DataInput(None, meas_shf, mesh=mesh)
dfm = DataInput(None, fm_qgeo, mesh=mesh)
#dsq = DataInput(None, sec_qgeo, mesh=mesh)