da_print_be_stats_h_regional.inc

1 subroutine da_print_be_stats_h_regional(outunit, variable, nk, scale_length)
2 
3    !-----------------------------------------------------------------------
4    ! Purpose: TBD
5    !-----------------------------------------------------------------------
6 
7    implicit none 
8 
9    integer, intent(inout)   :: outunit            ! Output file unit.
10    character*10, intent(in) :: variable           ! Variable name.
11    integer, intent(in)      :: nk                 ! Dimension of vertical index
12    real, intent(in)         :: scale_length(1:nk) ! Correlation scale lengths
13 
14    integer                  :: k                  ! Loop counter.
15 
16    write(unit=stdout,fmt='(3a,i5)') &
17       ' Scale length for variable ', trim(variable), ' in unit ', outunit
18 
19    do k = 1, nk
20      write(unit=outunit,fmt='(i4,1pe15.5)')k, scale_length(k)
21    end do
22 
23    outunit = outunit + 1
24    write(unit=stdout,fmt=*) ' '
25 
26 end subroutine da_print_be_stats_h_regional
27 
28