clen2thick.f File Reference

Go to the source code of this file.

Functions/Subroutines
real *8 function	clen2thick (z, cosz, leng)
real *8 function	clen2thickap (z, cosz, leng)
real *8 function	clen2thickex (z, cosz, leng, n)

Function/Subroutine Documentation

clen2thick()

real*8 function clen2thick	(	real*8	z,
		real*8	cosz,
		real*8	leng
	)

Definition at line 27 of file clen2thick.f.

!
!     z: real*8. input.  vertical height in m.
!  cosz: real*8. input.  cos of zenith angle at z.
!   leng: real*8. input.  length along cosz direction.
!   function value.      thickness of air in kg/m2.
!
      implicit none
#include  "Ztrack.h"
#include  "Ztrackp.h"
#include  "Ztrackv.h"
#include  "Zatmos.h"

      real*8 z, cosz, leng, clen2thickex, clen2thickap
      real*8  clen2thickta

      if(usetbl .and. z .lt. htop) then
         clen2thick = clen2thickta(z, leng)
      elseif(exactthick) then
         clen2thick = clen2thickex(z, cosz, leng, 10)
      elseif(abs(cosz) .gt. 0.6) then
         clen2thick = clen2thickap(z, cosz, leng)
      else
         clen2thick = clen2thickex(z, cosz, leng, 10)
      endif

clen2thickap()

real*8 function clen2thickap	(	real*8	z,
		real*8	cosz,
		real*8	leng
	)

Definition at line 80 of file clen2thick.f.

References cnewh().

!
!     z: real*8. input.  vertical height in m.
!  cosz: real*8. input.  cos of zenith angle at z.
!   leng: real*8. input.  length along cosz direction.
!   function value.      thickness of air in kg/m2.
!
      implicit none
!----      include 'Zearth.h'
#include  "Zearth.h"

      real*8 z, cosz, leng
!

      real*8 cosm, rs, re, cvh2thick, cnewcos, ze, t1, t2,
     *  cnewh



      rs = z + eradius
      re = cnewh(rs, cosz, leng)
      ze = re - eradius
      cosm = cnewcos(rs, cosz, leng/2)
      t1 = cvh2thick(z)
      t2 = cvh2thick(ze)
      clen2thickap = (t2 - t1)/cosm

Here is the call graph for this function:

clen2thickex()

real*8 function clen2thickex	(	real*8	z,
		real*8	cosz,
		real*8	leng,
		integer	n
	)

Definition at line 131 of file clen2thick.f.

References catmosrho(), k16pgaussleg(), and true.

!
!     z: real*8. input.  vertical height in m.
!  cosz: real*8. input.  cos of zenith angle at z.
!   leng: real*8. input.  length along cosz direction.
!     n:   integer.input.  how many points for Gauss-Legendre integ.
!   function value.      thickness of air in kg/m2.
!
      implicit none
!----      include 'Zearth.h'
#include  "Zearth.h"

      real*8 z, cosz, leng
      integer n
!
      real*8  seg/1000./, ans, a, b, inte
!
      external catmosrho
      real*8 catmosrho

      common /ccatmosrho/ coss, rs
      real*8 coss, rs
!/////
      real(8),parameter::eps=1.d-5
      real(8):: error
      integer icon
!///////

      inte =  0.
      b  =0.
      coss = cosz
      rs = z + eradius
      do while (.true.)
         a = b
         b = min(leng, a + seg)
!         call kdexpIntF(catmosrho, a, b, eps, ans, error, icon)
         call k16pgaussleg(catmosrho, a, b, n, ans)
         inte = inte + ans
!         write(*, *) 'a =', a, ' b=',b, ' ans=',ans
         if(b .eq. leng) goto 10
      enddo
 10   continue
      clen2thickex = inte

Here is the call graph for this function: