fix(mf6bmiUtil): use base model in get_grid_type_model (#2845)

#2705 switched PRT to ExplicitModelType, but did not accommodate this in the BMI implementation. Up to now get_grid_type_model only worked with NumericalModelType. Use BaseModelType instead.

Author: wpbonelli

doc/ReleaseNotes/develop.toml

@@ -82,4 +82,9 @@ description = "Previously, GWF model output files could only be referenced by a
 [[items]]
 section = "fixes"
 subsection = "model"
-description = "PRT was described as configuring a single release at the beginning of the simulation by default, but a release was in fact scheduled for the first time step of every stress period. Fix the default release behavior so a single release occurs at the beginning of the simulation."
+description = "PRT was described as configuring a single release at the beginning of the simulation by default, but a release was in fact scheduled for the first time step of every stress period. Fix the default release behavior so a single release occurs at the beginning of the simulation."
+
+[[items]]
+section = "fixes"
+subsection = "basic"
+description = "The BMI interface's get\\_grid\\_model\\_type function worked only for numerical model types such as GWF, GWT, and GWE, but would crash for PRT models. Fix the function's implementation to work with all model types."

srcbmi/mf6bmiUtil.f90

@@ -238,18 +238,19 @@ end function getSolution
   subroutine get_grid_type_model(model_name, grid_type_f)
     ! -- modules
     use ListsModule, only: basemodellist
-    use NumericalModelModule, only: NumericalModelType, GetNumericalModelFromList
+    use BaseModelModule, only: BaseModelType, GetBaseModelFromList
     ! -- dummy variables
     character(len=LENMODELNAME) :: model_name
     character(len=LENGRIDTYPE) :: grid_type_f
     ! -- local variables
     integer(I4B) :: i
-    class(NumericalModelType), pointer :: numericalModel
+    class(BaseModelType), pointer :: baseModel
 
     do i = 1, basemodellist%Count()
-      numericalModel => GetNumericalModelFromList(basemodellist, i)
-      if (numericalModel%name == model_name) then
-        call numericalModel%dis%get_dis_type(grid_type_f)
+      baseModel => GetBaseModelFromList(basemodellist, i)
+      if (baseModel%name == model_name) then
+        call baseModel%dis%get_dis_type(grid_type_f)
+        return
       end if
     end do
   end subroutine get_grid_type_model