diff options
Diffstat (limited to 'src/SourceData.F77')
-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, |