update(rhs, hcof, AdvancedInput): bug fixes for setting variable values for advanced inputs

* update rhs and hcof to copy values to pointer instead of overwriting the pointer
* add a check for AdvancedInput variables that do not have pointer support in xmipy

Author: jlarsen-usgs

modflowapi/extensions/data.py

@@ -558,10 +558,10 @@ def get_variable(self, name, model=None, package=None):
 
         try:
             values = self.mf6.get_value_ptr(var_addr)
+            self._ptrs[name.lower()] = values
         except xmipy.errors.InputError:
             values = self.mf6.get_value(var_addr)
 
-        self._ptrs[name.lower()] = values
         return values.copy()
 
     def set_variable(self, name, values, model=None, package=None):
@@ -601,7 +601,18 @@ def set_variable(self, name, values, model=None, package=None):
                 f"current shape={values.shape}, valid shape={values0.shape}"
             )
 
-        self._ptrs[name.lower()] = values
+        if name.lower() not in self._ptrs:
+            # this is a set value situation
+            self.mf6.set_value(
+                self.mf6.get_var_address(
+                    name.upper(),
+                    self.parent.model.name,
+                    self.parent.pkg_name
+                ),
+                values
+            )
+        else:
+            self._ptrs[name.lower()] = values
 
 
 class ScalarInput:

modflowapi/extensions/pakbase.py

@@ -341,7 +341,7 @@ def rhs(self, values):
         if self._rhs is None:
             return
 
-        self._rhs = values
+        self._rhs[:] = values[:]
 
     @property
     def hcof(self):
@@ -362,7 +362,7 @@ def hcof(self, values):
         if self._hcof is None:
             return
 
-        self.__hcof = values
+        self._hcof[:] = values[:]
 
 
 class ListPackage(PackageBase):