from src.utilities import DataInput, DataOutput
from dolfin import *
from pylab import sqrt, copy
from time import time
# 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)
import sys
src_directory = '../../../'
sys.path.append(src_directory)
from src.utilities import DataInput,DataOutput
from data.data_factory import DataFactory
from src.physics import VelocityBalance_2
from dolfin import Mesh, set_log_active
from plot.plothelp.plotting import plotIce
import os
set_log_active(True)
# collect the raw data :
searise = DataFactory.get_searise()
measure = DataFactory.get_gre_measures()
#bamber = DataFactory.get_bamber()
direc = os.path.dirname(os.path.realpath(__file__))
# load a mesh :
mesh = Mesh("../meshes/mesh.xml")
# create data objects to use with varglas :
dsr = DataInput(None, searise, mesh=mesh, create_proj=True)
#dbm = DataInput(None, bamber, mesh=mesh)
dms = DataInput(None, measure, mesh=mesh, create_proj=True, flip=True)
dms.change_projection(dsr)
dsr.set_data_min('H', 10.0, 10.0)
