diff options
author | ryoji <ryoji@b6f3ac56-194f-0410-8878-cdf6079d7f1b> | 2001-09-07 10:17:21 +0000 |
---|---|---|
committer | ryoji <ryoji@b6f3ac56-194f-0410-8878-cdf6079d7f1b> | 2001-09-07 10:17:21 +0000 |
commit | e75c73d4293614a07ee02daa7fb77269ea70bdcd (patch) | |
tree | 0eca6e9173b156b7fa1a7d1de81953a1e3d7e97c | |
parent | 4c4fa457b9edba2712c225a9df949eb75d977aaf (diff) |
put sergio's new initial data.
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/IDAxiOddBrillBH/trunk@33 b6f3ac56-194f-0410-8878-cdf6079d7f1b
-rw-r--r-- | param.ccl | 3 | ||||
-rw-r--r-- | src/IDAxiOddBrillBH.F | 10 |
2 files changed, 11 insertions, 2 deletions
@@ -58,6 +58,9 @@ BOOLEAN brandt_seidel "Brand-Seidel initial data" { } "no" +BOOLEAN sergio "Sergio initial data" +{ +} "no" BOOLEAN verbose "" { diff --git a/src/IDAxiOddBrillBH.F b/src/IDAxiOddBrillBH.F index f06a6d0..90d0cbe 100644 --- a/src/IDAxiOddBrillBH.F +++ b/src/IDAxiOddBrillBH.F @@ -53,9 +53,9 @@ c Perhaps this and others should go into cctk.h real*8 t1,t2,t3,t4 real*8 gtil,dngtil,dnngtil,dnnngtil,dnnnngtil,dnnnnngtil, $ dnnnnnngtil,dnnnnnnngtil - real*8 rhsmax,rmax,odd_get2d,adm,Jmom + real*8 rhsmax,rmax,odd_get2d,adm,Jmom,mass,a,R,rBL real*8,parameter :: rbh_tol = 1.0d-7,rbh_eps = 1.0d-10 - integer,parameter :: itmax = 50 + integer,parameter :: itmax = 100 real*8 pi,zero,one,two CCTK_REAL :: xmin,xmax,ymin,ymax,zmin,zmax integer :: ne, nq @@ -129,6 +129,9 @@ c #include "Development/IDAxiOddBrillBH/src/gauss.x" if(brandt_seidel == 1) then #include "Development/IDAxiOddBrillBH/src/ksq_bs.x" + else if (sergio == 1) then +#include "Development/IDAxiOddBrillBH/src/kerr.x" +#include "Development/IDAxiOddBrillBH/src/ksq_sergio.x" else #include "Development/IDAxiOddBrillBH/src/ksq_axi_rdbh.x" endif @@ -456,6 +459,9 @@ c Extrinsic Curvture:Cactus only needs physical extrinsic curvature #include "Development/IDAxiOddBrillBH/src/cgauss.x" if (brandt_seidel == 1) then #include "Development/IDAxiOddBrillBH/src/kij_bs.x" + else if (sergio == 1) then +#include "Development/IDAxiOddBrillBH/src/ckerr.x" +#include "Development/IDAxiOddBrillBH/src/kij_sergio.x" else #include "Development/IDAxiOddBrillBH/src/kij_axi_rdbh.x" endif |