feat(idm): integrate PRP package (#2444)

Author: mjreno

doc/mf6io/framework/processing_of_input.tex

@@ -66,7 +66,8 @@ \subsection{Supported Components}
 PRT/NAM & PRT name file \\
 PRT/DIS & DIS6 \\ 
 PRT/DISV & DISV6 \\
-PRT/MIP & MIP \\
+PRT/MIP & MIP6 \\
+PRT/PRP & PRP6 \\
 \hline
 EXG/GWFGWF & GWF6-GWF6 \\
 EXG/GWFGWT & GWF6-GWT6 \\

doc/mf6io/mf6ivar/dfn/prt-prp.dfn

@@ -15,6 +15,7 @@ name print_input
 type keyword
 reader urword
 optional true
+mf6internal iprpak
 longname print input to listing file
 description REPLACE print_input {'{#1}': 'all model stress package'}
 
@@ -23,6 +24,7 @@ name dev_exit_solve_method
 type integer
 reader urword
 optional true
+mf6internal iexmeth
 longname exit solve method
 description the method for iterative solution of particle exit location and time in the generalized Pollock's method.  0 default, 1 Brent, 2 Chandrupatla.  The default is Brent's method.
 
@@ -31,6 +33,7 @@ name exit_solve_tolerance
 type double precision
 reader urword
 optional true
+mf6internal extol
 longname exit solve tolerance
 description the convergence tolerance for iterative solution of particle exit location and time in the generalized Pollock's method.  A value of 0.00001 works well for many problems, but the value that strikes the best balance between accuracy and runtime is problem-dependent.
 default_value 1e-5
@@ -101,6 +104,7 @@ type record trackcsv fileout trackcsvfile
 shape
 reader urword
 tagged true
+mf6internal trackcsvfr
 optional true
 longname
 description
@@ -149,6 +153,7 @@ name stop_at_weak_sink
 type keyword
 reader urword
 optional true
+mf6internal istopweaksink
 longname stop at weak sink
 description is a text keyword to indicate that a particle is to terminate when it enters a cell that is a weak sink.  By default, particles are allowed to pass though cells that are weak sinks.
 
@@ -175,6 +180,7 @@ shape
 reader urword
 tagged true
 optional true
+mf6internal releasetr
 longname
 description
 removed 6.6.0
@@ -209,6 +215,7 @@ shape
 reader urword
 tagged true
 optional true
+mf6internal release_timesfr
 longname
 description
 removed 6.6.0
@@ -220,6 +227,7 @@ reader urword
 in_record true
 tagged true
 shape
+mf6internal release_timesfn
 longname
 description keyword indicating release times file name will follow
 removed 6.6.0
@@ -243,6 +251,7 @@ type string
 valid drop stop stay
 reader urword
 optional true
+mf6internal idrymeth
 longname what to do in dry-but-active cells
 description is a string indicating how particles should behave in dry-but-active cells (as can occur with the Newton formulation).  The value can be ``DROP'', ``STOP'', or ``STAY''.  The default is ``DROP'', which passes particles vertically and instantaneously to the water table. ``STOP'' causes particles to terminate. ``STAY'' causes particles to remain stationary but active.
 
@@ -260,6 +269,7 @@ name release_time_tolerance
 type double precision
 reader urword
 optional true
+mf6internal rttol
 longname release time coincidence tolerance
 description real number indicating the tolerance within which to consider consecutive release times coincident. Coincident release times will be merged into a single release time. The default is $\epsilon \times 10^{11}$, where $\epsilon$ is machine precision.
 
@@ -268,6 +278,7 @@ name release_time_frequency
 type double precision
 reader urword
 optional true
+mf6internal rtfreq
 longname release time frequency
 description real number indicating the time frequency at which to release particles. This option can be used to schedule releases at a regular interval for the duration of the simulation, starting at the simulation start time. The release schedule is the union of this option, the RELEASETIMES block, and PERIOD block RELEASESETTING selections. If none of these are provided, a single release time is configured at the beginning of the first time step of the simulation's first stress period.
 
@@ -277,6 +288,7 @@ type string
 valid none eager
 reader urword
 optional true
+mf6internal ichkmeth
 longname coordinate checking method
 description approach for verifying that release point coordinates are in the cell with the specified ID. By default, release point coordinates are checked at release time.
 default eager
@@ -378,6 +390,7 @@ name releasetimes
 type recarray time
 shape (nreleasetimes)
 reader urword
+optional true
 longname
 description
 

make/makefile

@@ -151,6 +151,7 @@ $(OBJDIR)/swf-chdidm.o \
 $(OBJDIR)/swf-cdbidm.o \
 $(OBJDIR)/sim-tdisidm.o \
 $(OBJDIR)/sim-namidm.o \
+$(OBJDIR)/prt-prpidm.o \
 $(OBJDIR)/prt-namidm.o \
 $(OBJDIR)/prt-mipidm.o \
 $(OBJDIR)/prt-disvidm.o \
@@ -519,6 +520,7 @@ $(OBJDIR)/NumericalSolution.o \
 $(OBJDIR)/MappedMemory.o \
 $(OBJDIR)/NCModel.o \
 $(OBJDIR)/Mf6FileStoInput.o \
+$(OBJDIR)/Mf6FileSetting.o \
 $(OBJDIR)/Mf6FileList.o \
 $(OBJDIR)/Mf6FileLayerArray.o \
 $(OBJDIR)/Mf6FileGridArray.o \

msvs/mf6core.vfproj

@@ -227,6 +227,7 @@
 		<File RelativePath="..\src\Idm\prt-disvidm.f90"/>
 		<File RelativePath="..\src\Idm\prt-mipidm.f90"/>
 		<File RelativePath="..\src\Idm\prt-namidm.f90"/>
+		<File RelativePath="..\src\Idm\prt-prpidm.f90"/>
 		<File RelativePath="..\src\Idm\sim-namidm.f90"/>
 		<File RelativePath="..\src\Idm\sim-tdisidm.f90"/>
 		<File RelativePath="..\src\Idm\swf-cdbidm.f90"/>
@@ -546,6 +547,7 @@
 		<File RelativePath="..\src\Utilities\Idm\mf6blockfile\Mf6FileGridArray.f90"/>
 		<File RelativePath="..\src\Utilities\Idm\mf6blockfile\Mf6FileLayerArray.f90"/>
 		<File RelativePath="..\src\Utilities\Idm\mf6blockfile\Mf6FileList.f90"/>
+		<File RelativePath="..\src\Utilities\Idm\mf6blockfile\Mf6FileSetting.f90"/>
 		<File RelativePath="..\src\Utilities\Idm\mf6blockfile\Mf6FileStoInput.f90"/>
 		<File RelativePath="..\src\Utilities\Idm\mf6blockfile\StructArray.f90"/>
 		<File RelativePath="..\src\Utilities\Idm\mf6blockfile\StructVector.f90"/></Filter>