mod_atmos_phy_cp_driver Module Reference

module ATMOSPHERE / Physics Cumulus More...

Functions/Subroutines
subroutine, public	atmos_phy_cp_driver_setup
	Setup. More...
subroutine, public	atmos_phy_cp_driver_resume
	Redume. More...
subroutine, public	atmos_phy_cp_driver (update_flag)
	Driver. More...

Detailed Description

module ATMOSPHERE / Physics Cumulus

Description

Cumulus parameterization driver

Author

Team SCALE

History

• 2014-05-04 (H.Yashiro) [new]

NAMELIST

• No namelist group

History Output

name	description	unit	variable
CBMFX	cloud base mass flux	kg/m2/s	MFLX_cloudbase

Function/Subroutine Documentation

atmos_phy_cp_driver_setup()

subroutine, public mod_atmos_phy_cp_driver::atmos_phy_cp_driver_setup

(

)

Setup.

Definition at line 52 of file mod_atmos_phy_cp_driver.f90.

References mod_atmos_admin::atmos_phy_cp_type, mod_atmos_admin::atmos_sw_phy_cp, scale_stdio::io_fid_log, and scale_stdio::io_l.

Referenced by mod_atmos_driver::atmos_driver_setup().

!    use scale_atmos_phy_cp, only: &
!       ATMOS_PHY_CP_setup
    use mod_atmos_admin, only: &
       atmos_phy_cp_type, &
       atmos_sw_phy_cp
    implicit none
    !---------------------------------------------------------------------------

    if( io_l ) write(io_fid_log,*)
    if( io_l ) write(io_fid_log,*) '++++++ Module[DRIVER] / Categ[ATMOS PHY_CP] / Origin[SCALE-RM]'

    if ( atmos_sw_phy_cp ) then

       ! setup library component
       !call ATMOS_PHY_CP_setup( ATMOS_PHY_CP_TYPE )

    else
       if( io_l ) write(io_fid_log,*) '*** this component is never called.'
    end if

    return

Here is the caller graph for this function:

atmos_phy_cp_driver_resume()

subroutine, public mod_atmos_phy_cp_driver::atmos_phy_cp_driver_resume

(

)

Redume.

Definition at line 78 of file mod_atmos_phy_cp_driver.f90.

References atmos_phy_cp_driver(), mod_atmos_admin::atmos_sw_phy_cp, scale_prof::prof_rapend(), and scale_prof::prof_rapstart().

Referenced by mod_atmos_driver::atmos_driver_resume2().

    use mod_atmos_admin, only: &
       atmos_sw_phy_cp
    implicit none

    if ( atmos_sw_phy_cp ) then

       ! run once (only for the diagnostic value)
       call prof_rapstart('ATM_Cumulus', 1)
       call atmos_phy_cp_driver( update_flag = .true. )
       call prof_rapend  ('ATM_Cumulus', 1)

    endif

    return

Here is the call graph for this function:

Here is the caller graph for this function:

atmos_phy_cp_driver()

subroutine, public mod_atmos_phy_cp_driver::atmos_phy_cp_driver

(

logical, intent(in)

update_flag

)

Driver.

Definition at line 97 of file mod_atmos_phy_cp_driver.f90.

References scale_tracer::aq_name, mod_atmos_phy_cp_vars::atmos_phy_cp_mflx_cloudbase, mod_atmos_phy_cp_vars::atmos_phy_cp_momx_t, mod_atmos_phy_cp_vars::atmos_phy_cp_momy_t, mod_atmos_phy_cp_vars::atmos_phy_cp_momz_t, mod_atmos_phy_cp_vars::atmos_phy_cp_rhoq_t, mod_atmos_phy_cp_vars::atmos_phy_cp_rhot_t, scale_grid_index::ie, scale_grid_index::is, scale_grid_index::je, scale_grid_index::js, scale_grid_index::ke, scale_grid_index::ks, mod_atmos_vars::momx_tp, mod_atmos_vars::momy_tp, mod_atmos_vars::momz_tp, scale_tracer::qa, mod_atmos_vars::rhoq_tp, mod_atmos_vars::rhot_tp, and scale_rm_statistics::statistics_checktotal.

Referenced by mod_atmos_driver::atmos_driver(), and atmos_phy_cp_driver_resume().

!     use scale_time, only: &
!        dt_CP => TIME_DTSEC_ATMOS_PHY_CP
    use scale_rm_statistics, only: &
       statistics_checktotal, &
       stat_total
    use scale_history, only: &
       hist_in
!    use scale_atmos_phy_cp, only: &
!       ATMOS_PHY_CP
    use mod_atmos_vars, only: &
!        DENS,              &
!        MOMZ,              &
!        MOMX,              &
!        MOMY,              &
!        RHOT,              &
!        QTRC,              &
       momz_t => momz_tp, &
       momx_t => momx_tp, &
       momy_t => momy_tp, &
       rhot_t => rhot_tp, &
       rhoq_t => rhoq_tp
    use mod_atmos_phy_cp_vars, only: &
       momz_t_cp      => atmos_phy_cp_momz_t,        &
       momx_t_cp      => atmos_phy_cp_momx_t,        &
       momy_t_cp      => atmos_phy_cp_momy_t,        &
       rhot_t_cp      => atmos_phy_cp_rhot_t,        &
       rhoq_t_cp      => atmos_phy_cp_rhoq_t,        &
       mflx_cloudbase => atmos_phy_cp_mflx_cloudbase
    implicit none

    logical, intent(in) :: update_flag

    real(RP) :: total ! dummy

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

    if ( update_flag ) then

!       call ATMOS_PHY_CP( DENS,          & ! [IN]
!                          MOMZ,          & ! [IN]
!                          MOMX,          & ! [IN]
!                          MOMY,          & ! [IN]
!                          RHOT,          & ! [IN]
!                          QTRC,          & ! [IN]
!                          MOMZ_t_CP,     & ! [INOUT]
!                          MOMX_t_CP,     & ! [INOUT]
!                          MOMY_t_CP,     & ! [INOUT]
!                          RHOT_t_CP,     & ! [INOUT]
!                          RHOQ_t_CP,     & ! [INOUT]
!                          MFLX_cloudbase ) ! [INOUT]

       ! tentative!
!OCL XFILL
       momz_t_cp(:,:,:)    = 0.0_rp
!OCL XFILL
       momx_t_cp(:,:,:)    = 0.0_rp
!OCL XFILL
       momy_t_cp(:,:,:)    = 0.0_rp
!OCL XFILL
       rhot_t_cp(:,:,:)    = 0.0_rp
!OCL XFILL
       rhoq_t_cp(:,:,:,:)  = 0.0_rp
!OCL XFILL
       mflx_cloudbase(:,:) = 0.0_rp

       call hist_in( mflx_cloudbase(:,:), 'CBMFX', 'cloud base mass flux', 'kg/m2/s' )
    endif

    !$omp parallel do private(i,j,k) OMP_SCHEDULE_ collapse(2)
    do j  = js, je
    do i  = is, ie
    do k  = ks, ke
       momz_t(k,i,j) = momz_t(k,i,j) + momz_t_cp(k,i,j)
       momx_t(k,i,j) = momx_t(k,i,j) + momx_t_cp(k,i,j)
       momy_t(k,i,j) = momy_t(k,i,j) + momy_t_cp(k,i,j)
       rhot_t(k,i,j) = rhot_t(k,i,j) + rhot_t_cp(k,i,j)
    enddo
    enddo
    enddo

    !$omp parallel do private(i,j,k) OMP_SCHEDULE_ collapse(3)
    do iq = 1, qa
    do j  = js, je
    do i  = is, ie
    do k  = ks, ke
       rhoq_t(k,i,j,iq) = rhoq_t(k,i,j,iq) + rhoq_t_cp(k,i,j,iq)
    enddo
    enddo
    enddo
    enddo

    if ( statistics_checktotal ) then
       call stat_total( total, momz_t_cp(:,:,:), 'MOMZ_t_CP' )
       call stat_total( total, momx_t_cp(:,:,:), 'MOMX_t_CP' )
       call stat_total( total, momy_t_cp(:,:,:), 'MOMY_t_CP' )
       call stat_total( total, rhot_t_cp(:,:,:), 'RHOT_t_CP' )

       do iq = 1, qa
          call stat_total( total, rhoq_t_cp(:,:,:,iq), trim(aq_name(iq))//'_t_CP' )
       enddo
    endif

    return

Here is the caller graph for this function: