diff options
| author | lanfer <lanfer@41e88fdd-2190-4c69-9c84-4659c8cf322e> | 1999-11-01 13:57:13 +0000 |
|---|---|---|
| committer | lanfer <lanfer@41e88fdd-2190-4c69-9c84-4659c8cf322e> | 1999-11-01 13:57:13 +0000 |
| commit | 7a48676e5e26928d4e16418f8cfe5b463c8a78ed (patch) | |
| tree | da6d67169669f240497dcee6fde2fc2bfecf7ce3 | |
| parent | 9314db503bc9c969300c807f2017dac4bd95dca5 (diff) | |
took out debug statement, using ell database, still is different then analytic data
git-svn-id: http://svn.cactuscode.org/arrangements/CactusWave/IDScalarWaveElliptic/trunk@14 41e88fdd-2190-4c69-9c84-4659c8cf322e
| -rw-r--r-- | src/SourceData.F77 | 8 |
1 files changed, 6 insertions, 2 deletions
diff --git a/src/SourceData.F77 b/src/SourceData.F77 index 43bfe34..14bcc70 100644 --- a/src/SourceData.F77 +++ b/src/SourceData.F77 @@ -41,10 +41,8 @@ c where Q is the total charge and R is the sphere radius Mcoeff(i,j,k) = 0.0d0 - print *,r(i,j,k),radius if (r(i,j,k) <= radius) then Ncoeff(i,j,k) = charge_factor - print *,"Here" else Ncoeff(i,j,k) = 0d0 end if @@ -69,6 +67,12 @@ c where Q is the total charge and R is the sphere radius c Call elliptic solver to fill out phi +c call Ell_SetStringKey(ierr, "yes", "EllLinFlat::Bbd::Robin") +c write (*,*) "SourceData: SetString ",ierr + + call Ell_SetRealKey(ierr, 1.0d0, "EllLinFlat::Bnd::Robin::V0") + call Ell_SetIntKey (ierr, 1, "EllLinFlat::Bnd::Robin::V1") + call Ell_LinFlatSolver(ierr,cctkGH, & iphi, & iMcoeff, iNcoeff, |