Skip to content

Commit aa5a942

Browse files
committed
bug fix
1 parent a44e099 commit aa5a942

File tree

1 file changed

+5
-5
lines changed

1 file changed

+5
-5
lines changed

classFiles/SVEclass.m

Lines changed: 5 additions & 5 deletions
Original file line numberDiff line numberDiff line change
@@ -825,9 +825,9 @@ function TransPostProcessor(obj)
825825
function [spatialDirection,slicedPlane] = plotTransientFront(obj)
826826
% [spatialDirection,slicedPlane] = plotTransientFront():
827827
% Plots transient front for all time steps.
828-
load(['TopologyBundle_',num2str(obj.nx),'_',num2str(obj.realizationNumber),'.mat'],'NodeCoords')
829-
load(['TopologyBundle_',num2str(obj.nx),'_',num2str(obj.realizationNumber),'.mat'],'meshProperties')
830-
load([obj.transientName,'.mat']);
828+
load([obj.path2Realization,'TopologyBundle_',num2str(obj.nx),'_',num2str(obj.realizationNumber),'.mat'],'NodeCoords')
829+
load([obj.path2Realization,'TopologyBundle_',num2str(obj.nx),'_',num2str(obj.realizationNumber),'.mat'],'meshProperties')
830+
load([obj.path2Realization,obj.transientName,'.mat']);
831831

832832
[ff nts] = size(a_store);
833833
slicedPlane.mean = zeros(meshProperties.nx,1);
@@ -836,13 +836,13 @@ function TransPostProcessor(obj)
836836

837837
for iTime=1:nts
838838
for i=0:meshProperties.nx-1
839-
[nodeSlice foo] = find(abs(NodeCoords(:,2)-i*meshProperties.dx)<10*eps);
839+
[nodeSlice foo] = find(abs(NodeCoords(:,2)-i*meshProperties.dx)<100*eps);
840840
slicedPlane.mean(i+1) = mean(a_store(nodeSlice,iTime));
841841
slicedPlane.std(i+1) = std(a_store(nodeSlice,iTime));
842842
end
843843
hold on
844844
plot(spatialDirection,slicedPlane.mean)
845-
xlabel('spatial direction [cm]')
845+
xlabel('spatial direction')
846846
end
847847
end
848848
% linear elasticity

0 commit comments

Comments
 (0)