Skip to content

Commit

Permalink
Update body and wake FoR multiplication
Browse files Browse the repository at this point in the history
silly mistake really - position in Gframe of multibody is offset by cga_0 not cga(t)
  • Loading branch information
kccwing authored Dec 9, 2024
1 parent 69a4cdf commit 0a79ca4
Showing 1 changed file with 2 additions and 2 deletions.
4 changes: 2 additions & 2 deletions sharpy/postproc/aerogridplot.py
Original file line number Diff line number Diff line change
Expand Up @@ -208,7 +208,7 @@ def plot_body(self):
# cga = algebra.euler2rot([0, self.data.trimmed_values[0], 0])
# coords[counter, :] += np.dot(cga, struct_tstep.for_pos[0:3])
# except AttributeError:
coords[counter, :] += np.dot(struct_tstep.cga(), struct_tstep.for_pos[0:3])
coords[counter, :] += np.dot(self.data.structure.timestep_info[0].cga(), struct_tstep.for_pos[0:3])
if self.settings['include_forward_motion']:
coords[counter, 0] -= self.settings['dt']*self.ts*self.settings['u_inf']

Expand Down Expand Up @@ -331,7 +331,7 @@ def plot_wake(self):
# cga = algebra.euler2rot([0, self.data.trimmed_values[0], 0])
# coords[counter, :] += np.dot(cga, self.data.structure.timestep_info[self.ts].for_pos[0:3])
# except AttributeError:
coords[counter, :] += np.dot(self.data.structure.timestep_info[self.ts].cga(), self.data.structure.timestep_info[self.ts].for_pos[0:3])
coords[counter, :] += np.dot(self.data.structure.timestep_info[0].cga(), self.data.structure.timestep_info[self.ts].for_pos[0:3])
if self.settings['include_forward_motion']:
coords[counter, 0] -= self.settings['dt']*self.ts*self.settings['u_inf']

Expand Down

0 comments on commit 0a79ca4

Please sign in to comment.