NASA Ocean Color

ocssw V2022

       subroutine watlev
 c
 c     compute the water leaving radiance
 c
 c**********************************************************************
 c
 c     include the common blocks
 c
       implicit real*8 (a-h,o-z)
       include 'afrt_rt2.cmn'
 c
 c**********************************************************************
 c     compute the downward diffused flux
       call fluxlvl(fio,fzz,0)
 c     write(*,*)'downward flux above ocean surface',fzz
       sum = 0.0d0
       sumtt=0.0d0
       sumd=0.0d0
       sumu=0.0d0
       do j = 1,(nx-1)
          psum = 0.0
          do k = 1,2
             psum=(fio(k,j,1)+fio(k,j,jpart))+psum
             do i=2,nophi/2
                psum=psum+2.*fio(k,j,i)
             enddo
          enddo
          sum=sum+psum*dcmusq(j)*(1.0d0-rfair(j))*ddphi
       enddo
       sumtt=(eo(1)+eo(2))*cosmu(kkx)*(1.0d0-rfair(kkx))
       sumd=sum+sumtt
 c     compute the upwelling flux and radiance just below the 
 c     ocean surface
       sumu=alw*sumd
       radu=sumu/pi
       albtdr(ksza)=sumtt
       albtdf(ksza)=sum
       albtdt(ksza)=sumd
       albtrf(ksza)=sumu
 c     write(*,150)sum,sumd,alw,sumu,radu
 150   format('sum,sumd,alw,sumu,radu'/1p5e11.3)
 c     compute the radiance just above the ocean surface
       do i=1,(nx-1)
          ii=nmum1-i+1
          do j=1,jpart
             radxi(ii,j)=((1.0d0-rfwat(i))*radu)/(xrw**2+xiw*2)
             radxi(i,j)=0.0d0
          enddo
 c     write(*,175)i,ii,radxi(ii,1)
 175   format('radiance just above the surface',2i4,1x,1pe11.3)
       enddo
 c
       return
       end
 c**********************************************************************

fluxlvl

subroutine fluxlvl(buft, sumg, iflag)

Definition: fluxlvl.f:2

#define pi

Definition: vincenty.c:23

watlev

subroutine watlev

Definition: watlev.f:2

fio

Definition: RsViirs.h:71