watlev.f

Go to the documentation of this file.

      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 'common_all.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,jjjj
         psum = 0.0
         do k = 1,2
            psum=(fio(k,j,1)+fio(k,j,jpart))+psum
            do i=2,noph2
               psum=psum+2.*fio(k,j,i)
            enddo
         enddo
         sum=sum+psum*dmus2(j)*(1.0d0-rfair(j))*ddphi
      enddo
      sumtt=(eo(1)+eo(2))*bmu(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,jjjj
         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**********************************************************************