aboutsummaryrefslogtreecommitdiff
diff options
context:
space:
mode:
authorschnetter <schnetter@850bcc8b-0e4f-0410-8c26-8d28fbf1eda9>2004-12-29 22:46:54 +0000
committerschnetter <schnetter@850bcc8b-0e4f-0410-8c26-8d28fbf1eda9>2004-12-29 22:46:54 +0000
commit7e72d1c33bc8aeed446ce1a0247189f9e53773b9 (patch)
treebabe65b90bedf4fe58cb64104f8330bbb1358b01
parentee2bb42cece14c3a96f872d1644d34fb836aa999 (diff)
Pass dt in addition to dx to the dissipating function
git-svn-id: http://svn.cactuscode.org/arrangements/CactusNumerical/Dissipation/trunk@6 850bcc8b-0e4f-0410-8c26-8d28fbf1eda9
-rw-r--r--src/apply_dissipation.F774
-rw-r--r--src/dissipation.c6
2 files changed, 6 insertions, 4 deletions
diff --git a/src/apply_dissipation.F77 b/src/apply_dissipation.F77
index 38d55f7..b5d6939 100644
--- a/src/apply_dissipation.F77
+++ b/src/apply_dissipation.F77
@@ -2,12 +2,12 @@ c $Header$
#include "cctk.h"
- subroutine apply_dissipation (var, rhs, ni, nj, nk, dx, order, epsdis)
+ subroutine apply_dissipation (var, rhs, ni, nj, nk, dx, dt, order, epsdis)
implicit none
integer ni, nj, nk
CCTK_REAL var(ni,nj,nk), rhs(ni,nj,nk)
- CCTK_REAL dx(3)
+ CCTK_REAL dx(3), dt
CCTK_INT order
CCTK_REAL epsdis
diff --git a/src/dissipation.c b/src/dissipation.c
index fd59e31..44a8eb2 100644
--- a/src/dissipation.c
+++ b/src/dissipation.c
@@ -14,6 +14,7 @@ CCTK_FNAME(apply_dissipation) (CCTK_REAL const * const var,
int const * const nj,
int const * const nk,
CCTK_REAL const * const dx,
+ CCTK_REAL const * const dt,
CCTK_INT const * const order,
CCTK_REAL const * const epsdis);
@@ -40,7 +41,7 @@ apply (int const varindex, char const * const optstring, void * const arg)
cGroup vardata, rhsdata;
CCTK_REAL const * varptr;
CCTK_REAL * rhsptr;
- CCTK_REAL dx[3];
+ CCTK_REAL dt, dx[3];
int n;
int d;
int ierr;
@@ -53,6 +54,7 @@ apply (int const varindex, char const * const optstring, void * const arg)
}
}
+ dt = CCTK_DELTA_TIME;
for (d=0; d<3; ++d) {
dx[d] = CCTK_DELTA_SPACE(d);
}
@@ -104,5 +106,5 @@ apply (int const varindex, char const * const optstring, void * const arg)
CCTK_FNAME(apply_dissipation)
(varptr, rhsptr, &cctk_lsh[0], &cctk_lsh[1], &cctk_lsh[2],
- dx, &order, &epsdis);
+ dx, &dt, &order, &epsdis);
}