Read data and append it to array, best practice?

[interkosmos]

Fortran 2003 also introduced the intrinsic move_alloc() routine that can be used to enlarge arrays.