da_status_rad.inc

1 subroutine da_status_rad(iv,str,line)
2 
3    implicit none
4 
5    type(ob_type),    intent(in) :: iv
6    character(len=*), intent(in) :: str
7    integer,          intent(in) :: line
8 
9    integer :: i
10    integer :: isum(10)
11    real    :: rsum(40)
12 
13    if (trace_use) call da_trace_entry("da_status_rad")
14 
15    isum(:) = 0
16    rsum(:) = 0.0
17 
18    isum(3) = -1
19    rsum(3) = -1.0
20    rsum(6) = -1.0
21    rsum(7) = -1.0
22 
23    write (unit=stdout,fmt='(A,I6)') str,line
24    do i=1,rtminit_nsensor
25       write (unit=stdout,fmt='(A,I2)') "For instrument ",i 
26       isum(1) = sum(iv%instid(i)%loc_i(:))
27       isum(2) = sum(iv%instid(i)%loc_j(:))
28       ! isum(3) = sum(iv%instid(i)%loc_k(:))
29       isum(4) = sum(iv%instid(i)%scanpos(:))
30       isum(5) = sum(iv%instid(i)%scanline(:))
31       isum(6) = sum(iv%instid(i)%ifgat(:))
32       isum(7) = sum(iv%instid(i)%landsea_mask(:))
33       isum(8) = sum(iv%instid(i)%cloud_flag(:,:))
34       isum(9) = sum(iv%instid(i)%isflg(:))
35 
36       rsum(1)  = sum(iv%instid(i)%loc_dx(:))
37       rsum(2)  = sum(iv%instid(i)%loc_dy(:))
38       ! rsum(3)  = sum(iv%instid(i)%loc_dz(:))
39       rsum(4)  = sum(iv%instid(i)%loc_dxm(:))
40       rsum(5)  = sum(iv%instid(i)%loc_dym(:))
41       ! rsum(6)  = sum(iv%instid(i)%loc_dzm(:))
42       ! rsum(7)  = sum(iv%instid(i)%zk(:,:))
43       rsum(11) = sum(iv%instid(i)%tb_xb(:,:))
44       rsum(12) = sum(iv%instid(i)%satzen(:))
45       rsum(13) = sum(iv%instid(i)%satazi(:))
46       rsum(14) = sum(iv%instid(i)%solzen(:))
47       rsum(15) = sum(iv%instid(i)%solazi(:))
48       rsum(16) = sum(iv%instid(i)%emiss(:,:))
49       rsum(17) = sum(iv%instid(i)%pm(:,:))
50       rsum(18) = sum(iv%instid(i)%tm(:,:))
51       rsum(19) = sum(iv%instid(i)%qm(:,:))
52       rsum(20) = sum(iv%instid(i)%qrn(:,:))
53       rsum(21) = sum(iv%instid(i)%qcw(:,:))
54       rsum(22) = sum(iv%instid(i)%qci(:,:))
55       rsum(23) = sum(iv%instid(i)%qsn(:,:))
56       rsum(24) = sum(iv%instid(i)%qgr(:,:))
57       rsum(25) = sum(iv%instid(i)%ps(:))
58       rsum(26) = sum(iv%instid(i)%ts(:))
59       rsum(27) = sum(iv%instid(i)%t2m(:))
60       rsum(28) = sum(iv%instid(i)%mr2m(:))
61       rsum(29) = sum(iv%instid(i)%u10(:))
62       rsum(30) = sum(iv%instid(i)%v10(:))
63       rsum(31) = sum(iv%instid(i)%clwp(:))
64       rsum(32) = sum(iv%instid(i)%smois(:))
65       rsum(33) = sum(iv%instid(i)%tslb(:))
66       rsum(34) = sum(iv%instid(i)%snowh(:))
67       rsum(35) = sum(iv%instid(i)%elevation(:))
68       rsum(36) = sum(iv%instid(i)%soiltyp(:))
69       rsum(37) = sum(iv%instid(i)%vegtyp(:))
70       rsum(38) = sum(iv%instid(i)%vegfra(:))
71 
72       isum(10) = sum(iv%instid(i)%tb_qc(:,:))
73       rsum(39) = sum(iv%instid(i)%tb_inv(:,:))
74       rsum(40) = sum(iv%instid(i)%tb_error(:,:))
75 
76       write(unit=stdout,fmt='(A,10I7)')         "1:10 ",isum(1:10)
77       write(unit=stdout,fmt='(A,10(1X,E12.6))') "1:5  ",rsum(1:5)
78       write(unit=stdout,fmt='(A,10(1X,E12.6))') "6:10 ",rsum(6:10)
79       write(unit=stdout,fmt='(A,10(1X,E12.6))') "11:15",rsum(11:15)
80       write(unit=stdout,fmt='(A,10(1X,E12.6))') "16:20",rsum(16:20)
81       write(unit=stdout,fmt='(A,10(1X,E12.6))') "21:25",rsum(21:25)
82       write(unit=stdout,fmt='(A,10(1X,E12.6))') "26:30",rsum(26:30)
83       write(unit=stdout,fmt='(A,10(1X,E12.6))') "31:35",rsum(31:35)
84       write(unit=stdout,fmt='(A,10(1X,E12.6))') "36:40",rsum(36:40)
85    end do
86 
87    if (trace_use) call da_trace_exit("da_status_rad")
88 
89 end subroutine da_status_rad
90 
91