diff options
| -rw-r--r-- | src/RadiationBoundary.c | 8 |
1 files changed, 2 insertions, 6 deletions
diff --git a/src/RadiationBoundary.c b/src/RadiationBoundary.c index 2ba4262..2288921 100644 --- a/src/RadiationBoundary.c +++ b/src/RadiationBoundary.c @@ -996,12 +996,10 @@ void CCTK_FCALL CCTK_FNAME (BndRadiativeVN) CCTK_REAL _radius0_inv = 1/_radius[_0], \ _radius1_inv = 1/_radius[_1]; \ \ - \ - if (radpower > 0) \ + if (radpower >= 0) \ { \ CCTK_REAL H; \ \ - \ H = 0.25 * radpower * dxyz[0] \ * (_xyz[_0]*SQR (_radius0_inv) + _xyz[_1]*SQR (_radius1_inv)); \ H = (1 + H) / (1 - H); \ @@ -1108,12 +1106,10 @@ void CCTK_FCALL CCTK_FNAME (BndRadiativeVN) CCTK_REAL _radius0_inv = 1/_radius[_0], \ _radius1_inv = 1/_radius[_1]; \ \ - \ - if (radpower > 0) \ + if (radpower >= 0) \ { \ CCTK_REAL H; \ \ - \ H = 0.25 * radpower * dxyz[0] \ * (_xyz[_0]*SQR (_radius0_inv) + _xyz[_1]*SQR (_radius1_inv)); \ H = (1 - H) / (1 + H); \ |