var/da/da_radar/da_radzicevar_pkx.inc

   1   subroutine da_radzicevar_pkx(ice_abc,pxabk_all,pkx)
   2   !----------------------------------------------------------------------
   3   ! purpose: calculate Pkx
   4   !----------------------------------------------------------------------
   5   implicit none
   6   real :: ice_abc(3),pxabk_all(3)
   7   real :: pkx
   8   integer :: i
   9
  10   pkx=0
  11   do i=1,3
  12     if(i<3) then
  13       pkx=pkx+ice_abc(i)*pxabk_all(i)
  14     else
  15       pkx=pkx+2.*ice_abc(i)*pxabk_all(i)
  16     endif
  17   enddo
  18
  19   end subroutine da_radzicevar_pkx