FORCING: Make the eddy orientation random in the 2D case.

Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
This commit is contained in:
Grzegorz Kowal 2017-08-11 18:18:43 -03:00
parent df7cc0f89f
commit e0e327f465

View File

@ -279,7 +279,7 @@ module forcing
end do
ap(:) = amp * ap(:) / tmp / del * dt
#else /* NDIMS == 3 */
ap(:) = (/ 0.0d+00, 0.0d+00, amp / del /) * dt
ap(:) = sign(1.0d+00, randomn()) * (/ 0.0d+00, 0.0d+00, amp / del /) * dt
#endif /* NDIMS == 3 */
! iterate over data blocks and add forcing components