Re: [FLASH-USERS] A possible problem in output.F90

From: Tomasz Plewa <tomek@flash.uchicago.edu>
Date: Tue Oct 26 2004 - 09:12:41 CDT

Thanks, Claudio. This is a good idea.

Tomek

--
On Tue, Oct 26, 2004 at 12:24:38PM +0100, Claudio Dalla Vecchia wrote:
> 
> Hi, I am running FLASH 2.4 on SUN and I spotted a problem in output.F90.
> Around line 269,
> 
> tpnext = tinitial + tplot + int((time-tinitial)/tplot)*tplot
>                             ^^^^^^^^^^^^^^^^^^^^^^^^^^
> 
> the INT function could give an integer overflow if the argument is too 
> large.  Notice that can happen in other places of the same routine.  
> Using the AINT function should help:
> 
> tpnext = tinitial + tplot + aint((time-tinitial)/tplot)*tplot
> 
> Cheers,
> Claudio
> 
> -- 
> ***********************************************************************
> Claudio Dalla Vecchia        e-mail: claudio.dalla-vecchia@durham.ac.uk
> Institute for Computational     web: http://star-www.dur.ac.uk/~caius/
> Cosmology, South Road           tel: +44 (0)191 33 43787     linux user
> Durham DH1 3LE - UK             fax: +44 (0)191 33 43645        #275369
> -----------------------------------------------------------------------
>   "Politicians hide themeselves away / They only started the war
>  Why should they go out to fight / They leave that all to the poor"
>                                           ("War Pigs" - Black Sabbath)  
> ***********************************************************************
-- 
Tue, 09:12 CDT (14:12 GMT), Oct-26-2004
_______________________________________________________________________________
   Tomasz Plewa                                      www:   flash.uchicago.edu
   Computational Physics and Validation Group        email: tomek@uchicago.edu
   The ASC FLASH Center, The University of Chicago   phone: 773.834.3227
   5640 South Ellis, RI 475, Chicago, IL 60637       fax:   773.834.3230
_______________________________________________________________________________
Received on Tue Oct 26 09:12:53 2004

This archive was generated by hypermail 2.1.8 : Thu Aug 31 2006 - 21:20:48 CDT