var/da/da_vtox_transforms/da_get_spoles.inc

   1 subroutine da_get_spoles(x,              &
   2           ids, ide, jds, jde, &
   3           ims, ime, jms, jme, kms, kme,  &
   4           its, ite, jts, jte, kts, kte  )
   5
   6    !---------------------------------------------------------------------------
   7    ! Purpose: Treatment for Scalar field at Poles
   8    !---------------------------------------------------------------------------
   9
  10    implicit none
  11
  12
  13    integer, intent(in)    :: ids, ide, jds, jde
  14    integer, intent(in)    :: ims, ime, jms, jme, kms, kme
  15    integer, intent(in)    :: its, ite, jts, jte, kts, kte
  16    real,    intent(inout) :: x(ims:ime,jms:jme,kms:kme)
  17
  18    integer                :: k
  19    real                   :: tmpvar,tmps,tmp_s
  20
  21    if (trace_use) call da_trace_entry("da_get_spoles")
  22
  23    tmpvar      = 1.0/real(ide-ids+1)
  24
  25    do k = kts, kte
  26       tmps =0.0  ; tmp_s =0.0
  27       if (jts == jds)  tmp_s = tmpvar*sum(x(its:ite,jts+1,k))
  28
  29       tmps = wrf_dm_sum_real( tmp_s)
  30       if (jts == jds) x(its:ite,jts,k) = tmps
  31
  32       tmps =0.0  ; tmp_s =0.0
  33       if (jte == jde) tmp_s = tmpvar*sum(x(its:ite,jte-1,k))
  34
  35       tmps = wrf_dm_sum_real( tmp_s)
  36       if (jte == jde) x(its:ite,jte,k) = tmps
  37    end do
  38
  39    if (trace_use) call da_trace_exit("da_get_spoles")
  40
  41 end subroutine da_get_spoles
  42
  43