var/da/da_par_util/da_alloc_and_copy_be_arrays.inc

   1 subroutine da_alloc_and_copy_be_arrays (vg , v, jts,jte, kts,kte)
   2
   3    !-----------------------------------------------------------------------
   4    ! Purpose: TBD
   5    !-----------------------------------------------------------------------
   6
   7    implicit none
   8
   9    type (be_subtype), intent(inout) :: vg       ! Global backgrd error struct.
  10    type (be_subtype), intent(inout) :: v        ! Local backgrd error struct.
  11    integer, intent(in)              :: jts,jte  ! Tile dimension.
  12    integer, intent(in)              :: kts,kte  ! Tile dimension.
  13
  14    ! Allocate local-grid structure.
  15    v % mz = vg % mz
  16    if (v % mz > 0) then
  17       allocate  (v % val(jts:jte,1:v % mz))
  18       allocate  (v % evec(jts:jte,kts:kte,1:v % mz))
  19       allocate  (v % rf_alpha(1:v % mz))
  20    end if
  21
  22    ! Make local copies from global arrays.
  23    if (v % mz > 0) then
  24        v % val(jts:jte,1:v % mz) = vg % val(jts:jte,1:v % mz)
  25        v % evec(jts:jte,kts:kte,1:v % mz) = vg % evec(jts:jte,kts:kte,1:v % mz)
  26        v % rf_alpha(1:v % mz) = vg % rf_alpha(1:v % mz)
  27    end if
  28
  29    ! Deallocate global arrays.
  30    if (v % mz > 0) then
  31       deallocate  (vg % val )
  32       deallocate  (vg % evec)
  33       deallocate  (vg % rf_alpha)
  34    end if
  35
  36 end subroutine da_alloc_and_copy_be_arrays
  37
  38