diff options
| author | rhaas <rhaas@ac85fae7-cede-4708-beff-ae01c7fa1c26> | 2013-01-11 15:04:15 +0000 |
|---|---|---|
| committer | rhaas <rhaas@ac85fae7-cede-4708-beff-ae01c7fa1c26> | 2013-01-11 15:04:15 +0000 |
| commit | fd05b50e5ccf463bc86ffcd19c520754e5bc03e7 (patch) | |
| tree | ff04fa3b3e3afc6461b0977149c47143ed714d73 /src | |
| parent | 7c39d4964e8b048681a402347b2dc7c4d57f8fee (diff) | |
GRHydro_InitData: add parameter to offset radial origin
r_code = r_KS - r0 with r0=1M by default. This avoids placing the
central singularity on the grid.
From: Roland Haas <roland.haas@physics.gatech.edu>
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/GRHydro_InitData/trunk@196 ac85fae7-cede-4708-beff-ae01c7fa1c26
Diffstat (limited to 'src')
| -rw-r--r-- | src/GRHydro_BondiM.c | 15 |
1 files changed, 15 insertions, 0 deletions
diff --git a/src/GRHydro_BondiM.c b/src/GRHydro_BondiM.c index ef7c1c5..f947abb 100644 --- a/src/GRHydro_BondiM.c +++ b/src/GRHydro_BondiM.c @@ -1097,6 +1097,21 @@ void GRHydro_BondiM(CCTK_ARGUMENTS) xpos[YY] = y[i] ; xpos[ZZ] = z[i] ; + if (bondi_radial_offset > 0.) { + double rspher_orig = sqrt(SQR(xpos[XX])+SQR(xpos[YY])+SQR(xpos[ZZ])); + double rspher_new = rspher_orig + bondi_radial_offset; + if(rspher_orig < SMALL_BONDI) { + xpos[XX] = rspher_new; + xpos[YY] = 0.; + xpos[ZZ] = 0.; + } else { + for(int n = XX ; n <= ZZ ; n++) { + xpos[n] *= rspher_new/rspher_orig; + } + } + } + + switch( coord_type ) { case COORD_BOYERLINDQUIST : |