NASA Ocean Color

ocssw V2022

 module viirs_ler_luts
  
 !private
  
 public    ::  load_viirs_ler_luts
 !public    ::  unload_viirs_ler_luts
  
 real, dimension(20800)    ::  logi0
 real, dimension(20800)    ::  z1i0
 real, dimension(20800)    ::  z2i0
 real, dimension(20800)    ::  ti0
 real, dimension(260)      ::  sb
 real, dimension(160)      ::  li0r
 real, dimension(160)      ::  z1i0r
 real, dimension(160)      ::  z2i0r
 real, dimension(160)      ::  ti0r
 real, dimension(2)        ::  sbr
  
 real, dimension(10)       ::  xxzlog = (/0.0000000, 0.00977964, 0.0395086, 0.0904221, &
 &                                        0.164818, 0.266515, 0.401776, 0.581261,      &
 &                                        0.824689, 1.17436/)
 !     -- these numbers correspond to satellite zenith angle
 !     -- node points of 0.,16.,30.,40.,48.,54.,58.,60. degrees
 real, dimension(8)       ::  xxlog = (/0.000000,0.0395086,0.143841,0.266515,0.401776, &
 &                                       0.531394,0.635031,0.693147/)
  
 real, dimension(10)       ::  xzlog
 real, dimension(8)        ::  xlog
 real, dimension(4,7)      ::  densol
 real, dimension(4,5)      ::  denscn
 real, dimension(4)        ::  cthet0
 real, dimension(4)        ::  ctheta
 real, dimension(16)       ::  cofs
 integer                   ::  indsol
 integer                   ::  indscn
 integer                   ::  iofset
 real                      ::  p1
 real                      ::  pr
 real, dimension(10)       ::  dum10
  
 common /lpoly/ xzlog, xlog, densol, denscn, cthet0, ctheta, cofs, &
 &              indsol, indscn, iofset, p1, pr, dum10
               
 real :: xlat
 real :: xlong
 real :: sza
 real :: xthet
 real :: xphi
 real :: cphi
 real :: c2phi
 real :: pteran
 real :: xnvalm
 integer ::  isnow
 integer ::  partial
 real :: so2ind
 real :: resn
 real :: sens
 real :: rsens
 real :: ozbst
 real :: ref
 real :: estozn
 real :: ozcld
 real :: pcloud
 real :: prfrac
 real :: clfrac
 real :: rayval
 real :: r412
 real :: r470
 real :: sfref412
 real :: sfref470
 real :: sfref650
 real :: qdif412
 real :: qdif470
 real :: qdif650
 real :: stdv
  
 common /sample/ xlat,xlong,sza,xthet,xphi,cphi,c2phi,pteran,    &
 & xnvalm(6),isnow,so2ind,resn(5),sens(5),                  &
 & rsens,ozbst,ref,estozn,ozcld,pcloud,prfrac,clfrac,partial, &
 & rayval(5),r412,r470,sfref412,sfref470,sfref650,qdif412,       &
 & qdif470,qdif650, stdv
         
 real, dimension(26) :: realbuf
 common /bufout/ realbuf
                  
 common /table/  logi0, z1i0, z2i0, ti0, sb, li0r, z1i0r, z2i0r, ti0r, sbr
  
 contains
  
 subroutine load_viirs_ler_luts(lut_ler_file, status)
   implicit none
   
   character(len=255), intent(in)    ::  lut_ler_file
   
   real        ::  xdenom
   integer     ::  i, j, k
   integer     ::  status
   
   open(1000, file=trim(lut_ler_file), form="unformatted", convert="big_endian", &
   &       status="old", iostat=status)
   if (status /= 0) then
     print *, "ERROR: Failed to open LER LUT: ", status
     return
   end if  
   
   read (1000, iostat=status) logi0, z1i0, z2i0, ti0, sb
   read (1000, iostat=status) li0r, z1i0r, z2i0r, ti0r, sbr
   if (status /= 0) then
     if (status > 0) then
       print *, "ERROR: Failed to read LER LUT file: ", status
       status = -1
       return
     end if
   end if
   
   close(1000)
   
   xzlog(:) = xxzlog(:)
   xlog(:) = xxlog(:)
   
 ! -- calculate values needed in table interpolation
   do j = 1, 7
     do k = j, j + 3
       xdenom = 1.0
       do i = j, j + 3
         if (i .ne. k) xdenom = xdenom * (xzlog(k) - xzlog(i))
       end do
       densol(k - j + 1, j) = xdenom
     end do
   end do
   do j = 1, 5
     do k = j, j + 3
       xdenom = 1.0
       do i = j, j + 3
         if (i .ne. k) xdenom = xdenom * (xlog(k) - xlog(i))
       end do
       denscn(k - j + 1, j) = xdenom
     end do
   end do
       
   return
 end subroutine load_viirs_ler_luts
  
 integer function wrap_total() result(status)
   implicit none
   
   status = 0
   return
   
 end function wrap_total
  
 end module viirs_ler_luts

viirs_ler_luts::sfref650

real sfref650

Definition: viirs_ler_luts.f95:71

viirs_ler_luts::sbr

real, dimension(2) sbr

Definition: viirs_ler_luts.f95:17

viirs_ler_luts::li0r

real, dimension(160) li0r

Definition: viirs_ler_luts.f95:13

viirs_ler_luts::rayval

real rayval

Definition: viirs_ler_luts.f95:66

viirs_ler_luts::densol

real, dimension(4, 7) densol

Definition: viirs_ler_luts.f95:29

viirs_ler_luts::z1i0r

real, dimension(160) z1i0r

Definition: viirs_ler_luts.f95:14

viirs_ler_luts::xxlog

real, dimension(8) xxlog

Definition: viirs_ler_luts.f95:24

viirs_ler_luts::ozbst

real ozbst

Definition: viirs_ler_luts.f95:59

viirs_ler_luts::c2phi

real c2phi

Definition: viirs_ler_luts.f95:50

viirs_ler_luts::stdv

real stdv

Definition: viirs_ler_luts.f95:75

viirs_ler_luts::estozn

real estozn

Definition: viirs_ler_luts.f95:61

viirs_ler_luts::resn

real resn

Definition: viirs_ler_luts.f95:56

table

an array had not been initialized Several spelling and grammar corrections were which is read from the appropriate MCF the above metadata values were hard coded A problem calculating the average background DN for SWIR bands when the moon is in the space view port was corrected The new algorithm used to calculate the average background DN for all reflective bands when the moon is in the space view port is now the same as the algorithm employed by the thermal bands For non SWIR changes in the averages are typically less than Also for non SWIR the black body DNs remain a backup in case the SV DNs are not available For SWIR the changes in computed averages were larger because the old which used the black body suffered from contamination by the micron leak As a consequence of the if SV DNs are not available for the SWIR the EV pixels will not be the granule time is used to identify the appropriate tables within the set given for one LUT table(e.g. m1) to be used for interpolation. The table values will be linearly interpolated using the tables corresponding to the node times bracketing the granule time. If the granule time falls before the time of the first node or after the time of the last node

viirs_ler_luts::ref

real ref

Definition: viirs_ler_luts.f95:60

viirs_ler_luts::r412

real r412

Definition: viirs_ler_luts.f95:67

viirs_ler_luts::prfrac

real prfrac

Definition: viirs_ler_luts.f95:64

viirs_ler_luts::pr

real pr

Definition: viirs_ler_luts.f95:38