<HTML> <BODY BGCOLOR=#ccccdd LINK=#0000aa VLINK=#0000ff ALINK=#ff0000 ><BASE TARGET="bottom_target"><PRE>
<A NAME='DA_INNERPROD'><A href='../../html_code/etkf/da_innerprod.inc.html#DA_INNERPROD' TARGET='top_target'><IMG SRC="../../gif/bar_red.gif" border=0></A>
subroutine da_innerprod(mata,matc,ni,nj) 1,2
!-----------------------------------------------------------------------
! Purpose: TBD
!-----------------------------------------------------------------------
implicit none
integer, intent(in) :: ni, nj
real, intent(in) :: mata(ni,nj)
real*8, intent(out) :: matc(nj,nj)
integer :: i1, i2, k ! Loop counters.
if (trace_use) call da_trace_entry
("da_innerprod")
matc = 0.0
do i1=1,nj
do i2=1,nj
do k=1,ni
matc(i1,i2) = matc(i1,i2) + mata(k,i1)*mata(k,i2)
end do
end do
end do
if (trace_use) call da_trace_exit("da_innerprod")
end subroutine da_innerprod