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