Fortran 2003 has square bracket syntax for array concatenation, Intel fortran compiler supports it too. I wrote a simple code here for matrix concatenation:
prog
Fill it using 1-D arrays then reshape your mat3.