From f5fda658d49efce252918ad5bfbc85c1aed8edb6 Mon Sep 17 00:00:00 2001 From: Ian Hinder Date: Wed, 2 May 2012 18:06:55 +0200 Subject: ML_BSSN_Helper: Correct background value of phi in NewRad boundary condition Previously the conformalMethod parameter was not taken into account. This will change numerical results. --- m/prototype/ML_BSSN_Helper/param.ccl | 1 + m/prototype/ML_BSSN_Helper/src/NewRad.c | 2 +- 2 files changed, 2 insertions(+), 1 deletion(-) (limited to 'm') diff --git a/m/prototype/ML_BSSN_Helper/param.ccl b/m/prototype/ML_BSSN_Helper/param.ccl index f9267d0..d0efc9d 100644 --- a/m/prototype/ML_BSSN_Helper/param.ccl +++ b/m/prototype/ML_BSSN_Helper/param.ccl @@ -21,6 +21,7 @@ USES CCTK_INT rhs_timelevels USES KEYWORD calculate_ADMBase_variables_at USES CCTK_REAL harmonicF +USES CCTK_REAL conformalMethod EXTENDS KEYWORD my_initial_boundary_condition { diff --git a/m/prototype/ML_BSSN_Helper/src/NewRad.c b/m/prototype/ML_BSSN_Helper/src/NewRad.c index c5529a1..f5a00e7 100644 --- a/m/prototype/ML_BSSN_Helper/src/NewRad.c +++ b/m/prototype/ML_BSSN_Helper/src/NewRad.c @@ -24,7 +24,7 @@ ML_BSSN_NewRad (CCTK_ARGUMENTS) CCTK_REAL const v0 = sqrt (harmonicF); - newrad (cctkGH, phi , phirhs , 0.0, v0 ); + newrad (cctkGH, phi , phirhs , conformalMethod ? 1 : 0, v0 ); newrad (cctkGH, gt11 , gt11rhs , 1.0, 1.0); newrad (cctkGH, gt12 , gt12rhs , 0.0, 1.0); -- cgit v1.2.3