module atmosphere / physics / radiation / offline More...

Functions/Subroutines
subroutine, public	atmos_phy_rd_offline_setup
	Setup. More...

subroutine, public	atmos_phy_rd_offline_flux (KA, KS, KE, IA, IS, IE, JA, JS, JE, time_now, flux_rad, SFLX_rad_dn)
	Radiation main. More...

Detailed Description

module atmosphere / physics / radiation / offline

Description: Atmospheric radiation transfer process for offline usage (input radiation flux from the file)

Author: Team SCALE

NAMELIST

PARAM_ATMOS_PHY_RD_OFFLINE

name	type	default value	comment
ATMOS_PHY_RD_OFFLINE_BASENAME	character(len=H_LONG)	''
ATMOS_PHY_RD_OFFLINE_BASENAME_ADD_NUM	logical	.false.
ATMOS_PHY_RD_OFFLINE_NUMBER_OF_FILES	integer	1
ATMOS_PHY_RD_OFFLINE_AXISTYPE	character(len=H_SHORT)	'XYZ'
ATMOS_PHY_RD_OFFLINE_ENABLE_PERIODIC_YEAR	logical	.false.
ATMOS_PHY_RD_OFFLINE_ENABLE_PERIODIC_MONTH	logical	.false.
ATMOS_PHY_RD_OFFLINE_ENABLE_PERIODIC_DAY	logical	.false.
ATMOS_PHY_RD_OFFLINE_STEP_FIXED	integer	0
ATMOS_PHY_RD_OFFLINE_DEFVAL	real(RP)		= UNDEF
ATMOS_PHY_RD_OFFLINE_CHECK_COORDINATES	logical	.true.
ATMOS_PHY_RD_OFFLINE_STEP_LIMIT	integer	0
ATMOS_PHY_RD_OFFLINE_DIFFUSE_RATE	real(RP)	0.5_RP
ATMOS_PHY_RD_OFFLINE_NIR_RATE	real(RP)	0.5_RP

History Output: No history output

Function/Subroutine Documentation

◆ atmos_phy_rd_offline_setup()

subroutine, public scale_atmos_phy_rd_offline::atmos_phy_rd_offline_setup

Setup.

Definition at line 58 of file scale_atmos_phy_rd_offline.F90.

     use scale_prc, only: &
        prc_abort
     use scale_file_external_input, only: &
        file_external_input_regist
     use scale_const, only: &
        undef => const_undef
     implicit none
  
     character(len=H_SHORT) :: vars_3d   (num_vars_3d)
     character(len=H_SHORT) :: vars_2d   (num_vars_2d)
     character(len=H_SHORT) :: vars_2d_op(num_vars_2d_op)
  
     data vars_3d    / 'RFLX_LW_up', 'RFLX_LW_dn', 'RFLX_SW_up', 'RFLX_SW_dn' /
     data vars_2d    / 'SFLX_LW_up', 'SFLX_LW_dn', 'SFLX_SW_up', 'SFLX_SW_dn' /
     data vars_2d_op / 'SFLX_NIR_dn_dir', 'SFLX_NIR_dn_dif', 'SFLX_VIS_dn_dir', 'SFLX_VIS_dn_dif' /
  
     character(len=H_LONG)  :: ATMOS_PHY_RD_offline_basename              = ''
     logical                :: ATMOS_PHY_RD_offline_basename_add_num      = .false.
     integer                :: ATMOS_PHY_RD_offline_number_of_files       = 1
     character(len=H_SHORT) :: ATMOS_PHY_RD_offline_axistype              = 'XYZ'
     logical                :: ATMOS_PHY_RD_offline_enable_periodic_year  = .false.
     logical                :: ATMOS_PHY_RD_offline_enable_periodic_month = .false.
     logical                :: ATMOS_PHY_RD_offline_enable_periodic_day   = .false.
     integer                :: ATMOS_PHY_RD_offline_step_fixed            = 0
     real(RP)               :: ATMOS_PHY_RD_offline_defval               ! = UNDEF
     logical                :: ATMOS_PHY_RD_offline_check_coordinates     = .true.
     integer                :: ATMOS_PHY_RD_offline_step_limit            = 0
  
     namelist / param_atmos_phy_rd_offline / &
        atmos_phy_rd_offline_basename,              &
        atmos_phy_rd_offline_basename_add_num,      &
        atmos_phy_rd_offline_number_of_files,       &
        atmos_phy_rd_offline_axistype,              &
        atmos_phy_rd_offline_enable_periodic_year,  &
        atmos_phy_rd_offline_enable_periodic_month, &
        atmos_phy_rd_offline_enable_periodic_day,   &
        atmos_phy_rd_offline_step_fixed,            &
        atmos_phy_rd_offline_defval,                &
        atmos_phy_rd_offline_check_coordinates,     &
        atmos_phy_rd_offline_step_limit,            &
        atmos_phy_rd_offline_diffuse_rate,          &
        atmos_phy_rd_offline_nir_rate
  
     integer :: n, ierr
     !---------------------------------------------------------------------------
  
     log_newline
     log_info("ATMOS_PHY_RD_offline_setup",*) 'Setup'
     log_info("ATMOS_PHY_RD_offline_setup",*) 'Offline radiation process'
  
     atmos_phy_rd_offline_defval = undef
  
     !--- read namelist
     rewind(io_fid_conf)
     read(io_fid_conf,nml=param_atmos_phy_rd_offline,iostat=ierr)
     if( ierr < 0 ) then !--- missing
        log_info("ATMOS_PHY_RD_offline_setup",*) 'Not found namelist. Default used.'
     elseif( ierr > 0 ) then !--- fatal error
        log_error("ATMOS_PHY_RD_offline_setup",*) 'Not appropriate names in namelist PARAM_ATMOS_PHY_RD_OFFLINE. Check!'
        call prc_abort
     endif
     log_nml(param_atmos_phy_rd_offline)
  
     if ( atmos_phy_rd_offline_basename == '' ) then
        log_error("ATMOS_PHY_RD_offline_setup",*) 'ATMOS_PHY_RD_offline_basename is necessary'
        call prc_abort
     end if
  
     do n = 1, num_vars_3d
        call file_external_input_regist( atmos_phy_rd_offline_basename,              & ! [IN]
                                         atmos_phy_rd_offline_basename_add_num,      & ! [IN]
                                         atmos_phy_rd_offline_number_of_files,       & ! [IN]
                                         vars_3d(n),                                 & ! [IN]
                                         atmos_phy_rd_offline_axistype,              & ! [IN]
                                         atmos_phy_rd_offline_enable_periodic_year,  & ! [IN]
                                         atmos_phy_rd_offline_enable_periodic_month, & ! [IN]
                                         atmos_phy_rd_offline_enable_periodic_day,   & ! [IN]
                                         atmos_phy_rd_offline_step_fixed,            & ! [IN]
                                         atmos_phy_rd_offline_defval,                & ! [IN]
                                         check_coordinates = atmos_phy_rd_offline_check_coordinates, & ! [IN]
                                         step_limit        = atmos_phy_rd_offline_step_limit         ) ! [IN]
     end do
  
     do n = 1, num_vars_2d
        call file_external_input_regist( atmos_phy_rd_offline_basename,              & ! [IN]
                                         atmos_phy_rd_offline_basename_add_num,      & ! [IN]
                                         atmos_phy_rd_offline_number_of_files,       & ! [IN]
                                         vars_2d(n),                                 & ! [IN]
                                         'XY',                                       & ! [IN]
                                         atmos_phy_rd_offline_enable_periodic_year,  & ! [IN]
                                         atmos_phy_rd_offline_enable_periodic_month, & ! [IN]
                                         atmos_phy_rd_offline_enable_periodic_day,   & ! [IN]
                                         atmos_phy_rd_offline_step_fixed,            & ! [IN]
                                         atmos_phy_rd_offline_defval,                & ! [IN]
                                         check_coordinates = atmos_phy_rd_offline_check_coordinates, & ! [IN]
                                         step_limit        = atmos_phy_rd_offline_step_limit         ) ! [IN]
     end do
  
     do n = 1, num_vars_2d_op
        call file_external_input_regist( atmos_phy_rd_offline_basename,              & ! [IN]
                                         atmos_phy_rd_offline_basename_add_num,      & ! [IN]
                                         atmos_phy_rd_offline_number_of_files,       & ! [IN]
                                         vars_2d_op(n),                              & ! [IN]
                                         'XY',                                       & ! [IN]
                                         atmos_phy_rd_offline_enable_periodic_year,  & ! [IN]
                                         atmos_phy_rd_offline_enable_periodic_month, & ! [IN]
                                         atmos_phy_rd_offline_enable_periodic_day,   & ! [IN]
                                         atmos_phy_rd_offline_step_fixed,            & ! [IN]
                                         atmos_phy_rd_offline_defval,                & ! [IN]
                                         check_coordinates = atmos_phy_rd_offline_check_coordinates, & ! [IN]
                                         step_limit        = atmos_phy_rd_offline_step_limit,        & ! [IN]
                                         exist             = vars_2d_exist(n)                        ) ! [OUT]
        if ( vars_2d_exist(n) ) then
           log_info("ATMOS_PHY_RD_offline_setup",*) '', trim(vars_2d_op(n)), ' found.'
        else
           log_info("ATMOS_PHY_RD_offline_setup",*) '', trim(vars_2d_op(n)), ' not found.'
        end if
     end do
  
     return

References scale_const::const_undef, scale_io::io_fid_conf, and scale_prc::prc_abort().

Referenced by mod_atmos_phy_rd_driver::atmos_phy_rd_driver_setup().

Here is the call graph for this function:

Here is the caller graph for this function:

◆ atmos_phy_rd_offline_flux()

subroutine, public scale_atmos_phy_rd_offline::atmos_phy_rd_offline_flux	(	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(dp), intent(in)	time_now,
		real(rp), dimension (ka,ia,ja,2,2), intent(out)	flux_rad,
		real(rp), dimension(ia,ja,n_rad_dir,n_rad_rgn), intent(out)	SFLX_rad_dn
	)

Radiation main.

Definition at line 187 of file scale_atmos_phy_rd_offline.F90.

     use scale_prc, only: &
        prc_abort
     use scale_file_external_input, only: &
        file_external_input_update
     use scale_atmos_phy_rd_common, only: &
        i_sw, &
        i_lw, &
        i_dn, &
        i_up
     implicit none
  
     integer,  intent(in)  :: KA, KS, KE
     integer,  intent(in)  :: IA, IS, IE
     integer,  intent(in)  :: JA, JS, JE
     real(DP), intent(in)  :: time_now
     real(RP), intent(out) :: flux_rad   (KA,IA,JA,2,2)
     real(RP), intent(out) :: SFLX_rad_dn(IA,JA,N_RAD_DIR,N_RAD_RGN)
  
     real(RP) :: buffer(IA,JA)
     logical  :: error, error_sum, error_sflx
  
     integer  :: i, j
     !---------------------------------------------------------------------------
  
     log_progress(*) 'atmosphere / physics / radiation / offline'
  
     ! [note] external data input is now support only SCALE-RM
  
     error_sum = .false.
  
     ! 3D
     call file_external_input_update( 'RFLX_LW_up', time_now, flux_rad(:,:,:,i_lw,i_up), error )
     error_sum = ( error .OR. error_sum )
  
     call file_external_input_update( 'RFLX_LW_dn', time_now, flux_rad(:,:,:,i_lw,i_dn), error )
     error_sum = ( error .OR. error_sum )
  
     call file_external_input_update( 'RFLX_SW_up', time_now, flux_rad(:,:,:,i_sw,i_up), error )
     error_sum = ( error .OR. error_sum )
  
     call file_external_input_update( 'RFLX_SW_dn', time_now, flux_rad(:,:,:,i_sw,i_dn), error )
     error_sum = ( error .OR. error_sum )
  
  
     ! 2D
     call file_external_input_update( 'SFLX_LW_up', time_now, buffer(:,:), error )
     if ( error ) then
        error_sum = .true.
     else
        !$omp parallel do default(none) OMP_SCHEDULE_ &
        !$omp private(i,j) &
        !$omp shared(KS,IS,IE,JS,JE) &
        !$omp shared(flux_rad,buffer)
        !$acc kernels
        do j = js, je
        do i = is, ie
           flux_rad(ks-1,i,j,i_lw,i_up) = buffer(i,j)
        end do
        end do
        !$acc end kernels
     end if
  
     call file_external_input_update( 'SFLX_LW_dn', time_now, buffer(:,:), error )
     if ( error ) then
        error_sum = .true.
     else
        !$omp parallel do default(none) OMP_SCHEDULE_ &
        !$omp private(i,j) &
        !$omp shared(KS,IS,IE,JS,JE) &
        !$omp shared(flux_rad,buffer)
        !$acc kernels
        do j = js, je
        do i = is, ie
           flux_rad(ks-1,i,j,i_lw,i_dn) = buffer(i,j)
        end do
        end do
        !$acc end kernels
     end if
  
     call file_external_input_update( 'SFLX_SW_up', time_now, buffer(:,:), error )
     if ( error ) then
        error_sum = .true.
     else
        !$omp parallel do default(none) OMP_SCHEDULE_ &
        !$omp private(i,j) &
        !$omp shared(KS,IS,IE,JS,JE) &
        !$omp shared(flux_rad,buffer)
        !$acc kernels
        do j = js, je
        do i = is, ie
           flux_rad(ks-1,i,j,i_sw,i_up) = buffer(i,j)
        end do
        end do
        !$acc end kernels
     end if
  
     call file_external_input_update( 'SFLX_SW_dn', time_now, buffer(:,:), error )
     if ( error ) then
        error_sum = .true.
     else
        !$omp parallel do default(none) OMP_SCHEDULE_ &
        !$omp private(i,j) &
        !$omp shared(KS,IS,IE,JS,JE) &
        !$omp shared(flux_rad,buffer)
        !$acc kernels
        do j = js, je
        do i = is, ie
           flux_rad(ks-1,i,j,i_sw,i_dn) = buffer(i,j)
        end do
        end do
        !$acc end kernels
     end if
  
     !$omp parallel do default(none) OMP_SCHEDULE_ &
     !$omp private(i,j) &
     !$omp shared(KS,IS,IE,JS,JE,ATMOS_PHY_RD_offline_diffuse_rate) &
     !$omp shared(SFLX_rad_dn,flux_rad)
     !$acc kernels
     do j = js, je
     do i = is, ie
        sflx_rad_dn(i,j,i_r_direct ,i_r_ir) = 0.0_rp
        sflx_rad_dn(i,j,i_r_diffuse,i_r_ir) = flux_rad(ks-1,i,j,i_lw,i_dn)
     end do
     end do
     !$acc end kernels
  
     ! 2D optional
  
     error_sflx = .false.
  
     if ( vars_2d_exist(1) ) then
        call file_external_input_update( 'SFLX_NIR_dn_dir', time_now, sflx_rad_dn(:,:,i_r_direct,i_r_vis), error )
        error_sum = ( error .OR. error_sum )
     else
        error_sflx = .true.
     endif
  
     if ( vars_2d_exist(2) ) then
        call file_external_input_update( 'SFLX_NIR_dn_dif', time_now, sflx_rad_dn(:,:,i_r_direct,i_r_vis), error )
        error_sum = ( error .OR. error_sum )
     else
        error_sflx = .true.
     endif
  
     if ( vars_2d_exist(3) ) then
        call file_external_input_update( 'SFLX_VIS_dn_dir', time_now, sflx_rad_dn(:,:,i_r_direct,i_r_vis), error )
        error_sum = ( error .OR. error_sum )
     else
        error_sflx = .true.
     endif
  
     if ( vars_2d_exist(4) ) then
        call file_external_input_update( 'SFLX_VIS_dn_dif', time_now, sflx_rad_dn(:,:,i_r_direct,i_r_vis), error )
        error_sum = ( error .OR. error_sum )
     else
        error_sflx = .true.
     endif
  
     if ( error_sflx ) then ! reconstruct from lowermost SW flux
        !$omp parallel do default(none) OMP_SCHEDULE_ &
        !$omp private(i,j) &
        !$omp shared(KS,IS,IE,JS,JE,ATMOS_PHY_RD_offline_diffuse_rate,ATMOS_PHY_RD_offline_NIR_rate) &
        !$omp shared(SFLX_rad_dn,flux_rad)
        !$acc kernels
        do j = js, je
        do i = is, ie
           sflx_rad_dn(i,j,i_r_direct ,i_r_nir) = ( 1.0_rp-atmos_phy_rd_offline_diffuse_rate ) &
                                                * (        atmos_phy_rd_offline_nir_rate     ) * flux_rad(ks-1,i,j,i_sw,i_dn)
           sflx_rad_dn(i,j,i_r_diffuse,i_r_nir) = (        atmos_phy_rd_offline_diffuse_rate ) &
                                                * (        atmos_phy_rd_offline_nir_rate     ) * flux_rad(ks-1,i,j,i_sw,i_dn)
           sflx_rad_dn(i,j,i_r_direct ,i_r_vis) = ( 1.0_rp-atmos_phy_rd_offline_diffuse_rate ) &
                                                * ( 1.0_rp-atmos_phy_rd_offline_nir_rate     ) * flux_rad(ks-1,i,j,i_sw,i_dn)
           sflx_rad_dn(i,j,i_r_diffuse,i_r_vis) = (        atmos_phy_rd_offline_diffuse_rate ) &
                                                * ( 1.0_rp-atmos_phy_rd_offline_nir_rate     ) * flux_rad(ks-1,i,j,i_sw,i_dn)
        enddo
        enddo
        !$acc end kernels
     endif
  
     if ( error_sum ) then
        log_error("ATMOS_PHY_RD_offline_flux",*) 'Requested data is not found!'
        call prc_abort
     endif
  
     return

References scale_atmos_phy_rd_common::i_dn, scale_atmos_phy_rd_common::i_lw, scale_cpl_sfc_index::i_r_diffuse, scale_cpl_sfc_index::i_r_direct, scale_cpl_sfc_index::i_r_ir, scale_cpl_sfc_index::i_r_nir, scale_cpl_sfc_index::i_r_vis, scale_atmos_phy_rd_common::i_sw, scale_atmos_phy_rd_common::i_up, and scale_prc::prc_abort().

Referenced by mod_atmos_phy_rd_driver::atmos_phy_rd_driver_calc_tendency().