diff options
Diffstat (limited to 'src/WaveBinary.F77')
-rw-r--r-- | src/WaveBinary.F77 | 6 |
1 files changed, 4 insertions, 2 deletions
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 |