diff options
author | miguel <miguel@a678b1cf-93e1-4b43-a69d-d43939e66649> | 2001-04-12 10:51:59 +0000 |
---|---|---|
committer | miguel <miguel@a678b1cf-93e1-4b43-a69d-d43939e66649> | 2001-04-12 10:51:59 +0000 |
commit | 6fd09e4d86b3a13278c83874f976c3ba2aa7cd42 (patch) | |
tree | e8575326bc17006dcdcfccb0c4b3851131a34b04 /src/brilldata.F | |
parent | 65c6f4ada5ed70fb476c58de5866b98eb7c911e8 (diff) |
Cleaning up.
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/IDBrillData/trunk@36 a678b1cf-93e1-4b43-a69d-d43939e66649
Diffstat (limited to 'src/brilldata.F')
-rw-r--r-- | src/brilldata.F | 25 |
1 files changed, 12 insertions, 13 deletions
diff --git a/src/brilldata.F b/src/brilldata.F index 24a745e..253ab15 100644 --- a/src/brilldata.F +++ b/src/brilldata.F @@ -1,12 +1,11 @@ + #include "cctk.h" #include "cctk_Parameters.h" #include "cctk_Arguments.h" - subroutine brilldata(CCTK_FARGUMENTS) +#include "CactusEinstein/Einstein/src/Einstein.h" -c Author: Carsten Gundlach. -c -c Driver routine for calculating Brill wave initial data. + subroutine brilldata(CCTK_FARGUMENTS) implicit none @@ -14,12 +13,11 @@ c Driver routine for calculating Brill wave initial data. DECLARE_CCTK_PARAMETERS DECLARE_CCTK_FUNCTIONS -c Declare arrays for the linear elliptic solver call: - - integer Mlin_index,Nsrc_index,field_index,ierr - integer metpsi_index(7) + integer Mlin_index,Nsrc_index,field_index,ierr + integer metpsi_index(7) + integer ierr - CCTK_REAL AbsTol(3), RelTol(3) + CCTK_REAL AbsTol(3),RelTol(3) c Set up background metric and coefficients for linear solve. @@ -41,8 +39,8 @@ c conformal factor. First some preliminaries. call CCTK_VarIndex(metpsi_index(7),"einstein::psi") call CCTK_VarIndex(field_index,"IDBrillData::brillpsi") - call CCTK_VarIndex(Mlin_index,"IDBrillData::Mlinear") - call CCTK_VarIndex(Nsrc_index,"IDBrillData::Nsource") + call CCTK_VarIndex(Mlin_index, "IDBrillData::Mlinear") + call CCTK_VarIndex(Nsrc_index, "IDBrillData::Nsource") AbsTol(1)= brill_thresh AbsTol(2)= -1 @@ -83,9 +81,10 @@ c Elliptic solver. . field_index,Mlin_index,Nsrc_index,AbsTol,RelTol,"bam") end if -c Synchronize conformal factor. +c Synchronization and symmetry boundaries. - call CCTK_SyncGroup(cctkGH,"einstein::confac") + call CCTK_SyncGroup(cctkGH,"IDBrillData::brillconf") + call CartSymGN(ierr,cctkGH,"IDBrillData::brillconf") c Reconstruct physical metric. |