On Tue, Jul 27, 2010 at 08:11:39AM -0400, Jeff Squyres wrote:
> On Jul 26, 2010, at 11:06 PM, Hugo Gagnon wrote:
> > 8 integer, parameter :: dp = kind(1.d0)
> > 9 real(kind=dp) :: inside(5), outside(5)
> I'm not a fortran expert -- is kind(1.d0) really double precision? According to http://gcc.gnu.org/onlinedocs/gcc-3.4.6/g77/Kind-Notation.html, kind(2) is double precision (but that's for a different compiler, and I don't quite grok the ".d0" notation).
kind (x) has type default integer and value equal to the kind
type parameter value of x.
p. 161 from Metcalf et al (2007) Fortran 95/2003 explained.
% cat tmp.f90
integer, parameter :: dp = kind(1.d0)
real(kind=dp) :: inside(5), outside(5)
end program z
% g95 -L/usr/local/lib tmp.f90
Kind 8 is (on most arch) 8-byte real, i.e. typically
Room 2.6, Queen's Building
Mech Eng Dept
University Walk, Bristol BS8 1TR, UK
Tel: +44 (0)117 331 5944
Fax: +44 (0)117 929 4423