var/da/da_radiance/da_read_biascoef.inc

   1 subroutine da_read_biascoef(string,nchan,nscan,nband,npred,global, &
   2       scanbias,scanbias_b,coef,coef0, vstd_dep)
   3
   4    integer,           intent(in)    :: nchan,nscan,nband,npred
   5    logical,           intent(in)    :: global
   6    character(len=20), intent(in)    :: string
   7    real,              intent(inout) :: scanbias(nchan,nscan)
   8    real,              intent(inout) :: scanbias_b(nchan,nscan,nband)
   9    real,              intent(inout) :: coef(nchan,npred),coef0(nchan)
  10    real,              intent(inout) :: vstd_dep(nchan)
  11
  12    integer :: iunit,j,i, ii,jj, iost
  13    integer :: nobs(nchan)
  14    real    :: vmean_abs(nchan)
  15    real    :: vstd_abs(nchan)
  16    real    :: vmean_dep(nchan)
  17
  18    character(len=80) :: filename
  19
  20    call da_trace_entry("da_read_biascoef")
  21
  22    call da_get_unit(iunit)
  23    filename='biascorr/'//trim(adjustl(string))//'.bcor'
  24    open(unit=iunit,file=filename, form='formatted',iostat = iost, status='old')
  25    if (iost /= 0) then
  26       message(1)="Cannot open radiance biascorr file "//adjustl(filename)
  27       call da_error(__FILE__,__LINE__,message(1:1))
  28    end if
  29
  30    ! read (iunit,'(4i6)') nchan,nscan,nband,npred
  31    read (iunit,'(4i6)')
  32    do i=1, nchan
  33       read (iunit,'(i5,i10,4F8.2)') ii,nobs(i),vmean_abs(i),vstd_abs(i), vmean_dep(i),vstd_dep(i)
  34    end do
  35
  36    do i=1, nchan
  37       read (iunit,'(i5,5F12.5)') ii,(coef(i,j),j=1,npred),coef0(i)
  38    end do
  39
  40    read (iunit,*)
  41    read (iunit,*)
  42    if (global) then   ! global coefs not available now, use regional one
  43       do j=1, nchan
  44          read(iunit,'(i5,90F7.2)') jj, scanbias(j,1:nscan)
  45       end do
  46    else
  47       do j=1, nchan
  48          read(iunit,'(i5,90F7.2)') jj, scanbias(j,1:nscan)
  49       end do
  50    end if
  51
  52    close(iunit)
  53    call da_free_unit(iunit)
  54
  55    call da_trace_exit("da_read_biascoef")
  56
  57 end subroutine da_read_biascoef
  58
  59