From c5d41fcaf08d57ae14a565654aa830be6da6d9d4 Mon Sep 17 00:00:00 2001 From: tradke Date: Wed, 8 Mar 2000 10:28:08 +0000 Subject: 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 --- param.ccl | 8 ++++---- src/WaveBinary.F77 | 6 ++++-- 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 -- cgit v1.2.3