diff options
| author | rhaas <rhaas@ac85fae7-cede-4708-beff-ae01c7fa1c26> | 2013-03-28 01:48:53 +0000 |
|---|---|---|
| committer | rhaas <rhaas@ac85fae7-cede-4708-beff-ae01c7fa1c26> | 2013-03-28 01:48:53 +0000 |
| commit | 7d5bbf2886014e63d0011f31d5a643d919f58f47 (patch) | |
| tree | 0d192b695a0915e37f9cb49392fe896f5ce768ad /src | |
| parent | fd05b50e5ccf463bc86ffcd19c520754e5bc03e7 (diff) | |
GRHydro: add parameters to set Bondi solution only in certain range.
From: Roland Haas <rhaas@tapir.caltech.edu>
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/GRHydro_InitData/trunk@197 ac85fae7-cede-4708-beff-ae01c7fa1c26
Diffstat (limited to 'src')
| -rw-r--r-- | src/GRHydro_BondiM.c | 17 |
1 files changed, 17 insertions, 0 deletions
diff --git a/src/GRHydro_BondiM.c b/src/GRHydro_BondiM.c index f947abb..ca693ad 100644 --- a/src/GRHydro_BondiM.c +++ b/src/GRHydro_BondiM.c @@ -84,6 +84,7 @@ //some math helpers #define SQR(x) ((x)*(x)) +#define CUBE(x) ((x)*(x)*(x)) static CCTK_REAL Mdot, rs, vs_sq, vs, cs_sq, cs, rhos, hs, K, Qdot, gamma_eos, r_sol; @@ -902,6 +903,19 @@ static void calc_vel_bondi( CCTK_REAL vtmp, CCTK_REAL x[NDIM], CCTK_REAL x_sphe ***********************************************************************************/ void GRHydro_BondiM(CCTK_ARGUMENTS) { + GRHydro_BondiM_Internal(CCTK_PASS_CTOC, 1e100, -1e100); +} + +void GRHydro_BondiM_Range(CCTK_ARGUMENTS) +{ + DECLARE_CCTK_PARAMETERS; + GRHydro_BondiM_Internal(CCTK_PASS_CTOC, bondi_freeze_inner_radius, + bondi_freeze_outer_radius); +} + + +void GRHydro_BondiM_Internal(CCTK_ARGUMENTS, CCTK_REAL range_min, CCTK_REAL range_max) +{ DECLARE_CCTK_ARGUMENTS; DECLARE_CCTK_PARAMETERS; CCTK_REAL det; @@ -1134,6 +1148,9 @@ void GRHydro_BondiM(CCTK_ARGUMENTS) rspher = x_spher[RR]; + if(rspher > range_min && rspher < range_max) + continue; + /* Find nearest point in the Bondi solution : */ j = (int) ( 0.5 + (log10(rspher) - logrmin) / dlogr ) ; |