aboutsummaryrefslogtreecommitdiff
diff options
context:
space:
mode:
authortradke <tradke@f5a6acaf-da7d-456b-b0a8-35edbc60b392>2000-03-08 10:28:08 +0000
committertradke <tradke@f5a6acaf-da7d-456b-b0a8-35edbc60b392>2000-03-08 10:28:08 +0000
commitc5d41fcaf08d57ae14a565654aa830be6da6d9d4 (patch)
treebb2a4aa803b2d887f2d5736231a3c49ab45f596b
parent45d6be576b3e146a447193355a29b1dc5dc1c632 (diff)
Made parameters binary_size, binary_omega, binary_charge, and binary_radius
steerable. Thomas git-svn-id: http://svn.cactuscode.org/arrangements/CactusWave/IDScalarWave/trunk@32 f5a6acaf-da7d-456b-b0a8-35edbc60b392
-rw-r--r--param.ccl8
-rw-r--r--src/WaveBinary.F776
2 files changed, 8 insertions, 6 deletions
diff --git a/param.ccl b/param.ccl
index 85f45c8..6df8a22 100644
--- a/param.ccl
+++ b/param.ccl
@@ -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