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    if (trace_use) call da_trace_entry("da_print_be_stats_h_regional")
17 
18    write(unit=stdout,fmt='(3a,i5)') &
19       ' Scale length for variable ', trim(variable), ' in unit ', outunit
20 
21    do k = 1, nk
22      write(unit=outunit,fmt='(i4,1pe15.5)')k, scale_length(k)
23    end do
24 
25    outunit = outunit + 1
26    write(unit=stdout,fmt=*) ' '
27 
28    if (trace_use) call da_trace_exit("da_print_be_stats_h_regional")
29 
30 end subroutine da_print_be_stats_h_regional
31 
32