Be aware of Fortran's array indexing:
value = array(I,J)
The next location in array's memory is array(I+1,J).
IOW the left most index is stride-1, the second (from left) index is stride-size(left most index).
This indexing order is reverse of what you find in C/C++.