def test_get_modellayers_mf2005():
import flopy
modelname = 'test_mf2005'
ml = flopy.modflow.Modflow(modelname, exe_name='mf2005')
# Model domain and grid definition
Lx = 300000.
Ly = 400000.
ztop = 50.
zbot = -150.
nlay = 4
nrow = 40
ncol = 30
delr = Lx / ncol
delc = Ly / nrow
delv = (ztop - zbot) / nlay
botm = np.linspace(ztop, zbot, nlay + 1)
# Create the discretization object
dis = flopy.modflow.ModflowDis(ml,
nlay,
nrow,
ncol,
delr=delr,
delc=delc,
top=ztop,
botm=botm[1:],
xul=0,
yul=700000)
dino_gw = ttf.test_obscollection_dinozip_gw()
modellayers = dino_gw.gwobs.get_modellayers(ml)
return modellayers
def test_get_nearest_point():
# check two of the same observation collections
# every point must find itself as the nearest point
dino_gw = ttf.test_obscollection_dinozip_gw()
fl = ttf.test_obscollection_fieldlogger()
dino_gw[['nearest point',
'distance nearest point']] = dino_gw.geo.get_nearest_point(fl)
return dino_gw
def test_obscollection_dino_to_mapgraph():
try:
from art_tools import hpd_extension
gw = ttf.test_obscollection_dinozip_gw()
gw.art.plot_mapgraphs(plot_ylim='min_dy')
return
except ModuleNotFoundError as e:
print(e)
return
def test_obscollection_dino_to_imap():
dino_gw = ttf.test_obscollection_dinozip_gw()
dino_gw.geo.set_lat_lon(verbose=True)
dino_gw.plots.interactive_map(plot_dir,
plot_columns=['stand_m_tov_nap'],
fname='imap.html',
legend_name='grondwater DINO',
add_filter_to_legend=True,
hoover_names=['gws'],
zoom_start=9,
verbose=True)
return
def test_get_nearest_polygon():
import geopandas as gpd
from shapely.geometry import Polygon
# check two of the same observation collections
# every point must find itself as the nearest point
dino_gw = ttf.test_obscollection_dinozip_gw()
extent = dino_gw.geo.get_extent()
polygon1 = Polygon(((extent[0], extent[2]), (extent[0], extent[3]),
(extent[1], extent[3]), (extent[1], extent[2]),
(extent[0], extent[2])))
polygon2 = Polygon(((-extent[0], -extent[2]), (-extent[0], -extent[3]),
(-extent[1], -extent[3]), (-extent[1], -extent[2]),
(-extent[0], -extent[2])))
gdf = gpd.GeoDataFrame({'name': [1, 2], 'geometry': [polygon1, polygon2]})
dino_gw[['nearest polygon', 'distance nearest polygon'
]] = dino_gw.geo.get_nearest_polygon(gdf)
assert (dino_gw['nearest polygon'] == 0.0).all()
assert (dino_gw['distance nearest polygon'] == 0.0).all()
return dino_gw[['nearest polygon', 'distance nearest polygon']]
def test_get_modellayers_mf6_structured():
import flopy
# Create the Flopy simulation object
model_name = 'test_mf6_structured'
sim = flopy.mf6.MFSimulation(sim_name=model_name,
exe_name='mf6',
version='mf6')
# Create the Flopy groundwater flow (gwf) model object
model_nam_file = '{}.nam'.format(model_name)
gwf = flopy.mf6.ModflowGwf(sim,
modelname=model_name,
model_nam_file=model_nam_file)
Lx = 300000.
Ly = 400000.
ztop = 50.
zbot = -150.
nlay = 4
nrow = 40
ncol = 30
delr = Lx / ncol
delc = Ly / nrow
delv = (ztop - zbot) / nlay
botm = np.linspace(ztop, zbot, nlay + 1)
dis = flopy.mf6.ModflowGwfdis(gwf,
xorigin=0,
yorigin=300000,
nlay=4,
nrow=nrow,
ncol=ncol,
delr=Lx / ncol,
delc=Ly / nrow,
top=ztop,
botm=botm[1:])
dino_gw = ttf.test_obscollection_dinozip_gw()
modellayers = dino_gw.gwobs.get_modellayers(gwf)
return modellayers
def test_within_extent():
dino_gw = ttf.test_obscollection_dinozip_gw()
extent = [210350, 213300, 473300, 474000]
dino_gw.geo.within_extent(extent, inplace=True)
assert dino_gw.shape[0] == 4
return dino_gw
def test_set_filter_num():
dino_gw = ttf.test_obscollection_dinozip_gw()
dino_gw.gwobs.set_filter_num(if_exists='replace')
return dino_gw
def test_get_regis_layer():
dino_gw = ttf.test_obscollection_dinozip_gw()
return dino_gw.gwobs.get_regis_layers()
