dir_b = sys.argv[1] + '/0' # directory to save # set the output directory : out_dir = dir_b + str(i) + '/' set_log_active(True) #set_log_level(PROGRESS) thklim = 200.0 # collect the raw data : searise = DataFactory.get_searise(thklim = thklim) bamber = DataFactory.get_bamber(thklim = thklim) fm_qgeo = DataFactory.get_gre_qgeo_fox_maule() #measure = DataFactory.get_gre_measures() rignot = DataFactory.get_gre_rignot() # define the mesh : mesh = MeshFactory.get_greenland_coarse() # 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) dfm = DataInput(None, fm_qgeo, mesh=mesh) drg = DataInput(None, rignot, mesh=mesh) # change the projection of the measures data to fit with other data : #dms.change_projection(dsr) drg.change_projection(dsr)