Permalink
Browse files

Fixes made to Fargo (didn't update B3i in 2d xy shearing sheet) and R…

…emap_exchange_ix1 (j loop limits in Step 5e were wrong)
  • Loading branch information...
mkunz
mkunz committed Jun 9, 2015
1 parent 75ead55 commit 3fcd9e930b2ed4a1349991070c3cb55ec54fb279
Showing with 3 additions and 1 deletion.
  1. +2 −0 src/bvals_shear.c
  2. +1 −1 src/particles/exchange.c
View
@@ -5304,6 +5304,8 @@ void Fargo(DomainS *pD)
pG->B2i[k][j][i] += (pG->dx2/pG->dx1)*
(ri[i+1]*FargoFlx[k][i+1][jj].U[NFARGO-1] - ri[i]*FargoFlx[k][i][jj].U[NFARGO-1]);
#endif
pG->B3i[k][j][i] +=
(FargoFlx[k][i][jj+1].U[NFARGO-2] - FargoFlx[k][i][jj].U[NFARGO-2]);
}
pG->B1i[k][j][ie+1] -=
(FargoFlx[k][ie+1][jj+1].U[NFARGO-1]-FargoFlx[k][ie+1][jj].U[NFARGO-1]);
View
@@ -936,7 +936,7 @@ void Remap_exchange_ix1(DomainS *pD)
for (k=kb; k<=kt; k++) {
for(i=0; i<NExc+NOfst; i++){
for (j=js; j<=je-joverlap; j++) {
for (j=js+joverlap; j<=je; j++) {
for (n=0; n<NVar; n++) {
*(pSnd++) = TempZns[k][i][j].U[n];
}

0 comments on commit 3fcd9e9

Please sign in to comment.