from varglas.meshing import MeshGenerator, MeshRefiner
from varglas.io import DataInput
from varglas.data.data_factory import DataFactory
from pylab import *
from scipy.interpolate import interp2d
#===============================================================================
# data preparation :
out_dir = 'dump/meshes/'
# collect the raw data :
bamber = DataFactory.get_bamber()
rignot = DataFactory.get_gre_rignot()
searise = DataFactory.get_searise()
gimp = DataFactory.get_gre_gimp()
# create data objects to use with varglas :
dbm = DataInput(bamber, gen_space=False)
drg = DataInput(rignot, gen_space=False)
dsr = DataInput(searise, gen_space=False)
dgm = DataInput(gimp, gen_space=False)
drg.set_data_max('U_ob', boundary=1000.0, val=1000.0)
drg.set_data_min('U_ob', boundary=0.0, val=0.0)
#===============================================================================
# form field from which to refine :
# antarctica :: Info : 4449632 vertices 25902918 elements
drg.data['ref'] = (0.02 + 1 / (1 + drg.data['U_ob'])) * 50000
## plot to check :
