mod_atmos_phy_ch_driver.f90

Go to the documentation of this file.

!-------------------------------------------------------------------------------
!-------------------------------------------------------------------------------
#include "inc_openmp.h"
module mod_atmos_phy_ch_driver
  !-----------------------------------------------------------------------------
  !
  !++ used modules
  !
  use scale_precision
  use scale_stdio
  use scale_prof
  use scale_grid_index
  use scale_tracer
  !-----------------------------------------------------------------------------
  implicit none
  private
  !-----------------------------------------------------------------------------
  !
  !++ Public procedure
  !
  public :: atmos_phy_ch_driver_setup
  public :: atmos_phy_ch_driver_resume
  public :: atmos_phy_ch_driver

  !-----------------------------------------------------------------------------
  !
  !++ Public parameters & variables
  !
  !-----------------------------------------------------------------------------
  !
  !++ Private procedure
  !
  !-----------------------------------------------------------------------------
  !
  !++ Private parameters & variables
  !
  !-----------------------------------------------------------------------------
contains
  !-----------------------------------------------------------------------------
  subroutine atmos_phy_ch_driver_setup
!    use scale_atmos_phy_ch, only: &
!       ATMOS_PHY_CH_setup
    use mod_atmos_admin, only: &
       atmos_phy_ch_type, &
       atmos_sw_phy_ch
    implicit none
    !---------------------------------------------------------------------------

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

    if ( atmos_sw_phy_ch ) then

       ! setup library component
       !call ATMOS_PHY_CH_setup( ATMOS_PHY_CH_TYPE )

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

    return
  end subroutine atmos_phy_ch_driver_setup

  !-----------------------------------------------------------------------------
  subroutine atmos_phy_ch_driver_resume
    use mod_atmos_admin, only: &
       atmos_sw_phy_ch
    implicit none

    if ( atmos_sw_phy_ch ) then

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

    end if

    return
  end subroutine atmos_phy_ch_driver_resume

  !-----------------------------------------------------------------------------
  subroutine atmos_phy_ch_driver( update_flag )
!     use scale_time, only: &
!        dt_CH => TIME_DTSEC_ATMOS_PHY_CH
    use scale_rm_statistics, only: &
       statistics_checktotal, &
       stat_total
    use scale_history, only: &
       hist_in
!    use scale_atmos_phy_ch, only: &
!       ATMOS_PHY_CH
    use mod_atmos_vars, only: &
!        DENS,              &
!        MOMZ,              &
!        MOMX,              &
!        MOMY,              &
!        RHOT,              &
!        QTRC,              &
       rhoq_t => rhoq_tp
    use mod_atmos_phy_ch_vars, only: &
       rhoq_t_ch => atmos_phy_ch_rhoq_t, &
       o3        => atmos_phy_ch_o3
    implicit none

    logical, intent(in) :: update_flag

    real(RP) :: total ! dummy

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

    if ( update_flag ) then

!       call ATMOS_PHY_CH( DENS,          & ! [IN]
!                          MOMZ,          & ! [IN]
!                          MOMX,          & ! [IN]
!                          MOMY,          & ! [IN]
!                          RHOT,          & ! [IN]
!                          QTRC,          & ! [IN]
!                          RHOQ_t_CH,     & ! [INOUT]
!                          O3             ) ! [INOUT]

       ! tentative!
!OCL XFILL
       rhoq_t_ch(:,:,:,:) = 0.0_rp
!OCL XFILL
       o3(:,:,:)   = 0.0_rp

       call hist_in( o3(:,:,:), 'Ozone', 'Ozone', 'PPM' )
    endif

    !$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_ch(k,i,j,iq)
    enddo
    enddo
    enddo
    enddo

    if ( statistics_checktotal ) then
       do iq = 1, qa
          call stat_total( total, rhoq_t_ch(:,:,:,iq), trim(aq_name(iq))//'_t_CH' )
       enddo
    endif

    return
  end subroutine atmos_phy_ch_driver

end module mod_atmos_phy_ch_driver