mod_atmos_driver Module Reference

module ATMOSPHERE driver More...

Functions/Subroutines
subroutine, public	atmos_driver_tracer_setup
	Tracer setup. More...
subroutine, public	atmos_driver_setup (init)
	Setup. More...
subroutine, public	atmos_driver_calc_tendency (force)
	Calculation tendency. More...
subroutine, public	atmos_driver_calc_tendency_from_sflux (force)
	Calculation tendency from surface flux with coupler. More...
subroutine, public	atmos_driver_update
	advance atmospheric state More...
subroutine, public	atmos_driver_finalize (init)
	Finalize. More...
subroutine, public	atmos_surface_get
	Get surface boundary condition. More...
subroutine, public	atmos_surface_set (countup)
	Set surface boundary condition. More...
subroutine	atmos_driver_boundary_update

Detailed Description

module ATMOSPHERE driver

Description

Atmosphere module

Author

Team SCALE

Function/Subroutine Documentation

atmos_driver_tracer_setup()

subroutine, public mod_atmos_driver::atmos_driver_tracer_setup

Tracer setup.

Definition at line 54 of file mod_atmos_driver.F90.

    use mod_atmos_phy_mp_driver, only: &
       atmos_phy_mp_driver_tracer_setup, &
       atmos_phy_mp_user_qhyd2qtrc,      &
       atmos_phy_mp_driver_qhyd2qtrc_onlyqv
    use mod_atmos_phy_mp_vars, only: &
       qa_mp, &
       qs_mp, &
       qe_mp
    use mod_atmos_phy_ae_driver, only: &
       atmos_phy_ae_driver_tracer_setup
    use mod_atmos_phy_ch_driver, only: &
       atmos_phy_ch_driver_tracer_setup
    use mod_atmos_phy_tb_driver, only: &
       atmos_phy_tb_driver_tracer_setup
    use mod_atmos_phy_bl_driver, only: &
       atmos_phy_bl_driver_tracer_setup
    use scale_atmos_hydrometeor, only: &
       atmos_hydrometeor_regist, &
       atmos_hydrometeor_dry
    use mod_atmos_admin, only: &
       atmos_use_qv
    use mod_atmos_phy_lt_driver, only: &
       atmos_phy_lt_driver_tracer_setup
    implicit none
 
    !---------------------------------------------------------------------------
 
    log_newline
    log_info("ATMOS_driver_tracer_setup",*) 'Setup'
 
 
    call atmos_phy_mp_driver_tracer_setup
    call atmos_phy_ch_driver_tracer_setup
    call atmos_phy_ae_driver_tracer_setup
    call atmos_phy_tb_driver_tracer_setup
    call atmos_phy_bl_driver_tracer_setup
    call atmos_phy_lt_driver_tracer_setup
 
    if ( atmos_hydrometeor_dry .and. atmos_use_qv ) then
       log_info("ATMOS_driver_tracer_setup",*) 'Regist QV'
       call atmos_hydrometeor_regist( 0, 0,                                    & ! (in)
            (/'QV'/),                                                          & ! (in)
            (/'Ratio of Water Vapor mass to total mass (Specific humidity)'/), & ! (in)
            (/'kg/kg'/),                                                       & ! (in)
            qs_mp                                                              ) ! (out)
       qa_mp = 1
       qe_mp = qs_mp
       atmos_phy_mp_user_qhyd2qtrc => atmos_phy_mp_driver_qhyd2qtrc_onlyqv
 
    end if
 
    return

References scale_atmos_hydrometeor::atmos_hydrometeor_dry, scale_atmos_hydrometeor::atmos_hydrometeor_regist(), mod_atmos_phy_ae_driver::atmos_phy_ae_driver_tracer_setup(), mod_atmos_phy_bl_driver::atmos_phy_bl_driver_tracer_setup(), mod_atmos_phy_ch_driver::atmos_phy_ch_driver_tracer_setup(), mod_atmos_phy_lt_driver::atmos_phy_lt_driver_tracer_setup(), mod_atmos_phy_mp_driver::atmos_phy_mp_driver_qhyd2qtrc_onlyqv(), mod_atmos_phy_mp_driver::atmos_phy_mp_driver_tracer_setup(), mod_atmos_phy_mp_driver::atmos_phy_mp_user_qhyd2qtrc, mod_atmos_phy_tb_driver::atmos_phy_tb_driver_tracer_setup(), mod_atmos_admin::atmos_use_qv, mod_atmos_phy_mp_vars::qa_mp, mod_atmos_phy_mp_vars::qe_mp, and mod_atmos_phy_mp_vars::qs_mp.

Referenced by mod_rm_driver::rm_driver(), and mod_rm_prep::rm_prep().

Here is the call graph for this function:

Here is the caller graph for this function:

atmos_driver_setup()

subroutine, public mod_atmos_driver::atmos_driver_setup

(

logical, intent(in), optional

init

)

Setup.

Definition at line 111 of file mod_atmos_driver.F90.

    use scale_time, only: &
       time_nowdate
    use scale_atmos_solarins, only: &
       atmos_solarins_setup
    use scale_atmos_refstate, only: &
       atmos_refstate_setup
    use mod_atmos_bnd_driver, only: &
       atmos_boundary_driver_setup
    use mod_atmos_dyn_driver, only: &
       atmos_dyn_driver_setup
    use mod_atmos_phy_mp_driver, only: &
       atmos_phy_mp_driver_setup
    use mod_atmos_phy_mp_vars, only: &
       qa_mp, &
       qs_mp, &
       qe_mp
    use mod_atmos_phy_ch_driver, only: &
       atmos_phy_ch_driver_setup
    use mod_atmos_phy_ae_driver, only: &
       atmos_phy_ae_driver_setup
    use mod_atmos_phy_rd_driver, only: &
       atmos_phy_rd_driver_setup
    use mod_atmos_phy_sf_driver, only: &
       atmos_phy_sf_driver_setup
    use mod_atmos_phy_tb_driver, only: &
       atmos_phy_tb_driver_setup
    use mod_atmos_phy_bl_driver, only: &
       atmos_phy_bl_driver_setup
    use mod_atmos_phy_cp_driver, only: &
       atmos_phy_cp_driver_setup
    use mod_atmos_phy_lt_driver, only: &
       atmos_phy_lt_driver_setup
    use scale_atmos_grid_cartesc, only: &
       cz => atmos_grid_cartesc_cz, &
       fz => atmos_grid_cartesc_fz
    use scale_atmos_grid_cartesc_real, only: &
       base_lon => atmos_grid_cartesc_real_basepoint_lon, &
       base_lat => atmos_grid_cartesc_real_basepoint_lat, &
       real_cz  => atmos_grid_cartesc_real_cz, &
       real_fz  => atmos_grid_cartesc_real_fz, &
       real_phi => atmos_grid_cartesc_real_phi
    implicit none
    logical, intent(in), optional :: init
    logical :: not_init
    !---------------------------------------------------------------------------
 
    not_init = .true.
    if ( present(init) ) not_init = .not. init
 
    log_newline
    log_info("ATMOS_driver_setup",*) 'Setup'
 
    log_newline
    log_info("ATMOS_driver_setup",*) 'Setup each atmospheric components ...'
 
    !--- setup solar insolation
    call atmos_solarins_setup( base_lon, base_lat, time_nowdate(1) )
 
    call prof_rapstart('ATM_Refstate', 2)
    call atmos_refstate_setup( ka, ks, ke, ia, isb, ieb, ja, jsb, jeb, &
                               cz(:), fz(:), real_cz(:,:,:), real_fz(:,:,:), real_phi(:,:,:) )
    call prof_rapend  ('ATM_Refstate', 2)
 
    call prof_rapstart('ATM_Boundary', 2)
    if ( not_init ) call atmos_boundary_driver_setup
    call prof_rapend  ('ATM_Boundary', 2)
 
    ! setup each components
    if ( not_init ) call atmos_dyn_driver_setup
    if ( not_init ) call atmos_phy_lt_driver_setup
                    call atmos_phy_mp_driver_setup
    if ( not_init ) call atmos_phy_ae_driver_setup
    if ( not_init ) call atmos_phy_ch_driver_setup
    if ( not_init ) call atmos_phy_rd_driver_setup
    if ( not_init ) call atmos_phy_sf_driver_setup
    if ( not_init ) call atmos_phy_tb_driver_setup
                    call atmos_phy_bl_driver_setup
    if ( not_init ) call atmos_phy_cp_driver_setup
 
    log_newline
    log_info("ATMOS_driver_setup",*) 'Finish setup of each atmospheric components.'
 
    return

References mod_atmos_bnd_driver::atmos_boundary_driver_setup(), mod_atmos_dyn_driver::atmos_dyn_driver_setup(), scale_atmos_grid_cartesc::atmos_grid_cartesc_cz, scale_atmos_grid_cartesc::atmos_grid_cartesc_fz, scale_atmos_grid_cartesc_real::atmos_grid_cartesc_real_basepoint_lat, scale_atmos_grid_cartesc_real::atmos_grid_cartesc_real_basepoint_lon, scale_atmos_grid_cartesc_real::atmos_grid_cartesc_real_cz, scale_atmos_grid_cartesc_real::atmos_grid_cartesc_real_fz, scale_atmos_grid_cartesc_real::atmos_grid_cartesc_real_phi, mod_atmos_phy_ae_driver::atmos_phy_ae_driver_setup(), mod_atmos_phy_bl_driver::atmos_phy_bl_driver_setup(), mod_atmos_phy_ch_driver::atmos_phy_ch_driver_setup(), mod_atmos_phy_cp_driver::atmos_phy_cp_driver_setup(), mod_atmos_phy_lt_driver::atmos_phy_lt_driver_setup(), mod_atmos_phy_mp_driver::atmos_phy_mp_driver_setup(), mod_atmos_phy_rd_driver::atmos_phy_rd_driver_setup(), mod_atmos_phy_sf_driver::atmos_phy_sf_driver_setup(), mod_atmos_phy_tb_driver::atmos_phy_tb_driver_setup(), scale_atmos_refstate::atmos_refstate_setup(), scale_atmos_solarins::atmos_solarins_setup(), scale_atmos_grid_cartesc_index::ia, scale_atmos_grid_cartesc_index::ieb, scale_atmos_grid_cartesc_index::isb, scale_atmos_grid_cartesc_index::ja, scale_atmos_grid_cartesc_index::jeb, scale_atmos_grid_cartesc_index::jsb, scale_atmos_grid_cartesc_index::ka, scale_atmos_grid_cartesc_index::ke, scale_atmos_grid_cartesc_index::ks, scale_prof::prof_rapend(), scale_prof::prof_rapstart(), mod_atmos_phy_mp_vars::qa_mp, mod_atmos_phy_mp_vars::qe_mp, mod_atmos_phy_mp_vars::qs_mp, and scale_time::time_nowdate.

Referenced by mod_rm_driver::rm_driver(), and mod_rm_prep::rm_prep().

Here is the call graph for this function:

Here is the caller graph for this function:

atmos_driver_calc_tendency()

subroutine, public mod_atmos_driver::atmos_driver_calc_tendency

(

logical, intent(in)

force

)

Calculation tendency.

Definition at line 199 of file mod_atmos_driver.F90.

    use mod_atmos_vars, only: &
       dens_tp, &
       momz_tp, &
       rhou_tp, &
       rhov_tp, &
       rhot_tp, &
       rhoh_p,  &
       rhoq_tp, &
       momx_tp, &
       momy_tp
    use mod_atmos_phy_mp_driver, only: &
       atmos_phy_mp_driver_calc_tendency
    use mod_atmos_phy_ch_driver, only: &
       atmos_phy_ch_driver_calc_tendency
    use mod_atmos_phy_ae_driver, only: &
       atmos_phy_ae_driver_calc_tendency
    use mod_atmos_phy_rd_driver, only: &
       atmos_phy_rd_driver_calc_tendency
    use mod_atmos_phy_sf_driver, only: &
       atmos_phy_sf_driver_calc_tendency
    use mod_atmos_phy_tb_driver, only: &
       atmos_phy_tb_driver_calc_tendency
    use mod_atmos_phy_cp_driver, only: &
       atmos_phy_cp_driver_calc_tendency
    use mod_atmos_phy_bl_driver, only: &
       atmos_phy_bl_driver_calc_tendency
    use mod_admin_time, only: &
       do_phy_mp => time_doatmos_phy_mp, &
       do_phy_ae => time_doatmos_phy_ae, &
       do_phy_ch => time_doatmos_phy_ch, &
       do_phy_rd => time_doatmos_phy_rd, &
       do_phy_sf => time_doatmos_phy_sf, &
       do_phy_tb => time_doatmos_phy_tb, &
       do_phy_bl => time_doatmos_phy_bl, &
       do_phy_cp => time_doatmos_phy_cp, &
       do_phy_lt => time_doatmos_phy_lt
    use mod_atmos_admin, only: &
       atmos_sw_phy_mp, &
       atmos_sw_phy_ae, &
       atmos_sw_phy_ch, &
       atmos_sw_phy_rd, &
       atmos_sw_phy_sf, &
       atmos_sw_phy_tb, &
       atmos_sw_phy_bl, &
       atmos_sw_phy_cp
    use mod_cpl_admin, only: &
       cpl_sw
    implicit none
    logical, intent(in) :: force
    !---------------------------------------------------------------------------
 
    !########## Get Surface Boundary from coupler ##########
    call atmos_surface_get
 
    !########## calculate tendency ##########
    ! reset tendencies
    !$omp parallel workshare
    !$acc kernels
!OCL XFILL
    dens_tp(:,:,:)   = 0.0_rp
    !$acc end kernels
    !$acc kernels
!OCL XFILL
    momz_tp(:,:,:)   = 0.0_rp
    !$acc end kernels
    !$acc kernels
!OCL XFILL
    rhou_tp(:,:,:)   = 0.0_rp
    !$acc end kernels
    !$acc kernels
!OCL XFILL
    rhov_tp(:,:,:)   = 0.0_rp
    !$acc end kernels
    !$acc kernels
!OCL XFILL
    rhot_tp(:,:,:)   = 0.0_rp
    !$acc end kernels
    !$acc kernels
!OCL XFILL
    rhoh_p(:,:,:)   = 0.0_rp
    !$acc end kernels
    !$acc kernels
!OCL XFILL
    rhoq_tp(:,:,:,:) = 0.0_rp
    !$acc end kernels
    !$acc kernels
!OCL XFILL
    momx_tp(:,:,:)   = 0.0_rp
    !$acc end kernels
    !$acc kernels
!OCL XFILL
    momy_tp(:,:,:)   = 0.0_rp
    !$acc end kernels
    !$omp end parallel workshare
 
    ! Microphysics
    if ( atmos_sw_phy_mp ) then
       call prof_rapstart('ATM_Microphysics', 1)
       call atmos_phy_mp_driver_calc_tendency( update_flag = do_phy_mp .or. force )
       call prof_rapend  ('ATM_Microphysics', 1)
    endif
    ! Aerosol
    if ( atmos_sw_phy_ae ) then
       call prof_rapstart('ATM_Aerosol', 1)
       call atmos_phy_ae_driver_calc_tendency( update_flag = do_phy_ae .or. force )
       call prof_rapend  ('ATM_Aerosol', 1)
    endif
    ! Chemistry
    if ( atmos_sw_phy_ch ) then
       call prof_rapstart('ATM_Chemistry', 1)
       call atmos_phy_ch_driver_calc_tendency( update_flag = do_phy_ch .or. force )
       call prof_rapend  ('ATM_Chemistry', 1)
    endif
    ! Radiation
    if ( atmos_sw_phy_rd ) then
       call prof_rapstart('ATM_Radiation', 1)
       call atmos_phy_rd_driver_calc_tendency( update_flag = do_phy_rd .or. force )
       call prof_rapend  ('ATM_Radiation', 1)
    endif
    ! Turbulence
    if ( atmos_sw_phy_tb ) then
       call prof_rapstart('ATM_Turbulence', 1)
       call atmos_phy_tb_driver_calc_tendency( update_flag = do_phy_tb .or. force )
       call prof_rapend  ('ATM_Turbulence', 1)
    endif
    ! Cumulus
    if ( atmos_sw_phy_cp ) then
       call prof_rapstart('ATM_Cumulus', 1)
       call atmos_phy_cp_driver_calc_tendency( update_flag = do_phy_cp .or. force )
       call prof_rapend  ('ATM_Cumulus', 1)
    endif
    if ( .not. cpl_sw ) then
       ! Surface Flux
       if ( atmos_sw_phy_sf ) then
          call prof_rapstart('ATM_SurfaceFlux', 1)
          call atmos_phy_sf_driver_calc_tendency( update_flag = do_phy_sf .or. force )
          call prof_rapend  ('ATM_SurfaceFlux', 1)
       endif
       ! Planetary Boundary layer
       if ( atmos_sw_phy_bl ) then
          call prof_rapstart('ATM_PBL', 1)
          call atmos_phy_bl_driver_calc_tendency( update_flag = do_phy_bl .or. force )
          call prof_rapend  ('ATM_PBL', 1)
       endif
    end if
 
    !########## Set Surface Boundary Condition ##########
    call atmos_surface_set( countup = .true. )
 
    return

References mod_atmos_phy_ae_driver::atmos_phy_ae_driver_calc_tendency(), mod_atmos_phy_bl_driver::atmos_phy_bl_driver_calc_tendency(), mod_atmos_phy_ch_driver::atmos_phy_ch_driver_calc_tendency(), mod_atmos_phy_cp_driver::atmos_phy_cp_driver_calc_tendency(), mod_atmos_phy_mp_driver::atmos_phy_mp_driver_calc_tendency(), mod_atmos_phy_rd_driver::atmos_phy_rd_driver_calc_tendency(), mod_atmos_phy_sf_driver::atmos_phy_sf_driver_calc_tendency(), mod_atmos_phy_tb_driver::atmos_phy_tb_driver_calc_tendency(), atmos_surface_get(), atmos_surface_set(), mod_atmos_admin::atmos_sw_phy_ae, mod_atmos_admin::atmos_sw_phy_bl, mod_atmos_admin::atmos_sw_phy_ch, mod_atmos_admin::atmos_sw_phy_cp, mod_atmos_admin::atmos_sw_phy_mp, mod_atmos_admin::atmos_sw_phy_rd, mod_atmos_admin::atmos_sw_phy_sf, mod_atmos_admin::atmos_sw_phy_tb, mod_cpl_admin::cpl_sw, mod_atmos_vars::dens_tp, mod_atmos_vars::momx_tp, mod_atmos_vars::momy_tp, mod_atmos_vars::momz_tp, scale_prof::prof_rapend(), scale_prof::prof_rapstart(), mod_atmos_vars::rhoh_p, mod_atmos_vars::rhoq_tp, mod_atmos_vars::rhot_tp, mod_atmos_vars::rhou_tp, mod_atmos_vars::rhov_tp, mod_admin_time::time_doatmos_phy_ae, mod_admin_time::time_doatmos_phy_bl, mod_admin_time::time_doatmos_phy_ch, mod_admin_time::time_doatmos_phy_cp, mod_admin_time::time_doatmos_phy_lt, mod_admin_time::time_doatmos_phy_mp, mod_admin_time::time_doatmos_phy_rd, mod_admin_time::time_doatmos_phy_sf, and mod_admin_time::time_doatmos_phy_tb.

Referenced by mod_rm_driver::restart_read(), and mod_rm_driver::rm_driver().

Here is the call graph for this function:

Here is the caller graph for this function:

atmos_driver_calc_tendency_from_sflux()

subroutine, public mod_atmos_driver::atmos_driver_calc_tendency_from_sflux

(

logical, intent(in)

force

)

Calculation tendency from surface flux with coupler.

Definition at line 354 of file mod_atmos_driver.F90.

    use mod_atmos_phy_sf_driver, only: &
       atmos_phy_sf_driver_calc_tendency
    use mod_atmos_phy_bl_driver, only: &
       atmos_phy_bl_driver_calc_tendency
    use mod_cpl_admin, only: &
       cpl_sw
    use mod_atmos_admin, only: &
       atmos_sw_phy_sf, &
       atmos_sw_phy_bl
    use mod_admin_time, only: &
       do_phy_sf => time_doatmos_phy_sf, &
       do_phy_bl => time_doatmos_phy_bl
    implicit none
    logical, intent(in) :: force
    !---------------------------------------------------------------------------
 
    if ( cpl_sw ) then
 
       !########## Get Surface Boundary Condition ##########
       call atmos_surface_get
 
       ! Surface Flux
       if ( atmos_sw_phy_sf ) then
          call prof_rapstart('ATM_SurfaceFlux', 1)
          call atmos_phy_sf_driver_calc_tendency( update_flag = do_phy_sf .or. force )
          call prof_rapend  ('ATM_SurfaceFlux', 1)
       endif
 
       ! Planetary Boundary layer
       if ( atmos_sw_phy_bl ) then
          call prof_rapstart('ATM_PBL', 1)
          call atmos_phy_bl_driver_calc_tendency( update_flag = do_phy_bl .or. force )
          call prof_rapend  ('ATM_PBL', 1)
       endif
 
    end if
 
    return

References mod_atmos_phy_bl_driver::atmos_phy_bl_driver_calc_tendency(), mod_atmos_phy_sf_driver::atmos_phy_sf_driver_calc_tendency(), atmos_surface_get(), mod_atmos_admin::atmos_sw_phy_bl, mod_atmos_admin::atmos_sw_phy_sf, mod_cpl_admin::cpl_sw, scale_prof::prof_rapend(), scale_prof::prof_rapstart(), mod_admin_time::time_doatmos_phy_bl, and mod_admin_time::time_doatmos_phy_sf.

Referenced by mod_rm_driver::restart_read(), and mod_rm_driver::rm_driver().

Here is the call graph for this function:

Here is the caller graph for this function:

atmos_driver_update()

subroutine, public mod_atmos_driver::atmos_driver_update

advance atmospheric state

Definition at line 397 of file mod_atmos_driver.F90.

    use mod_atmos_admin, only: &
       atmos_sw_dyn,    &
       atmos_sw_phy_mp, &
       atmos_sw_phy_ae, &
       atmos_sw_phy_lt
    use mod_admin_time, only: &
       do_dyn    => time_doatmos_dyn,    &
       do_phy_mp => time_doatmos_phy_mp, &
       do_phy_ae => time_doatmos_phy_ae
    use scale_atmos_refstate, only: &
       atmos_refstate_update
    use mod_atmos_vars, only: &
       atmos_vars_calc_diagnostics,&
       atmos_vars_fillhalo,        &
       dens,                       &
       temp,                       &
       pres,                       &
       pott,                       &
       qv
    use mod_atmos_bnd_driver, only: &
       atmos_boundary_driver_update, &
       atmos_boundary_driver_send,   &
       atmos_boundary_update_flag
    use mod_atmos_dyn_driver, only: &
       atmos_dyn_driver
    use mod_atmos_phy_mp_driver, only: &
       atmos_phy_mp_driver_adjustment
    use mod_atmos_phy_ae_driver, only: &
       atmos_phy_ae_driver_adjustment
    use mod_atmos_phy_lt_driver, only: &
       atmos_phy_lt_driver_adjustment
    use scale_atmos_grid_cartesc, only: &
       cz   => atmos_grid_cartesc_cz,  &
       fz   => atmos_grid_cartesc_fz,  &
       fdz  => atmos_grid_cartesc_fdz, &
       rcdz => atmos_grid_cartesc_rcdz
    use scale_atmos_grid_cartesc_real, only: &
       real_cz  => atmos_grid_cartesc_real_cz, &
       real_fz  => atmos_grid_cartesc_real_fz, &
       real_phi => atmos_grid_cartesc_real_phi, &
       area     => atmos_grid_cartesc_real_area
    use scale_time, only: &
       time_nowdaysec, &
       time_dtsec
    implicit none
 
    !---------------------------------------------------------------------------
 
    !########## Dynamics ##########
    if ( atmos_sw_dyn ) then
       if ( atmos_boundary_update_flag ) then
          call prof_rapstart('ATM_Boundary', 2)
          call atmos_boundary_driver_update( time_nowdaysec - time_dtsec * 0.5_dp )
          call prof_rapend  ('ATM_Boundary', 2)
          call atmos_vars_fillhalo
       endif
       call prof_rapstart('ATM_Dynamics', 1)
       call atmos_dyn_driver( do_dyn )
       call prof_rapend  ('ATM_Dynamics', 1)
    endif
 
 
    !########## Calculate diagnostic variables ##########
    call atmos_vars_calc_diagnostics
 
 
    !########## Adjustment ##########
    ! Microphysics
    if ( atmos_sw_phy_mp ) then
       call prof_rapstart('ATM_Microphysics', 1)
       call atmos_phy_mp_driver_adjustment
       call prof_rapend  ('ATM_Microphysics', 1)
       call atmos_vars_fillhalo
       call atmos_vars_calc_diagnostics
    endif
    ! Aerosol
    if ( atmos_sw_phy_ae ) then
       call prof_rapstart('ATM_Aerosol', 1)
       call atmos_phy_ae_driver_adjustment
       call prof_rapend  ('ATM_Aerosol', 1)
       call atmos_vars_fillhalo
       call atmos_vars_calc_diagnostics
    endif
    ! Lightning
    if ( atmos_sw_phy_lt ) then
       call prof_rapstart('ATM_Lightning', 1)
       call atmos_phy_lt_driver_adjustment
       call prof_rapend  ('ATM_Lightning', 1)
       call atmos_vars_fillhalo
       ! calc_diagnostics is not necessary
    end if
 
 
    !########## Lateral/Top Boundary Condition ###########
    if ( atmos_boundary_update_flag ) then
       call prof_rapstart('ATM_Boundary', 2)
       call atmos_boundary_driver_update( time_nowdaysec )
       call prof_rapend  ('ATM_Boundary', 2)
       call atmos_vars_fillhalo
       call atmos_vars_calc_diagnostics
    endif
 
 
    !########## Send Lateral/Top Boundary Condition (Online nesting) ###########
    if ( atmos_boundary_update_flag ) then
       call prof_rapstart('ATM_Boundary', 2)
       call atmos_boundary_driver_send
       call prof_rapend  ('ATM_Boundary', 2)
    endif
 
 
    !########## Reference State ###########
    call prof_rapstart('ATM_Refstate', 2)
    call atmos_refstate_update( ka, ks, ke, ia, is, ie, isb, ieb, ja, js, je, jsb, jeb, &
                                dens(:,:,:), pott(:,:,:), temp(:,:,:), pres(:,:,:), qv(:,:,:), & ! [IN]
                                cz(:), fz(:), fdz(:), rcdz(:),                                 & ! [IN]
                                real_cz(:,:,:), real_fz(:,:,:), real_phi(:,:,:), area(:,:),    & ! [IN]
                                time_nowdaysec                                                 ) ! [IN]
    call prof_rapend  ('ATM_Refstate', 2)
 
    return

References mod_atmos_bnd_driver::atmos_boundary_driver_send(), mod_atmos_bnd_driver::atmos_boundary_driver_update(), mod_atmos_bnd_driver::atmos_boundary_update_flag, mod_atmos_dyn_driver::atmos_dyn_driver(), scale_atmos_grid_cartesc::atmos_grid_cartesc_cz, scale_atmos_grid_cartesc::atmos_grid_cartesc_fdz, scale_atmos_grid_cartesc::atmos_grid_cartesc_fz, scale_atmos_grid_cartesc::atmos_grid_cartesc_rcdz, scale_atmos_grid_cartesc_real::atmos_grid_cartesc_real_area, scale_atmos_grid_cartesc_real::atmos_grid_cartesc_real_cz, scale_atmos_grid_cartesc_real::atmos_grid_cartesc_real_fz, scale_atmos_grid_cartesc_real::atmos_grid_cartesc_real_phi, mod_atmos_phy_ae_driver::atmos_phy_ae_driver_adjustment(), mod_atmos_phy_lt_driver::atmos_phy_lt_driver_adjustment(), mod_atmos_phy_mp_driver::atmos_phy_mp_driver_adjustment(), scale_atmos_refstate::atmos_refstate_update(), mod_atmos_admin::atmos_sw_dyn, mod_atmos_admin::atmos_sw_phy_ae, mod_atmos_admin::atmos_sw_phy_lt, mod_atmos_admin::atmos_sw_phy_mp, mod_atmos_vars::atmos_vars_calc_diagnostics(), mod_atmos_vars::atmos_vars_fillhalo(), mod_atmos_vars::dens, scale_atmos_grid_cartesc_index::ia, scale_atmos_grid_cartesc_index::ie, scale_atmos_grid_cartesc_index::ieb, scale_atmos_grid_cartesc_index::is, scale_atmos_grid_cartesc_index::isb, scale_atmos_grid_cartesc_index::ja, scale_atmos_grid_cartesc_index::je, scale_atmos_grid_cartesc_index::jeb, scale_atmos_grid_cartesc_index::js, scale_atmos_grid_cartesc_index::jsb, scale_atmos_grid_cartesc_index::ka, scale_atmos_grid_cartesc_index::ke, scale_atmos_grid_cartesc_index::ks, mod_atmos_vars::pott, mod_atmos_vars::pres, scale_prof::prof_rapend(), scale_prof::prof_rapstart(), mod_atmos_vars::qv, mod_atmos_vars::temp, mod_admin_time::time_doatmos_dyn, mod_admin_time::time_doatmos_phy_ae, mod_admin_time::time_doatmos_phy_mp, scale_time::time_dtsec, and scale_time::time_nowdaysec.

Referenced by mod_rm_driver::rm_driver().

Here is the call graph for this function:

Here is the caller graph for this function:

atmos_driver_finalize()

subroutine, public mod_atmos_driver::atmos_driver_finalize

(

logical, intent(in), optional

init

)

Finalize.

Definition at line 523 of file mod_atmos_driver.F90.

    use scale_atmos_refstate, only: &
       atmos_refstate_finalize
    use mod_atmos_bnd_driver, only: &
       atmos_boundary_driver_finalize
    use mod_atmos_dyn_driver, only: &
       atmos_dyn_driver_finalize
    use mod_atmos_phy_lt_driver, only: &
       atmos_phy_lt_driver_finalize
    use mod_atmos_phy_mp_driver, only: &
       atmos_phy_mp_driver_finalize
    use mod_atmos_phy_ae_driver, only: &
       atmos_phy_ae_driver_finalize
    use mod_atmos_phy_ch_driver, only: &
       atmos_phy_ch_driver_finalize
    use mod_atmos_phy_rd_driver, only: &
       atmos_phy_rd_driver_finalize
    use mod_atmos_phy_tb_driver, only: &
       atmos_phy_tb_driver_finalize
    use mod_atmos_phy_bl_driver, only: &
       atmos_phy_bl_driver_finalize
    use mod_atmos_phy_cp_driver, only: &
       atmos_phy_cp_driver_finalize
    implicit none
    logical, intent(in), optional :: init
    logical :: not_init
    !---------------------------------------------------------------------------
 
    log_newline
    log_info("ATMOS_driver_finalize",*) 'Finalize'
 
    not_init = .true.
    if ( present(init) ) not_init = .not. init
 
    if ( not_init ) call atmos_dyn_driver_finalize
    if ( not_init ) call atmos_phy_lt_driver_finalize
                    call atmos_phy_mp_driver_finalize
    if ( not_init ) call atmos_phy_ae_driver_finalize
    if ( not_init ) call atmos_phy_ch_driver_finalize
    if ( not_init ) call atmos_phy_rd_driver_finalize
    !if ( not_init ) call ATMOS_PHY_SF_driver_finalize
    if ( not_init ) call atmos_phy_tb_driver_finalize
                    call atmos_phy_bl_driver_finalize
    if ( not_init )call atmos_phy_cp_driver_finalize
 
    if ( not_init ) call atmos_boundary_driver_finalize
 
    call atmos_refstate_finalize
 
    return

References mod_atmos_bnd_driver::atmos_boundary_driver_finalize(), mod_atmos_dyn_driver::atmos_dyn_driver_finalize(), mod_atmos_phy_ae_driver::atmos_phy_ae_driver_finalize(), mod_atmos_phy_bl_driver::atmos_phy_bl_driver_finalize(), mod_atmos_phy_ch_driver::atmos_phy_ch_driver_finalize(), mod_atmos_phy_cp_driver::atmos_phy_cp_driver_finalize(), mod_atmos_phy_lt_driver::atmos_phy_lt_driver_finalize(), mod_atmos_phy_mp_driver::atmos_phy_mp_driver_finalize(), mod_atmos_phy_rd_driver::atmos_phy_rd_driver_finalize(), mod_atmos_phy_tb_driver::atmos_phy_tb_driver_finalize(), and scale_atmos_refstate::atmos_refstate_finalize().

Referenced by mod_rm_driver::rm_driver(), and mod_rm_prep::rm_prep().

Here is the call graph for this function:

Here is the caller graph for this function:

atmos_surface_get()

subroutine, public mod_atmos_driver::atmos_surface_get

Get surface boundary condition.

Definition at line 577 of file mod_atmos_driver.F90.

    use mod_atmos_phy_sf_vars, only: &
       sfc_temp   => atmos_phy_sf_sfc_temp,   &
       sfc_albedo => atmos_phy_sf_sfc_albedo, &
       sfc_z0m    => atmos_phy_sf_sfc_z0m,    &
       sfc_z0h    => atmos_phy_sf_sfc_z0h,    &
       sfc_z0e    => atmos_phy_sf_sfc_z0e,    &
       sflx_mw    => atmos_phy_sf_sflx_mw,    &
       sflx_mu    => atmos_phy_sf_sflx_mu,    &
       sflx_mv    => atmos_phy_sf_sflx_mv,    &
       sflx_sh    => atmos_phy_sf_sflx_sh,    &
       sflx_lh    => atmos_phy_sf_sflx_lh,    &
       sflx_shex  => atmos_phy_sf_sflx_shex,  &
       sflx_lhex  => atmos_phy_sf_sflx_lhex,  &
       sflx_qvex  => atmos_phy_sf_sflx_qvex,  &
       sflx_gh    => atmos_phy_sf_sflx_gh,    &
       sflx_qtrc  => atmos_phy_sf_sflx_qtrc,  &
       sflx_engi  => atmos_phy_sf_sflx_engi,  &
       u10        => atmos_phy_sf_u10,        &
       v10        => atmos_phy_sf_v10,        &
       t2         => atmos_phy_sf_t2,         &
       q2         => atmos_phy_sf_q2
    use mod_cpl_admin, only: &
       cpl_sw
    use mod_cpl_vars, only: &
       cpl_getsfc_atm
    implicit none
    !---------------------------------------------------------------------------
 
    call prof_rapstart('ATM_SfcExch', 2)
 
    if ( cpl_sw ) then
       call cpl_getsfc_atm( sfc_temp(:,:),     & ! [OUT]
                            sfc_albedo(:,:,:,:), & ! [OUT]
                            sfc_z0m(:,:),     & ! [OUT]
                            sfc_z0h(:,:),     & ! [OUT]
                            sfc_z0e(:,:),     & ! [OUT]
                            sflx_mw(:,:),     & ! [OUT]
                            sflx_mu(:,:),     & ! [OUT]
                            sflx_mv(:,:),     & ! [OUT]
                            sflx_sh(:,:),     & ! [OUT]
                            sflx_lh(:,:),     & ! [OUT]
                            sflx_shex(:,:),     & ! [OUT]
                            sflx_lhex(:,:),     & ! [OUT]
                            sflx_qvex(:,:),     & ! [OUT]
                            sflx_gh(:,:),     & ! [OUT]
                            sflx_qtrc(:,:,:),   & ! [OUT]
                            sflx_engi(:,:),     & ! [OUT]
                            u10(:,:),     & ! [OUT]
                            v10(:,:),     & ! [OUT]
                            t2(:,:),     & ! [OUT]
                            q2(:,:)      ) ! [OUT]
    endif
 
    call prof_rapend  ('ATM_SfcExch', 2)
 
    return

References mod_atmos_phy_sf_vars::atmos_phy_sf_q2, mod_atmos_phy_sf_vars::atmos_phy_sf_sfc_albedo, mod_atmos_phy_sf_vars::atmos_phy_sf_sfc_temp, mod_atmos_phy_sf_vars::atmos_phy_sf_sfc_z0e, mod_atmos_phy_sf_vars::atmos_phy_sf_sfc_z0h, mod_atmos_phy_sf_vars::atmos_phy_sf_sfc_z0m, mod_atmos_phy_sf_vars::atmos_phy_sf_sflx_engi, mod_atmos_phy_sf_vars::atmos_phy_sf_sflx_gh, mod_atmos_phy_sf_vars::atmos_phy_sf_sflx_lh, mod_atmos_phy_sf_vars::atmos_phy_sf_sflx_lhex, mod_atmos_phy_sf_vars::atmos_phy_sf_sflx_mu, mod_atmos_phy_sf_vars::atmos_phy_sf_sflx_mv, mod_atmos_phy_sf_vars::atmos_phy_sf_sflx_mw, mod_atmos_phy_sf_vars::atmos_phy_sf_sflx_qtrc, mod_atmos_phy_sf_vars::atmos_phy_sf_sflx_qvex, mod_atmos_phy_sf_vars::atmos_phy_sf_sflx_sh, mod_atmos_phy_sf_vars::atmos_phy_sf_sflx_shex, mod_atmos_phy_sf_vars::atmos_phy_sf_t2, mod_atmos_phy_sf_vars::atmos_phy_sf_u10, mod_atmos_phy_sf_vars::atmos_phy_sf_v10, mod_cpl_vars::cpl_getsfc_atm(), mod_cpl_admin::cpl_sw, scale_prof::prof_rapend(), and scale_prof::prof_rapstart().

Referenced by atmos_driver_calc_tendency(), atmos_driver_calc_tendency_from_sflux(), and mod_rm_prep::rm_prep().

Here is the call graph for this function:

Here is the caller graph for this function:

atmos_surface_set()

subroutine, public mod_atmos_driver::atmos_surface_set

(

logical, intent(in)

countup

)

Set surface boundary condition.

Definition at line 638 of file mod_atmos_driver.F90.

    use scale_atmos_grid_cartesc_real, only: &
       real_fz => atmos_grid_cartesc_real_fz
    use scale_atmos_bottom, only: &
       bottom_estimate => atmos_bottom_estimate
    use mod_atmos_vars, only: &
       dens, &
       qv,   &
       temp, &
       pres, &
       w,    &
       u,    &
       v,    &
       prec, &
       prec_engi
    use mod_atmos_phy_sf_vars, only: &
       sfc_temp  => atmos_phy_sf_sfc_temp
    use mod_atmos_phy_mp_vars, only: &
       sflx_rain_mp => atmos_phy_mp_sflx_rain, &
       sflx_snow_mp => atmos_phy_mp_sflx_snow, &
       sflx_engi_mp => atmos_phy_mp_sflx_engi
    use mod_atmos_phy_cp_vars, only: &
       sflx_rain_cp => atmos_phy_cp_sflx_rain, &
       sflx_snow_cp => atmos_phy_cp_sflx_snow, &
       sflx_engi_cp => atmos_phy_cp_sflx_engi
    use mod_atmos_phy_rd_vars, only: &
       sflx_rad_dn => atmos_phy_rd_sflx_down, &
       cossza      => atmos_phy_rd_cossza
    use mod_atmos_phy_bl_vars, only: &
       atm_pbl => atmos_phy_bl_zi
    use mod_cpl_admin, only: &
       cpl_sw
    use mod_cpl_vars, only: &
       cpl_putatm
    implicit none
 
    ! arguments
    logical, intent(in) :: countup
 
    ! works
    real(RP) :: SFC_DENS(IA,JA)
    real(RP) :: SFC_PRES(IA,JA)
 
    real(RP) :: TEMP1(IA,JA)
    real(RP) :: PRES1(IA,JA)
    real(RP) :: W1   (IA,JA)
    real(RP) :: U1   (IA,JA)
    real(RP) :: V1   (IA,JA)
    real(RP) :: DENS1(IA,JA)
    real(RP) :: QV1  (IA,JA)
 
    integer  :: i,j
    !---------------------------------------------------------------------------
 
    call prof_rapstart('ATM_SfcExch', 2)
 
    ! sum of rainfall from mp and cp
    !$omp parallel do private(i,j) OMP_SCHEDULE_
    !$acc kernels
    do j = js, je
    do i = is, ie
       prec(i,j) = sflx_rain_mp(i,j) + sflx_rain_cp(i,j) + sflx_snow_mp(i,j) + sflx_snow_cp(i,j)
       prec_engi(i,j) = sflx_engi_mp(i,j) + sflx_engi_cp(i,j)
    enddo
    enddo
    !$acc end kernels
 
    if ( cpl_sw ) then
 
       !$acc data create(SFC_DENS,SFC_PRES,TEMP1,PRES1,W1,U1,V1,DENS1,QV1)
 
       ! planetary boundary layer
       call bottom_estimate( ka, ks, ke, ia, is, ie, ja, js, je, &
                             dens(:,:,:), pres(:,:,:), qv(:,:,:), & ! [IN]
                             sfc_temp(:,:),                       & ! [IN]
                             real_fz(:,:,:),                      & ! [IN]
                             sfc_dens(:,:), sfc_pres(:,:)         ) ! [OUT]
 
       !$omp parallel do
       !$acc kernels
       do j = js, je
       do i = is, ie
          temp1(i,j) = temp(ks,i,j)
          pres1(i,j) = pres(ks,i,j)
          w1(i,j) = w(ks,i,j)
          u1(i,j) = u(ks,i,j)
          v1(i,j) = v(ks,i,j)
          dens1(i,j) = dens(ks,i,j)
          qv1(i,j) = qv(ks,i,j)
       end do
       end do
       !$acc end kernels
 
       call cpl_putatm( temp1(:,:),     & ! [IN]
                        pres1(:,:),     & ! [IN]
                        w1(:,:),     & ! [IN]
                        u1(:,:),     & ! [IN]
                        v1(:,:),     & ! [IN]
                        dens1(:,:),     & ! [IN]
                        qv1(:,:),     & ! [IN]
                        atm_pbl(:,:),     & ! [IN]
                        sfc_dens(:,:),     & ! [IN]
                        sfc_pres(:,:),     & ! [IN]
                        sflx_rad_dn(:,:,:,:), & ! [IN]
                        cossza(:,:),     & ! [IN]
                        prec(:,:),     & ! [IN]
                        prec_engi(:,:),     & ! [IN]
                        countup               ) ! [IN]
 
       !$acc end data
 
    endif
 
    call prof_rapend  ('ATM_SfcExch', 2)
 
    return

References scale_atmos_bottom::atmos_bottom_estimate(), scale_atmos_grid_cartesc_real::atmos_grid_cartesc_real_fz, mod_atmos_phy_bl_vars::atmos_phy_bl_zi, mod_atmos_phy_cp_vars::atmos_phy_cp_sflx_engi, mod_atmos_phy_cp_vars::atmos_phy_cp_sflx_rain, mod_atmos_phy_cp_vars::atmos_phy_cp_sflx_snow, mod_atmos_phy_mp_vars::atmos_phy_mp_sflx_engi, mod_atmos_phy_mp_vars::atmos_phy_mp_sflx_rain, mod_atmos_phy_mp_vars::atmos_phy_mp_sflx_snow, mod_atmos_phy_rd_vars::atmos_phy_rd_cossza, mod_atmos_phy_rd_vars::atmos_phy_rd_sflx_down, mod_atmos_phy_sf_vars::atmos_phy_sf_sfc_temp, mod_cpl_vars::cpl_putatm(), mod_cpl_admin::cpl_sw, mod_atmos_vars::dens, scale_atmos_grid_cartesc_index::ia, scale_atmos_grid_cartesc_index::ie, scale_atmos_grid_cartesc_index::is, scale_atmos_grid_cartesc_index::ja, scale_atmos_grid_cartesc_index::je, scale_atmos_grid_cartesc_index::js, scale_atmos_grid_cartesc_index::ka, scale_atmos_grid_cartesc_index::ke, scale_atmos_grid_cartesc_index::ks, mod_atmos_vars::prec, mod_atmos_vars::prec_engi, mod_atmos_vars::pres, scale_prof::prof_rapend(), scale_prof::prof_rapstart(), mod_atmos_vars::qv, mod_atmos_vars::temp, mod_atmos_vars::u, mod_atmos_vars::v, and mod_atmos_vars::w.

Referenced by atmos_driver_calc_tendency(), and mod_rm_driver::restart_read().

Here is the call graph for this function:

Here is the caller graph for this function:

atmos_driver_boundary_update()

subroutine mod_atmos_driver::atmos_driver_boundary_update

Definition at line 756 of file mod_atmos_driver.F90.