var/da/da_radar/da_radzicevar_pxabk.inc

   1   subroutine da_radzicevar_pxabk(alpha_rxa,alpha_rxb,k,n,pxabk)
   2   !----------------------------------------------------------------------
   3   ! purpose: calculate PxA, PxB, PxC for Pxk
   4   !----------------------------------------------------------------------
   5   implicit none
   6   integer :: n
   7   integer :: k
   8   real :: alpha_rxa(n),alpha_rxb(n)
   9   real :: pxabk
  10   integer :: i,j
  11
  12   pxabk=0
  13   do j=0,n-1
  14     do i=0,n-1
  15       if(i+j==k) then
  16         pxabk=pxabk+alpha_rxa(i+1)*alpha_rxb(j+1)
  17       endif
  18     enddo
  19   enddo
  20
  21   end subroutine da_radzicevar_pxabk