scale_atmos_thermodyn Module Reference

module ATMOSPHERE / Thermodynamics More...

Functions/Subroutines
subroutine, public	atmos_thermodyn_setup
	Setup. More...
subroutine, public	atmos_thermodyn_tempre (temp, pres, Ein, dens, qdry, q, CVq, Rq)
subroutine, public	atmos_thermodyn_tempre2 (temp, pres, dens, pott, qdry, q, CPq, Rq)

Detailed Description

module ATMOSPHERE / Thermodynamics

Description

Thermodynamics module

Author

Team SCALE

History

• 2011-10-24 (T.Seiki) [new] Import from NICAM
• 2012-02-10 (H.Yashiro) [mod] Reconstruction
• 2012-03-23 (H.Yashiro) [mod] Explicit index parameter inclusion
• 2012-12-22 (S.Nishizawa) [mod] Use thermodyn macro set

Function/Subroutine Documentation

atmos_thermodyn_setup()

subroutine, public scale_atmos_thermodyn::atmos_thermodyn_setup

(

)

Setup.

Definition at line 122 of file scale_atmos_sub_thermodyn.F90.

References scale_grid_index::ieb, scale_stdio::io_fid_log, scale_stdio::io_l, scale_grid_index::isb, scale_grid_index::jeb, scale_grid_index::jsb, scale_grid_index::ke, scale_grid_index::ks, and scale_tracer::qa.

Referenced by mod_rm_driver::scalerm(), and mod_rm_prep::scalerm_prep().

    implicit none
    !---------------------------------------------------------------------------

    if( io_l ) write(io_fid_log,*)
    if( io_l ) write(io_fid_log,*) '++++++ Module[THERMODYN] / Categ[ATMOS SHARE] / Origin[SCALElib]'
    if( io_l ) write(io_fid_log,*) '*** No namelists.'

    return

Here is the caller graph for this function:

atmos_thermodyn_tempre()

subroutine, public scale_atmos_thermodyn::atmos_thermodyn_tempre	(	real(rp), dimension(ka,ia,ja), intent(out)	temp,
		real(rp), dimension(ka,ia,ja), intent(out)	pres,
		real(rp), dimension (ka,ia,ja), intent(in)	Ein,
		real(rp), dimension(ka,ia,ja), intent(in)	dens,
		real(rp), dimension(ka,ia,ja), intent(in)	qdry,
		real(rp), dimension (ka,ia,ja,qa), intent(in)	q,
		real(rp), dimension (qa), intent(in)	CVq,
		real(rp), dimension (qa), intent(in)	Rq
	)

Definition at line 951 of file scale_atmos_sub_thermodyn.F90.

References scale_grid_index::ieb, scale_grid_index::isb, scale_grid_index::jeb, scale_grid_index::jsb, scale_grid_index::ke, scale_grid_index::ks, and scale_tracer::qa.

    implicit none

    real(RP), intent(out) :: temp(KA,IA,JA)    ! temperature
    real(RP), intent(out) :: pres(KA,IA,JA)    ! pressure
    real(RP), intent(in)  :: Ein (KA,IA,JA)    ! internal energy
    real(RP), intent(in)  :: dens(KA,IA,JA)    ! density
    real(RP), intent(in)  :: qdry(KA,IA,JA)    ! dry concentration
    real(RP), intent(in)  :: q   (KA,IA,JA,QA) ! water concentration
    real(RP), intent(in)  :: CVq (QA)          ! specific heat
    real(RP), intent(in)  :: Rq  (QA)          ! gas constant

    real(RP) :: cv, Rmoist

    integer  :: i, j, k, iqw
    !---------------------------------------------------------------------------

    !$omp parallel do private(i,j,k,iqw,cv,Rmoist) OMP_SCHEDULE_ collapse(2)
    do j = jsb, jeb
    do i = isb, ieb
    do k = ks, ke
       cv     = qdry(k,i,j) * cvdry
       rmoist = qdry(k,i,j) * rdry
       do iqw =1, qa
          cv     = cv     + q(k,i,j,iqw) * cvq(iqw)
          rmoist = rmoist + q(k,i,j,iqw) * rq(iqw)
       enddo

       temp(k,i,j) = ein(k,i,j) / cv
       pres(k,i,j) = dens(k,i,j) * rmoist * temp(k,i,j)
    enddo
    enddo
    enddo

    return

atmos_thermodyn_tempre2()

subroutine, public scale_atmos_thermodyn::atmos_thermodyn_tempre2	(	real(rp), dimension(ka,ia,ja), intent(out)	temp,
		real(rp), dimension(ka,ia,ja), intent(out)	pres,
		real(rp), dimension(ka,ia,ja), intent(in)	dens,
		real(rp), dimension(ka,ia,ja), intent(in)	pott,
		real(rp), dimension(ka,ia,ja), intent(in)	qdry,
		real(rp), dimension (ka,ia,ja,qa), intent(in)	q,
		real(rp), dimension (qa), intent(in)	CPq,
		real(rp), dimension (qa), intent(in)	Rq
	)

Definition at line 992 of file scale_atmos_sub_thermodyn.F90.

References scale_grid_index::ieb, scale_grid_index::isb, scale_grid_index::jeb, scale_grid_index::jsb, scale_grid_index::ke, scale_grid_index::ks, and scale_tracer::qa.

    implicit none

    real(RP), intent(out) :: temp(KA,IA,JA)    ! temperature
    real(RP), intent(out) :: pres(KA,IA,JA)    ! pressure
    real(RP), intent(in)  :: dens(KA,IA,JA)    ! density
    real(RP), intent(in)  :: pott(KA,IA,JA)    ! potential temperature
    real(RP), intent(in)  :: qdry(KA,IA,JA)    ! dry concentration
    real(RP), intent(in)  :: q   (KA,IA,JA,QA) ! water concentration
    real(RP), intent(in)  :: CPq (QA)          ! specific heat
    real(RP), intent(in)  :: Rq  (QA)          ! gas constant

    real(RP) :: Rmoist, cp

    integer  :: i, j, k, iqw
    !---------------------------------------------------------------------------

    !$omp parallel do private(i,j,k,iqw,cp,Rmoist) OMP_SCHEDULE_ collapse(2)
    do j = jsb, jeb
    do i = isb, ieb
    do k = ks, ke
       cp     = qdry(k,i,j) * cpdry
       rmoist = qdry(k,i,j) * rdry
       do iqw = 1, qa
          cp     = cp     + q(k,i,j,iqw) * cpq(iqw)
          rmoist = rmoist + q(k,i,j,iqw) * rq(iqw)
       enddo
       calc_pre(pres(k,i,j), dens(k,i,j), pott(k,i,j), rmoist, cp, pre00)

       temp(k,i,j) = pres(k,i,j) / ( dens(k,i,j) * rmoist )
    enddo
    enddo
    enddo

    return