On Thu, 24 Nov 2005, John Reid wrote:
> What has happened is the default real is 8 bytes and you have declared a
> 'half-length' variable r4. Adding the statement
>
> write(*,*) r4, huge(1.), huge(r4)
>
> gives the output
>
> INF 1.7976931348623157E+308 3.4028235E+38
OK, I see the point. Giving HUGE an argument with the kind type
for single precision generates the corresponding value. That,
exactly, was what I didn't recognise.
Thanks a lot, John.
Roland
+----------------------------------------------------------------------+
| Roland Schilling Home-Office: +49(89)32929-670 |
| Max-Planck-Institut für Gravitationsphysik |
| -- Albert-Einstein-Institut -- |
| D-85748 Garching E-mail: [log in to unmask] |
| Germany http://www.geo600.uni-hannover.de |
+----------------------------------------------------------------------+
|