catmosutil.f File Reference

Go to the source code of this file.

Functions/Subroutines
real *8 function	cnewcos (H, cost, L)
real *8 function	cnewsin (H, cost, L)
real *8 function	cnewh (H, cost, L)
subroutine	cnewcossin (h1, cos1, leng, h2, cos2, sin2)
real *8 function	clenbetween2h (h1, h2, cost)
subroutine	clenbetw2h (h1, h2, cost, leng, icon)

Function/Subroutine Documentation

clenbetw2h()

subroutine clenbetw2h	(	real*8	h1,
		real*8	h2,
		real*8	cost,
		real*8	leng,
		integer	icon
	)

Definition at line 176 of file catmosutil.f.

References d, and d0.

Referenced by csetmintime(), and cthick2len().

       implicit none
       real*8 h1, h2, cost
       real*8 leng ! output
       integer icon ! output

       real*8 sint, costp, sintp
       character*120 text

       sint = sqrt(1.d0 - cost**2)
       sintp = h1*sint /h2
       icon = 0
       if(sintp .le. 1.0d0) then
          costp = sqrt(1.d0 - sintp**2)*sign(1.d0, cost)
          leng = h1* cost - h2 * costp
       else
          if(abs(1.d0-sintp**2) .lt. 1.d-6) then
             costp = 0.
             leng = h1* cost 
          else
             icon = 1
          endif 
       endif

Here is the caller graph for this function:

clenbetween2h()

real*8 function clenbetween2h	(	real*8	h1,
		real*8	h2,
		real*8	cost
	)

Definition at line 138 of file catmosutil.f.

References cerrormsg(), d, and d0.

       implicit none
       real*8 h1, h2, cost

       real*8 sint, costp, sintp
       character*120 text

       sint = sqrt(1.d0 - cost**2)
       sintp = h1*sint /h2
       if(sintp .le. 1.0d0) then

          costp = sqrt(1.d0 - sintp**2)*sign(1.d0, cost)
       else
          if(abs(1.d0-sintp**2) .lt. 1.d-6) then
             costp = 0.
          else
             write(text, *) 'h1, h2, cost=', h1, h2, cost, 
     *          ' sintp=',sintp
             call cerrormsg(text, 1)
             call
     *       cerrormsg('h1,h2,cost invalid at clenbetwee2h', 0)
          endif 
       endif
!
       clenbetween2h = h1* cost - h2 * costp

Here is the call graph for this function:

cnewcos()

real*8 function cnewcos	(	real*8	H,
		real*8	cost,
		real*8	L
	)

Definition at line 70 of file catmosutil.f.

References d0.

! 
      implicit none
      real*8 h, cost, l
!
      real*8 eps/1.d-8/,  tmp
      
      tmp = l/h
      if(tmp .lt. eps) then
          cnewcos = (cost - tmp) / 
     *      ( ( tmp * (1.-cost**2)/2 -cost)*tmp +1.)
      elseif(abs(cost) .ne. 1.d0) then
            cnewcos = (cost - tmp)/ sqrt( (tmp - cost*2)*tmp +1.)
      else
            cnewcos = cost
      endif

cnewcossin()

subroutine cnewcossin	(	real*8	h1,
		real*8	cos1,
		real*8	leng,
		real*8	h2,
		real*8	cos2,
		real*8	sin2
	)

Definition at line 112 of file catmosutil.f.

References d0.

Referenced by clen2thickt(), and ct2lt().

      implicit none
      real*8  h1  !  input.  radial distance from the earth center
      real*8  cos1  ! input.  cos of zenith angle at h1
      real*8 leng    ! input.  length in m along cos1 from h1
      real*8  h2    ! output.  radial  distance from the earth center
      real*8  cos2  ! outpu.  cos of zenith angle at h2.
      real*8  sin2  ! output.  sin of //

      real*8 sin1, cnewh

!        h1 cos1 - h2 cos2 = leng
!        h1 sin1 = h2 sin2
      
      sin1 = sqrt(1.d0 - cos1**2)
      h2 = cnewh(h1, cos1, leng)
      sin2 = h1 * sin1/h2 
      cos2 = (h1*cos1 -leng)/h2

Here is the caller graph for this function:

cnewh()

real*8 function cnewh	(	real*8	H,
		real*8	cost,
		real*8	L
	)

Definition at line 98 of file catmosutil.f.

References d0.

Referenced by clen2thickap().

      implicit none
      real*8 h, cost, l, tmp

      real*8 eps/1.d-8/

      tmp = l/h
      if(tmp .lt. eps) then
         cnewh = h * ( (tmp * (1.d0 -cost**2)/2 - cost )* tmp + 1.d0)
      else
         cnewh =h* sqrt( ( tmp - cost*2)*tmp + 1.d0 )
      endif

Here is the caller graph for this function:

cnewsin()

real*8 function cnewsin	(	real*8	H,
		real*8	cost,
		real*8	L
	)

Definition at line 88 of file catmosutil.f.

References d0.

      implicit none
      real*8 h, cost, l
!
      real*8 cnewh,  sint
      
      sint = sqrt(1.d0 - cost**2)
      cnewsin = h * sint/ cnewh(h, cost, l)