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 :
