wawi 0.0.11__py3-none-any.whl → 0.0.12__py3-none-any.whl

examples/3 Software interfacing/Abaqus model export/bergsoysund-export.py

@@ -42,25 +42,29 @@ with open('element.json', 'w') as f:
 phi_full_disp = wawi.ext.abq.get_nodal_phi(step_obj, nodes, flatten_components=True, field_outputs=['UT', 'UR'])
 
 #%% Get pontoon data
-node_matrix_pontoons = wawi.ext.abq.get_element_matrices(region_hydro, obj=part)
+pontoon_sets = ['P1', 'P2', 'P3', 'P4', 'P5', 'P6', 'P7']
+pontoon_types = ['p17', 'p26', 'p345', 'p345', 'p345', 'p26', 'p17']
+rotations = np.array([13.944407004001892, 9.2962713360012632, 4.6481356680006325, 0.0,
+                     -4.6481356680006334, -9.2962713360012632, -13.944407004001883])
+
+node_matrix_pontoons = np.vstack([wawi.ext.abq.get_element_matrices(db.rootAssembly.nodeSets[pset], obj=part) for pset in pontoon_sets])
 node_labels_pontoons = node_matrix_pontoons[:,0]
-nodes_pontoons = wawi.ext.abq.create_list_of_nodes(part, node_labels_pontoons)
-phi_h = wawi.ext.abq.get_nodal_phi(step_obj, nodes_pontoons, flatten_components=True, field_outputs=['UT', 'UR'])
 
 # Export pontoon.json
-pontoon_types = ['ptype_1']*7   
-rotations = -np.array([13.944407004001892, 9.2962713360012632, 4.6481356680006325, 0.0,
-                     -4.6481356680006334, -9.2962713360012632, -13.944407004001883])
-
 pontoon_data = OrderedDict()
 
 for ix, node in enumerate(node_labels_pontoons):
-    key = 'P'+str (ix+1)
+    key = pontoon_sets[ix]
     pontoon_data[key] = dict(coordinates=node_matrix_pontoons[ix, 1:].tolist(),
                             node=node,
                             rotation=rotations[ix],
                             pontoon_type=pontoon_types[ix])
 
+
+nodes_pontoons = wawi.ext.abq.create_list_of_nodes(part, node_labels_pontoons)
+phi_h = wawi.ext.abq.get_nodal_phi(step_obj, nodes_pontoons, flatten_components=True, field_outputs=['UT', 'UR'])
+
+
 with open('pontoon.json', 'w') as f:
     json.dump(pontoon_data, f)
 

examples/3 Software interfacing/Abaqus model export/hardanger-export.py

@@ -0,0 +1,47 @@
+from sys import path
+from collections import OrderedDict
+
+path.append('C:/Users/knutankv/git-repos/wawi/')   # this is an easy quick fix to enable importing wawi package in Abaqus environment
+savefolder = 'C:/Temp/'
+
+# NOTE: Only linear *beam* elements are supported in relevant sets!!
+
+import wawi.ext.abq
+import json
+import numpy as np
+
+#%% Get database object (ODB) 
+db = wawi.ext.abq.get_db('odb')
+
+#%% Definitions
+frequency_step = 'Step-11'  # define set name to extract modal properties from
+part = db.rootAssembly.instances['PART-1-1']    # define part 
+step_obj = db.steps[frequency_step]             
+
+part.ElementSet('ALL', part.elements)   # get all elements (for full element and node matrix)
+
+#%% Grab regions
+region_full = part.elementSets['ALL']
+
+#%% Get modal parameters
+fn, m = wawi.ext.abq.get_modal_parameters(frequency_step)
+
+#%% Get wind elements and mode shapes
+node_matrix, element_matrix = wawi.ext.abq.get_element_matrices(region_full, obj=part)
+node_labels = node_matrix[:,0]
+nodes = wawi.ext.abq.create_list_of_nodes(part, node_labels)
+
+# Export element definitions as json
+el_data = dict(node_matrix=node_matrix.tolist(), 
+               element_matrix=element_matrix.tolist())
+
+with open('element.json', 'w') as f:
+    json.dump(el_data, f)
+
+phi_full_disp = wawi.ext.abq.get_nodal_phi(step_obj, nodes, flatten_components=True, field_outputs=['U', 'UR'])
+
+## ------------------- EXPORT MODES ----------------
+modal_data = dict(omega_n=(fn*2*np.pi).tolist(), m=m.tolist(), phi=dict(full=phi_full_disp.tolist()))
+
+with open('modal.json', 'w') as f:
+    json.dump(modal_data, f)

tests/test_IABSE_step11a.py

@@ -3,10 +3,10 @@ import pytest
 import numpy as np
 from math import isclose
 
-# use the local wawi (Github folder) instead of the installed version (remove this when using the installed version)
-import sys
-import os
-sys.path.insert(0, os.path.abspath('C:\\Users\\aksef\\Documents\\GitHub\\wawi'))
+# use the local wawi 
+#import sys
+#import os
+#sys.path.insert(0, r'C:\Users\aksef\OneDrive - Multiconsult\project\wawi\wawi_keep')
 
 # import functions
 from wawi.io import import_folder
@@ -108,7 +108,6 @@ def test_mean_speed_45():
     # PSD response
     S_exp = np.transpose(np.array( [ [0.0167, -0.007j], [0.007j, 0.00293] ] ))
 
-
     model = import_folder(model_folder)
     model.aero.windstate = iabse_11a_windstate(V)
 
@@ -121,11 +120,11 @@ def test_mean_speed_45():
     omega = np.array([0.001, f])*np.pi*2
     HH = eval_3d_fun(model.get_frf_fun(opt = 1), omega)
 
-    model.aero.sections['girder0'].Admittance = lambda fred: np.full((4, 3), davenport(fred))
+    model.aero.sections['girder0'].admittance = lambda fred: np.full((4, 3), davenport(fred))
     model.run_freqsim(omega,
                     include_selfexcited=['aero'], 
                     include_action=['aero'],
-                    print_progress=False, merge_aero_sections=True)
+                    print_progress=False)
 
     global_dof_ix = np.array([2,3])
     S = model.get_result_psd(key='full', 
@@ -170,12 +169,14 @@ def test_response_spectra(Vbuff, RS_vert_exp, RS_tors_exp):
     model.aero.windstate = iabse_11a_windstate(Vbuff)
     model.modal_dry.xi0 = .3e-2
     model.aero.sections['girder0'].ADs = ADs(**flatplate_ads())
-    model.aero.sections['girder0'].Admittance = lambda fred: np.full((4, 3), davenport(fred))
+    model.aero.sections['girder0'].admittance = lambda fred: np.full((4, 3), davenport(fred))
     
     model.run_freqsim(omega_axis,
                     include_selfexcited=['aero'], 
                     include_action=['aero'],
-                    print_progress=False, merge_aero_sections=True)
+                    print_progress=False, 
+                    #merge_aero_sections=True
+                    )
 
     global_dof_ix = np.array([2,3])
     S = model.get_result_psd(key='full', 
@@ -204,12 +205,12 @@ def test_RMS_response(Vbuff, RMS_vert_exp, RMS_tors_exp):
     model.aero.windstate = iabse_11a_windstate(Vbuff)
     model.modal_dry.xi0 = .3e-2
     model.aero.sections['girder0'].ADs = ADs(**flatplate_ads())
-    model.aero.sections['girder0'].Admittance = lambda fred: np.full((4, 3), davenport(fred))
+    model.aero.sections['girder0'].admittance = lambda fred: np.full((4, 3), davenport(fred))
 
     model.run_freqsim(omega_axis,
                     include_selfexcited=['aero'], 
                     include_action=['aero'],
-                    print_progress=False, merge_aero_sections=True)
+                    print_progress=False)
     
     stds = model.get_result_std(key = 'full')
 

tests/test_IABSE_step11c.py

@@ -5,9 +5,9 @@ from math import isclose
 import dill
 
 # use the local wawi (Github folder) instead of the installed version (remove this when using the installed version)
-import sys
-import os
-sys.path.insert(0, os.path.abspath('C:\\Users\\aksef\\Documents\\GitHub\\wawi'))
+#import sys
+#import os
+#sys.path.insert(0, os.path.abspath('C:\\Users\\aksef\\Documents\\GitHub\\wawi'))
 
 # import functions
 from wawi.io import import_folder
@@ -166,11 +166,11 @@ def test_mean_speed_45():
     omega = np.array([0.001, f])*np.pi*2
     HH = eval_3d_fun(model.get_frf_fun(opt = 1), omega)
 
-    model.aero.sections['girder0'].Admittance = lambda fred: np.full((4, 3), davenport(fred))
+    model.aero.sections['girder0'].admittance = lambda fred: np.full((4, 3), davenport(fred))
     model.run_freqsim(omega,
                     include_selfexcited=['aero'], 
                     include_action=['aero'],
-                    print_progress=False, merge_aero_sections=True)
+                    print_progress=False)
 
     global_dof_ix = np.array([1,2,3])
     S = model.get_result_psd(key='full', 
@@ -225,13 +225,13 @@ def test_response_spectra(Vbuff, RS_horz_exp, RS_vert_exp, RS_tors_exp):
         AD_funs = dill.load(file)
     AD_s = AD_dict(AD_funs)
     model.aero.sections['girder0'].ADs = ADs(**AD_s)
-    model.aero.sections['girder0'].Admittance = lambda fred: np.full((4, 3), davenport(fred))
+    model.aero.sections['girder0'].admittance = lambda fred: np.full((4, 3), davenport(fred))
 
     # run analysis
     model.run_freqsim(omega_axis,
                       include_selfexcited=['aero'], 
                       include_action=['aero'],
-                      print_progress=False, merge_aero_sections=True)
+                      print_progress=False)
     
     # extract PSD
     global_dof_ix = np.array([1,2,3])
@@ -265,13 +265,13 @@ def test_RMS_response(Vbuff, RMS_horz_exp, RMS_vert_exp, RMS_tors_exp):
         AD_funs = dill.load(file)
     AD_s = AD_dict(AD_funs)
     model.aero.sections['girder0'].ADs = ADs(**AD_s)
-    model.aero.sections['girder0'].Admittance = lambda fred: np.full((4, 3), davenport(fred))
+    model.aero.sections['girder0'].admittance = lambda fred: np.full((4, 3), davenport(fred))
     
     # run analysis
     model.run_freqsim(omega_axis,
                     include_selfexcited=['aero'], 
                     include_action=['aero'],
-                    print_progress=False, merge_aero_sections=True)
+                    print_progress=False)
     
     # RMS responses
     stds = model.get_result_std(key = 'full')

tests/test_IABSE_step2a.py

@@ -4,9 +4,9 @@ from math import isclose
 import dill
 
 # use the local wawi (Github folder) instead of the installed version (remove this when using the installed version)
-import sys
-import os
-sys.path.insert(0, os.path.abspath('C:\\Users\\aksef\\Documents\\GitHub\\wawi'))
+#import sys
+#import os
+#sys.path.insert(0, os.path.abspath('C:\\Users\\aksef\\Documents\\GitHub\\wawi'))
 
 from wawi.io import import_folder
 from wawi.model import Windstate
@@ -184,14 +184,14 @@ def test_RS_spectra():
     model.aero.windstate = iabse_2a_windstate(V)
     # admittance
     for key in model.aero.sections:
-        model.aero.sections[key].Admittance = lambda fred: np.full((4, 3), davenport(fred))
+        model.aero.sections[key].admittance = lambda fred: np.full((4, 3), davenport(fred))
 
     # run analysis
     model.run_freqsim(omega,
                         include_selfexcited=['aero'], 
                         include_action=['aero'],
-                        print_progress=False, merge_aero_sections=True,
-                                            ensure_at_peaks=False) # keep the given frequency axis
+                        print_progress=False,
+                        ensure_at_peaks=False) # keep the given frequency axis
 
 
     # extract PSD - at the midpsan
@@ -237,7 +237,7 @@ def test_RMS_response(Vbuff, RMS_horz_exp1, RMS_vert_exp1, RMS_tors_exp1, RMS_ho
     model.aero.windstate = iabse_2a_windstate(Vbuff)
     # admittance
     for key in model.aero.sections:
-        model.aero.sections[key].Admittance = lambda fred: np.full((4, 3), davenport(fred))
+        model.aero.sections[key].admittance = lambda fred: np.full((4, 3), davenport(fred))
 
     #model.run_eig(w_initial=model.modal_dry.omega_n.tolist(), freq_kind=True, itmax=100) # alters integration 
 
@@ -245,7 +245,7 @@ def test_RMS_response(Vbuff, RMS_horz_exp1, RMS_vert_exp1, RMS_tors_exp1, RMS_ho
     model.run_freqsim(omega,
                         include_selfexcited=['aero'], 
                         include_action=['aero'],
-                        print_progress=False, merge_aero_sections=True) 
+                        print_progress=False) 
     # RMS responses
     stds = model.get_result_std(key = 'full')
     # global dofs

wawi/__init__.py

@@ -3,7 +3,7 @@
 '''
 __pdoc__ = {'wawi.ext.abq': False}
 
-__version__ = "0.0.11"
+__version__ = "0.0.12"
 
 # Other packages
 import numpy as np

wawi/modal.py

@@ -192,7 +192,7 @@ def iteig(K, C, M, omega=None, omega_ref=0, input_functions=True, itmax=None, to
         q[:, m] = phi[:,0]
 
         if print_progress:
-            pp(m+2,2*ndofs, sym='>', postfix=' finished with iterative modal analysis.')
+            pp(m+2,2*ndofs+1, sym='>', postfix=' finished with iterative modal analysis.')
 
     if print_progress:
         print(' ')
@@ -323,7 +323,7 @@ def iteig_freq(K, C, M, omega=None, itmax=15, reference_omega=0, input_functions
         q[:, m] = qi[:, m] 
 
         if print_progress:
-            pp(m+2,2*ndofs, sym='>', postfix=' finished with iterative modal analysis.')
+            pp(m+2 ,2*ndofs+1, sym='>', postfix=' finished with iterative modal analysis.')
 
     if print_progress:
         print(' ')

wawi/model/_model.py

@@ -16,6 +16,8 @@ from wawi.structural import freqsim
 from wawi.tools import print_progress as pp
 from wawi.wave import stochastic_linearize, waveaction_fft, harmonic_linearize
 from wawi.wind import windaction, windaction_static
+# from wawi.io import import_folder
+
 import dill
 
 
@@ -266,6 +268,17 @@ class Model:
             
         self.Cquad_lin = 0.0
     
+    # Not used due to circular referencing.
+    # @staticmethod
+    # def from_folder(*args, **kwargs):
+    #     '''
+    #     Alternative constructor. Passes all inputs to `wawi.io.import_folder`.
+    #     '''
+
+    #     model = import_folder(*args, **kwargs)
+        
+    #     return model
+    
 
     @staticmethod    #alternative constructor
     def from_wwi(path):
@@ -568,7 +581,7 @@ class Model:
             pl.add_point_labels(np.vstack([wind_origin]), [f'U={self.aero.windstate.U0:.2f} m/s from heading {self.aero.windstate.direction:.1f} deg (CW)'])
             
         # Plot wave arrow        
-        if plot_wave_direction and self.hydro.seastate is not None:
+        if plot_wave_direction and self.hydro is not None and self.hydro.seastate is not None:
             if self.hydro.seastate.homogeneous:
                 if wave_origin=='center':
                     wave_origin = origin*1
@@ -1335,4 +1348,10 @@ class Model:
         return ix
     
     def get_node_ixs(self, nodelabels):
-        return [self.get_node_ix(nodelabel) for nodelabel in nodelabels]
+        return [self.get_node_ix(nodelabel) for nodelabel in nodelabels]
+    
+    def get_aerosection_phi_and_x(self, key):
+        ixs = self.aero.phi_ixs[key]
+        x = np.vstack([node.coordinates for node in self.aero.eldef[key].nodes])
+
+        return ixs, x

{wawi-0.0.11.dist-info → wawi-0.0.12.dist-info}/METADATA

@@ -1,6 +1,6 @@
 Metadata-Version: 2.4
 Name: wawi
-Version: 0.0.11
+Version: 0.0.12
 Summary: WAve and WInd response prediction
 Author-email: "Knut A. Kvåle" <knut.a.kvale@ntnu.no>, Ole Øiseth <ole.oiseth@ntnu.no>, Aksel Fenerci <aksel.fenerci@ntnu.no>, Øivind Wiig Petersen <oyvind.w.petersen@ntnu.no>
 License: MIT License

{wawi-0.0.11.dist-info → wawi-0.0.12.dist-info}/RECORD

@@ -1,15 +1,15 @@
-examples/3 Software interfacing/Abaqus model export/bergsoysund-export.py,sha256=-1jrDNmevOcB1IV0iLOiW0z-QiNkOv8GkGk-7K8CcK0,2606
-examples/3 Software interfacing/Abaqus model export/test.py,sha256=kAmol3XZQ2jkO2qujEQrqx38Q24LbzrYZzqE5ng5-E8,2188
-tests/test_IABSE_step11a.py,sha256=guRsaP9gpERvSp0Ui3hO-WH764BcIe6Qx8vC4QVTAoU,8013
-tests/test_IABSE_step11c.py,sha256=U_NcL6BAmnQG7MSjkr8KYk9hH8ZPdZgg17pM9ugUBzo,10412
-tests/test_IABSE_step2a.py,sha256=L5AyTOJZ_sPmcmW_4h-sQfRwvVW71LZkXc1yz33DABM,10203
+examples/3 Software interfacing/Abaqus model export/bergsoysund-export.py,sha256=3HSV009z5oqIUgM6B7xucUpGeh2f8SzC9FSiyP0zd7A,2760
+examples/3 Software interfacing/Abaqus model export/hardanger-export.py,sha256=YiqS-TxnkydwVHE7k-57nP6CDsxqSXgWR3f8UytAubY,1603
+tests/test_IABSE_step11a.py,sha256=zTWOxR78F0PNhNqiqodEcL65vdJXyMP78nkJumhy-RA,7915
+tests/test_IABSE_step11c.py,sha256=TfxPlP70dhujKlSYbA-Gj4U-EJmYZ3jt6yjvX9wSXcA,10337
+tests/test_IABSE_step2a.py,sha256=9pAd1em2IEY4k5N-jAigV3VGup0QsbUxnwSiziuMbUg,10134
 tests/test_wind.py,sha256=r4rx6f8Tn-u4fn4gGPBMrL_JJQ2SVHGQiQ9sQuMQCPo,1707
-wawi/__init__.py,sha256=rH8H6W3AN1iZb8rU9K_8k-mycRSvQNxFq97Biy0-hzw,183
+wawi/__init__.py,sha256=7vUvWFZzW67YG3yHPiBO_bsphhLLgEF4N7wdUkJeEq8,183
 wawi/fe.py,sha256=22QKI1GlfsG7o_TpFXaKJfzmbO2_2zdIMaoJmaIZdmY,4001
 wawi/general.py,sha256=xHRoDkcchrL6Y7pTUqGsjBHh8YBLvX9dYcNXXCjQap8,13787
 wawi/identification.py,sha256=bVB6EVRR6J39OO9ckuzNJ6f0FwIo4cLqfYgrsIN89TE,1748
 wawi/io.py,sha256=cr8XLF7bwpTCxpZ84YjOyWjc5uuu-JjGv8vY9UKH-HI,25327
-wawi/modal.py,sha256=WjNGFsk0C3tYRy18Q9WNRCatmGJtq1JSv0WrkGV02Eo,20247
+wawi/modal.py,sha256=L1Rs_YVEADOmsQL-hJ43eb8Nt5VpYuB1VHOjYK2O9X0,20252
 wawi/plot.py,sha256=jllJcjZxTBqjzBoT4k9jLXVUnie8oqNr8371IJvCd3c,19791
 wawi/prob.py,sha256=0nCdKdwkNf4M6sHyCZuYlt06gD0NmqRNfl4KesgySWA,215
 wawi/random.py,sha256=MHPpyTlRJSJFkCmeTAmw4Q5K1BPoFVb0Nxg0jDhkuIM,871
@@ -29,10 +29,10 @@ wawi/model/__init__.py,sha256=u6B-dugP76LBkOUVMs0KAoQ81PeRHwFL0M8MbNeAaJA,218
 wawi/model/_aero.py,sha256=tvWMjMoVi9ZVBc3NZ_s-wSofbEk1cc5jf5s-Swv2RxQ,9337
 wawi/model/_dry.py,sha256=KpmFlSinoY6DrSyc3Y0M8w1-cCC7VCFep-uzcqZsHz4,3940
 wawi/model/_hydro.py,sha256=1tUj19OgvEjItNtBxXuPHf0zWtKv_ioTrOcgl9Ov3dg,26960
-wawi/model/_model.py,sha256=34qblCqGpIKwbsa38CJ3a8KbHYyX8LDWOFQhdRwK98Q,52504
+wawi/model/_model.py,sha256=ynZkyoU00crTAlu7sN0CmJsn-J8vwK7BDbMCh0pb-hI,53061
 wawi/model/_screening.py,sha256=NRYkKq928z2lqMSUTpbQLls04td_9R_4dhkjU3Gv1oQ,3716
-wawi-0.0.11.dist-info/licenses/LICENSE,sha256=bH1aWhrNbbPLrYnVFRaoYYzcUr-figHjry-kGB7Tc54,1076
-wawi-0.0.11.dist-info/METADATA,sha256=TPaFfAkztMCRaw4F8U5JT6XqOK2IjY3orqmBmguLoaY,6530
-wawi-0.0.11.dist-info/WHEEL,sha256=CmyFI0kx5cdEMTLiONQRbGQwjIoR1aIYB7eCAQ4KPJ0,91
-wawi-0.0.11.dist-info/top_level.txt,sha256=Nk5G_ZwgZRCb9ZMWZdr1M3QIskX6kCnlqeMl67N3zg8,20
-wawi-0.0.11.dist-info/RECORD,,
+wawi-0.0.12.dist-info/licenses/LICENSE,sha256=bH1aWhrNbbPLrYnVFRaoYYzcUr-figHjry-kGB7Tc54,1076
+wawi-0.0.12.dist-info/METADATA,sha256=lXTUPijzmyd6QOJ2LmCPJaF7H_NWkysiOZkJkhlKqlU,6530
+wawi-0.0.12.dist-info/WHEEL,sha256=CmyFI0kx5cdEMTLiONQRbGQwjIoR1aIYB7eCAQ4KPJ0,91
+wawi-0.0.12.dist-info/top_level.txt,sha256=Nk5G_ZwgZRCb9ZMWZdr1M3QIskX6kCnlqeMl67N3zg8,20
+wawi-0.0.12.dist-info/RECORD,,

examples/3 Software interfacing/Abaqus model export/test.py

@@ -1,67 +0,0 @@
-from sys import path
-from collections import OrderedDict
-
-path.append('C:/Users/knutankv/git-repos/wawi/')   # this is an easy quick fix to enable importing wawi package in Abaqus environment
-savefolder = 'C:/Temp/'
-
-import wawi.ext.abq
-import json
-import numpy as np
-
-#%% Get (O) database object
-db = wawi.ext.abq.get_db('odb')
-
-#%% Definitions
-frequency_step = 'Step-6'
-part = db.rootAssembly.instances['BRIDGE-1']
-step_obj = db.steps[frequency_step]
-
-if 'ALL' not in part.elementSets:   #CREATE SET OF ALL ELEMENTS IN PART
-    part.ElementSet('ALL', part.elements)
-
-#%% Grab regions
-region_full = part.elementSets['ALL']
-region_hydro = part.nodeSets['SPRING']
-
-#%% Get modal parameters
-fn, m = wawi.ext.abq.get_modal_parameters(frequency_step)
-
-#%% Get wind elements and mode shapes
-node_matrix, element_matrix = wawi.ext.abq.get_element_matrices(region_full, obj=part)
-node_labels = node_matrix[:,0]
-
-# Export element definitions as json
-el_data = dict(node_matrix=node_matrix.tolist(), element_matrix=element_matrix.tolist())
-
-with open('element.json', 'w') as f:
-    json.dump(el_data, f)
-
-phi_full_disp = wawi.ext.abq.get_nodal_phi(step_obj, node_labels, flatten_components=True)
-
-#%% Get pontoon data
-node_matrix_pontoons = wawi.ext.abq.get_element_matrices(region_hydro, obj=part, sort_nodes_fun=None)
-node_labels = node_matrix_pontoons[:,0]
-phi_h = wawi.ext.abq.get_nodal_phi(step_obj, node_labels, flatten_components=True)
-
-# Export pontoon.json
-pontoon_types = ['ptype_1']*48
-rotations = np.zeros(48)
-
-pontoon_data = OrderedDict()
-
-for ix, node in enumerate(node_labels):
-    key = 'P'+str (ix+1)
-    pontoon_data[key] = dict(coordinates=node_matrix_pontoons[ix, 1:].tolist(),
-                            node=node,
-                            rotation=rotations[ix],
-                            pontoon_type=pontoon_types[ix])
-
-with open('pontoon.json', 'w') as f:
-    json.dump(pontoon_data, f)
-
-## ------------------- EXPORT MODES ----------------
-modal_data = dict(omega_n=(fn*2*np.pi).tolist(), m=m.tolist(), phi=dict(full=phi_full_disp.tolist(),
-                                     hydro=phi_h.tolist()))
-
-with open('modal.json', 'w') as f:
-    json.dump(modal_data, f)