Adding Resetting Array Origin to Mid-Platform:

This will help with TurbSim and FFarm integration as farm needs to have an origin at the center of the farm.
The resetArrayCenter can accept centering the array based on all platforms or only those carrying FOWTs.
The function will change bathymetries, boundaries, platform locations, and everything associated with the platforms.

Author: Yuksel-Rudy

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