from dolfin import *
from pylab import sqrt, copy
from time import time
# make the directory if needed :
out_dir = './test_results/'
set_log_active(True)
#set_log_level(PROGRESS)
thklim = 200.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()
# define the meshes :
mesh = Mesh('meshes/mesh_high_new.xml')
flat_mesh = Mesh('meshes/mesh_high_new.xml')
#mesh = Mesh('meshes/mesh_low.xml')
#flat_mesh = Mesh('meshes/mesh_low.xml')
mesh.coordinates()[:, 2] /= 100000.0
flat_mesh.coordinates()[:, 2] /= 100000.0
# create data objects to use with varglas :
dsr = DataInput(None, searise, mesh=mesh)
dbm = DataInput(None, bamber, mesh=mesh)
sys.path.append(src_directory)
from src.utilities import DataInput,DataOutput
from data.data_factory import DataFactory
from meshes.mesh_factory import MeshFactory
from src.physics import VelocityBalance_2
from dolfin import Mesh, set_log_active
set_log_active(True)
thklim = 20.0
# collect the raw data :
searise = DataFactory.get_searise()
measure = DataFactory.get_gre_measures()
meas_shf = DataFactory.get_shift_gre_measures()
bamber = DataFactory.get_bamber(thklim = thklim)
# load a mesh :
mesh = Mesh("results/meshes/refined_mesh.xml")
# 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)
dms.change_projection(dsr)
H = dbm.get_projection("H_n")
S = dbm.get_projection("h_n")
