[FLASH-BUGS] Bug in the visc_explicit.F90 routine

[Tomasz Plewa]

You are right, Salvatore. This will be fixed in the next code release.

Tomek
--
On Wed, Apr 09, 2003 at 12:00:23PM +0200, Salvatore Orlando wrote:
> 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

--