var/da/da_biascorr_airmass/da_bias_scan.f90

   1   program da_bias_scan
   2
   3   use rad_bias, only : long, bias, jpchan, jpscan, da_read_biasprep, print_bias, &
   4      qc_amsua, qc_amsub, qc_ssmis, jband, boff, bdiv
   5
   6 ! PURPOSE.
   7 ! -------
   8 ! TO CALCULATE SCAN BIAS CORRECTIONS.
   9
  10 ! REFERENCE.
  11 ! ----------
  12 ! The ECMWF Tovs bias Correction, Harris and Kelly..
  13
  14   implicit none
  15
  16 ! Local Variables
  17
  18   type (bias) :: tovs
  19
  20   integer :: nobs(JPCHAN)
  21   integer :: nscan(JPSCAN), nscanch(JPCHAN,JPSCAN)
  22   integer :: nscanchb(JPCHAN,JPSCAN,JBAND)
  23
  24   real (kind=long) :: vmean_dep(JPCHAN), vmean_abs(JPCHAN)
  25   real (kind=long) :: vstd_dep(JPCHAN), vstd_abs(JPCHAN)
  26   real (kind=long) :: vec_dep(JPCHAN), vec_abs(JPCHAN)
  27
  28   real (kind=long) :: vmn(JPCHAN), vstd(JPCHAN)
  29   real (kind=long) :: vmnsc(JPCHAN,JPSCAN), vstdsc(JPCHAN,JPSCAN), vmnrl(JPCHAN,JPSCAN)
  30
  31   real (kind=long) :: vmnscb(JPCHAN,JPSCAN,JBAND), vstdscb(JPCHAN,JPSCAN,JBAND)
  32   real (kind=long) :: vmnrlb(JPCHAN,JPSCAN,JBAND)
  33
  34 ! smoothed values
  35   integer :: nscanchbt(JPCHAN,JPSCAN,JBAND)
  36   real (kind=long) :: vmnscbt(JPCHAN,JPSCAN,JBAND), vstdscbt(JPCHAN,JPSCAN,JBAND)
  37   real (kind=long) :: vmnrlbt(JPCHAN,JPSCAN,JBAND)
  38
  39   ! integer :: ibin
  40   integer :: ierr, iscanm
  41   integer :: j, jv, jscan, js, k , l
  42   ! integer :: path
  43   ! integer :: cloud_flag
  44
  45   integer :: iband
  46   ! integer :: ia
  47   ! integer :: iab
  48   ! integer :: aoff
  49   ! integer :: avb
  50
  51   ! CHARACTER(LEN=2) :: cchan
  52
  53   real (kind=long) :: vmn0, vmn0b(JBAND)
  54   ! real (kind=long) :: VDHL_WINdoW, VDLL_WINdoW
  55
  56   logical :: global, smoothing
  57
  58 ! definition of defaut namelist values
  59 !-----------------------------------------
  60   integer :: sband =1
  61   integer :: eband =18
  62   integer :: kscan =90
  63   integer :: ICOUNT = 0                              ! connter
  64   real (kind=long)   :: fac = 3.0                    ! Sdev factor
  65   namelist /INPUTS/ kscan,fac,global,sband,eband,smoothing
  66
  67 !------------------------------------------------------------------------------
  68 !       1.   SETUP.
  69 !            -----
  70
  71   read (5,INPUTS,end=100)
  72
  73   100 continue
  74
  75   write (6,INPUTS)
  76
  77   OPEN(unit=10,FORM='UNformatTED')       ! Open input file  : Input data from select
  78   OPEN(unit=11,FORM='UNformatTED')       ! Open output file : Scan mean core arrays
  79   OPEN(unit=12,FORM='UNformatTED')       ! Open output file : Scan bias coef
  80
  81
  82 !------------------------------------------------------------------------------
  83 !       2.   READ IN DATA, Q.C., CALC MEANS AND VARIANCES.
  84 !            ---- -- ----- ----- ---- ----- --- ---------
  85
  86   vmean_dep = 0.0                                          ! Clear matrices
  87   vstd_dep  = 0.0
  88   vmean_abs = 0.0                                          ! Clear matrices
  89   vstd_abs  = 0.0
  90   nobs      = 0
  91
  92   read (unit=10,end=265)  tovs%nchan, tovs%npred    ! Read in data
  93   rewind(unit=10)
  94
  95   allocate(tovs%tb(tovs%nchan))
  96   allocate(tovs%omb(tovs%nchan))
  97   allocate(tovs%bias(tovs%nchan))
  98   allocate(tovs%qc_flag(tovs%nchan))
  99   allocate(tovs%cloud_flag(tovs%nchan))
 100   allocate(tovs%pred(tovs%npred))
 101
 102 loop1:&
 103   do
 104       icount=icount+1
 105
 106 !
 107 !  2.1 read in data
 108 !
 109       call da_read_biasprep(tovs,10,ierr)
 110       if (ierr == 0) then      ! not end
 111           continue
 112       elseif (ierr == 1) then  ! end
 113           exit
 114       else                     ! error
 115           stop 'read error in da_scan_bias'
 116       end if
 117
 118       if (icount < 2) call print_bias(tovs)
 119
 120       vec_dep(1:tovs%nchan) = tovs%omb(1:tovs%nchan) ! obs minus guess
 121       vec_abs(1:tovs%nchan) = tovs%tb(1:tovs%nchan)  ! obs
 122
 123 !
 124 !   2.2 QC: extrme values/extrme departure/window channel
 125 !
 126       if (tovs%sensor_id == 3) then
 127            call qc_amsua(tovs)
 128       elseif(tovs%sensor_id == 4) then
 129            call qc_amsub(tovs)
 130       elseif(tovs%sensor_id == 15) then  ! mhs
 131            call qc_amsub(tovs)
 132       elseif(tovs%sensor_id == 10) then  ! ssmis
 133            call qc_ssmis(tovs)
 134       end if
 135
 136     do j=1,tovs%nchan
 137 !-------------------------------
 138 !  2.3 compute the statistics
 139 !-------------------------------
 140       if ( tovs%qc_flag(j) == 1 ) then
 141         nobs(j) = nobs(j) + 1
 142         vmean_dep(j) = vmean_dep(j) + vec_dep(j)
 143         vstd_dep(j) = vstd_dep(j) + vec_dep(j)*vec_dep(j)
 144         vmean_abs(j) = vmean_abs(j) + vec_abs(j)
 145         vstd_abs(j) = vstd_abs(j) + vec_abs(j)*vec_abs(j)
 146       end if
 147     end do
 148
 149   end do loop1
 150
 151   265 continue            ! MEANS AND STANDARD DEVIATIONS
 152
 153   write (6,270) icount,nobs(1:tovs%nchan)
 154   270 format (/1X,'NUMBER OF DATA Total/ACCEPTED'/1X,i10/1X,15I10)
 155
 156 !---------------------------------
 157 !  2.8 mean and std for channels
 158 !---------------------------------
 159   where (nobs(:) /= 0)
 160     vmean_dep(:) = vmean_dep(:)/nobs(:)
 161      vstd_dep(:) = vstd_dep(:)/nobs(:) - vmean_dep(:)**2
 162      vstd_dep(:) = SQRT(MAX(0.0_8,vstd_dep(:)))
 163     vmean_abs(:) = vmean_abs(:)/nobs(:)
 164      vstd_abs(:) = vstd_abs(:)/nobs(:) - vmean_abs(:)**2
 165      vstd_abs(:) = SQRT(MAX(0.0_8,vstd_abs(:)))
 166   end where
 167
 168   vmn(:) = vmean_dep(:)   ! used by outlier check later
 169   vstd(:) = vstd_dep(:)
 170
 171   write (6,288)
 172   288 format (/1X,'FIRST PASS: MEANS AND STANDARD DEVIATIONS')
 173
 174   do j=1,tovs%nchan
 175     jv = j
 176     write (6,289) jv, nobs(j), vmean_abs(j),vstd_abs(j), vmean_dep(j),vstd_dep(j)
 177   end do
 178   289 format (1X,I5,I10,4F15.2)
 179
 180 !-------------------------------------------------------------------------------
 181 !     3.   SECOND PASS THROUGH DATA, EXTRA Q.C.
 182 !          ------ ---- ------- ----- ----- ----
 183
 184   rewind(unit=10)
 185
 186   vmean_dep = 0.0                                          ! Clear matrices
 187   vstd_dep = 0.0
 188   vmean_abs = 0.0                                          ! Clear matrices
 189   vstd_abs = 0.0
 190   vmnsc = 0.0
 191   vstdsc = 0.0
 192   vmnscb = 0.0
 193   vstdscb = 0.0
 194   nobs = 0
 195   nscan = 0
 196   nscanch = 0
 197   nscanchb = 0
 198
 199 loop2:&
 200   do
 201
 202     call da_read_biasprep(tovs,10,ierr)
 203     if (ierr == 0) then      ! not end
 204          continue
 205     elseif (ierr == 1) then  ! end
 206          exit
 207     else                     ! error
 208          stop 'read error in da_scan_bias'
 209     end if
 210
 211       vec_dep(1:tovs%nchan) = tovs%omb(1:tovs%nchan)
 212       vec_abs(1:tovs%nchan) = tovs%tb(1:tovs%nchan)
 213
 214       if (tovs%sensor_id == 3) then
 215            call qc_amsua(tovs)
 216       elseif(tovs%sensor_id == 4) then
 217            call qc_amsub(tovs)
 218       elseif(tovs%sensor_id == 15) then ! mhs
 219            call qc_amsub(tovs)
 220       elseif(tovs%sensor_id == 10) then ! ssmis
 221            call qc_ssmis(tovs)
 222       end if
 223
 224 !---------------------------------------------------------------------------------
 225 ! 3.6 Reject outliers : using fac and stdv calculated in the first pass : loop 1
 226 !---------------------------------------------------------------------------------
 227     do j=1, tovs%nchan
 228       if ( (abs(vec_dep(j)-vmn(j)) > (vstd(j)*FAC)) ) then
 229           tovs%qc_flag(j) = -1
 230       end if
 231     end do
 232
 233 ! Good data : count and use in the ststistics
 234 !----------------------------------------------------
 235
 236     jscan = tovs%scanpos
 237     nscan(jscan) = nscan(jscan) + 1
 238
 239     iband = FLOOR(tovs%lat/BDIV) + BOFF   !! latitude band
 240
 241     do j=1, tovs%nchan
 242       if ( tovs%qc_flag(j) == 1 ) then
 243
 244 !  statistics for channel
 245 !--------------------------
 246         nobs(j) = nobs(j) + 1
 247         vmean_dep(j) = vmean_dep(j) + vec_dep(j)
 248          vstd_dep(j) = vstd_dep(j) + vec_dep(j)*vec_dep(j)
 249         vmean_abs(j) = vmean_abs(j) + vec_abs(j)
 250          vstd_abs(j) = vstd_abs(j) + vec_abs(j)*vec_abs(j)
 251
 252 !  statistics for channel+scan position
 253 !----------------------------------------
 254         nscanch(j,jscan) = nscanch(j,jscan) + 1
 255           vmnsc(j,jscan) = vmnsc(j,jscan) + vec_dep(j)
 256          vstdsc(j,jscan) = vstdsc(j,jscan) + vec_dep(j)*vec_dep(j)
 257
 258 !  statistics for channel + scan position + lat band
 259 !----------------------------------------------------------
 260         nscanchb(j,jscan,iband) =  nscanchb(j,jscan,iband) + 1
 261           vmnscb(j,jscan,iband) = vmnscb(j,jscan,iband) + vec_dep(j)
 262          vstdscb(j,jscan,iband) = vstdscb(j,jscan,iband) + vec_dep(j)*vec_dep(j)
 263
 264       end if
 265     end do
 266
 267   end do loop2
 268
 269 ! Write scan 'core' arrays NOTE: not divided by nobs
 270 !------------------------------------------------------
 271   write (11) nobs(:)         ! statistics for channel
 272   write (11) vmean_dep(:)
 273   write (11) vstd_dep(:)
 274   write (11) vmean_abs(:)
 275   write (11) vstd_abs(:)
 276
 277   write (11) nscanch(:,:)    ! statistics for channel + scan position
 278   write (11) vmnsc(:,:)
 279   write (11) vstdsc(:,:)
 280
 281   write (11) nscanchb(:,:,:) ! statistics for channel + scan position + lat band
 282   write (11) vmnscb(:,:,:)
 283   write (11) vstdscb(:,:,:)
 284
 285 ! Write out copious amounts of output
 286 !----------------------------------------
 287   write (6,270) icount,nobs(1:tovs%nchan)
 288
 289   where (nobs(:) /= 0)     ! for channels
 290     vmean_dep(:) = vmean_dep(:)/nobs(:)
 291     vstd_dep(:) = vstd_dep(:)/nobs(:) - vmean_dep(:)**2
 292     vstd_dep(:) = SQRT(MAX(0.0_8,vstd_dep(:)))
 293     vmean_abs(:) = vmean_abs(:)/nobs(:)
 294     vstd_abs(:) = vstd_abs(:)/nobs(:) - vmean_abs(:)**2
 295     vstd_abs(:) = SQRT(MAX(0.0_8,vstd_abs(:)))
 296   end where
 297
 298   do j=1, tovs%nchan
 299     where (nscanch(j,:) /= 0)  ! for channels + scan
 300       vmnsc(j,:) = vmnsc(j,:)/nscanch(j,:)
 301       vstdsc(j,:) = vstdsc(j,:)/nscanch(j,:) - vmnsc(j,:)**2
 302       vstdsc(j,:) = SQRT(MAX(0.0_8,vstdsc(j,:)))
 303     end where
 304     where (nscanchb(j,:,:) /= 0) ! for channels + scan + band
 305       vmnscb(j,:,:) = vmnscb(j,:,:)/nscanchb(j,:,:)
 306       vstdscb(j,:,:) = vstdscb(j,:,:)/nscanchb(j,:,:) - vmnscb(j,:,:)**2
 307       vstdscb(j,:,:) = SQRT(MAX(0.0_8,vstdscb(j,:,:)))
 308     end where
 309
 310 !  4.1 compute central scan position mean
 311 !--------------------------------------------
 312     if (MOD(KSCAN,2) == 0 ) then   ! even scan number
 313       iscanm = KSCAN/2     !! middle scan position
 314
 315       if( nscanch(j,iscanm) .ne. 0 .and. nscanch(j,iscanm+1) .ne.0 )then
 316          vmn0 = 0.5*(vmnsc(j,iscanm)+vmnsc(j,iscanm+1))
 317          vmn0b(:) = 0.5*(vmnscb(j,iscanm,:)+vmnscb(j,iscanm+1,:))
 318       end if
 319
 320       if( nscanch(j,iscanm) .eq. 0 .and. nscanch(j,iscanm+1) .ne.0 )then
 321          vmn0 = vmnsc(j,iscanm+1)
 322          vmn0b(:) = vmnscb(j,iscanm+1,:)
 323       end if
 324
 325       if( nscanch(j,iscanm) .ne. 0 .and. nscanch(j,iscanm+1) .eq.0 )then
 326          vmn0 = vmnsc(j,iscanm)
 327          vmn0b(:) = vmnscb(j,iscanm,:)
 328       end if
 329
 330     ELSE
 331       iscanm = kscan/2 + 1
 332       vmn0 = vmnsc(j,iscanm)
 333       vmn0b(:) = vmnscb(j,iscanm,:)
 334     end if
 335
 336 !  4.2 compute relative bias
 337 !------------------------------------
 338    do k=1,KSCAN
 339    if( nscanch(j,k) .ne. 0 )then
 340     vmnrl(j,k)=vmnsc(j,k) - vmn0
 341     do l=1, JBAND
 342       vmnrlb(j,K,l) = vmnscb(j,K,l) - vmn0b(l)
 343     end do
 344    end if
 345    end do
 346
 347   end do
 348
 349 ! prinit output
 350
 351 !------------------
 352   write (6,388)
 353   388 format (/1X,'SECOND PASS: MEANS AND STANDARD DEVIATIONS')
 354   do j=1, tovs%nchan
 355     jv = j
 356     write (6,289) jv, nobs(j), vmean_abs(j), vstd_abs(j), vmean_dep(j), vstd_dep(j)
 357   end do
 358
 359 !-------------------
 360   write (6,370) (js,js=1,KSCAN)
 361   370 format (/5X,'NUMBER AT EACH SCAN POSITION'/5X,90I7)
 362
 363   do j=1, tovs%nchan
 364     jv = j
 365     write (6,371) jv, nscanch(j,1:KSCAN)
 366   end do
 367   371 format(1X,I4,90I7)
 368
 369   write (6,391) (js,js=1,KSCAN)
 370   391 format (/1X,'BIASES FOR EACH SCAN ANGLE'/4X,90I7)
 371   do j=1, tovs%nchan
 372     jv = j
 373     write (6,393) jv, (vmnsc(j,js),js=1,KSCAN)
 374   end do
 375   393 format (1X,I3,90F7.2)
 376
 377   write (6,394) (js,js=1,KSCAN)
 378   394 format (/1X,'STD DEV FOR EACH SCAN ANGLE'/4X,90I7)
 379   do j=1, tovs%nchan
 380      jv = j
 381     write (6,396) jv, (vstdsc(j,js),js=1,KSCAN)
 382   end do
 383   396  format (1X,I3,90F7.2)
 384
 385 !----------------------
 386   if (global) then
 387     write (6,'(1x,a17,2i4)') 'SCAN COEFFICIENTS', jband,kscan
 388     do j=1,tovs%nchan
 389       jv = j
 390       do iband=1, jband
 391          write (6,'(2i3,90f7.2)') jv, iband, vmnrlb(j,1:KSCAN,iband)
 392       end do
 393     end do
 394
 395     if (smoothing) then
 396
 397      close(11)
 398
 399      OPEN(11,FORM='UNformatTED')
 400
 401      read (11) nobs(:)
 402      read (11) vmean_dep(:)
 403      read (11) vstd_dep(:)
 404      read (11) vmean_abs(:)
 405      read (11) vstd_abs(:)
 406
 407      read (11) nscanch(:,:)
 408      read (11) vmnsc(:,:)
 409      read (11) vstdsc(:,:)
 410
 411      read (11) nscanchb(:,:,:)
 412      read (11) vmnscb(:,:,:)
 413      read (11) vstdscb(:,:,:)
 414
 415 ! Perform smoothing on banded corrections.
 416 ! for absolute bias/stdev
 417
 418      vmnscbt = 0.0
 419      vstdscbt= 0.0
 420      nscanchbt = 0
 421
 422      do j=1, tovs%nchan
 423       do iband=1, JBAND
 424
 425 !      if (iband <= AVBAND) then
 426 !        aoff = AVBAND - iband + 1
 427 !        avb = AVBAND
 428 !      ELSEif ((iband < 2*AVBAND+1) .AND. (iband > AVBAND))  then
 429 !        aoff = 0
 430 !        avb = 2*AVBAND+1 - iband
 431 !      ELSEif ((iband >= 2*AVBAND+1) .AND. (JBAND-iband+1 >= 2*AVBAND+1)) then
 432 !        aoff = 0
 433 !        avb = 0
 434 !      ELSEif ((JBAND-iband+1 < 2*AVBAND+1) .AND. (JBAND-iband+1 > AVBAND)) then
 435 !        aoff = 0
 436 !        avb = 2*AVBAND+1 - (JBAND-iband+1)
 437 !      ELSEif (JBAND-iband+1 <= AVBAND) then
 438 !        aoff = -(AVBAND - (JBAND-iband+1) + 1)
 439 !        avb = AVBAND
 440 !      end if
 441
 442 !      do ia=-avb, avb
 443 !        iab = iband + ia + aoff
 444 !        vmnscbt(j,1:KSCAN,iband) = vmnscbt(j,1:KSCAN,iband) + vmnscb(j,1:KSCAN,IAb)
 445 !        vstdscbt(j,1:KSCAN,iband) = vstdscbt(j,1:KSCAN,iband) + vstdscb(j,1:KSCAN,IAb)
 446 !        nscanchbt(j,1:KSCAN,iband) = nscanchbt(j,1:KSCAN,iband) + nscanchb(j,1:KSCAN,IAb)
 447 !      end do
 448         vmnscbt(j,1:KSCAN,iband) =  vmnscb(j,1:KSCAN,iband)
 449         vstdscbt(j,1:KSCAN,iband) = vstdscb(j,1:KSCAN,iband)
 450         nscanchbt(j,1:KSCAN,iband) = nscanchb(j,1:KSCAN,iband)
 451       end do
 452     end do
 453
 454     do j=1, tovs%nchan
 455       where (nscanchbt(j,1:KSCAN,:) /= 0)
 456        vmnscbt(j,1:KSCAN,:) = vmnscbt(j,1:KSCAN,:)/nscanchbt(j,1:KSCAN,:)
 457        vstdscbt(j,1:KSCAN,:) = vstdscbt(j,1:KSCAN,:)/nscanchbt(j,1:KSCAN,:) - vmnscbt(j,1:KSCAN,:)**2
 458        vstdscbt(j,1:KSCAN,:) = SQRT(MAX(0.0_8,vstdscbt(j,1:KSCAN,:)))
 459       end where
 460
 461 ! get bias at nadir
 462     if (MOD(KSCAN,2) == 0 ) then  ! even scan number
 463       iscanm = KSCAN/2
 464
 465       if( nscanch(j,iscanm) .ne. 0 .and. nscanch(j,iscanm+1) .ne.0 )then
 466          vmn0b(:) = 0.5*(vmnscbt(j,iscanm,:)+vmnscbt(j,iscanm+1,:))
 467       end if
 468
 469       if( nscanch(j,iscanm) .eq. 0 .and. nscanch(j,iscanm+1) .ne.0 )then
 470          vmn0b(:) = vmnscbt(j,iscanm+1,:)
 471       end if
 472
 473       if( nscanch(j,iscanm) .ne. 0 .and. nscanch(j,iscanm+1) .eq.0 )then
 474          vmn0b(:) = vmnscbt(j,iscanm,:)
 475       end if
 476
 477     ELSE
 478       iscanm = kscan/2 + 1
 479       vmn0b(:) = vmnscbt(j,iscanm,:)
 480     end if
 481
 482 !  4.2 compute relative bias
 483 !------------------------------------
 484    do k=1,KSCAN
 485    if( nscanch(j,k) .ne. 0 )then
 486     do l=1, JBAND
 487       vmnrlb(j,k,l) = vmnscbt(j,k,l) - vmn0b(l)
 488     end do
 489    end if
 490    end do
 491
 492   end do
 493
 494 ! Perform smoothing on banded corrections.
 495 ! for relative bias
 496
 497 !  do j=1, tovs%nchan
 498 !    vmnrlbt(j,1:KSCAN,1) = 0.50*vmnrlb(j,1:KSCAN,1) &
 499 !                         + 0.50*vmnrlb(j,1:KSCAN,2)
 500 !    do iband=2, JBAND-1
 501 !      vmnrlbt(j,1:KSCAN,iband) = 0.25*vmnrlb(j,1:KSCAN,iband-1) &
 502 !                               + 0.50*vmnrlb(j,1:KSCAN,iband) &
 503 !                               + 0.25*vmnrlb(j,1:KSCAN,iband+1)
 504 !    end do
 505 !    vmnrlbt(j,1:KSCAN,JBAND) = 0.50*vmnrlb(j,1:KSCAN,JBAND-1) &
 506 !                             + 0.50*vmnrlb(j,1:KSCAN,JBAND)
 507 !  end do
 508
 509  if ( (eband-sband+1) >= 3 ) then
 510   do j=1, tovs%nchan
 511     vmnrlbt(j,1:KSCAN,sband) = 0.50*vmnrlb(j,1:KSCAN,sband) &
 512                          + 0.50*vmnrlb(j,1:KSCAN,sband+1)
 513     do iband=sband, eband-1
 514       vmnrlbt(j,1:KSCAN,iband) = 0.25*vmnrlb(j,1:KSCAN,iband-1) &
 515                                + 0.50*vmnrlb(j,1:KSCAN,iband) &
 516                                + 0.25*vmnrlb(j,1:KSCAN,iband+1)
 517     end do
 518     vmnrlbt(j,1:KSCAN,eband) = 0.50*vmnrlb(j,1:KSCAN,eband-1) &
 519                              + 0.50*vmnrlb(j,1:KSCAN,eband)
 520   end do
 521  elseif ( (eband-sband+1) == 2 ) then
 522     vmnrlbt(j,1:KSCAN,sband) = 0.50*vmnrlb(j,1:KSCAN,sband) &
 523                          + 0.50*vmnrlb(j,1:KSCAN,sband+1)
 524     vmnrlbt(j,1:KSCAN,eband) = 0.50*vmnrlb(j,1:KSCAN,eband-1) &
 525                              + 0.50*vmnrlb(j,1:KSCAN,eband)
 526  elseif ( (eband-sband+1) == 1 ) then
 527     vmnrlbt(j,1:KSCAN,sband) = vmnrlb(j,1:KSCAN,sband)
 528     vmnrlbt(j,1:KSCAN,eband) = vmnrlb(j,1:KSCAN,eband)
 529  end if
 530
 531 ! output relative bias
 532   do j=1, tovs%nchan
 533      jv = j
 534     do iband=1,JBAND
 535       write (12) jv, vmnrlbt(j,1:KSCAN,iband)
 536     end do
 537   end do
 538
 539 !----------------------
 540     write (6,'(1x,a30,2i4)') 'SMOOTHED SCAN COEFFICIENTS', jband,kscan
 541     do j=1,tovs%nchan
 542       jv = j
 543       do iband=1, jband
 544          write (6,'(2i3,90f7.2)') jv, iband, vmnrlbt(j,1:KSCAN,iband)
 545       end do
 546     end do
 547
 548  end if   ! end if smoothing
 549
 550  else  ! regional
 551 !---------------------
 552   write (6,397) (js,js=1,KSCAN)
 553   397 format (/1X,'RELATIVE BIASES FOR EACH SCAN ANGLE'/4X,90I7)
 554   do j=1, tovs%nchan
 555       jv = j
 556     write (6,399) jv, vmnrl(j,1:KSCAN)
 557   end do
 558   399 format (1X,I3,90F7.2)
 559
 560   do j=1, tovs%nchan
 561      jv = j
 562      write (12) jv, vmnrl(j,1:KSCAN)
 563   end do
 564
 565 !----------------------
 566  end if   ! global
 567
 568    deallocate(tovs%tb)
 569    deallocate(tovs%omb)
 570    deallocate(tovs%bias)
 571    deallocate(tovs%qc_flag)
 572    deallocate(tovs%cloud_flag)
 573    deallocate(tovs%pred)
 574
 575   close(unit=10)
 576   close(unit=11)
 577   close(unit=12)
 578
 579   end program da_bias_scan