Iea22

src/pynumad/analysis/beam_utils.py

@@ -358,8 +358,11 @@ def get_bd_source_radial_stations(bd_sum_file_name):
 
     with open(bd_sum_file_name) as blade_yaml:
         data = yaml.load(blade_yaml,Loader=yaml.FullLoader)
-    node_location=list(np.array(data['Init_Nodes_E1'])[:,2]/np.array(data['Init_Nodes_E1'])[-1,2])
-    return node_location
+    if 'trap' in data['Quadrature_method'].lower():
+        return None #The the source locations are target locations
+    else:
+        node_location=list(np.array(data['Init_Nodes_E1'])[:,2]/np.array(data['Init_Nodes_E1'])[-1,2])
+        return node_location
 
 # def get_bd_source_radial_stations(qp_location,blade_length,nqp): #Adds root and top to qp data
 

src/pynumad/analysis/make_models.py

@@ -223,11 +223,13 @@ def write_sierra_sd_model(template_file,wt_name,station_list,blade,materials_use
 
         blockLines+=f'block {material.name}\n'
         blockLines+=f'    material {material.name}\n'
+        blockLines+=f'    coordinate from_transfer\n'
         blockLines+=f'end \n'
         blockLines+='\n\n' 
 
     copy_and_replace(template_file, sierra_file_name,
         {
+            'ROOT_STATION': 'station'+str(station_list[0]).zfill(3),
             'BLADE_MATERIALS': materialLines,
             'IN_MESH':wt_name+'.g',
             'BLADE_BLOCKS': blockLines,