da_print_stats_rad.inc

1 subroutine da_print_stats_rad( stats_unit, nchan, rad )
2 !------------------------------------------------------------------------------
3 !  Purpose: print out statistics of omb, oma for radiance data.
4 !
5 !  METHOD:  print out average, rms, minimum, maximum of iv, re
6 !
7 !  HISTORY: 08/2005 - Creation            Zhiquan Liu
8 !
9 !------------------------------------------------------------------------------
10 
11    implicit none
12 
13    integer,           intent(in)    :: stats_unit, nchan
14    type (stats_rad_type), intent(in)    :: rad
15    
16    integer    :: k, n, nmin, nmax
17    integer    :: used_nchan 
18 
19    if (trace_use) call da_trace_entry("da_print_stats_rad")
20 
21     used_nchan = 0
22     do k=1, nchan                  !! loop for channels
23        if(rad%ichan(k)%num > 0) used_nchan = used_nchan + 1
24     end do
25     write(unit=stats_unit, fmt='((a,i5))')  ' used_nchan: ', used_nchan 
26 
27     write(unit=stats_unit, fmt='(6a)') &
28         ' Channel ', &
29         ' num  ', &
30         ' ave  ', &
31         ' rms  ', &
32         ' min  ', &
33         ' max  '
34 
35     do k=1, nchan                  !! loop for channels
36 
37     if(rad%ichan(k)%num > 1) then
38       n    = rad%ichan(k)%num
39       nmin = rad%ichan(k)%minimum%n
40       nmax = rad%ichan(k)%maximum%n
41      
42       write(unit=stats_unit, fmt='((i3,i7,4f8.2))') &
43         k, rad%ichan(k)%num, rad%ichan(k)%ave/real(n), sqrt(rad%ichan(k)%rms/real(n)), &
44         rad%ichan(k)%minimum%value, rad%ichan(k)%maximum%value
45     end if
46     end do
47 
48    if (trace_use) call da_trace_exit("da_print_stats_rad")
49 
50 
51 end subroutine da_print_stats_rad
52