Merge pull request #15 from Yuksel-Rudy/main

Adding Resetting Array Origin to Mid-Platform:

Author: lsirkis

famodel/mooring/mooring.py

@@ -198,7 +198,7 @@ def reposition(self, r_center=None, heading=None, project=None,
         r_center : list or nested list
             The x, y, z coordinates of the platform(s) (undisplaced) [m]. If shared mooring, must be a list of lists, with the
             coordinates for each platform connection. In this case, end A platform connection is the first entry and end B 
-            platform connection is the second entry.
+            platform connection is the second entry. If not given, r_center will be populated from what self is attached to.
         heading : float
             The absolute heading compass direction of the mooring line from end B
             [deg or rad] depending on degrees parameter (True or False). Must account for the platform heading as well.
@@ -227,11 +227,15 @@ def reposition(self, r_center=None, heading=None, project=None,
         # heading 2D unit vector
         u = np.array([np.cos(phi), np.sin(phi)])
 
-        if self.shared == 1:
-            r_centerB = np.array(r_center)[1]
-            r_centerA = np.array(r_center)[0]
+        if np.any(r_center):
+            if self.shared == 1:
+                r_centerA = np.array(r_center)[0]
+                r_centerB = np.array(r_center)[1]
+            else:
+                r_centerB = np.array(r_center)
         else:
-            r_centerB = np.array(r_center)
+            r_centerA = self.attached_to[0].r
+            r_centerB = self.attached_to[1].r
 
         # create fairlead radius list for end A and end B if needed
         if not rad_fair:

famodel/platform/platform.py

@@ -64,8 +64,10 @@ def __init__(self, id, r=[0,0,0], heading=0, mooring_headings=[60,180,300],rFair
     def setPosition(self, r, heading=None, degrees=False,project=None):
         '''
         Set the position/orientation of the platform as well as the associated
-        anchor points.
+        anchor points. 
         
+        "Note: must only be used for a platform that's only attached with anchored lines"
+
         Parameters
         ----------
         r : list

famodel/project.py

@@ -117,7 +117,8 @@ def __init__(self, lon=0, lat=0, file=None, depth=202,raft=1):
         # MoorPy system associated with the project
         self.ms  = None
 
-        
+        # RAFTDict associated with the project
+        self.RAFTDict = None
 
         # ----- if an input file has been passed, load it -----
         if file:
@@ -788,15 +789,17 @@ def loadDesign(self, d, raft=1):
             # create a name for the raft model
             RAFTDict['name'] = 'Project_Array'
             RAFTDict['type'] = 'input file for RAFT'
-    
+
+            self.RAFTDict = deepcopy(RAFTDict)
+
             # create RAFT model if necessary components exist
             if 'platforms' in RAFTDict or 'platform' in RAFTDict:
 
                     self.getRAFT(RAFTDict,pristine=1)
-                    
+                
+            
 
 
-        
 
 
     # ----- Site conditions processing functions -----
@@ -4114,20 +4117,21 @@ def extractFarmInfo(self, cmax=5, fmax=10/6, Cmeander=1.9):
         i = 0
         yaw_init = np.zeros((1, len(self.platformList.items())))
         for _, pf in self.platformList.items():
-            x, y, z   = pf.body.r6[0], pf.body.r6[1], pf.body.r6[2]
-            phi_deg       = np.degrees(pf.phi)  # float((90 - np.degrees(pf.phi)) % 360)  # Converting FAD's rotational convention (0deg N, +ve CW) into FF's rotational convention (0deg E, +ve CCW)
-            phi_deg       = (phi_deg + 180) % 360 - 180  # Shift range to -180 to 180
-            for att in pf.attachments.values():
-                if isinstance(att['obj'],Turbine):
-                    D    = 240   # att['obj'].D         (assuming 15MW)
-                    zhub = att['obj'].dd['hHub']
+            if pf.entity=='FOWT':
+                x, y, z   = pf.body.r6[0], pf.body.r6[1], pf.body.r6[2]
+                phi_deg       = np.degrees(pf.phi)  # float((90 - np.degrees(pf.phi)) % 360)  # Converting FAD's rotational convention (0deg N, +ve CW) into FF's rotational convention (0deg E, +ve CCW)
+                phi_deg       = (phi_deg + 180) % 360 - 180  # Shift range to -180 to 180
+                for att in pf.attachments.values():
+                    if isinstance(att['obj'],Turbine):
+                        D    = 240   # att['obj'].D         (assuming 15MW)
+                        zhub = att['obj'].dd['hHub']
 
-            wts[i] = {
-                'x': x, 'y': y, 'z': z, 'phi_deg': phi_deg, 'D': D, 'zhub': zhub, 
-                'cmax': cmax, 'fmax': fmax, 'Cmeander': Cmeander
-                }
-            yaw_init[0, i] = -phi_deg
-            i += 1
+                wts[i] = {
+                    'x': x, 'y': y, 'z': z, 'phi_deg': phi_deg, 'D': D, 'zhub': zhub, 
+                    'cmax': cmax, 'fmax': fmax, 'Cmeander': Cmeander
+                    }
+                yaw_init[0, i] = -phi_deg
+                i += 1
 
         # store farm-level wind turbine information
         self.wts = wts
@@ -4206,6 +4210,48 @@ def FFarmCompatibleMDOutput(self, filename, unrotateTurbines=True, renameBody=Tr
             with open(filename, 'w') as f:
                 f.writelines(newLines)
 
+    def resetArrayCenter(self, FOWTOnly=True):
+        '''
+        Function to reset array center such that the farm origin is the mid-point between all FOWT platforms:
+
+        Parameters
+        ----------
+        FOWTOnly : bool
+            find the center between only FOWT-entity platforms if True.
+        '''       
+        xCenter = np.mean([p.r[0] for p in self.platformList.values() if p.entity=='FOWT' or not FOWTOnly])
+        yCenter = np.mean([p.r[1] for p in self.platformList.values() if p.entity=='FOWT' or not FOWTOnly])          
+        delta   = np.array([xCenter, yCenter])
+        
+        # Change boundaries
+        self.boundary -= delta
+
+        # Change bathymetry
+        self.grid_x -= delta[0]
+        self.grid_y -= delta[1]
+        
+        # Change all platform locations and associated anchors/moorings/cables
+        for pf in self.platformList.values():
+            pf.r[:2]-=delta
+            for i, att in enumerate(pf.attachments.values()):
+                obj = att['obj']
+                if isinstance(obj, Mooring):           
+                    obj.reposition(project=self)
+
+                if isinstance(obj, Cable):
+                    obj.reposiiton(project=self)
+
+        # Change RAFTDict if available.
+        if self.RAFTDict:
+            x_idx = self.RAFTDict['array']['keys'].index('x_location')
+            y_idx = self.RAFTDict['array']['keys'].index('y_location')            
+            for i in range(len(self.platformList.values())):
+                self.RAFTDict['array']['data'][i][x_idx] -= delta[0]
+                self.RAFTDict['array']['data'][i][y_idx] -= delta[1]
+
+            if 'platforms' in self.RAFTDict or 'platform' in self.RAFTDict:
+                    self.getRAFT(self.RAFTDict,pristine=1)
+        
 
     def updateFailureProbability(self):
         '''
@@ -4219,7 +4265,7 @@ def updateFailureProbability(self):
         None.
 
         '''
-                
+        pass