Small fixes unload and addCablesConnections; update examples readme

Author: lsirkis

examples/README.md

@@ -7,7 +7,17 @@ The examples are described below:
 
 ## Numbered Examples
 
-The numbered example folders (01-08) contain small isolated examples showing ontologies and driver files for basic functionalities. Within each folder, a number of sub-functionalities are shown individually. Currently, only the 01_Visualization folder is developed, as this is a work in progress. 
+The numbered example folders (01-08) contain small isolated examples showing ontologies and driver files for basic functionalities. Within each folder, a number of sub-functionalities are shown individually. Currently, only the 01-03 example folders are developed, as this is a work in progress. 
+The numbered example folders are as follows:
+- 01_Visualization : 2D and 3D plotting of various components and array information
+- 02_Mooring_Analysis_MoorPy : Run and analyze mooring and cable simulations
+- 03_Frequency_Domain_Analysis_RAFT : Run and analyze floating platform/ floating wind simulations
+- 04_Geography : Load and set up geographical information such as soils, bathymetry, lease coordinates, etc
+- 05_Anchor_Capacity_and_Sizing : Calculate anchor capacity and size anchors to meet safety requirements
+- 06_Wakes_Analysis_FLORIS : Run and analyze wind turbine wake simulations
+- 07_Cost_Analysis : Calculate total cost of array components, output spreadsheet breakdown of costs
+- 08_Cable_Rerouting : Reroute intra-array cables around mooring lines and anchors
+- 09_Failure_Analysis : Run a failure analysis of the array
 
 ## FAModel project from a yaml file example
 This example shows how to make an FAModel project from a yaml file. The sample shows various possibilities for mooring and anchoring including shared moorings and shared anchors. It also shows some functionalities of FAModel including:

famodel/project.py

@@ -845,7 +845,6 @@ def loadSite(self, site):
                 self.grid_x = np.array(site['bathymetry']['x'])
                 self.grid_y = np.array(site['bathymetry']['y'])
                 self.grid_depth = np.array(site['bathymetry']['depths'])
-                # breakpoint()
                 # self.setGrid(xs,ys)
             else:
                 # assume a flat bathymetry
@@ -1000,8 +999,6 @@ def seabedIntersect(self, r, u):
 
             r_i = r_i + u*fac  # updated seabed crossing estimate 
 
-            if any(np.isnan(r_i)):
-                breakpoint()
 
         return r_i
 
@@ -1734,7 +1731,8 @@ def cableDesignInterpolation(self,depth,cables):
     def addCablesConnections(self,connDict,cableType_def='dynamic_cable_66',oss=False,
                              substation_r=[None],ss_id=200,id_method='location',
                              keep_old_cables=False, connect_ss=True, 
-                             cableConfig=None, configType=0,heading_buffer=30):
+                             cableConfig=None, configType=0,heading_buffer=30,
+                             route_anchors=True):
         '''Adds cables and connects them to existing platforms/substations based on info in connDict
         Designed to work with cable optimization output designed by Michael Biglu
 
@@ -1868,7 +1866,7 @@ def addCablesConnections(self,connDict,cableType_def='dynamic_cable_66',oss=Fals
 
 
             # create cable object
-            cab = Cable(cableType_def+str(i+lcab),d=dd)
+            cab = Cable('cable'+str(i+lcab),d=dd)
             self.cableList[cab.id] = cab
 
             # update upstream turbines property
@@ -1930,7 +1928,10 @@ def addCablesConnections(self,connDict,cableType_def='dynamic_cable_66',oss=Fals
                 # update routing
                 cab.subcomponents[cs].updateRouting(coords) # also updates static and general cable lengths
 
-        route_around_anchors(self)
+
+        if route_anchors:
+            route_around_anchors(self)
+
 
 
     def updatePositions(self):
@@ -3832,21 +3833,27 @@ def unload(self,file='project.yaml'):
         mapAnchNames = {}
         mscs = {}
         arrayMoor = []
+        pf_types = []
         for i,anch in enumerate(self.anchorList.values()):  
             newanch = True
             name = anch.dd['name'] if 'name' in anch.dd else str(len(anchConfigs))
             if len(anch.attachments)>1:
                 # shared anchor, add to arrayAnch list
                 arrayAnch.append([anch.id,name,anch.r[0],anch.r[1],0])
+            # add mass if available
+            aad = deepcopy(anch.dd['design'])
+            if anch.mass is not None and anch.mass>0:
+                aad['mass'] = anch.mass
+            if 'type' in anch.dd:
+                aad['type'] = anch.dd['type']
             if anchConfigs:
-                ac = [an for an,ad in anchConfigs.items() if ad==anch.dd['design']]
+                ac = [an for an,ad in anchConfigs.items() if ad==aad]
                 if len(ac)>0:
                     newanch = False
                     name = ac[0] # reset name to matching config name
             if newanch:
-                anchConfigs[name] = dict(anch.dd['design'])
-                if anch.mass is not None and anch.mass>0: 
-                    anchConfigs[name]['mass'] = anch.mass
+                anchConfigs[name] = aad
+
             mapAnchNames[anch.id] = name
 
         # build out platform info
@@ -3933,6 +3940,20 @@ def unload(self,file='project.yaml'):
                     mscs[mname] = msys
             else:
                 mname = 0
+                
+            if not 'type' in pf.dd:
+                pf_type_info = [pf.rFair, pf.zFair, pf.entity]
+                if not pf_type_info in pf_types:
+                    pf_types.append(pf_type_info)
+                    if not self.platformTypes:
+                        self.platformTypes = []
+                    pf.dd['type'] = len(self.platformTypes)
+                    self.platformTypes.append({'rFair': pf.rFair,
+                                                         'zFair': pf.zFair,
+                                                         'type': pf.entity})
+                else:
+                    tt = [n for n,ps in enumerate(pf_types) if ps==pf_type_info]
+                    pf.dd['type'] = tt[0]