Skip to content

Commit

Permalink
Fix TSoutput
Browse files Browse the repository at this point in the history
  • Loading branch information
CyprienBosserelle committed Apr 16, 2021
1 parent 22037d0 commit c37c79e
Showing 1 changed file with 3 additions and 3 deletions.
6 changes: 3 additions & 3 deletions src/Mainloop.cu
Original file line number Diff line number Diff line change
Expand Up @@ -178,15 +178,15 @@ template <class T> void pointoutputstep(Param XParam, Loop<T> &XLoop, Model<T> X
//
Pointout stepread;

stepread.time = XParam.totaltime;
stepread.time = XLoop.totaltime;
stepread.zs = 0.0;// That is a bit useless
stepread.h = 0.0;
stepread.u = 0.0;
stepread.v = 0.0;
XLoop.TSAllout[o].push_back(stepread);


storeTSout << <gridDim, blockDim, 0 >> > (XParam,(int)XParam.TSnodesout.size(), o, XLoop.nTSsteps, XParam.TSnodesout[o].block, XParam.TSnodesout[o].i, XParam.TSnodesout[o].j, XModel.bndblk.Tsout, XModel_g.evolv, XModel_g.TSstore);
storeTSout << <gridDim, blockDim, 0 >> > (XParam,(int)XParam.TSnodesout.size(), o, XLoop.nTSsteps, XParam.TSnodesout[o].block, XParam.TSnodesout[o].i, XParam.TSnodesout[o].j, XModel_g.bndblk.Tsout, XModel_g.evolv, XModel_g.TSstore);
}
CUDA_CHECK(cudaDeviceSynchronize());
}
Expand All @@ -199,7 +199,7 @@ template <class T> void pointoutputstep(Param XParam, Loop<T> &XLoop, Model<T> X

int i = memloc(XParam.halowidth, XParam.blkmemwidth, XParam.TSnodesout[o].i, XParam.TSnodesout[o].j, XParam.TSnodesout[o].block);

stepread.time = XParam.totaltime;
stepread.time = XLoop.totaltime;
stepread.zs = XModel.evolv.zs[i];
stepread.h = XModel.evolv.h[i];;
stepread.u = XModel.evolv.u[i];;
Expand Down

0 comments on commit c37c79e

Please sign in to comment.