diff options
-rw-r--r-- | param.ccl | 8 | ||||
-rw-r--r-- | src/WaveBinary.F77 | 6 |
2 files changed, 8 insertions, 6 deletions
@@ -36,22 +36,22 @@ KEYWORD binary_verbose "Rotating binary source verbose" "no" :: "no output" } "no" -REAL binary_size "Radial extension of the binary source" +REAL binary_size "Radial extension of the binary source" STEERABLE = ALWAYS { 0.0: :: "Some positive value" } 0.5 -REAL binary_omega "Frequency of the circular binary orbit" +REAL binary_omega "Frequency of the circular binary orbit" STEERABLE = ALWAYS { 0.0: :: "Some positive value" } 2.0 -REAL binary_charge "Charge of source" +REAL binary_charge "Charge of source" STEERABLE = ALWAYS { : :: "No restriction" } 0.1 -REAL binary_radius "Radius of the circular binary orbit" +REAL binary_radius "Radius of the circular binary orbit" STEERABLE = ALWAYS { 0.0: :: "Some positive value" } 2.0 diff --git a/src/WaveBinary.F77 b/src/WaveBinary.F77 index 26f296f..b457d1a 100644 --- a/src/WaveBinary.F77 +++ b/src/WaveBinary.F77 @@ -147,10 +147,12 @@ c Things to do on first iteration SAVE firstcall,charge_factor - if (firstcall.eq.1) then +c Because binary_charge and binary_size are a steerable parameters now, +c charge_factor needs to be recomputed at every iteration. +c if (firstcall.eq.1) then charge_factor = 3.0d0*binary_charge/ $ (4.0d0*3.1415*binary_size*binary_size*binary_size) - end if +c end if |