Here is a little MATLAB code and an accompanying fortran snippet to help with the translation:
First the MATLAB code:
F2=[1:13*76*2];
fid=fopen('test1','w');
fprintf(fid,'%16.8f %16.8f %16.8f %16.8f \n',F2);
fclose(fid)
f2=reshape(F2,[13,76,2]);
f2(1:13,1,1)
f2(1,1:10,1)...