cstdatmos0-multi-seg.f File Reference

Go to the source code of this file.

Functions/Subroutines
subroutine	cstdatmos0
real *8 function	cvh2den (z)
real *8 function	cvh2thick (z)
real *8 function	cvthick2h (t)
real *8 function	cvh2denp (z)
real *8 function	cvh2den2p (z)

Function/Subroutine Documentation

cstdatmos0()

subroutine cstdatmos0

(

)

Definition at line 38 of file cstdatmos0-multi-seg.f.

References copenf(), cskipcomment(), d, d0, d3, false, fd1, ha, hl, and true.

Referenced by cvh2den(), cvh2den2p(), cvh2denp(), cvh2thick(), and cvthick2h().

      implicit none
#include  "Zmanager.h"
#include  "Zearth.h"
#include  "Zstdatmos.h"

      integer icon,  ios, i
      external blkstdatmos
c----      include 'Zstdatmosf.h'     !  rho=rho(h) function
#include  "Zstdatmosf.h"

      if(first) then
         call copenf(tempdev, "stdatmos2%d", icon)
         if(icon .ne. 0) stop 9999
         call cskipcomment(tempdev, icon)
         if(icon .ne. 0) stop 9999

         nodes = 0
         do while( .true. )
            read(tempdev, *, iostat=ios)
     *           znode(nodes+1), tnode(nodes+1), pnode(nodes+1), 
     *           rho(nodes+1),  alfa(nodes+1), d0(nodes+1),
     *           dsum(nodes+1),  scaleh(nodes+1)
            if(ios .ne. 0) goto 10
            if(nodes .ge. maxnodes) then
               write(*,*) 'numbr of nodes for atmosphere > ', maxnodes
               stop 9999
            endif
            nodes = nodes + 1
!           write(0,*) znode(nodes), alfa(nodes)
         enddo
 10      continue
!        write(*, *) " # of nodal points =", nodes

         do i=2, nodes
            if(alfa(i-1) .ne. 0.) then
               fd1i(i) = fd1(znode(i), alfa(i-1), znode(i-1), 
     *                   scaleh(i-1))
               rhop(i) =  rho(i-1) *(-1.d0 -1.d0/alfa(i-1)) *
     *              alfa(i-1)/scaleh(i-1) 
               pwp(i) =-2.d0-1.d0/alfa(i-1)
            else
               fd0i(i) = fd0(znode(i), znode(i-1), scaleh(i-1) )
            endif
         enddo
         fd3 = ((ha-znode(2))/hl)**(pw+1.d0)
         first = .false.
      endif

Here is the call graph for this function:

Here is the caller graph for this function:

cvh2den()

real*8 function cvh2den

(

real*8

z

)

Definition at line 100 of file cstdatmos0-multi-seg.f.

References cstdatmos0(), cvh2den(), d0, ha, and hl.

!     *********************************   
!         gives density of air at a given vertical height.
!      h: vertical height in m
!      function value:  density in kg/m3.
!
      implicit none
c----      include 'Zearth.h'
#include  "Zearth.h"
c----      include 'Zstdatmos.h'
#include  "Zstdatmos.h"
!
      integer i

      real*8 z

      if(first) call cstdatmos0

!          check if z < 11.2km
      if(z .lt. znode(2)) then
         cvh2den = rho00 * ( (ha - z)/hl )**pw
      else
         do i = 3, nodes
            if(z .lt. znode(i) .or. i .eq. nodes ) then
               if(alfa(i-1) .ne. 0.)then
                  if(z .lt. znode(i-1) ) then
                     write(0,*) ' z=',z, ' znode =', znode(i-1)
                  endif
                  cvh2den = 
     *             rho(i-1) * 
     *             (1.d0 + alfa(i-1)* 
     *             (z-znode(i-1))/scaleh(i-1) )**(-1.0-1.d0/alfa(i-1)) 
                   goto 10
               else
                  cvh2den =
     *            rho(i-1) * exp(- (z-znode(i-1))/scaleh(i-1))         
                  goto 10
               endif
            endif
         enddo
 10      continue
      endif

Here is the call graph for this function:

cvh2den2p()

real*8 function cvh2den2p

(

real*8

z

)

Definition at line 326 of file cstdatmos0-multi-seg.f.

References cstdatmos0(), cvh2den2p(), d0, ha, and hl.

!     *********************************   
!         gives double derivative of the density of
!     air at a given vertical height.
!      h: vertical height in m
!      function value: d (drho/dz)/dz density in kg/m5
!
      implicit none
c----      include 'Zearth.h'
#include  "Zearth.h"
c----      include 'Zstdatmos.h'
#include  "Zstdatmos.h"
!
      integer i

      real*8 z

      if(first) call cstdatmos0

!          check if z < 11.2km
      if(z .lt. znode(2)) then
         cvh2den2p = rho00 * pw*(pw-1.d0)/hl/hl
     *          * ( (ha - z)/hl )**(pw-2.d0)
      else
         do i = 3, nodes
            if(z .lt. znode(i) .or. i .eq. nodes ) then
               if(alfa(i-1) .ne. 0.)then
                  cvh2den2p = 
     *               rhop(i) * pwp(i) *alfa(i-1)/scaleh(i-1)*
     *             (1.d0 + alfa(i-1)* 
     *             (z-znode(i-1))/scaleh(i-1) )**(pwp(i)-1.d0)
                   goto 10
               else
                  cvh2den2p =
     *            rho(i-1) * exp(- (z-znode(i-1))/scaleh(i-1))         
     *            /scaleh(i-1)**2
                  goto 10
               endif
            endif
         enddo
 10      continue
      endif

Here is the call graph for this function:

cvh2denp()

real*8 function cvh2denp

(

real*8

z

)

Definition at line 280 of file cstdatmos0-multi-seg.f.

References cstdatmos0(), cvh2denp(), d0, ha, and hl.

!     *********************************   
!         gives derivative of the density of
!     air at a given vertical height.
!      h: vertical height in m
!      function value: drho/dz density in kg/m4
!
      implicit none
c----      include 'Zearth.h'
#include  "Zearth.h"
c----      include 'Zstdatmos.h'
#include  "Zstdatmos.h"
!
      integer i

      real*8 z

      if(first) call cstdatmos0

!          check if z < 11.2km
      if(z .lt. znode(2)) then
         cvh2denp =- rho00 * pw/hl* ( (ha - z)/hl )**(pw-1.d0)
      else
         do i = 3, nodes
            if(z .lt. znode(i) .or. i .eq. nodes ) then
               if(alfa(i-1) .ne. 0.)then
                  cvh2denp = 
!       *             rho(i-1) *(-1.d0 -1.d0/alfa(i-1)) *
!       *             alfa(i-1)/scaleh(i-1) *     = rhop(i)
     *               rhop(i) *
     *             (1.d0 + alfa(i-1)* 
!     *             (z-znode(i-1))/scaleh(i-1) )**(-2.d0-1.d0/alfa(i-1)) 
     *             (z-znode(i-1))/scaleh(i-1) )**pwp(i)
                   goto 10
               else
                  cvh2denp =
     *            -rho(i-1) * exp(- (z-znode(i-1))/scaleh(i-1))         
     *            /scaleh(i-1)
                  goto 10
               endif
            endif
         enddo
 10      continue
      endif

Here is the call graph for this function:

cvh2thick()

real*8 function cvh2thick

(

real*8

z

)

Definition at line 144 of file cstdatmos0-multi-seg.f.

References cstdatmos0(), cvh2thick(), d0, fd1, ha, and hl.

!         vertical height to thickness of air above that point.
!      z:  input. vertical height in m
!  function value. output.  thickness in kg/m2.
!
!
      implicit none
c----      include 'Zearth.h'
#include  "Zearth.h"
c----      include 'Zstdatmos.h'
#include  "Zstdatmos.h"
!
      integer i

      real*8 z

c----      include 'Zstdatmosf.h'
#include  "Zstdatmosf.h"


      if(first) call cstdatmos0

!          check if z < 11.2km
      if(z .lt. znode(2)) then
           cvh2thick =hl* rho00 /(pw+1.d0)* (
     *         ( (ha-z)/hl)**(pw+1.d0) -
     *         fd3 )
     *      +  dsum(2)
!              where fd3 is
!     *           ((ha-znode(2))/hl)**(pw+1.d0) )
      else
!          The gramage between  given heights, z1 and z2  is by
!
!              d = d0 *(fd(z1) - fd(z2))  where
!         
!            fd(z) = (1+ a(z-z0)/H(z0))**(-1/a)                   (a != 0)
!
!                  =  exp(-(z-z0)/H )                             (a = 0)
         do i = 3, nodes-1
            if(z .lt. znode(i) .or. i .eq. nodes-1 ) then
               if(alfa(i-1) .ne. 0.)then
                   cvh2thick =  d0(i-1) * (
     *              fd1( z, alfa(i-1), znode(i-1), scaleh(i-1) )
     *              - fd1i(i)
     *             )   + dsum(i) 

!                      where fd1i(i) is
!     *              fd1(znode(i), alfa(i-1),
!     *              znode(i-1), scaleh(i-1) )

                   goto 10
               else
                  cvh2thick = d0(i-1) * (
     *             fd0(z, znode(i-1), scaleh(i-1) )
     *             - fd0i(i)               
     *             ) + dsum(i)
!                      where fd0i is
!     *            fd0(znode(i), znode(i-1), scaleh(i-1) )
!     *            )
                  goto 10
               endif
            endif
         enddo
 10      continue
      endif

Here is the call graph for this function:

cvthick2h()

real*8 function cvthick2h

(

real*8

t

)

Definition at line 211 of file cstdatmos0-multi-seg.f.

References cstdatmos0(), cvthick2h(), d0, ha, and hl.

!         vertical thickness to vertical height conversion
!     t:  input. vertical tikness in kg/m2
!  function value. output.  height in m.
!
!
      implicit none
c----      include 'Zearth.h'
#include  "Zearth.h"
c----      include 'Zstdatmos.h'
#include  "Zstdatmos.h"
!
      integer i

      real*8  t,  temp


      if(first) call cstdatmos0

!          check if z < 11.2km
      if(t .gt. dsum(2)) then
!           t =hl* rho00 /(pw+1.d0)* (
!     *         ( (ha-z)/hl)**(pw+1.d0) -
!     *           ((ha-znode(2))/hl)**(pw+1.d0) )
!     *      +  dsum(2)
!           solve above eq.
            temp =( t - dsum(2))/hl/rho00 * (pw+1.d0)
     *          +  fd3
            cvthick2h =ha -  hl* temp**(1.d0/(pw+1.d0))
      else
!          The gramage between  given heights, z1 and z2  is by
!
!              d = d0 *(fd(z1) - fd(z2))  where
!         
!            fd(z) = (1+ a(z-z0)/H(z0))**(-1/a)                   (a != 0)
!
!                  =  exp(-(z-z0)/H )                             (a = 0)
         do i = 3, nodes-1
            if(t .gt. dsum(i) .or. i .eq. nodes-1 ) then
               if(alfa(i-1) .ne. 0.)then
!                   cvh2thick =  d0(i-1) * (
!     *              fd1( z, alfa(i-1), znode(i-1), scaleh(i-1) )
!     *             -  fd1(znode(i), alfa(i-1),
!     *              znode(i-1), scaleh(i-1) )
!     *             )   + dsum(i) 
!                    solve above eq.
                   temp =( t -dsum(i))/d0(i-1)
     *              +  fd1i(i)
                   cvthick2h = (temp**(-alfa(i-1)) -1.d0)*
     *               scaleh(i-1)/alfa(i-1) + znode(i-1)
                   goto 10
               else
!                  cvh2thick = d0(i-1) * (
!     *             fd0(z, znode(i-1), scaleh(i-1) )
!     *            -fd0(znode(i), znode(i-1), scaleh(i-1) )
!     *            )  + dsum(i)
!                    solve above eq.
                   temp = (t -dsum(i)) /d0(i-1) +
     *             + fd0i(i)
!                   temp =  exp(-(z-z0)/H )    
                   cvthick2h =znode(i-1)- log(temp)*scaleh(i-1)
                  goto 10
               endif
            endif
         enddo
 10      continue
      endif

Here is the call graph for this function: