Hi I see that multidimensional arrays are working fine. But I have a weird datatype double precision, dimension(:,:), pointer :: coefs and if I try and print state%instrument%ils(num)%coefs(1,1) I get incorrect values. (I also tried *(state%instrument%ils(num)%coefs(1,1)) but that didn't work either. Are fortran pointers supported? Mathew