Dear all,
I think I have found a small bug in the
source/hydro/explicit/split/ppm/diffuse/visc_explicit.F90 routine (FLASH
2.2).
In the z-sweep section, see LINES 356 and 357 (see below)
!-----------------------------------------------------------------------------
! z-sweep
!-----------------------------------------------------------------------------
elseif (xyzswp .EQ. sweep_z .AND. ndim .EQ. 3 &
& .and. diffuse_visc) then
! compute dx for this block -- in the z direction
call dBaseGetCoords(iznl, iYcoord, block_no, xl)
<------------ 356
call dBaseGetCoords(iznr, iYcoord, block_no, xr)
<------------ 357
*******************************************************************************
They should be changed with
call dBaseGetCoords(iznl, iZcoord, block_no, xl)
call dBaseGetCoords(iznr, iZcoord, block_no, xr)
Best wishes
Salvatore
-- I N A F Istituto Nazionale di Astrofisica Osservatorio Astronomico di Palermo "G.S. Vaiana" Piazza del Parlamento 1, I-90134 Palermo, Italy Tel. +39-91-233 423, Fax. +39-91-233 444 Web http://www.astropa.unipa.it/~orlando
This archive was generated by hypermail 2b30 : Wed Apr 09 2003 - 05:00:42 CDT