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()