var/da/da_gen_be/da_print_be_stats_h_global.inc

   1 subroutine da_print_be_stats_h_global(outunit, variable, k, max_wavenumber, &
   2   total_power)
   3
   4    !-----------------------------------------------------------------------
   5    ! Purpose: TBD
   6    !-----------------------------------------------------------------------
   7
   8    implicit none
   9
  10    integer,      intent(inout) :: outunit        ! Output file unit.
  11    character*10, intent(in)    :: variable       ! Variable name.
  12    integer,      intent(in)    :: k              ! Vertical index.
  13    integer,      intent(in)    :: max_wavenumber ! Smallest scale required (ni/2-1)
  14
  15    ! Total Power spectrum (averaged over time/members)
  16    real,         intent(in) :: total_power(0:max_wavenumber)
  17
  18    integer :: n                     ! Loop counter.
  19    real    :: accum_power, variance ! Accumulated power, variance.
  20
  21    if (trace_use) call da_trace_entry("da_print_be_stats_h_global")
  22
  23    accum_power = 0.0
  24    variance = sum(total_power(0:max_wavenumber))
  25
  26    write(unit=stdout,fmt='(3a,i5,a,i5)')' Power spectra for variable ', trim(variable), &
  27                           ' and level ', k, ' in unit ', outunit
  28
  29    do n = 0, max_wavenumber
  30       accum_power = accum_power + total_power(n)
  31       write(unit=outunit,fmt='(2i4,2f18.5,f8.3)')k, n, total_power(n), accum_power, &
  32                                         accum_power / variance
  33    end do
  34
  35    outunit = outunit + 1
  36    write(unit=stdout,fmt=*) ''
  37
  38    if (trace_use) call da_trace_exit("da_print_be_stats_h_global")
  39
  40 end subroutine da_print_be_stats_h_global
  41
  42