file = outputfilename, form = 'unformatted', action='write', iostat = istat ) do l1 = 1, size(store_huge, dim=2) write(file_id, IOSTAT=istat) store_huge(:,l1) if( istat /= 0 ) write(*,*) 'write1 l1, istat', l1, istat enddo !write(file_id, IOSTAT=istat) store_...