Hi Bertrand, I tried your code on the follwing compilers: NAG 4.2, Intel 8.1 and 9.0, IBM xlf, Sun 6.2, Lahey-Fujitsu 6.2d and GNU g95. All are giving the same result, namely lbound(typ%in), ubound(typ%in) 1 5 lbound(pt_1), ubound(pt_1) 1 5 except for the GNU g95 which gives lbound(typ%in), ubound(typ%in) 1 5 lbound(pt_1), ubound(pt_1) 0 4 Which compilers did you try, and with which results? 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 | +----------------------------------------------------------------------+