scale_atmos_diagnostic Module Reference

module atmosphere / diagnostic More...

Functions/Subroutines
subroutine, public	atmos_diagnostic_get_therm_rhot (KA, KS, KE, IA, IS, IE, JA, JS, JE, DENS, RHOT, Rtot, CVtot, CPtot, POTT, TEMP, PRES, EXNER)
	ATMOS_DIAGNOSTIC_get_therm_rhot potential temperature, temperature, pressure. More...
subroutine, public	atmos_diagnostic_get_therm_rhoe (KA, KS, KE, IA, IS, IE, JA, JS, JE, DENS, RHOE, Rtot, CPtot, CVtot, TEMP, POTT, PRES, EXNER)
	ATMOS_DIAGNOSTIC_get_therm_rhoe potential temperature, temperature, pressure. More...
subroutine, public	atmos_diagnostic_get_phyd (KA, KS, KE, IA, IS, IE, JA, JS, JE, DENS, PRES, CZ, FZ, PHYD, PHYDH)
	ATMOS_DIAGNOSTIC_get_phyd hydrostatic pressure. More...
subroutine, public	atmos_diagnostic_get_n2 (KA, KS, KE, IA, IS, IE, JA, JS, JE, POTT, Rtot, CZ, N2)
	ATMOS_DIAGNOSTIC_get_n2 N^2. More...
subroutine, public	atmos_diagnostic_get_potv (KA, KS, KE, IA, IS, IE, JA, JS, JE, POTT, Rtot, POTV)
	ATMOS_DIAGNOSTIC_get_potv virtual potential temperature. More...
subroutine, public	atmos_diagnostic_get_teml (KA, KS, KE, IA, IS, IE, JA, JS, JE, TEMP, LHV, LHS, QC, QI, CPtot, TEML)
	ATMOS_DIAGNOSTIC_get_teml liqued water temperature. More...

Detailed Description

module atmosphere / diagnostic

Description

Calculate diagnostic variables

Author

Team SCALE

Function/Subroutine Documentation

atmos_diagnostic_get_therm_rhot()

subroutine, public scale_atmos_diagnostic::atmos_diagnostic_get_therm_rhot	(	integer, intent(in)	KA,
		integer, intent(in)	KS,
		integer, intent(in)	KE,
		integer, intent(in)	IA,
		integer, intent(in)	IS,
		integer, intent(in)	IE,
		integer, intent(in)	JA,
		integer, intent(in)	JS,
		integer, intent(in)	JE,
		real(rp), dimension (ka,ia,ja), intent(in)	DENS,
		real(rp), dimension (ka,ia,ja), intent(in)	RHOT,
		real(rp), dimension (ka,ia,ja), intent(in)	Rtot,
		real(rp), dimension(ka,ia,ja), intent(in)	CVtot,
		real(rp), dimension(ka,ia,ja), intent(in)	CPtot,
		real(rp), dimension (ka,ia,ja), intent(out)	POTT,
		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(out)	EXNER
	)

ATMOS_DIAGNOSTIC_get_therm_rhot potential temperature, temperature, pressure.

Definition at line 59 of file scale_atmos_diagnostic.F90.

Referenced by mod_atmos_vars::atmos_vars_calc_diagnostics().

    use scale_atmos_thermodyn, only: &
       thermodyn_rhot2temp_pres => atmos_thermodyn_rhot2temp_pres
    integer,  intent(in)  :: ka, ks, ke
    integer,  intent(in)  :: ia, is, ie
    integer,  intent(in)  :: ja, js, je

    real(RP), intent(in)  :: dens (ka,ia,ja)
    real(RP), intent(in)  :: rhot (ka,ia,ja)
    real(RP), intent(in)  :: rtot (ka,ia,ja)
    real(RP), intent(in)  :: cvtot(ka,ia,ja)
    real(RP), intent(in)  :: cptot(ka,ia,ja)

    real(RP), intent(out) :: pott (ka,ia,ja)
    real(RP), intent(out) :: temp (ka,ia,ja)
    real(RP), intent(out) :: pres (ka,ia,ja)
    real(RP), intent(out) :: exner(ka,ia,ja)

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

    call thermodyn_rhot2temp_pres( ka, ks, ke, ia, is, ie, ja, js, je, &
                                   dens(:,:,:), rhot(:,:,:),                & ! (in)
                                   rtot(:,:,:), cvtot(:,:,:), cptot(:,:,:), & ! (in)
                                   temp(:,:,:), pres(:,:,:)                 ) ! (out)


!OCL XFILL
    !$omp parallel do default(none) OMP_SCHEDULE_ collapse(2) &
    !$omp private(i,j,k) &
    !$omp shared(POTT,EXNER,RHOT,DENS,TEMP) &
    !$omp shared(KS,KE,IS,IE,JS,JE)
    do j = js, je
    do i = is, ie
    do k = ks, ke
       pott(k,i,j) = rhot(k,i,j) / dens(k,i,j)
       exner(k,i,j) = temp(k,i,j) / pott(k,i,j)
    enddo
    enddo
    enddo

    return

Here is the caller graph for this function:

atmos_diagnostic_get_therm_rhoe()

subroutine, public scale_atmos_diagnostic::atmos_diagnostic_get_therm_rhoe	(	integer, intent(in)	KA,
		integer, intent(in)	KS,
		integer, intent(in)	KE,
		integer, intent(in)	IA,
		integer, intent(in)	IS,
		integer, intent(in)	IE,
		integer, intent(in)	JA,
		integer, intent(in)	JS,
		integer, intent(in)	JE,
		real(rp), dimension (ka,ia,ja), intent(in)	DENS,
		real(rp), dimension (ka,ia,ja), intent(in)	RHOE,
		real(rp), dimension (ka,ia,ja), intent(in)	Rtot,
		real(rp), dimension(ka,ia,ja), intent(in)	CPtot,
		real(rp), dimension(ka,ia,ja), intent(in)	CVtot,
		real(rp), dimension (ka,ia,ja), intent(out)	TEMP,
		real(rp), dimension (ka,ia,ja), intent(out)	POTT,
		real(rp), dimension (ka,ia,ja), intent(out)	PRES,
		real(rp), dimension(ka,ia,ja), intent(out)	EXNER
	)

ATMOS_DIAGNOSTIC_get_therm_rhoe potential temperature, temperature, pressure.

Definition at line 110 of file scale_atmos_diagnostic.F90.

References scale_const::const_pre00.

    use scale_const, only: &
       pre00 => const_pre00
    implicit none
    integer, intent(in) :: ka, ks, ke
    integer, intent(in) :: ia, is, ie
    integer, intent(in) :: ja, js, je

    real(RP), intent(in) :: dens (ka,ia,ja)
    real(RP), intent(in) :: rhoe (ka,ia,ja)
    real(RP), intent(in) :: rtot (ka,ia,ja)
    real(RP), intent(in) :: cptot(ka,ia,ja)
    real(RP), intent(in) :: cvtot(ka,ia,ja)

    real(RP), intent(out) :: pott (ka,ia,ja)
    real(RP), intent(out) :: temp (ka,ia,ja)
    real(RP), intent(out) :: pres (ka,ia,ja)
    real(RP), intent(out) :: exner(ka,ia,ja)

    integer :: i, j, k

    !$omp parallel do
    do j = js, je
    do i = is, ie
    do k = ks, ke
       temp(k,i,j) = rhoe(k,i,j) / ( cvtot(k,i,j) * dens(k,i,j) )
       pres(k,i,j) = dens(k,i,j) * rtot(k,i,j) * temp(k,i,j)
       exner(k,i,j) = ( pres(k,i,j) / pre00 )**( rtot(k,i,j) / cptot(k,i,j) )
       pott(k,i,j) = temp(k,i,j) / exner(k,i,j)
    end do
    end do
    end do

    return

atmos_diagnostic_get_phyd()

subroutine, public scale_atmos_diagnostic::atmos_diagnostic_get_phyd	(	integer, intent(in)	KA,
		integer, intent(in)	KS,
		integer, intent(in)	KE,
		integer, intent(in)	IA,
		integer, intent(in)	IS,
		integer, intent(in)	IE,
		integer, intent(in)	JA,
		integer, intent(in)	JS,
		integer, intent(in)	JE,
		real(rp), dimension ( ka,ia,ja), intent(in)	DENS,
		real(rp), dimension ( ka,ia,ja), intent(in)	PRES,
		real(rp), dimension ( ka,ia,ja), intent(in)	CZ,
		real(rp), dimension (0:ka,ia,ja), intent(in)	FZ,
		real(rp), dimension ( ka,ia,ja), intent(out)	PHYD,
		real(rp), dimension(0:ka,ia,ja), intent(out)	PHYDH
	)

ATMOS_DIAGNOSTIC_get_phyd hydrostatic pressure.

Definition at line 157 of file scale_atmos_diagnostic.F90.

References scale_const::const_grav.

Referenced by mod_atmos_vars::atmos_vars_calc_diagnostics().

    use scale_const, only: &
       grav => const_grav
    implicit none

    integer,  intent(in)  :: ka, ks, ke
    integer,  intent(in)  :: ia, is, ie
    integer,  intent(in)  :: ja, js, je
    real(RP), intent(in)  :: dens (  ka,ia,ja)
    real(RP), intent(in)  :: pres (  ka,ia,ja)
    real(RP), intent(in)  :: cz   (  ka,ia,ja)
    real(RP), intent(in)  :: fz   (0:ka,ia,ja)
    real(RP), intent(out) :: phyd (  ka,ia,ja)
    real(RP), intent(out) :: phydh(0:ka,ia,ja)

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

    !$omp parallel do default(none) OMP_SCHEDULE_ &
    !$omp private(i,j,k) &
    !$omp shared(PHYD,PHYDH,DENS,PRES,CZ,FZ,GRAV) &
    !$omp shared(KS,KE,IS,IE,JS,JE)
    do j = js, je
    do i = is, ie
       phydh(ke,i,j) = pres(ke,i,j) - dens(ke,i,j) * grav * ( fz(ke,i,j) - cz(ke,i,j) )
       do k = ke, ks, -1
          phydh(k-1,i,j) = phydh(k,i,j) + dens(k,i,j) * grav * ( fz(k,i,j) - fz(k-1,i,j) )
!          PHYD (k  ,i,j) = 0.5_RP * ( PHYDH(k,i,j) + PHYDH(k-1,i,j) )
          phyd(k  ,i,j) = sqrt(phydh(k,i,j)) * sqrt(phydh(k-1,i,j))
       end do
    enddo
    enddo

    return

Here is the caller graph for this function:

atmos_diagnostic_get_n2()

subroutine, public scale_atmos_diagnostic::atmos_diagnostic_get_n2	(	integer, intent(in)	KA,
		integer, intent(in)	KS,
		integer, intent(in)	KE,
		integer, intent(in)	IA,
		integer, intent(in)	IS,
		integer, intent(in)	IE,
		integer, intent(in)	JA,
		integer, intent(in)	JS,
		integer, intent(in)	JE,
		real(rp), dimension(ka,ia,ja), intent(in)	POTT,
		real(rp), dimension(ka,ia,ja), intent(in)	Rtot,
		real(rp), dimension(ka,ia,ja), intent(in)	CZ,
		real(rp), dimension (ka,ia,ja), intent(out)	N2
	)

ATMOS_DIAGNOSTIC_get_n2 N^2.

Definition at line 204 of file scale_atmos_diagnostic.F90.

References scale_const::const_grav.

Referenced by mod_atmos_vars::atmos_vars_get_diagnostic_3d().

    use scale_const, only: &
       grav => const_grav
    integer, intent(in) :: ka, ks, ke
    integer, intent(in) :: ia, is, ie
    integer, intent(in) :: ja, js, je

    real(RP), intent(in)  :: pott(ka,ia,ja)
    real(RP), intent(in)  :: rtot(ka,ia,ja)

    real(RP), intent(in)  :: cz(ka,ia,ja)

    real(RP), intent(out) :: n2  (ka,ia,ja)

    real(RP) :: rpt(ka)

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

    !$omp parallel do default(none) OMP_SCHEDULE_ collapse(2) &
    !$omp private(i,j,k) &
    !$omp private(RPT) &
    !$omp shared(N2,POTT,Rtot,CZ,GRAV) &
    !$omp shared(KS,KE,IS,IE,JS,JE)
    do j = js, je
    do i = is, ie
       do k = ks, ke
          rpt(k) = rtot(k,i,j) * pott(k,i,j)
       end do

       n2(ks,i,j) = grav * ( rpt(ks+1) - rpt(ks) ) / ( ( cz(ks+1,i,j) - cz(ks,i,j) ) * rpt(ks) )
       do k = ks+1,ke-1
          n2(k,i,j) = grav * ( rpt(k+1) - rpt(k-1) ) / ( ( cz(k+1,i,j) - cz(k-1,i,j) ) * rpt(k ) )
       end do
       n2(ke,i,j) = grav * ( rpt(ke) - rpt(ke-1) ) / ( ( cz(ke,i,j) - cz(ke-1,i,j) ) * rpt(ke) )
    end do
    end do

    return

Here is the caller graph for this function:

atmos_diagnostic_get_potv()

subroutine, public scale_atmos_diagnostic::atmos_diagnostic_get_potv	(	integer, intent(in)	KA,
		integer, intent(in)	KS,
		integer, intent(in)	KE,
		integer, intent(in)	IA,
		integer, intent(in)	IS,
		integer, intent(in)	IE,
		integer, intent(in)	JA,
		integer, intent(in)	JS,
		integer, intent(in)	JE,
		real(rp), dimension(ka,ia,ja), intent(in)	POTT,
		real(rp), dimension(ka,ia,ja), intent(in)	Rtot,
		real(rp), dimension(ka,ia,ja), intent(out)	POTV
	)

ATMOS_DIAGNOSTIC_get_potv virtual potential temperature.

Definition at line 255 of file scale_atmos_diagnostic.F90.

References scale_const::const_rdry.

Referenced by mod_atmos_vars::atmos_vars_get_diagnostic_3d().

    use scale_const, only: &
       rdry => const_rdry
    integer, intent(in) :: ka, ks, ke
    integer, intent(in) :: ia, is, ie
    integer, intent(in) :: ja, js, je

    real(RP), intent(in)  :: pott(ka,ia,ja)
    real(RP), intent(in)  :: rtot(ka,ia,ja)

    real(RP), intent(out) :: potv(ka,ia,ja)

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

!OCL XFILL
    !$omp parallel do default(none) OMP_SCHEDULE_ collapse(2) &
    !$omp private(i,j,k) &
    !$omp shared(POTV,POTT,Rtot,Rdry) &
    !$omp shared(KS,KE,IS,IE,JS,JE)
    do j = js, je
    do i = is, ie
    do k = ks, ke
       potv(k,i,j) = pott(k,i,j) * rtot(k,i,j) / rdry
    end do
    end do
    end do

    return

Here is the caller graph for this function:

atmos_diagnostic_get_teml()

subroutine, public scale_atmos_diagnostic::atmos_diagnostic_get_teml	(	integer, intent(in)	KA,
		integer, intent(in)	KS,
		integer, intent(in)	KE,
		integer, intent(in)	IA,
		integer, intent(in)	IS,
		integer, intent(in)	IE,
		integer, intent(in)	JA,
		integer, intent(in)	JS,
		integer, intent(in)	JE,
		real(rp), dimension (ka,ia,ja), intent(in)	TEMP,
		real(rp), dimension (ka,ia,ja), intent(in)	LHV,
		real(rp), dimension (ka,ia,ja), intent(in)	LHS,
		real(rp), dimension (ka,ia,ja), intent(in)	QC,
		real(rp), dimension (ka,ia,ja), intent(in)	QI,
		real(rp), dimension(ka,ia,ja), intent(in)	CPtot,
		real(rp), dimension(ka,ia,ja), intent(out)	TEML
	)

ATMOS_DIAGNOSTIC_get_teml liqued water temperature.

Definition at line 298 of file scale_atmos_diagnostic.F90.

Referenced by mod_atmos_vars::atmos_vars_get_diagnostic_3d().

    integer, intent(in) :: ka, ks, ke
    integer, intent(in) :: ia, is, ie
    integer, intent(in) :: ja, js, je

    real(RP), intent(in)  :: temp (ka,ia,ja)
    real(RP), intent(in)  :: lhv  (ka,ia,ja)
    real(RP), intent(in)  :: lhs  (ka,ia,ja)
    real(RP), intent(in)  :: qc   (ka,ia,ja)
    real(RP), intent(in)  :: qi   (ka,ia,ja)
    real(RP), intent(in)  :: cptot(ka,ia,ja)

    real(RP), intent(out) :: teml(ka,ia,ja)

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

!OCL XFILL
    !$omp parallel do default(none) OMP_SCHEDULE_ collapse(2) &
    !$omp private(i,j,k) &
    !$omp shared(TEML,TEMP,LHV,LHS,QC,QI,CPtot) &
    !$omp shared(KS,KE,IS,IE,JS,JE)
    do j = js, je
    do i = is, ie
    do k = ks, ke
       teml(k,i,j) = temp(k,i,j) &
                   - ( lhv(k,i,j) * qc(k,i,j) + lhs(k,i,j) * qi(k,i,j) ) / cptot(k,i,j)
    end do
    end do
    end do

    return

Here is the caller graph for this function: