cNRLAtmos.f File Reference

Go to the source code of this file.

Modules
module	nrl_atmos

Functions/Subroutines
real(8) function	cvh2den (vh)
real(8) function	cvh2temp (vh)
real(8) function	cvh2denp (vh)
real(8) function	cvh2den2p (vh)
subroutine	cnrldataread (devno, filename)
subroutine	cnrlheaderread (io)
subroutine	cnrlheaderwrite (io)
subroutine	cnrlheaderw0 (io)
subroutine	cnrllatlongcheck (lat, long)
subroutine	cnrldatawrite (io, filename)
subroutine	cnrlgendata (lat, long, period)
subroutine	cnrlalloc
subroutine	cnrldatamanip
real(8) function	cthick2h (t)
real(8) function	cthick2den (t)
real(8) function	cvh2thick (vh)
real(8) function	cvh2scaleh (vh)

Variables
integer, save	nsize
integer, save	nsizem
logical, save	first =.true.
real(4), save	glat
real(4), save	glong
integer, save	gday1
integer, save	gday2
integer, save	ghour1
integer, save	ghour2
real(8)	drhodh2
real(4), save	heightx =499d3
real(4), save	thickx
real(4), save	a =3.6399d-09
real(4), save	b =63000.0
real(4), save	al
real(4), save	bl
real(8), parameter	heightspace =1000.d3
real(4), dimension(:), allocatable, save	h
real(4), dimension(:), allocatable, save	den
real(4), dimension(:), allocatable, save	thick
real(4), dimension(:), allocatable, save	temp
real(4), dimension(:,:), allocatable, save	h2dencoef
real(4), dimension(:,:), allocatable, save	h2thickcoef
real(4), dimension(:,:), allocatable, save	h2tempcoef
real(4), dimension(:,:), allocatable, save	thick2hcoef
real(4), dimension(:,:), allocatable, save	thick2dencoef

Function/Subroutine Documentation

cnrlalloc()

subroutine cnrlalloc

(

)

Definition at line 512 of file cNRLAtmos.f.

References nrl_atmos::den, nrl_atmos::h, nrl_atmos::h2dencoef, nrl_atmos::h2tempcoef, nrl_atmos::h2thickcoef, nrl_atmos::nsize, nrl_atmos::nsizem, nrl_atmos::temp, nrl_atmos::thick, nrl_atmos::thick2dencoef, and nrl_atmos::thick2hcoef.

Referenced by cnrldataread(), and cnrlgendata().

      use nrl_atmos
      implicit none
      if( allocated( h ) ) return !   !!!!
      allocate( h(nsize) )
      allocate( den(nsize) )
      allocate( temp(nsize) )
!!        allocate( denp(nsize) )
      allocate( thick(nsize) )
      allocate( h2dencoef(nsizem,3) )
      allocate( h2tempcoef(nsizem,3) )
!!        allocate( h2denpCoef(nsizem,3) )
      allocate( h2thickcoef(nsizem,3) )
      allocate( thick2dencoef(nsizem,3) )
      allocate( thick2hcoef(nsizem,3) )

Here is the caller graph for this function:

cnrldatamanip()

subroutine cnrldatamanip

(

)

Definition at line 529 of file cNRLAtmos.f.

References nrl_atmos::a, nrl_atmos::al, nrl_atmos::b, nrl_atmos::bl, nrl_atmos::den, e3, nrl_atmos::h, nrl_atmos::h2dencoef, nrl_atmos::h2tempcoef, nrl_atmos::h2thickcoef, nrl_atmos::heightx, kscsplcoef(), kscsplinteg(), nrl_atmos::nsize, nrl_atmos::nsizem, nrl_atmos::temp, nrl_atmos::thick, nrl_atmos::thick2dencoef, nrl_atmos::thick2hcoef, and nrl_atmos::thickx.

Referenced by cnrlgendata().

      use nrl_atmos
      implicit none
      integer:: i
      
      real(4)::dum

      call kscsplcoef(h, den, nsize, h2dencoef, nsizem)


      do i = nsize/2, nsize
         if( h(i) > heightx-50.0e3 ) exit
      enddo
      b = (h(nsize)-h(i)) /log(den(i)/den(nsize))
      a = den(i)*exp(h(i)/b)
      write(0,*) ' a,b=',a,b

      bl = (h(2)-h(1))/log(den(1)/den(2)) 
      al = den(1)*exp(h(1)/bl)
!         integral of den at each segment
      thick(:) = 0.
      do i =  1, nsizem
         call kscsplinteg(h, den, nsize, h2dencoef, nsizem,  
     *        h(i), h(i+1), thick(i) )
      enddo
           !  thickness above h> h(nsize) 
           !  ~  3.e-3  thickness at h(nsize) kg/m^2.  very rough
      thickx =  a*b*exp(-h(nsize)/b)
      thick(nsize) = thickx
!       make  cummulative 
      do i = nsizem, 1, -1
         thick(i) = thick(i) + thick(i+1)
      enddo
      
      call kscsplcoef(h, thick, nsize, h2thickcoef, nsizem)
      call kscsplcoef(thick, den, nsize, thick2dencoef, nsizem)
      call kscsplcoef(thick, h, nsize, thick2hcoef, nsizem)
      call kscsplcoef(h, temp, nsize, h2tempcoef, nsizem)

Here is the call graph for this function:

Here is the caller graph for this function:

cnrldataread()

subroutine cnrldataread	(	integer, intent(in)	devno,
		character(*), intent(in)	filename
	)

Definition at line 140 of file cNRLAtmos.f.

References cnrlalloc(), cnrlheaderread(), copenf(), cskipcomment(), nrl_atmos::den, nrl_atmos::h, nrl_atmos::nsize, nrl_atmos::nsizem, nrl_atmos::temp, and true.

      use nrl_atmos
      implicit none
      integer,intent(in):: devno ! temporary  logical dev. # of NRL data 
                          ! container
      character(*),intent(in)::filename ! file name of the container
                 ! file should contain height (m) vs density (kg/m3) table
                 ! first  h is -400 m.  step h should be ~ 500 m.  (100 m is
                 ! too small and spline interpolation becomes NG for the 
                 ! 2nd derivative 500 m gives as good as 100 m in the stable
                 ! region. 
      character(len=64):: line  ! 1 line should be <= 64 cha.

      integer:: i, icon
      call copenf(devno, filename, icon)
      if(icon /= 0 ) then
         write(0,*) ' file=',filename
         write(0,*) ' cannot be opened for cNRLdataRead'
         stop
      endif
      call cnrlheaderread(devno)

      call cskipcomment(devno, icon)
!             find # of lines
      if(icon /= 0 ) then
         write(0,*) ' filename=',trim(filename)
         write(0,*) ' has no #------------- line'
         stop
      endif
      nsize=0
      do while ( .true. )
         read(devno,'(a)', end=100)  line
         nsize = nsize + 1
      enddo
 100  continue

      rewind devno

      nsizem = nsize - 1
      write(0,*) ' nsize=',nsize
      call cnrlalloc            ! allocate arrays

      call cskipcomment(devno, icon)
      do i = 1, nsize
         read(devno, *)  h(i), den(i), temp(i)
      enddo
      close(devno)
      write(0,*) ' NRL atmospheric data has been read from'
      write(0,*)  filename

Here is the call graph for this function:

cnrldatawrite()

subroutine cnrldatawrite	(	integer, intent(in)	io,
		character(*), intent(in)	filename
	)

Definition at line 329 of file cNRLAtmos.f.

References cnrlheaderwrite(), copenfw(), nrl_atmos::den, nrl_atmos::h, nrl_atmos::nsize, and nrl_atmos::temp.

      use nrl_atmos
      implicit none
      integer,intent(in):: io  ! write dev #. If this is 6,
             ! filename is not referred.
      character(*),intent(in):: filename  ! filename path

      integer:: i, icon
      
      if( io /=  6) then 
         call copenfw(io, filename, icon)
         if( icon /= 0 ) then
            write(0,*) ' file=',trim(filename)
            write(0,*) 'cannot be opened for cNRLdataWrite'
            stop
         endif
      endif
      call cnrlheaderwrite(io)
      do i = 1, nsize
         write(io, '(i7, 1p, E12.5, 0p, f8.1)') 
     *     int( h(i) ), den(i), temp(i)
      enddo
      if( io /= 6 ) then
         close(io)
      endif

Here is the call graph for this function:

cnrlgendata()

subroutine cnrlgendata	(	real(8), intent(in)	lat,
		real(8), intent(in)	long,
		integer, dimension(4), intent(in)	period
	)

Definition at line 356 of file cNRLAtmos.f.

References cnrlalloc(), cnrldatamanip(), nrl_atmos::den, nrl_atmos::gday1, nrl_atmos::gday2, nrl_atmos::ghour1, nrl_atmos::ghour2, nrl_atmos::glat, nrl_atmos::glong, nrl_atmos::h, nrl_atmos::heightx, nrl_atmos::nsize, nrl_atmos::nsizem, and nrl_atmos::temp.

      use nrl_atmos
!        get the  average of height vs air density table of the atmosphere 
!        at a given place during the given period specified by
!        starting time and ending time.  
      implicit none
      real(8),intent(in):: lat  ! latitute in deg. of the place
      real(8),intent(in):: long ! longitude in deg. of the place
!                  time info below is for the local time of the place
      integer,intent(in):: period(4) !  period(1):  starting day (Jan.1 is 1)  
                                   !                Dec. 31 is 365 (366 for leap
                                   !                year. but at present, year
                                   !                is not considered. so 365 is
                                   !                better?)    
                                   !    period(2):  endign day 1 ~ 365
                                   !    period(3):  starting hours
                                   !                0,24 is midnight. 12 is noon.
                                   !    period(4)   same for   ending hour
                                   !   
!          Let dayi=period(i) (i=1,2), hourj-2=period(j) (j=3,4)
!          If day1 > day2, the period is understood as straddling Dec. to Jan.
!          If hour1 > hour2, the period is understood as straddling midnight
!           
!          The data is generated by taking the average of data for sampled
!          times in the period:
!          a)  Samples days from day1 to day2 (7 day step; day1 and day2 is always
!              included. If the last but one day is close to day2 (< nearday) 
!              such day is skipped. 
!          b)  For each sampled day, get data for sampled hours between 
!          　hour1 and hour2. (4 hour step; hour1 and hour2 are always included.
!            If the last but one hour is close to hour2 (< nearhour) such  hour
!            is skipped.
!          The user can generate data, e.g, 
!          the average during day time of  the winter season (say, from Dec to Feb),
!          or night time of the same period.  Also it is possible to take the 
!          average of  whole days during the summer season etc.
!     As to the height, data is generated  from lowh to heightx, 
!         above heightx extrapolation by exp formula  is used.
!         Also below lowh, extrapolation is made but it must not be large.
 
      real(8),parameter:: lowh=-400d0  ! generate data from this height (m a.s.l)
      real(8),parameter:: step=500d0   ! height step m. this is enough      
      real(8),parameter:: dayinsec = 24*3600
      integer,parameter:: daystep  = 7  !7 day step in whole days
      integer,parameter:: nearday  = 3   !
      integer,parameter:: hourstep = 4   ! 4 hour step in each  day
      integer,parameter:: nearhour = 2
      integer,save::iday, day, tempday
      integer,save::hour, temphour
      real(8):: dalti

      real(4):: totalden(9), tempera(2)

      integer::nsum
!                usec: UT in sec in the day
!                stl:  apparent local soloar time in the day in hours
!               alti: atltitude in km
      real(4),save::usec, stl, alti
      integer::i

      nsize = (heightx-lowh)/step+1   
      nsizem = nsize - 1

      call cnrlalloc  ! allocate arrays

      den(:) = 0.     ! clear rho and temperature area
      temp(:) = 0.

      glat = lat       ! save lat, long
      glong = long

      gday1 = period(1)
      gday2 = period(2)
      ghour1 = period(3)
      ghour2 = period(4)

      nsum = 0
      day=gday1

      if(day > gday2) then
         day = day-366
      endif

      do while ( day <=   gday2 )
         hour = ghour1
         if( hour > ghour2) then
            hour = hour - 24
         endif

         do while ( hour <= ghour2 )
            if( day > 0 ) then
               tempday = day
            elseif( day < 0 ) then
               tempday = day+366
            elseif( day == 0 ) then
               tempday = 1
            endif
            
            if( hour > 0 ) then
               temphour = hour
            else
               temphour = hour + 24
            endif
 !             local apparent solar time in hours
            stl = temphour
             !  STL=UT(sec)/3600+GLONG/15   hours
            usec = (3600*stl - long*240)  ! UT in sec
            iday = tempday      ! day; year may be added; if y2015, 15*1000 may be
                           ! added
            dalti = lowh      ! lowest height
            nsum = nsum + 1
                ! for given place and time, compute 
                ! density at all heights
            do  i=1,  nsize
               alti = dalti/1000.  ! to km
!                         day  UT(sec) km   deg   deg   hours
               call gtd7(iday, usec, alti, glat, glong, stl,
     &              150., 150., 4., 48, totalden, tempera)
               den(i) = den(i) + totalden(6)   ! g/cm3
               temp(i) = temp(i) +  tempera(2) ! K
               dalti = dalti + step
            enddo

            if( hour == ghour2 ) exit

            hour = hour + hourstep
            if( ghour2 - hour <= nearhour) then
               hour = ghour2
            endif
         enddo

         if( day == gday2 ) exit

         day  = day + daystep
         if( gday2 -  day <= nearday ) then
            day = gday2
         endif
      enddo
!            get average
      den(:) = den(:)*1000./nsum  ! kg/m3
      temp(:) = temp(:) / nsum  ! K

!        make height table
      dalti = lowh
      do i=1, nsize
         h(i) = dalti      ! m
         dalti = dalti + step
      enddo
!           generate other data.
      call cnrldatamanip
      write(0,*) 'NRL atmospheric data has been made for'
      write(0,*) 'Latitute=',glat, ' Londitude=',glong, ' deg'
      write(0,*) 'for period=',period(:)

Here is the call graph for this function:

cnrlheaderread()

subroutine cnrlheaderread

(

integer, intent(in)

io

)

Definition at line 191 of file cNRLAtmos.f.

References nrl_atmos::gday1, nrl_atmos::gday2, nrl_atmos::ghour1, nrl_atmos::ghour2, nrl_atmos::glat, and nrl_atmos::glong.

Referenced by cnrldataread().

      use nrl_atmos
      implicit none
      integer,intent(in):: io  ! read dev #

      character(len=10)::term1, term2
      character(len=20)::term3

      read(io,*) term1, term2, term3
      if( term2 /= "NRL") then
         write(0,*)
     *    'NRL data is needed since ATMOSPHERE 3 is specified'
         write(0,*)
     *    'in Zcondc.h '    
         write(0,*)
     *    'For the NRL atmosphere data: 1st line must be like'
         write(0,*) '# NRL  atmos...  but we have'
         write(0,'(a," ", a," ", a)') 
     *     trim(term1), trim(term2), trim(term3)
         stop
      endif
      read(io,*) term1, term2, glat
      if( term2 /= "lat") then
         write(0,*)
     *   'Atmosphere data for NRL: 2nd line must be like'
         write(0,*) '# lat  32.0'
         stop
      endif
      read(io,*) term1, term2, glong
      if( term2 /= "long") then
         write(0,*)
     *   'Atmosphere data for NRL: 3rd line must be like'
         write(0,*) '# long  132.0'
         stop
      endif
      read(io,*) term1, term2, gday1
      if( term2 /= "day1") then
         write(0,*)
     *   'Atmosphere data for NRL: 4th line must be like'
         write(0,*) '# day1  125'
         stop
      endif

      read(io,*) term1, term2, gday2
      if( term2 /= "day2") then
         write(0,*)
     *   'Atmosphere data for NRL: 6th line must be like'
         write(0,*) '# day2  135'
         stop
      endif

      read(io,*) term1, term2, ghour1
      if( term2 /= "hour1") then
         write(0,*)
     *   'Atmosphere data for NRL: 5th line must be like'
         write(0,*) '# hour1  3'
         stop
      endif

      read(io,*) term1, term2, ghour2
      if( term2 /= "hour2") then
         write(0,*)
     *   'Atmosphere data for NRL: 7th line must be like'
         write(0,*) '# hour2  1.0'
         stop
      endif

Here is the caller graph for this function:

cnrlheaderw0()

subroutine cnrlheaderw0

(

integer, intent(in)

io

)

Definition at line 277 of file cNRLAtmos.f.

References nrl_atmos::gday1, nrl_atmos::gday2, nrl_atmos::ghour1, nrl_atmos::ghour2, nrl_atmos::glat, and nrl_atmos::glong.

Referenced by cnrlheaderwrite(), and main().

      use nrl_atmos
      implicit none
      
      integer,intent(in):: io    ! write dev #

      write(io,'(a)') "# NRL atmosphere: # 3"
      write(io,'(a,f10.1)')  "# lat ", glat
      write(io,'(a, f10.1)') "# long ", glong
      write(io,'(a, i3)') "# day1 ", gday1
      write(io,'(a, i3)') "# day2 ", gday2
      write(io,'(a, i3)') "# hour1 ", ghour1
      write(io,'(a, i3)') "# hour2 ", ghour2

Here is the caller graph for this function:

cnrlheaderwrite()

subroutine cnrlheaderwrite

(

integer, intent(in)

io

)

Definition at line 259 of file cNRLAtmos.f.

References cnrlheaderw0().

Referenced by cnrldatawrite().

      use nrl_atmos
      implicit none
      
      integer,intent(in):: io    ! write dev #

      call cnrlheaderw0(io)

      write(io,'(a)') 
     * "# 3 terms of the 1st 7 lines must be as above (except numbers)"
      write(io,'(a)')
     * "# any comment may be put after 3 terms. and lines below"
      write(io,'(a)')
     * "# until #-----------------  line "
      write(io,'(a)') "# h(m) den(kg/m3) T(K)"
      write(io,'(a)') "#-------------------------------"

Here is the call graph for this function:

Here is the caller graph for this function:

cnrllatlongcheck()

subroutine cnrllatlongcheck	(	real(8), intent(in)	lat,
		real(8), intent(in)	long
	)

Definition at line 292 of file cNRLAtmos.f.

References nrl_atmos::glat, and nrl_atmos::glong.

Referenced by ciniatmos().

!        should be called  when AtmosFile /= " "
!     for ATMOSPHERE 3 (NRL atmosphere)
!     if lat, long in AtmosFile differ from
!     this input, stop is made.
      use nrl_atmos
      implicit none
      real(8),intent(in):: lat   ! LatitOfSite in deg.
      real(8),intent(in):: long ! LongitOfSite in deg
! these two should be the same those in the Namlist input.
!     and if differ from glat glong  stop is made.
!     We may replace lat and long by glat glong to keep
!     consisency with Namelist input. But we worry about
!     doing so. (future program update may use LatitOfSite
!     and LongitOfSite before this rouine is called. etc).
!
      logical::warning
      warning =  abs(lat -glat) > 1.0
      warning= warning .or. abs(long -glong) > 2.0
      if( warning ) then
         write(0,*) 'ATMOSPHERE 3 (NRL atmosphere) is being used'
         write(0,*) 'and AtmosFile contains filename in which '
         write(0,*) 'latitude and longitude are gvien.'
         write(0,*) 'They differ from LatitOfSite and LongitOfSite'
         write(0,*) 'more than 1 or 2 deg.;diff must be < 1 or 2 deg'
         write(0,*) '              latitude         longitude '
         write(0,*) 'AtmosFile  ', glat, glong
         write(0,*) 'param      ', sngl(lat), sngl(long)
         write(0,*)
     *        'Give LatitOfSite and LongitOfSite values close to'
         write(0,*) 'those in AtmosFile; Or '
         write(0,*)
     *   'Consider using NRL_period without using AtmosFile' 
         stop
      endif

Here is the caller graph for this function:

cthick2den()

real(8) function cthick2den

(

real(8), intent(in)

t

)

Definition at line 591 of file cNRLAtmos.f.

References nrl_atmos::al, nrl_atmos::b, nrl_atmos::bl, nrl_atmos::den, nrl_atmos::h, kscsplintp(), nrl_atmos::nsize, nrl_atmos::nsizem, nrl_atmos::thick, nrl_atmos::thick2dencoef, and nrl_atmos::thickx.

Referenced by __atmosd2.f__(), __atmosd4.f__(), cthick2den(), and main().

      use nrl_atmos
      implicit none
      real(8),intent(in):: t    ! kg/m2
      
      real(8)::ans              !  kg/m3
      
      real(4):: sans, st
      if( t > thickx ) then
         if( t<= thick(1) ) then
            st = t
            call kscsplintp(thick, den, nsize, thick2dencoef, nsizem,
     *           st, sans)
            ans = sans
         else
            ans= (t-thick(1))/bl + al*exp(-h(1)/bl)
         endif
      else
       !   ans = a*exp( b*log(t/a/b)/b) =a t/(ab) = t/b
         ans = t/b
      endif

Here is the call graph for this function:

Here is the caller graph for this function:

cthick2h()

real(8) function cthick2h

(

real(8), intent(in)

t

)

Definition at line 569 of file cNRLAtmos.f.

References nrl_atmos::a, nrl_atmos::al, nrl_atmos::b, nrl_atmos::bl, nrl_atmos::h, kscsplintp(), nrl_atmos::nsize, nrl_atmos::nsizem, nrl_atmos::thick, nrl_atmos::thick2hcoef, and nrl_atmos::thickx.

      use nrl_atmos
      implicit none
      real(8),intent(in)::t     ! kg/m2
      real(8)::ans
      
      real(4):: sans, st
      
      if( t > thickx ) then
         if( t <= thick(1) ) then
            st = t 
            call kscsplintp(thick, h, nsize, thick2hcoef, nsizem, st, 
     *         sans)
            ans = sans
         else
            ans = -bl* log( (t-thick(1))/al/bl + exp(-h(1)/bl ))
         endif
      else
         ans = -b*log(t/a/b)
      endif

Here is the call graph for this function:

cvh2den()

real(8) function cvh2den

(

real(8), intent(in)

vh

)

Definition at line 39 of file cNRLAtmos.f.

References nrl_atmos::a, nrl_atmos::al, nrl_atmos::b, nrl_atmos::bl, nrl_atmos::den, nrl_atmos::h, nrl_atmos::h2dencoef, nrl_atmos::heightx, kscsplintp(), nrl_atmos::nsize, and nrl_atmos::nsizem.

      use nrl_atmos
      implicit none
      real(8),intent(in)::vh    !  m
      real(8)::ans
      
      real(4):: svh, sans
      
      svh = vh
      if( vh < heightx ) then
         if( vh >= h(1) ) then
            call kscsplintp(h, den, nsize, h2dencoef, nsizem, svh, sans)
            ans = sans
         else
            ans = al*exp(-vh/bl)
         endif
      else
         ans = a*exp(-vh/b)
      endif

Here is the call graph for this function:

cvh2den2p()

real(8) function cvh2den2p

(

real(8), intent(in)

vh

)

Definition at line 129 of file cNRLAtmos.f.

References nrl_atmos::drhodh2.

!                this must be called after cvh2denp
      use nrl_atmos
      implicit none
      real(8),intent(in)::vh    !  m
      real(8)::ans

      ans = drhodh2

cvh2denp()

real(8) function cvh2denp

(

real(8), intent(in)

vh

)

Definition at line 84 of file cNRLAtmos.f.

References nrl_atmos::a, nrl_atmos::al, nrl_atmos::b, nrl_atmos::bl, d3, nrl_atmos::drhodh2, nrl_atmos::h, nrl_atmos::h2dencoef, nrl_atmos::heightx, kscspldif(), nrl_atmos::nsize, and nrl_atmos::nsizem.

!               get drho/dh and its derivative 
      use nrl_atmos
      implicit none
      real(8),intent(in)::vh    !  m
      real(8)::ans
      
      real(4)::svh, sans, dum
      real(8):: rho 
      svh = vh
      if( vh < heightx ) then
         if( vh >= h(1) ) then
            call kscspldif(h, nsize, h2dencoef, nsizem, svh, sans, dum)
            drhodh2 = dum
            ans = sans
            if( vh > 65.0d3) then
!
!           2nd derivative (dum) is not stable at ~ 72 km or so.
!           and h> 200km 
!           This is rather due to original data. 
!!
!           It's very small (|d(drho/dh)/dh| ~ 10^-(3~4)|(drho/dh)|)
!           so we put it 0 over 65 km.  Also if we take small
!           h step (100 m ), result becomes worse than   500 m step
!           so we use 500 m step.  single precision and double 
!           precision does  not change the situaiton. So we use
!           single precision. Also taking log of density cannot
!           improve the situaiton.

!!!!             drhodh2 = 0.    ! at prsent don't do this
            endif
         else
            call kscspldif(h, nsize, h2dencoef, nsizem, h(1), sans, dum)
              ! not so good but don't care 
            ans = -al/bl*exp(-vh/bl) + al/bl*exp(-h(1)/bl) + sans
            drhodh2 = al/bl/bl*(exp(-vh/bl) - exp(-h(1)/bl))
     *             + dum
         endif
      else
         ans = -a/b*exp(-vh/b)
         drhodh2 = -ans/b
      endif

Here is the call graph for this function:

cvh2scaleh()

real(8) function cvh2scaleh

(

real(8), intent(in)

vh

)

Definition at line 635 of file cNRLAtmos.f.

Referenced by __atmosd2.f__(), cvh2scaleh(), and main().

      use nrl_atmos
      implicit none
      real(8),intent(in):: vh   ! m
      real(8)::ans              ! m 

      real(8),external:: cvh2den, cvh2denp
!
      ans = - cvh2den(vh)/ cvh2denp(vh)

Here is the caller graph for this function:

cvh2temp()

real(8) function cvh2temp

(

real(8), intent(in)

vh

)

Definition at line 61 of file cNRLAtmos.f.

References nrl_atmos::h, nrl_atmos::h2tempcoef, nrl_atmos::heightx, kscsplintp(), nrl_atmos::nsize, nrl_atmos::nsizem, and nrl_atmos::temp.

Referenced by __atmosd2.f__().

      use nrl_atmos
      implicit none
      real(8),intent(in):: vh  ! m
      real(8):: ans !  temperature Kelvin

      real(4):: svh, sans

      svh = vh
      if( vh < heightx ) then
         if( vh >= h(1) ) then
            call kscsplintp(h, temp, nsize, h2tempcoef, nsizem,
     *       svh, sans)
            ans = sans
         else
            ans = temp(1)-(vh-h(1))*0.3/100.
         endif
      else
         ans =  temp(nsize)
      endif

Here is the call graph for this function:

Here is the caller graph for this function:

cvh2thick()

real(8) function cvh2thick

(

real(8), intent(in)

vh

)

Definition at line 614 of file cNRLAtmos.f.

References nrl_atmos::a, nrl_atmos::al, nrl_atmos::b, nrl_atmos::bl, nrl_atmos::h, nrl_atmos::h2thickcoef, nrl_atmos::heightx, kscsplintp(), nrl_atmos::nsize, nrl_atmos::nsizem, and nrl_atmos::thick.

      use nrl_atmos
      implicit none
      real(8),intent(in):: vh   ! m
      real(8)::ans              !  kg/m2
      
      real(4):: sans, svh
      if( vh < heightx ) then
         if( vh >= h(1) ) then
            svh = vh
            call kscsplintp(h, thick, nsize, h2thickcoef, nsizem, 
     *           svh, sans)
            ans = sans
         else
            ans = al*bl*(exp(-vh/bl) -exp(-h(1)/bl)) + thick(1)
         endif
      else
         ans =  a*b*exp(-vh/b)
      endif

Here is the call graph for this function: