diff options
| -rw-r--r-- | src/SourceData.F | 3 |
1 files changed, 2 insertions, 1 deletions
diff --git a/src/SourceData.F b/src/SourceData.F index 2913efb..8bd4222 100644 --- a/src/SourceData.F +++ b/src/SourceData.F @@ -26,6 +26,7 @@ c where Q is the total charge and R is the sphere radius pi = 4.0*atan(1.0) charge_factor = 4.0d0*pi*charge*3.0d0/(4.0d0*pi*radius**3) + print *,charge_factor do k=1, cctk_lsh(3) do j=1, cctk_lsh(2) @@ -58,7 +59,7 @@ c where Q is the total charge and R is the sphere radius RelTol(3)=-1 c Call to elliptic solver will go here - call Ell_LinFlatSolver(cctkGH, + call Ell_LinFlatSolver(ierr,cctkGH, & iphi, & iMcoeff, iNcoeff, & AbsTol, RelTol, |