Update addprocess.py

Geodels · Mar 26, 2024 · 16fb569 · 16fb569
1 parent ced68db
commit 16fb569
Showing 1 changed file with 4 additions and 2 deletions.
diff --git a/gospl/tools/addprocess.py b/gospl/tools/addprocess.py
@@ -37,6 +37,9 @@ def __init__(self):
 
         if self.flexOn:
             self.localFlex = np.zeros(self.lpoints)
+            self.boundflex = self.boundCond.replace('0','2')
+            self.boundflex = self.boundflex.replace('1','0')
+            self.boundflex = self.boundflex.replace('2','1')
         if self.oroOn:
             self.oroEPS = np.finfo(float).eps
 
@@ -124,10 +127,9 @@ def applyFlexure(self):
                 regNewZ[self.regOnIDs] = newZ[self.regIDs[self.regOnIDs,0]]
             regNewZ = regNewZ.reshape(self.reg_ny,self.reg_nx) 
 
-            # Force boundary to be all reflective
             newh = flexure(regNewZ,regOldZ,self.reg_ny,self.reg_nx,
                            self.reg_yl,self.reg_xl,self.flex_rhos,
-                           self.flex_rhoa,self.flex_eet,0) #int(self.boundCond))
+                           self.flex_rhoa,self.flex_eet,int(self.boundflex))
             rflexTec = (newh - regNewZ).ravel()
 
             # Interpolate back to goSPL mesh