<p><b>duda</b> 2010-08-03 11:24:43 -0600 (Tue, 03 Aug 2010)</p><p>BRANCH COMMIT<br>
<br>
Correction to third-order flux operator from Bill.<br>
<br>
M src/core_nhyd_atmos/module_time_integration.F<br>
</p><hr noshade><pre><font color="gray">Modified: branches/atmos_nonhydrostatic/src/core_nhyd_atmos/module_time_integration.F
===================================================================
--- branches/atmos_nonhydrostatic/src/core_nhyd_atmos/module_time_integration.F 2010-08-02 22:04:09 UTC (rev 455)
+++ branches/atmos_nonhydrostatic/src/core_nhyd_atmos/module_time_integration.F 2010-08-03 17:24:43 UTC (rev 456)
@@ -1097,7 +1097,7 @@
flux3(q_im2, q_im1, q_i, q_ip1, ua, coef3) = &
flux4(q_im2, q_im1, q_i, q_ip1, ua) + &
- coef3*ua*((q_ip1 - q_im2)-3.*(q_i-q_im1))/12.0
+ coef3*abs(ua)*((q_ip1 - q_im2)-3.*(q_i-q_im1))/12.0
! coef_3rd_order = 0.
! if (config_scalar_adv_order == 3) coef_3rd_order = 1.0
@@ -1433,7 +1433,7 @@
flux3(q_im2, q_im1, q_i, q_ip1, ua, coef3) = &
flux4(q_im2, q_im1, q_i, q_ip1, ua) + &
- coef3*ua*((q_ip1 - q_im2)-3.*(q_i-q_im1))/12.0
+ coef3*abs(ua)*((q_ip1 - q_im2)-3.*(q_i-q_im1))/12.0
!------
@@ -1834,7 +1834,7 @@
flux3(q_im2, q_im1, q_i, q_ip1, ua, coef3) = &
flux4(q_im2, q_im1, q_i, q_ip1, ua) + &
- coef3*ua*((q_ip1 - q_im2)-3.*(q_i-q_im1))/12.0
+ coef3*abs(ua)*((q_ip1 - q_im2)-3.*(q_i-q_im1))/12.0
!-----------
</font>
</pre>