Issue 173 coordinate system integration

.travis.yml

@@ -9,20 +9,22 @@ install:
   - pip install scipy
   - pip install xlsxwriter
   - pip install xlrd
-  - pip install psycopg2-binary
+#  - pip install psycopg2-binaryyam
+  - pip install psycopg2==2.7.5
   - pip install sqlalchemy
   - pip install pytest
 
-
 env:
   global:
     - LANDBOSSE_INPUT_DIR: $TRAVIS_BUILD_DIR/project_input_template
     - LANDBOSSE_OUTPUT_DIR: $TRAVIS_BUILD_DIR/project_input_template
 
 # command to run tests
 script:
+  - pytest $TRAVIS_BUILD_DIR/landbosse/landbosse_api/API_test/API_test.py
   - python $TRAVIS_BUILD_DIR/main.py --validate
 
 branches:
   only:
+    - pip_installable
     - develop

MANIFEST.in

@@ -1,3 +1,6 @@
 include README.md
 include CHANGELOG.md
 include LICENSE.txt
+include landbosse/landbosse_api/project_list.xlsx
+include landbosse/landbosse_api/project_data/ge15_public.xlsx
+include landbosse/landbosse_api/project_data/az_rolling.srw

README.md

@@ -1,4 +1,6 @@
-[![Build Status](https://travis-ci.com/WISDEM/LandBOSSE.svg?branch=issue_136_add_tests)](https://travis-ci.com/WISDEM/LandBOSSE)
+
+<!--- [![Build Status](https://travis-ci.com/WISDEM/LandBOSSE.svg?branch=pip_installable)](https://travis-ci.com/WISDEM/LandBOSSE)  -->
+
 # LandBOSSE
 
 ## Welcome to LandBOSSE!

landbosse/excelio/XlsxReader.py

@@ -361,32 +361,63 @@ def create_master_input_dictionary(self, project_data_dataframes, project_parame
         # Incomplete project dict will hold the input dictionary
         # configurations.
         incomplete_input_dict = dict()
+        incomplete_input_dict['error'] = dict()
 
         # Read all project_data sheets.
         # The erection module takes in a bunch of keys and values under the
         # 'project_data' key in the incomplete_input_dict
 
         # Apply the labor multipliers
-        labor_cost_multiplier = project_parameters['Labor cost multiplier']
+        incomplete_input_dict['labor_cost_multiplier'] = project_parameters['Labor cost multiplier']
+        labor_cost_multiplier = incomplete_input_dict['labor_cost_multiplier']
         self.apply_labor_multiplier_to_project_data_dict(project_data_dataframes, labor_cost_multiplier)
+        # Apply the Collection mode
+        if 'Collection mode' not in project_parameters:
+            incomplete_input_dict['collection_mode'] = 'auto'
+            project_parameters['Collection mode'] = 'auto'
+        else:
+            incomplete_input_dict['collection_mode'] = project_parameters['Collection mode']
+
+
+        incomplete_input_dict['row_spacing_rotor_diameters'] = project_parameters['Row spacing (times rotor diameter)']
+        incomplete_input_dict['turbine_spacing_rotor_diameters'] = project_parameters['Turbine spacing (times rotor diameter)']
+
+        if project_parameters['Collection mode'] == 'manual':
+            incomplete_input_dict['collection_layout'] = project_data_dataframes['collection_layout']
+            erection_input_worksheets = [
+                'crane_specs',
+                'equip',
+                'crew',
+                'equip_price',
+                'crew_price',
+                'material_price',
+                'components',
+                'collection_layout'
+            ]
+        else:
+            erection_input_worksheets = [
+                'crane_specs',
+                'equip',
+                'crew',
+                'equip_price',
+                'crew_price',
+                'material_price',
+                'components'
+            ]
+        if project_parameters['Collection mode'] == 'manual' and (
+                len(project_data_dataframes['collection_layout']['Adjacency matrix']) != (
+                project_parameters['Number of turbines'] + 1)):
+            exit('ERROR: mismatch between # turbines and # turbine locations')
 
-        erection_input_worksheets = [
-            'crane_specs',
-            'equip',
-            'crew',
-            'equip_price',
-            'crew_price',
-            'material_price',
-            'components'
-        ]
 
         erection_project_data_dict = dict()
         for worksheet in erection_input_worksheets:
             erection_project_data_dict[worksheet] = project_data_dataframes[worksheet]
 
         # Add the erection project data to the incomplete_input_dict
         incomplete_input_dict['project_data'] = erection_project_data_dict
-
+
+
         # Get the first set of data
         incomplete_input_dict['rsmeans'] = project_data_dataframes['rsmeans']
         incomplete_input_dict['site_facility_building_area_df'] = project_data_dataframes['site_facility_building_area']
@@ -429,7 +460,6 @@ def create_master_input_dictionary(self, project_data_dataframes, project_parame
             project_parameters['Breakpoint between base and topping (percent)']
         incomplete_input_dict['fuel_usd_per_gal'] = project_parameters['Fuel cost USD per gal']
         incomplete_input_dict['rate_of_deliveries'] = project_parameters['Rate of deliveries (turbines per week)']
-        incomplete_input_dict['turbine_spacing_rotor_diameters'] = project_parameters['Turbine spacing (times rotor diameter)']
         incomplete_input_dict['depth'] = project_parameters['Foundation depth m']
         incomplete_input_dict['rated_thrust_N'] = project_parameters['Rated Thrust (N)']
         incomplete_input_dict['bearing_pressure_n_m2'] = project_parameters['Bearing Pressure (n/m2)']
@@ -440,6 +470,7 @@ def create_master_input_dictionary(self, project_data_dataframes, project_parame
             incomplete_input_dict['road_distributed_wind'] = True
         else:
             incomplete_input_dict['road_distributed_wind'] = False
+
         incomplete_input_dict['site_prep_area_m2'] = project_parameters['Site prep area for Distributed wind (m2)']
         incomplete_input_dict['road_length_adder_m'] = project_parameters['Road length adder (m)']
         incomplete_input_dict['fraction_new_roads'] = project_parameters['Percent of roads that will be constructed']
@@ -448,7 +479,9 @@ def create_master_input_dictionary(self, project_data_dataframes, project_parame
         incomplete_input_dict['plant_capacity_MW'] = project_parameters['Turbine rating MW'] * project_parameters['Number of turbines']
         incomplete_input_dict['row_spacing_rotor_diameters'] = project_parameters['Row spacing (times rotor diameter)']
         incomplete_input_dict['user_defined_distance_to_grid_connection'] = project_parameters['Flag for user-defined home run trench length (0 = no; 1 = yes)']
-        incomplete_input_dict['distance_to_grid_connection_km'] = project_parameters['Combined Homerun Trench Length to Substation (km)']
+        incomplete_input_dict['distance_to_interconnect_mi'] = project_parameters['Distance to interconnect (miles)']
+        incomplete_input_dict['distance_to_grid_connection_km'] = project_parameters[
+            'Combined Homerun Trench Length to Substation (km)']
         incomplete_input_dict['crew'] = incomplete_input_dict['project_data']['crew']
         incomplete_input_dict['crew_cost'] = incomplete_input_dict['project_data']['crew_price']
 
@@ -468,7 +501,8 @@ def create_master_input_dictionary(self, project_data_dataframes, project_parame
             'Combined Homerun Trench Length to Substation (km)']
 
         # Add inputs for transmission & Substation modules:
-        incomplete_input_dict['distance_to_interconnect_mi'] = project_parameters['Distance to interconnect (miles)']
+
+        #incomplete_input_dict['distance_to_interconnect_km'] = project_parameters['Combined Homerun Trench Length to Substation (km)']
         incomplete_input_dict['interconnect_voltage_kV'] = project_parameters['Interconnect Voltage (kV)']
         new_switchyard = True
         if project_parameters['New Switchyard (y/n)'] == 'y':
@@ -488,6 +522,7 @@ def create_master_input_dictionary(self, project_data_dataframes, project_parame
         incomplete_input_dict['overtime_multiplier'] = project_parameters['Overtime multiplier']
         incomplete_input_dict['allow_same_flag'] = True if project_parameters['Allow same flag'] == 'y' else False
 
+
         override_total_mgmt_cost_col_name = 'Override total management cost for distributed (0 does not override)'
         if override_total_mgmt_cost_col_name in project_parameters and project_parameters[override_total_mgmt_cost_col_name] > 0:
             incomplete_input_dict['override_total_management_cost'] = \
@@ -581,7 +616,6 @@ def apply_cost_and_scaling_modifications_to_project_parameters(self, project_par
         hub_height_m = project_parameters['Hub height m']
         flag_use_user_homerun = project_parameters['Flag for user-defined home run trench length (0 = no; 1 = yes)']
         nameplate = project_parameters['Turbine rating MW']
-
         distance_to_interconnect_mi = 0.0 if project_size_MW <= 20 else (0.009375 * project_size_MW + 0.625)
         interconnect_voltage_kV = 0.4398 * project_size_MW + 60.204
         new_switchyard_y_n = 'n' if project_size_MW <= 40 else 'y'

landbosse/landbosse_api/API_test/API_test.py

@@ -0,0 +1,132 @@
+import os
+from ..run import *
+from unittest import TestCase
+
+
+class TestLandBOSSE_API(TestCase):
+    def setUp(self):
+        self.file_path = os.path.dirname(__file__) + '/project_data/az_rolling.srw'
+        api_inputs = dict()
+        api_inputs['num_turbines'] = 100
+        api_inputs['project_id'] = 'foundation_validation_ge15'
+        self.results = run_landbosse(api_inputs)
+        print(self.results)
+
+    def test_total_BOS_cost(self):\
+        self.assertEqual(44208796.36106637, self.results['total_bos_cost'])
+
+    def test_total_management_cost(self):
+        delta = 0.002    # max % tolerance in expected and actual value
+        self.assertAlmostEqual(10570956.97,
+                               self.results['total_management_cost'],
+                               delta=int(delta * 10570956.97))
+
+        self.assertAlmostEqual(188371.90,
+                               self.results['insurance_usd'],
+                               delta=int(delta * 188371.90))
+
+        self.assertAlmostEqual(400723.14,
+                               self.results['construction_permitting_usd'],
+                               delta=int(delta * 400723.14))
+
+        self.assertAlmostEqual(1887763.492,
+                               self.results['project_management_usd'],
+                               delta=int(delta * 1887763.492))
+
+        self.assertAlmostEqual(336378.394,
+                               self.results['bonding_usd'],
+                               delta=int(delta * 336378.394))
+
+        self.assertAlmostEqual(4379646.689,
+                               self.results['markup_contingency_usd'],
+                               delta=int(delta * 4379646.689))
+
+        self.assertAlmostEqual(2071325,
+                               self.results['engineering_usd'],
+                               delta=int(delta * 2071325))
+
+        self.assertAlmostEqual(1306748.35,
+                               self.results['site_facility_usd'],
+                               delta=int(delta * 1306748.35))
+
+        total_management_cost = self.results['insurance_usd'] + \
+                                self.results['construction_permitting_usd'] + \
+                                self.results['project_management_usd'] + \
+                                self.results['bonding_usd'] + \
+                                self.results['markup_contingency_usd'] + \
+                                self.results['engineering_usd'] + \
+                                self.results['site_facility_usd']
+
+        self.assertAlmostEqual(self.results['total_management_cost'],
+                               total_management_cost,
+                               delta=int(delta * self.results['total_management_cost']))
+
+    def test_siteprep_cost(self):
+        self.assertEqual(self.results['total_sitepreparation_cost'], 3209538.616279147)
+        self.assertEqual(round(self.results['sitepreparation_equipment_rental_usd'], 2), 412013.84)
+        self.assertEqual(round(self.results['sitepreparation_labor_usd'], 4), 609463.1297)
+        self.assertEqual(round(self.results['sitepreparation_material_usd'], 4), 797751.5346)
+        self.assertEqual(round(self.results['sitepreparation_mobilization_usd'], 5), 83009.35197)
+        self.assertEqual(self.results['sitepreparation_other_usd'], 1307300.76)
+
+        total_siteprep_cost = self.results['sitepreparation_equipment_rental_usd'] + \
+                              self.results['sitepreparation_labor_usd'] + \
+                              self.results['sitepreparation_material_usd'] + \
+                              self.results['sitepreparation_mobilization_usd'] + \
+                              self.results['sitepreparation_other_usd']
+
+        self.assertEqual(self.results['total_sitepreparation_cost'], total_siteprep_cost)
+
+    def test_foundation_cost(self):
+        self.assertEqual(self.results['total_foundation_cost'], 10036157.011254452)
+        self.assertEqual(self.results['foundation_equipment_rental_usd'], 307553.56983444514)
+        self.assertEqual(round(self.results['foundation_labor_usd'], 3), 3677659.987)
+        self.assertEqual(round(self.results['foundation_material_usd'], 3), 5573031.216)
+        self.assertEqual(round(self.results['foundation_mobilization_usd'], 4), 477912.2386)
+
+        total_foundation_cost = self.results['foundation_equipment_rental_usd'] + \
+                                self.results['foundation_labor_usd'] + \
+                                self.results['foundation_material_usd'] + \
+                                self.results['foundation_mobilization_usd']
+
+        self.assertEqual(self.results['total_foundation_cost'], total_foundation_cost)
+
+    def test_total_erection_cost(self):
+        self.assertEqual(6811700.7960274285, self.results['total_erection_cost'])
+        self.assertEqual(748911.4927891536, self.results['erection_equipment_rental_usd'])
+        self.assertEqual(5134653.8032382745, self.results['erection_labor_usd'])
+        self.assertEqual(0, self.results['erection_material_usd'])
+        self.assertEqual(0, self.results['erection_other_usd'])
+        self.assertEqual(28825.5, self.results['erection_fuel_usd'])
+        self.assertEqual(899310, self.results['erection_mobilization_usd'])
+
+        total_erection_cost = self.results['erection_equipment_rental_usd'] + \
+                              self.results['erection_labor_usd'] + \
+                              self.results['erection_material_usd'] + \
+                              self.results['erection_other_usd'] + \
+                              self.results['erection_fuel_usd'] + \
+                              self.results['erection_mobilization_usd']
+
+        self.assertEqual(self.results['total_erection_cost'], total_erection_cost)
+
+    def test_total_gridconnection_cost(self):
+           self.assertEqual(4084775.152898776, self.results['total_gridconnection_cost'])
+
+    def test_total_collection_cost(self):
+           self.assertEqual(self.results['total_collection_cost'], 4869435.746468595)
+           self.assertEqual(self.results['collection_equipment_rental_usd'], 426614.6862960762)
+           self.assertEqual(self.results['collection_labor_usd'], 1402888.0267842393)
+           self.assertEqual(self.results['collection_material_usd'], 2808055.140699299)
+           self.assertEqual(self.results['collection_mobilization_usd'], 231877.89268898076)
+
+    def test_total_substation_cost(self):
+           self.assertEqual(4859182.072082901, self.results['total_substation_cost'])
+
+    def test_read_weather_data(self):
+        weather_data = read_weather_data(self.file_path)
+        self.assertEqual(len(weather_data), 8760)
+
+    def test_read_data(self):
+        input_dictionary = dict()
+        project_list = read_data(input_dictionary)
+        self.assertEqual(len(project_list.columns), 47)