anglw.f

Go to the documentation of this file.

      subroutine anglw
c
c     compute angles and their sin/cos in water
c
c**********************************************************************
c
c     include the common blocks
c
      implicit real*8 (a-h,o-z)
      include 'afrt_rt2.cmn'
c
c**********************************************************************
c
c     compute angles of refraction and sin/cos assuming flat ocean
c
      do i=1,nx
         cnuw(i)=dsqrt(1.0d0-cmu(i)**2)/xrw
         rmuw(i)=dasin(cnuw(i))/conv
         cmuw(i)=dsqrt(1.0d0-cnuw(i)**2)
      enddo
c     compute the center of the bin and its sin/cos
      do i=1,(nx-1)
         thew(i)=0.5d0*(rmuw(i)+rmuw(i+1))
         bmuw(i)=dcos(thew(i)*conv)
         bnuw(i)=dsqrt(1.0d0-bmuw(i)**2)
         dmuw(i)=cmuw(i)-cmuw(i+1)
         dmuws2(i)=0.5d0*(cmuw(i)**2-cmuw(i+1)**2)
         dnxy=xrw**2+xiw**2
         yrw=xrw/dnxy   
         yiw=xiw/dnxy
         call frs(the(i),xrw,xiw,rfair(i))
         call frs(thew(i),yrw,yiw,rfwat(i))
      enddo 
      return
      end
c**********************************************************************