How to impose a minimum layer thickness in ModelMuse?
Mohammad Hadi Sahil
Robert Lanning
IMPORT MODEL LAYERS AND ADJUST FOR VERTICAL SEPARATION
Online reference:
hhttps://youtu.be/2GOv4o--1Y8?t=2271
This procedure imports a SHP but doesn't assign it to a layer so it is “user defined”. It is then used to define the base of an aquifer as a “User Defined” shape. This gives more options in adjusting top and bottom formulas of a narrow layer if it overlaps another layer. This example adjusts bottom of a narrow second layer in a model (Layer2_bottom) because portions of its profile stick above the bottom of Layer 1 above it.
- Begin by creating an elevation points SHP file for Layer 2_bottom by importing CSV file into
QGIS (see Attachment E). This created file is named Layer2-A and is the elev point file for
Layer 2’s Base.
* In ModelMuse Import, Select Layer2-A.shp
- Set Values of cells by interpolation
- single multipart object
- visible but not selected
- Number of Z formulas = zero
- under DATA tab select Layer2_bottom / New data set / none for interpolation. OK
- all points show up as circles on grid (Object/show-hide and deselect)
* From DATA->EDIT DATA SETS->User Defined->Created from shapefile-> Layer2-A
- Set TYPE = Real
- ORIENTATION = 2D TOP
- INTERPOLATION = Natural Neighbors
(note: Nearest Point result very blocky profile & Fitted Surface exaggerates profile)
- APPLY
* From Data Sets->Required->Layer Definition
- Select desired Layer2_bottom
- Edit Formula. From right panel select Layer2-A, dbl click to insert into formula
- OK/APPLY – creates a smooth profile but with portions above top of Layer1_bottom
* DATA->EDIT DATA SETS->Required->Layer Definition-> Layer2-bottom
EDIT FORMULA use Functions->Logical to build "If ___ then ___"
*FORMAT: IF(state condition, if true do this, if false, do that)
-Use box on right to build formula by double-clicking selected components
- For this situation built this “CONDITION,IF,THEN” formula:
IF(Layer1_Bottom-3<Layer2-A,Layer1_Bottom-3,Layer2-A)
This means: Condition = Layer1 bottom elev minus 3 ft is less than Layer2 elev
If TRUE, make Layer2 base = Layer1 base minus 3 ft
IF FALSE, make Layer2 base = Layer2-A elevation
This sets minimum thickness of Layer 2 as 3 ft (adjust for your model as needed)
- OK. APPLY. This should fix all profile crossing issues
Robel
--
This group was created in 2004 by Mr. C. P. Kumar, Former Scientist 'G', National Institute of Hydrology, Roorkee. Please visit his webpage at https://www.angelfire.com/nh/cpkumar/
---
You received this message because you are subscribed to the Google Groups "MODFLOW Users Group" group.
To unsubscribe from this group and stop receiving emails from it, send an email to modflow+u...@googlegroups.com.
To view this discussion on the web visit https://groups.google.com/d/msgid/modflow/3cedb8a3-88d5-4ce0-89a8-4d5dcc1a601bn%40googlegroups.com.
Mohammad Hadi Sahil
Robert Lanning
HOW TO CREATE CONTOUR MAP IN QGIS FROM CSV ELEVATION FILE
Online reference:
https://www.youtube.com/watch?v=btKaqZ9YVR0
1. Need CSV file with point_ID, easting, northing, and elev
2. Open QGIS and install CONTOUR plugin if not installed
-Plugins (top menu row) ->Manage & Install Plugins
-Type CONTOUR in search row, when shows up, install
3. Go to LAYER tab at top -> Data Source Manager
- Select DELIMITED TEXT
- Browse in top box to find CSV file and select file
- confirm Geometry definitions, including CRS (whatever using for your model)
- Add and Close window
4. In left panel, right click on new file name
- Select EXPORT -> Save Features As (this saves CSV points as *.SHP)
- Use "3 dots" on right to select location, name file, and SAVE
(Hint: Create SHP subfolder to keep all Shapefiles in one place)
4A. If desire the point symbols and names to show up on contour map
- Select either original CSV file or SHP file created in #4
- Right click -> Properties. Select Symbol and adjust size and color of points
- Right click -> Properties. Select Select LABELS. From pull down select Single Label
- Adjust text, size, color, placement
5. Right click, Activity Table - Check correct data in columns
- if you put extra REMARKS column with text defining various data in rows, can use to filter if needed
6. Under VECTOR menu tab, select Contour
- verify Point Layer, Data Value (elev column you want)
- select Method contours (Fixed Elevation)
- Select contour interval
- select Max number contours (want 50 or more)
- modify number of decimels in lower section
- Add and Close
7. Smooth contours
- Select Processing toolbox "Gear" if not open on right
- Type in Smooth in search, select Vector Geometry Smooth, double click
- Change iterations to 10, leave other parameters
- Run and Close
- check for smoothing results
- If OK, delete the former layer (one prior to smoothing)
- Rename Smoothed file
- Select icon or right click and select Make Permanent, select location & name. Save
8. Adjust contours
- Highlight file in left window
- right click, select PROPERTIES
- In left menu select SYMBOLOGY
- can change color, width of contour lines
9. Add labels to contours
- right click on file, select PROPERTIES
- Select LABELS, in box pull down at top select SINGLE LABELS
- Adjust components (Text:type, size, color. Placement mode, repeating)
- Save
10. Create Map
- Go to Project-> Layout Manager
- Select EMPTY LAYOUT (under New from Template) - Select CREATE
- Name Print Layout
- Right click in open white space, select Page Properties, select paper size
- From Menu select ADD ITEM->Add Map. Click & drag cursor to fill space.
this inserts the contour map
- In right column select ITEM PROPERTIES and set desired scale (1:xxxx)
by "evening out" the number in the box; e.g. 509952807->500,000,000
- Select other options and adjust (FRAME, GRID - set interval for X,Y)
Can adjust grid lines color & form
- Select DRAW COORDINATES to label grid lines
Adjust fields for orientation, distance from frame, etc.
- Can add Scale Bar, North Arrow, etc.
ADD ITEM (top menu)-> Scale Bar. Use cursor to place
Same for adding North Arrow, Legend, etc.
11. Save QGIS Project
12. To create PDF
- Project -> Import/Export -> Select Export to PDF
- Select location and name file
Note: Can also save basic contour map in Step #9 using same process
Note: Original draft, SHP files can be deleted if desired
Mohammad Hadi Sahil
Mohammad Hadi Sahil
--
This group was created in 2004 by Mr. C. P. Kumar, Former Scientist 'G', National Institute of Hydrology, Roorkee. Please visit his webpage at https://www.angelfire.com/nh/cpkumar/
---
You received this message because you are subscribed to the Google Groups "MODFLOW Users Group" group.
To unsubscribe from this group and stop receiving emails from it, send an email to modflow+u...@googlegroups.com.
To view this discussion on the web visit https://groups.google.com/d/msgid/modflow/fb29736f-2bff-4609-b06b-bac1769e9291n%40googlegroups.com.