import varglas.model as model from varglas.data.data_factory import DataFactory from varglas.mesh.mesh_factory import MeshFactory from varglas.io import DataInput, DataOutput from varglas.helper import default_config from fenics import * from scipy.io import loadmat # get the input args : out_dir = 'dump/bed/09/bv_smb/' in_dir = 'dump/bed/09/' mesh = Mesh(in_dir + 'submesh.xdmf') Q = FunctionSpace(mesh, 'CG', 1) bamber = DataFactory.get_bamber() dbm = DataInput(bamber, mesh=mesh) S = Function(Q) B = Function(Q) File(in_dir + 'S_s.xml') >> S File(in_dir + 'B_s.xml') >> B adot_v = loadmat(in_dir + 'smb/smb.mat') dbm.data['smb'] = adot_v['map_data'] adot = dbm.get_expression('smb', near=True) config = default_config() config['output_path'] = out_dir
import varglas.physical_constants as pc import varglas.model as model from varglas.mesh.mesh_factory import MeshFactory from varglas.data.data_factory import DataFactory from varglas.io import DataInput from fenics import * out_dir = 'dump/vars_ant_spacing/' thklim = 1.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() rignot = DataFactory.get_gre_rignot() # define the mesh : mesh = Mesh('dump/meshes/gre_mesh_ant_spacing.xml') # create data objects to use with varglas : dsr = DataInput(searise, mesh=mesh) dbm = DataInput(bamber, mesh=mesh) dfm = DataInput(fm_qgeo, mesh=mesh) drg = DataInput(rignot, mesh=mesh) # change the projection of all data to Rignot projection : dsr.change_projection(drg) dbm.change_projection(drg) dfm.change_projection(drg) # get the expressions used by varglas :