mod_atmos_admin.F90

Go to the documentation of this file.

!-------------------------------------------------------------------------------
!-------------------------------------------------------------------------------
#include "scalelib.h"
module mod_atmos_admin
  !-----------------------------------------------------------------------------
  !
  !++ used modules
  !
  use scale_precision
  use scale_io
  use scale_prof
  !-----------------------------------------------------------------------------
  implicit none
  private
  !-----------------------------------------------------------------------------
  !
  !++ Public procedure
  !
  public :: atmos_admin_setup
  public :: atmos_admin_getscheme

  !-----------------------------------------------------------------------------
  !
  !++ Public parameters & variables
  !
  logical,                public :: atmos_do          = .true. ! main switch for the model

  character(len=H_SHORT), public :: atmos_dyn_type    = 'NONE'
  character(len=H_SHORT), public :: atmos_phy_mp_type = 'NONE'
  character(len=H_SHORT), public :: atmos_phy_ae_type = 'NONE'
  character(len=H_SHORT), public :: atmos_phy_ch_type = 'NONE'
  character(len=H_SHORT), public :: atmos_phy_rd_type = 'NONE'
  character(len=H_SHORT), public :: atmos_phy_sf_type = 'NONE'
  character(len=H_SHORT), public :: atmos_phy_tb_type = 'NONE'
  character(len=H_SHORT), public :: atmos_phy_bl_type = 'NONE'
  character(len=H_SHORT), public :: atmos_phy_cp_type = 'NONE'

  character(len=H_SHORT), public :: atmos_phy_precip_type = 'Upwind-Euler'

  logical,                public :: atmos_use_average = .false.

  logical,                public :: atmos_sw_dyn
  logical,                public :: atmos_sw_phy_mp
  logical,                public :: atmos_sw_phy_ae
  logical,                public :: atmos_sw_phy_ch
  logical,                public :: atmos_sw_phy_rd
  logical,                public :: atmos_sw_phy_sf
  logical,                public :: atmos_sw_phy_tb
  logical,                public :: atmos_sw_phy_bl
  logical,                public :: atmos_sw_phy_cp

  !-----------------------------------------------------------------------------
  !
  !++ Private procedure
  !
  !-----------------------------------------------------------------------------
  !
  !++ Private parameters & variables
  !
  !-----------------------------------------------------------------------------
contains
  !-----------------------------------------------------------------------------
  subroutine atmos_admin_setup
    use scale_prc, only: &
       prc_abort
    implicit none

    namelist / param_atmos / &
       atmos_do,          &
       atmos_dyn_type,    &
       atmos_phy_mp_type, &
       atmos_phy_ae_type, &
       atmos_phy_ch_type, &
       atmos_phy_rd_type, &
       atmos_phy_sf_type, &
       atmos_phy_tb_type, &
       atmos_phy_bl_type, &
       atmos_phy_cp_type, &
       atmos_phy_precip_type, &
       atmos_use_average

    integer :: ierr
    !---------------------------------------------------------------------------

    log_newline
    log_info("ATMOS_ADMIN_setup",*) 'Setup'

    !--- read namelist
    rewind(io_fid_conf)
    read(io_fid_conf,nml=param_atmos,iostat=ierr)
    if( ierr < 0 ) then !--- missing
       log_info("ATMOS_ADMIN_setup",*) 'Not found namelist. Default used.'
    elseif( ierr > 0 ) then !--- fatal error
       log_error("ATMOS_ADMIN_setup",*) 'Not appropriate names in namelist PARAM_ATMOS. Check!'
       call prc_abort
    endif
    log_nml(param_atmos)

    !-----< module component check >-----

    log_newline
    log_info("ATMOS_ADMIN_setup",*) 'Atmosphere model components '

    if ( atmos_do ) then
       log_info_cont(*) 'Atmosphere model       : ON'
    else
       log_info_cont(*) 'Atmosphere model       : OFF'
    endif

    log_info_cont(*) 'Dynamics...'

    if    ( atmos_dyn_type == 'OFF' ) then
       log_info_cont(*) '+ Dynamical core       : OFF'
       log_info_cont(*) '+ Advection            : OFF'
       atmos_sw_dyn = .false.
    elseif( atmos_dyn_type == 'NONE' ) then
       ! The advection is disbled
       ! The tendencies calculated by physical processed are added
       log_info_cont(*) '+ Dynamical core       : ON, ', trim(atmos_dyn_type)
       log_info_cont(*) '+ Advection            : OFF'
       atmos_sw_dyn = .true.
    else ! default
       log_info_cont(*) '+ Dynamical core       : ON, ', trim(atmos_dyn_type)
       log_info_cont(*) '+ Advection            : ON'
       atmos_sw_dyn = .true.
    endif

    log_info_cont(*) 'Physics...'

    if ( atmos_phy_mp_type /= 'OFF' .AND. atmos_phy_mp_type /= 'NONE' ) then
       log_info_cont(*) '+ Cloud Microphysics   : ON, ', trim(atmos_phy_mp_type)
       atmos_sw_phy_mp = .true.
    else
       log_info_cont(*) '+ Cloud Microphysics   : OFF'
       atmos_sw_phy_mp = .false.
    endif

    if ( atmos_phy_ae_type /= 'OFF' .AND. atmos_phy_ae_type /= 'NONE' ) then
       log_info_cont(*) '+ Aerosol Microphysics : ON, ', trim(atmos_phy_ae_type)
       atmos_sw_phy_ae = .true.
    else
       log_info_cont(*) '+ Aerosol Microphysics : OFF'
       atmos_sw_phy_ae = .false.
    endif

    if ( atmos_phy_ch_type /= 'OFF' .AND. atmos_phy_ch_type /= 'NONE' ) then
       log_info_cont(*) '+ Chemistry            : ON, ', trim(atmos_phy_ch_type)
       atmos_sw_phy_ch = .true.
    else
       log_info_cont(*) '+ Chemistry            : OFF'
       atmos_sw_phy_ch = .false.
    endif

    if ( atmos_phy_rd_type /= 'OFF' .AND. atmos_phy_rd_type /= 'NONE' ) then
       log_info_cont(*) '+ Radiative transfer   : ON, ', trim(atmos_phy_rd_type)
       atmos_sw_phy_rd = .true.
    else
       log_info_cont(*) '+ Radiative transfer   : OFF'
       atmos_sw_phy_rd = .false.
    endif

    if ( atmos_phy_sf_type /= 'OFF' .AND. atmos_phy_sf_type /= 'NONE' ) then
       log_info_cont(*) '+ Surface Flux         : ON, ', trim(atmos_phy_sf_type)
       atmos_sw_phy_sf = .true.
    else
       log_info_cont(*) '+ Surface Flux         : OFF'
       atmos_sw_phy_sf = .false.
    endif

    if ( atmos_phy_tb_type /= 'OFF' .AND. atmos_phy_tb_type /= 'NONE' ) then
       log_info_cont(*) '+ Sub-grid Turbulence  : ON, ', trim(atmos_phy_tb_type)
       atmos_sw_phy_tb = .true.
    else
       log_info_cont(*) '+ Sub-grid Turbulence  : OFF'
       atmos_sw_phy_tb = .false.
    endif

    if ( atmos_phy_bl_type /= 'OFF' .AND. atmos_phy_bl_type /= 'NONE' ) then
       log_info_cont(*) '+ PBL Turbulence       : ON, ', trim(atmos_phy_bl_type)
       atmos_sw_phy_bl = .true.
    else
       log_info_cont(*) '+ PBL Turbulence       : OFF'
       atmos_sw_phy_bl = .false.
    endif

    if ( atmos_phy_cp_type /= 'OFF' .AND. atmos_phy_cp_type /= 'NONE' ) then
       log_info_cont(*) '+ Convection Param.    : ON, ', trim(atmos_phy_cp_type)
       atmos_sw_phy_cp = .true.
    else
       log_info_cont(*) '+ Convection Param.    : OFF'
       atmos_sw_phy_cp = .false.
    endif

    if ( atmos_use_average ) then
       log_info_cont(*) '+ Use time-averaging value for physics? : YES'
    else
       log_info_cont(*) '+ Use time-averaging value for physics? : NO'
    endif

    return
  end subroutine atmos_admin_setup

  !-----------------------------------------------------------------------------
  subroutine atmos_admin_getscheme( &
       component_name, &
       scheme_name     )
    use scale_prc, only: &
       prc_abort
    implicit none

    character(len=*),       intent(in)  :: component_name
    character(len=H_SHORT), intent(out) :: scheme_name
    !---------------------------------------------------------------------------

    select case(component_name)
    case("DYN")
       scheme_name = atmos_dyn_type
    case("PHY_MP")
       scheme_name = atmos_phy_mp_type
    case("PHY_AE")
       scheme_name = atmos_phy_ae_type
    case("PHY_CH")
       scheme_name = atmos_phy_ch_type
    case("PHY_RD")
       scheme_name = atmos_phy_rd_type
    case("PHY_SF")
       scheme_name = atmos_phy_sf_type
    case("PHY_TB")
       scheme_name = atmos_phy_tb_type
    case("PHY_BL")
       scheme_name = atmos_phy_bl_type
    case("PHY_CP")
       scheme_name = atmos_phy_cp_type
    case default
       log_error("ATMOS_ADMIN_getscheme",*) 'Unsupported component_name. Check!', trim(component_name)
       call prc_abort
    end select

    return
  end subroutine atmos_admin_getscheme

end module mod_atmos_admin