mod_atmos_bnd_driver Module Reference

module ATMOSPHERE / Boundary treatment More...

Functions/Subroutines
subroutine, public	atmos_boundary_driver_setup
	Setup. More...
subroutine, public	atmos_boundary_driver_set
	set More...
subroutine	atmos_boundary_set_file
	Set boundary value for real case experiment. More...
subroutine	atmos_boundary_set_online
	Set boundary value for real case experiment [online daughter]. More...
subroutine, public	atmos_boundary_driver_finalize
	Finalize boundary value. More...
subroutine, public	atmos_boundary_driver_update (last_step)
	Update boundary value with a constant time boundary. More...
subroutine, public	atmos_boundary_driver_send
	Send data to child domain. More...
subroutine	set_boundary (use_velz)
subroutine	update_ref_index
	Update indices of array of boundary references. More...
subroutine	calc_mass (ref)
subroutine	set_offset

Variables
integer, public	bnd_qa
integer, allocatable, public	of
integer, allocatable, public	tracer
integer, public	at
integer, public	boundary
integer, dimension(:), allocatable, public	bnd_iq
integer, allocatable, public	index
real(rp), dimension(:,:,:), allocatable, public	atmos_boundary_dens
real(rp), allocatable, public	reference
real(rp), dimension(0-1), allocatable, public	dens
real(rp), dimension(:,:,:), allocatable, public	atmos_boundary_velz
real(rp), dimension(0-1), allocatable, public	velz
real(rp), dimension(:,:,:), allocatable, public	atmos_boundary_velx
real(rp), dimension(0-1), allocatable, public	velx
real(rp), dimension(:,:,:), allocatable, public	atmos_boundary_vely
real(rp), dimension(0-1), allocatable, public	vely
real(rp), dimension(:,:,:), allocatable, public	atmos_boundary_pott
real(rp), dimension(0-1), allocatable, public	pott
real(rp), dimension(:,:,:,:), allocatable, public	atmos_boundary_qtrc
real(rp), dimension(0-1), allocatable, public	qtrc
real(rp), dimension(:,:,:), allocatable, public	atmos_boundary_alpha_dens
real(rp), allocatable, public	damping
real(rp), allocatable, public	coefficient
real(rp), allocatable, public	for
real(rp), dimension(:,:,:), allocatable, public	atmos_boundary_alpha_velz
real(rp), dimension(:,:,:), allocatable, public	atmos_boundary_alpha_velx
real(rp), dimension(:,:,:), allocatable, public	atmos_boundary_alpha_vely
real(rp), dimension(:,:,:), allocatable, public	atmos_boundary_alpha_pott
real(rp), dimension(:,:,:,:), allocatable, public	atmos_boundary_alpha_qtrc
real(rp), dimension(:,:,:), allocatable, public	atmos_boundary_mflux_offset_x
real(rp), allocatable, public	mass
real(rp), allocatable, public	flux
real(rp), dimension(south, north), allocatable, public	offset
real(rp), dimension(:,:,:), allocatable, public	atmos_boundary_mflux_offset_y
real(rp), public	atmos_boundary_smoother_fact = 0.2_RP
real(rp), public	fact
real(rp), public	smoother
real(rp), public	to
logical, public	atmos_boundary_update_flag = .false.
logical, public	switch
logical, public	real
logical, public	case

Detailed Description

module ATMOSPHERE / Boundary treatment

Description

Boundary treatment of model domain Additional forcing, Sponge layer, rayleigh dumping

Author

Team SCALE

NAMELIST

• PARAM_ATMOS_BOUNDARY

  name	type	default value	comment
  ATMOS_BOUNDARY_TYPE	character(len=H_SHORT)	'NONE'
  ATMOS_BOUNDARY_IN_BASENAME	character(len=H_LONG)	''
  ATMOS_BOUNDARY_IN_CHECK_COORDINATES	logical	.true.
  ATMOS_BOUNDARY_IN_AGGREGATE	logical
  ATMOS_BOUNDARY_OUT_BASENAME	character(len=H_LONG)	''
  ATMOS_BOUNDARY_OUT_TITLE	character(len=H_MID)	'SCALE-RM BOUNDARY CONDITION'	title of the output file
  ATMOS_BOUNDARY_OUT_DTYPE	character(len=H_SHORT)	'DEFAULT'	REAL4 or REAL8
  ATMOS_BOUNDARY_OUT_AGGREGATE	logical
  ATMOS_BOUNDARY_USE_VELZ	logical	.false.	read from file?
  ATMOS_BOUNDARY_USE_VELX	logical	.false.	read from file?
  ATMOS_BOUNDARY_USE_VELY	logical	.false.	read from file?
  ATMOS_BOUNDARY_USE_PT	logical	.false.	read from file?
  ATMOS_BOUNDARY_USE_DENS	logical	.false.	read from file?
  ATMOS_BOUNDARY_USE_QV	logical	.false.	read from file?
  ATMOS_BOUNDARY_USE_QHYD	logical	.false.	read from file?
  ATMOS_BOUNDARY_USE_CHEM	logical	.false.	read from file?
  ATMOS_BOUNDARY_DENS_ADJUST	logical	.false.
  ATMOS_BOUNDARY_DENS_ADJUST_TAU	real(RP)	-1.0_RP
  ATMOS_BOUNDARY_VALUE_VELZ	real(RP)	0.0_RP	velocity w at boundary, 0 [m/s]
  ATMOS_BOUNDARY_VALUE_VELX	real(RP)	0.0_RP	velocity u at boundary, 0 [m/s]
  ATMOS_BOUNDARY_VALUE_VELY	real(RP)	0.0_RP	velocity v at boundary, 0 [m/s]
  ATMOS_BOUNDARY_VALUE_PT	real(RP)	300.0_RP	potential temp. at boundary, 300 [K]
  ATMOS_BOUNDARY_VALUE_QTRC	real(RP)	0.0_RP	tracer at boundary, 0 [kg/kg]
  ATMOS_BOUNDARY_ALPHAFACT_DENS	real(RP)	1.0_RP	alpha factor again default
  ATMOS_BOUNDARY_ALPHAFACT_VELZ	real(RP)	1.0_RP	alpha factor again default
  ATMOS_BOUNDARY_ALPHAFACT_VELX	real(RP)	1.0_RP	alpha factor again default
  ATMOS_BOUNDARY_ALPHAFACT_VELY	real(RP)	1.0_RP	alpha factor again default
  ATMOS_BOUNDARY_ALPHAFACT_PT	real(RP)	1.0_RP	alpha factor again default
  ATMOS_BOUNDARY_ALPHAFACT_QTRC	real(RP)	1.0_RP	alpha factor again default
  ATMOS_BOUNDARY_SMOOTHER_FACT	real(RP)	0.2_RP	> fact for smoother to damping
  ATMOS_BOUNDARY_FRACZ	real(RP)	1.0_RP	fraction of boundary region for dumping (z) (0-1)
  ATMOS_BOUNDARY_FRACX	real(RP)	1.0_RP	fraction of boundary region for dumping (x) (0-1)
  ATMOS_BOUNDARY_FRACY	real(RP)	1.0_RP	fraction of boundary region for dumping (y) (0-1)
  ATMOS_BOUNDARY_TAUZ	real(RP)		maximum value for damping tau (z) [s]
  ATMOS_BOUNDARY_TAUX	real(RP)		maximum value for damping tau (x) [s]
  ATMOS_BOUNDARY_TAUY	real(RP)		maximum value for damping tau (y) [s]
  ATMOS_BOUNDARY_UPDATE_DT	real(DP)	0.0_DP	inteval time of boudary data update [s]
  ATMOS_BOUNDARY_START_DATE	integer, dimension(6)	(/ -9999, 0, 0, 0, 0, 0 /)	boundary initial date
  ATMOS_BOUNDARY_LINEAR_V	logical	.false.	linear or non-linear profile of relax region
  ATMOS_BOUNDARY_LINEAR_H	logical	.true.	linear or non-linear profile of relax region
  ATMOS_BOUNDARY_EXP_H	real(RP)	2.0_RP	factor of non-linear profile of relax region
  ATMOS_BOUNDARY_INTERP_TYPE	character(len=H_SHORT)	'lerp_initpoint'	type of boundary interporation
  ATMOS_GRID_NUDGING_UV	logical	.false.	grid nudging
  ATMOS_GRID_NUDGING_PT	logical	.false.	grid nudging
  ATMOS_GRID_NUDGING_QV	logical	.false.	grid nudging
  ATMOS_GRID_NUDGING_TAU	real(RP)		Damping tau for grid nudging [s]

  
  

History Output

name	description	unit	variable
DENS_BND	Boundary Density	kg/m3	ATMOS_BOUNDARY_DENS
PT_BND	Boundary potential temperature	K	ATMOS_BOUNDARY_POTT
VELX_BND	Boundary velocity x-direction	m/s	ATMOS_BOUNDARY_VELX
VELY_BND	Boundary velocity y-direction	m/s	ATMOS_BOUNDARY_VELY
VELZ_BND	Boundary velocity z-direction	m/s	ATMOS_BOUNDARY_VELZ
{TRACER_NAME}_BND	{TRACER_NAME} in boundary; {TRACER_NAME} depends on the physics schemes, e.g., QV, QC, QR.	kg/kg	ATMOS_BOUNDARY_QTRC

Function/Subroutine Documentation

atmos_boundary_driver_setup()

subroutine, public mod_atmos_bnd_driver::atmos_boundary_driver_setup

Setup.

Definition at line 223 of file mod_atmos_bnd_driver.F90.

    use scale_prc, only: &
       prc_abort
    use scale_const, only: &
       undef => const_undef
    use scale_time, only: &
       dt => time_dtsec
    use scale_file, only: &
       file_aggregate
    use scale_comm_cartesc_nest, only: &
       online_boundary_diagqhyd, &
       online_boundary_use_qhyd, &
       use_nesting,         &
       offline,             &
       online_iam_parent,   &
       online_iam_daughter, &
       nestqa => comm_cartesc_nest_bnd_qa
    use scale_atmos_hydrometeor, only: &
       atmos_hydrometeor_dry, &
       i_qv
    use mod_atmos_phy_mp_vars, only: &
       qs_mp, &
       qe_mp
    use mod_atmos_phy_ch_vars, only: &
       qs_ch, &
       qe_ch
    use scale_atmos_grid_cartesc_real, only: &
       atmos_grid_cartesc_real_areazuy_x
    implicit none
 
    namelist / param_atmos_boundary / &
       atmos_boundary_type,           &
       atmos_boundary_in_basename,    &
       atmos_boundary_in_check_coordinates, &
       atmos_boundary_in_aggregate,   &
       atmos_boundary_out_basename,   &
       atmos_boundary_out_title,      &
       atmos_boundary_out_dtype,      &
       atmos_boundary_out_aggregate,  &
       atmos_boundary_use_velz,       &
       atmos_boundary_use_velx,       &
       atmos_boundary_use_vely,       &
       atmos_boundary_use_pt,         &
       atmos_boundary_use_dens,       &
       atmos_boundary_use_qv,         &
       atmos_boundary_use_qhyd,       &
       atmos_boundary_use_chem,       &
       atmos_boundary_dens_adjust,    &
       atmos_boundary_dens_adjust_tau, &
       atmos_boundary_value_velz,     &
       atmos_boundary_value_velx,     &
       atmos_boundary_value_vely,     &
       atmos_boundary_value_pt,       &
       atmos_boundary_value_qtrc,     &
       atmos_boundary_alphafact_dens, &
       atmos_boundary_alphafact_velz, &
       atmos_boundary_alphafact_velx, &
       atmos_boundary_alphafact_vely, &
       atmos_boundary_alphafact_pt,   &
       atmos_boundary_alphafact_qtrc, &
       atmos_boundary_smoother_fact,  &
       atmos_boundary_fracz,          &
       atmos_boundary_fracx,          &
       atmos_boundary_fracy,          &
       atmos_boundary_tauz,           &
       atmos_boundary_taux,           &
       atmos_boundary_tauy,           &
       atmos_boundary_update_dt,      &
       atmos_boundary_start_date,     &
       atmos_boundary_linear_v,       &
       atmos_boundary_linear_h,       &
       atmos_boundary_exp_h,          &
       atmos_boundary_interp_type,    &
       atmos_grid_nudging_uv,         &
       atmos_grid_nudging_pt,         &
       atmos_grid_nudging_qv,         &
       atmos_grid_nudging_tau
 
    integer :: k, i, j, iq
    integer :: ierr
    !---------------------------------------------------------------------------
 
    log_newline
    log_info("ATMOS_BOUNDARY_setup",*) 'Setup'
 
 
    atmos_boundary_tauz = dt * 10.0_rp
    atmos_boundary_taux = dt * 10.0_rp
    atmos_boundary_tauy = dt * 10.0_rp
 
    atmos_boundary_in_aggregate  = file_aggregate
    atmos_boundary_out_aggregate = file_aggregate
 
    atmos_grid_nudging_tau = 10.0_rp * 24.0_rp * 3600.0_rp   ! 10days [s]
 
    !--- read namelist
    rewind(io_fid_conf)
    read(io_fid_conf,nml=param_atmos_boundary,iostat=ierr)
    if( ierr < 0 ) then !--- missing
       log_info("ATMOS_BOUNDARY_setup",*) 'Not found namelist. Default used.'
    elseif( ierr > 0 ) then !--- fatal error
       log_error("ATMOS_BOUNDARY_setup",*) 'Not appropriate names in namelist PARAM_ATMOS_BOUNDARY. Check!'
       call prc_abort
    endif
    log_nml(param_atmos_boundary)
 
    ! setting switches
    if( .NOT. use_nesting ) then
       atmos_boundary_online = .false.
    else
       if( offline ) then
          atmos_boundary_online = .false.
       else
          atmos_boundary_online = .true.
       endif
    endif
    do_parent_process   = .false.
    do_daughter_process = .false.
    atmos_boundary_online_master = .false.
    if ( atmos_boundary_online ) then
       if ( online_iam_parent ) then
          do_parent_process = .true.
          if ( .NOT. online_iam_daughter ) then
             atmos_boundary_online_master = .true.
          endif
       endif
       if ( online_iam_daughter ) then
          do_daughter_process = .true.
          atmos_boundary_use_qhyd = online_boundary_use_qhyd
       endif
    endif
 
    allocate( bnd_iq(qa) )
    bnd_iq(:) = -1
    bnd_qa = 0
    if ( .not. atmos_hydrometeor_dry ) then
       bnd_qa = bnd_qa + 1
       bnd_iq(i_qv) = bnd_qa
       if( atmos_boundary_use_qhyd ) then
          do iq = qs_mp+1, qe_mp
             bnd_qa = bnd_qa + 1
             bnd_iq(iq) = bnd_qa
          end do
       end if
    end if
    if( atmos_boundary_use_chem ) then
       do iq = qs_ch, qe_ch
          bnd_qa = bnd_qa + 1
          bnd_iq(iq) = bnd_qa
       end do
    endif
 
    if ( atmos_boundary_dens_adjust_tau <= 0.0_rp ) then
       atmos_boundary_dens_adjust_tau = max( real(atmos_boundary_update_dt,kind=rp) / 6.0_rp, &
                                             atmos_boundary_taux, atmos_boundary_tauy )
    end if
 
 
 
    allocate( atmos_boundary_dens(ka,ia,ja) )
    allocate( atmos_boundary_velz(ka,ia,ja) )
    allocate( atmos_boundary_velx(ka,ia,ja) )
    allocate( atmos_boundary_vely(ka,ia,ja) )
    allocate( atmos_boundary_pott(ka,ia,ja) )
    allocate( atmos_boundary_qtrc(ka,ia,ja,bnd_qa) )
    atmos_boundary_dens(:,:,:)   = undef
    atmos_boundary_velz(:,:,:)   = undef
    atmos_boundary_velx(:,:,:)   = undef
    atmos_boundary_vely(:,:,:)   = undef
    atmos_boundary_pott(:,:,:)   = undef
    atmos_boundary_qtrc(:,:,:,:) = undef
 
    allocate( atmos_boundary_alpha_dens(ka,ia,ja) )
    allocate( atmos_boundary_alpha_velz(ka,ia,ja) )
    allocate( atmos_boundary_alpha_velx(ka,ia,ja) )
    allocate( atmos_boundary_alpha_vely(ka,ia,ja) )
    allocate( atmos_boundary_alpha_pott(ka,ia,ja) )
    allocate( atmos_boundary_alpha_qtrc(ka,ia,ja,bnd_qa) )
    atmos_boundary_alpha_dens(:,:,:)   = 0.0_rp
    atmos_boundary_alpha_velz(:,:,:)   = 0.0_rp
    atmos_boundary_alpha_velx(:,:,:)   = 0.0_rp
    atmos_boundary_alpha_vely(:,:,:)   = 0.0_rp
    atmos_boundary_alpha_pott(:,:,:)   = 0.0_rp
    atmos_boundary_alpha_qtrc(:,:,:,:) = 0.0_rp
 
    allocate( atmos_boundary_mflux_offset_x(ka,ja,2) )
    allocate( atmos_boundary_mflux_offset_y(ka,ia,2) )
    atmos_boundary_mflux_offset_x(:,:,:) = 0.0_rp
    atmos_boundary_mflux_offset_y(:,:,:) = 0.0_rp
 
    if ( atmos_boundary_type == 'REAL' .OR. do_daughter_process ) then
       l_bnd = .true.
    else
       l_bnd = .false.
    end if
 
    if ( l_bnd ) then
 
       select case(atmos_boundary_interp_type)
       case('same_parent')
          get_boundary => get_boundary_same_parent
       case('nearest_neighbor')
          get_boundary => get_boundary_nearest_neighbor
       case('lerp_initpoint')
          get_boundary => get_boundary_lerp_initpoint
       case('lerp_midpoint')
          get_boundary => get_boundary_lerp_midpoint
       case default
          log_error("ATMOS_BOUNDARY_setup",*) 'Wrong parameter in ATMOS_BOUNDARY_interp_TYPE. Check!'
          call prc_abort
       end select
 
       allocate( atmos_boundary_ref_dens(ka,ia,ja,ref_size) )
       allocate( atmos_boundary_ref_velz(ka,ia,ja,ref_size) )
       allocate( atmos_boundary_ref_velx(ka,ia,ja,ref_size) )
       allocate( atmos_boundary_ref_vely(ka,ia,ja,ref_size) )
       allocate( atmos_boundary_ref_pott(ka,ia,ja,ref_size) )
       allocate( atmos_boundary_ref_qtrc(ka,ia,ja,bnd_qa,ref_size) )
       atmos_boundary_ref_dens(:,:,:,:)   = undef
       atmos_boundary_ref_velz(:,:,:,:)   = undef
       atmos_boundary_ref_velx(:,:,:,:)   = undef
       atmos_boundary_ref_vely(:,:,:,:)   = undef
       atmos_boundary_ref_pott(:,:,:,:)   = undef
       atmos_boundary_ref_qtrc(:,:,:,:,:) = undef
 
       ! initialize boundary value (reading file or waiting parent domain)
       if ( do_daughter_process ) then
          call atmos_boundary_initialize_online
       else
          if ( atmos_boundary_in_basename /= '' ) then
             call atmos_boundary_initialize_file
          else
             log_error("ATMOS_BOUNDARY_setup",*) 'You need specify ATMOS_BOUNDARY_IN_BASENAME'
             call prc_abort
          endif
       endif
 
       call atmos_boundary_setalpha
 
       atmos_boundary_update_flag = .true.
 
       ! for mass flux offset
       allocate( areazuy_w(ka,ja), areazuy_e(ka,ja) )
       allocate( mflux_offset_x(ka,ja,2,2), mflux_offset_y(ka,ia,2,2) )
 
       !$omp parallel do
       do j = js, je
       do k = ks, ke
          areazuy_w(k,j) = atmos_grid_cartesc_real_areazuy_x(k,is-1,j)
          areazuy_e(k,j) = atmos_grid_cartesc_real_areazuy_x(k,ie  ,j)
       end do
       end do
       !$omp parallel do
       do j = js, je
       do k = ks, ke
          mflux_offset_x(k,j,:,:) = 0.0_rp
       end do
       end do
       !$omp parallel do
       do i = is, ie
       do k = ks, ke
          mflux_offset_y(k,i,:,:) = 0.0_rp
       end do
       end do
 
    elseif ( atmos_boundary_type == 'NONE' ) then
 
       atmos_boundary_update_flag = .false.
 
    elseif ( atmos_boundary_type == 'CONST' ) then
 
       call atmos_boundary_setalpha
 
       atmos_boundary_update_flag = .false.
 
    elseif ( atmos_boundary_type == 'INIT' ) then
 
       call atmos_boundary_setalpha
 
       atmos_boundary_update_flag = .false.
 
    elseif ( atmos_boundary_type == 'OFFLINE' ) then
 
       if ( atmos_boundary_in_basename /= '' ) then
          call atmos_boundary_read
       else
          log_error("ATMOS_BOUNDARY_setup",*) 'You need specify ATMOS_BOUNDARY_IN_BASENAME'
          call prc_abort
       endif
 
       atmos_boundary_update_flag = .false.
 
    else
       log_error("ATMOS_BOUNDARY_setup",*) 'unsupported ATMOS_BOUNDARY_TYPE. Check!', trim(atmos_boundary_type)
       call prc_abort
    endif
 
    if ( use_nesting ) atmos_boundary_update_flag = .true.
 
 
    !----- report data -----
    log_newline
    log_info("ATMOS_BOUNDARY_setup",*) 'Atmospheric boundary parameters '
    log_info_cont(*) 'Atmospheric boundary type                      : ', atmos_boundary_type
    log_newline
    log_info_cont(*) 'Is VELZ used in atmospheric boundary?          : ', atmos_boundary_use_velz
    log_info_cont(*) 'Is VELX used in atmospheric boundary?          : ', atmos_boundary_use_velx
    log_info_cont(*) 'Is VELY used in atmospheric boundary?          : ', atmos_boundary_use_vely
    log_info_cont(*) 'Is PT used in atmospheric boundary?            : ', atmos_boundary_use_pt
    log_info_cont(*) 'Is DENS used in atmospheric boundary?          : ', atmos_boundary_use_dens
    log_info_cont(*) 'Is QV   used in atmospheric boundary?          : ', atmos_boundary_use_qv
    log_info_cont(*) 'Is QHYD used in atmospheric boundary?          : ', atmos_boundary_use_qhyd
    log_info_cont(*) 'Is CHEM used in atmospheric boundary?          : ', atmos_boundary_use_chem
    log_newline
    log_info_cont(*) 'Atmospheric boundary VELZ values               : ', atmos_boundary_value_velz
    log_info_cont(*) 'Atmospheric boundary VELX values               : ', atmos_boundary_value_velx
    log_info_cont(*) 'Atmospheric boundary VELY values               : ', atmos_boundary_value_vely
    log_info_cont(*) 'Atmospheric boundary PT values                 : ', atmos_boundary_value_pt
    log_info_cont(*) 'Atmospheric boundary QTRC values               : ', atmos_boundary_value_qtrc
    log_newline
    log_info_cont(*) 'Atmospheric boundary smoother factor           : ', atmos_boundary_smoother_fact
    log_info_cont(*) 'Atmospheric boundary z-fraction                : ', atmos_boundary_fracz
    log_info_cont(*) 'Atmospheric boundary x-fraction                : ', atmos_boundary_fracx
    log_info_cont(*) 'Atmospheric boundary y-fraction                : ', atmos_boundary_fracy
    log_info_cont(*) 'Atmospheric boundary z-relaxation time         : ', atmos_boundary_tauz
    log_info_cont(*) 'Atmospheric boundary x-relaxation time         : ', atmos_boundary_taux
    log_info_cont(*) 'Atmospheric boundary y-relaxation time         : ', atmos_boundary_tauy
    log_newline
    log_info_cont(*) 'Atmospheric boundary update dt                 : ', atmos_boundary_update_dt
    log_info_cont(*) 'Atmospheric boundary start date                : ', atmos_boundary_start_date(:)
    log_newline
    log_info_cont(*) 'Linear profile in vertically relax region      : ', atmos_boundary_linear_v
    log_info_cont(*) 'Linear profile in horizontally relax region    : ', atmos_boundary_linear_h
    log_info_cont(*) 'Non-linear factor in horizontally relax region : ', atmos_boundary_exp_h
    log_newline
    log_info_cont(*) 'Online nesting for lateral boundary            : ', atmos_boundary_online
 
    log_info_cont(*) 'Does lateral boundary exist in this domain?    : ', l_bnd
    if ( l_bnd ) then
       log_info_cont(*) 'Lateral boundary interporation type            : ', trim(atmos_boundary_interp_type)
    endif
    log_newline
    log_info_cont(*) 'Is grid nudging used for VELX & VELY?          : ', atmos_grid_nudging_uv
    log_info_cont(*) 'Is grid nudging used for POTT?                 : ', atmos_grid_nudging_pt
    log_info_cont(*) 'Is grid nudging used for QV?                   : ', atmos_grid_nudging_qv
    log_info_cont(*) 'Relaxation time for grid nudging               : ', atmos_grid_nudging_tau
 
    log_info_cont(*) 'Density adjustment                             : ', atmos_boundary_dens_adjust
    if ( atmos_boundary_dens_adjust ) then
       log_info_cont(*) 'Density relaxation time                        : ', atmos_boundary_dens_adjust_tau
    end if
 
    if ( online_boundary_diagqhyd ) then
       allocate( q_work(ka,ia,ja,nestqa) )
    end if
 
    allocate( zero_x(ka,ja), zero_y(ka,ia) )
    !$omp parallel do
    do j = js, je
    do k = ks, ke
       zero_x(k,j) = 0.0_rp
    end do
    end do
    !$omp parallel do
    do i = is, ie
    do k = ks, ke
       zero_y(k,i) = 0.0_rp
    end do
    end do
 
 
 
    return

References atmos_boundary_alpha_dens, atmos_boundary_alpha_pott, atmos_boundary_alpha_qtrc, atmos_boundary_alpha_velx, atmos_boundary_alpha_vely, atmos_boundary_alpha_velz, atmos_boundary_dens, atmos_boundary_mflux_offset_x, atmos_boundary_mflux_offset_y, atmos_boundary_pott, atmos_boundary_qtrc, atmos_boundary_smoother_fact, atmos_boundary_update_flag, atmos_boundary_velx, atmos_boundary_vely, atmos_boundary_velz, scale_atmos_grid_cartesc_real::atmos_grid_cartesc_real_areazuy_x, scale_atmos_hydrometeor::atmos_hydrometeor_dry, bnd_iq, bnd_qa, scale_comm_cartesc_nest::comm_cartesc_nest_bnd_qa, scale_const::const_undef, scale_file::file_aggregate, scale_atmos_hydrometeor::i_qv, scale_atmos_grid_cartesc_index::ia, scale_atmos_grid_cartesc_index::ie, scale_io::io_fid_conf, 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_tracer::k, scale_atmos_grid_cartesc_index::ka, scale_atmos_grid_cartesc_index::ke, scale_atmos_grid_cartesc_index::ks, scale_comm_cartesc_nest::offline, scale_comm_cartesc_nest::online_boundary_diagqhyd, scale_comm_cartesc_nest::online_boundary_use_qhyd, scale_comm_cartesc_nest::online_iam_daughter, scale_comm_cartesc_nest::online_iam_parent, scale_prc::prc_abort(), scale_tracer::qa, mod_atmos_phy_ch_vars::qe_ch, mod_atmos_phy_mp_vars::qe_mp, mod_atmos_phy_ch_vars::qs_ch, mod_atmos_phy_mp_vars::qs_mp, real, scale_precision::rp, scale_time::time_dtsec, and scale_comm_cartesc_nest::use_nesting.

Referenced by mod_atmos_driver::atmos_driver_setup().

Here is the call graph for this function:

Here is the caller graph for this function:

atmos_boundary_driver_set()

subroutine, public mod_atmos_bnd_driver::atmos_boundary_driver_set

set

Definition at line 600 of file mod_atmos_bnd_driver.F90.

    use mod_atmos_vars, only: &
       dens, &
       momz, &
       momx, &
       momy, &
       rhot, &
       qtrc, &
       qv,   &
       qe
    use mod_atmos_phy_mp_vars, only: &
       qs_mp, &
       qe_mp
    implicit none
 
    if ( do_parent_process ) then !online [parent]
       call atmos_boundary_firstsend( &
            dens, momz, momx, momy, rhot, qtrc(:,:,:,qs_mp:qe_mp), qv, qe )
    end if
 
    if ( l_bnd ) then
 
       ! initialize boundary value (reading file or waiting parent domain)
       if ( do_daughter_process ) then
          call atmos_boundary_set_online
       else
          if ( atmos_boundary_in_basename /= '' ) then
             call atmos_boundary_set_file
          endif
       endif
 
    elseif ( atmos_boundary_type == 'CONST' ) then
 
       call atmos_boundary_generate
 
    elseif ( atmos_boundary_type == 'INIT' ) then
 
       call atmos_boundary_setinitval( dens, & ! [IN]
                                       momz, & ! [IN]
                                       momx, & ! [IN]
                                       momy, & ! [IN]
                                       rhot, & ! [IN]
                                       qtrc  ) ! [IN]
    endif
 
    if( atmos_boundary_out_basename /= '' ) then
       call atmos_boundary_write
    endif
 
    if ( atmos_boundary_update_flag ) then
 
       call history_bnd( &
            atmos_boundary_dens, &
            atmos_boundary_velz, &
            atmos_boundary_velx, &
            atmos_boundary_vely, &
            atmos_boundary_pott, &
            atmos_boundary_qtrc )
    end if
 
    return

References atmos_boundary_alpha_dens, atmos_boundary_alpha_pott, atmos_boundary_alpha_qtrc, atmos_boundary_alpha_velx, atmos_boundary_alpha_vely, atmos_boundary_alpha_velz, atmos_boundary_dens, atmos_boundary_pott, atmos_boundary_qtrc, atmos_boundary_set_file(), atmos_boundary_set_online(), atmos_boundary_update_flag, atmos_boundary_velx, atmos_boundary_vely, atmos_boundary_velz, scale_atmos_grid_cartesc::atmos_grid_cartesc_cbfx, scale_atmos_grid_cartesc::atmos_grid_cartesc_cbfy, scale_atmos_grid_cartesc::atmos_grid_cartesc_cbfz, scale_atmos_grid_cartesc::atmos_grid_cartesc_fbfx, scale_atmos_grid_cartesc::atmos_grid_cartesc_fbfy, scale_atmos_grid_cartesc::atmos_grid_cartesc_fbfz, scale_atmos_refstate::atmos_refstate_dens, bnd_iq, bnd_qa, scale_calendar::calendar_combine_daysec(), scale_calendar::calendar_date2char(), scale_calendar::calendar_date2daysec(), scale_const::const_eps, scale_const::const_pi, dens, scale_file_cartesc::file_cartesc_close(), scale_file_cartesc::file_cartesc_create(), scale_file_cartesc::file_cartesc_def_var(), scale_file_cartesc::file_cartesc_enddef(), scale_file_cartesc::file_cartesc_open(), scale_atmos_hydrometeor::i_qv, 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_tracer::k, scale_atmos_grid_cartesc_index::ka, scale_atmos_grid_cartesc_index::ke, scale_atmos_grid_cartesc_index::ks, mod_atmos_vars::momx, mod_atmos_vars::momy, mod_atmos_vars::momz, scale_comm_cartesc_nest::online_use_velz, scale_prc::prc_abort(), scale_tracer::qa, mod_atmos_vars::qe, mod_atmos_phy_mp_vars::qe_mp, mod_atmos_phy_mp_vars::qs_mp, qtrc, mod_atmos_vars::qv, mod_atmos_vars::rhot, scale_time::time_nowdate, scale_tracer::tracer_name, and scale_tracer::tracer_unit.

Referenced by mod_rm_driver::restart_read().

Here is the call graph for this function:

Here is the caller graph for this function:

atmos_boundary_set_file()

subroutine mod_atmos_bnd_driver::atmos_boundary_set_file

Set boundary value for real case experiment.

Definition at line 1412 of file mod_atmos_bnd_driver.F90.

    use scale_prc, only: &
       prc_abort
    use scale_const, only: &
       pi => const_pi
    use scale_time, only: &
       time_nowdate,      &
       time_dtsec
    use scale_calendar, only: &
       calendar_date2daysec,    &
       calendar_combine_daysec
    implicit none
 
    real(RP) :: bnd_DENS(KA,IA,JA)        ! damping coefficient for DENS (0-1)
    real(RP) :: bnd_VELZ(KA,IA,JA)        ! damping coefficient for VELZ (0-1)
    real(RP) :: bnd_VELX(KA,IA,JA)        ! damping coefficient for VELX (0-1)
    real(RP) :: bnd_VELY(KA,IA,JA)        ! damping coefficient for VELY (0-1)
    real(RP) :: bnd_POTT(KA,IA,JA)        ! damping coefficient for POTT (0-1)
    real(RP) :: bnd_QTRC(KA,IA,JA,BND_QA) ! damping coefficient for QTRC (0-1)
 
    integer  :: run_time_startdate(6)
    integer  :: run_time_startday
    real(DP) :: run_time_startsec
    real(DP) :: run_time_startms
    integer  :: run_time_offset_year
    real(DP) :: run_time_nowdaysec
 
    real(DP) :: boundary_diff_daysec
    real(RP) :: boundary_inc_offset
    integer  :: fillgaps_steps
 
    real(RP) :: total
 
    integer  :: i, j, k, n, iq
    !---------------------------------------------------------------------------
 
    if ( atmos_boundary_update_dt <= 0.0_dp ) then
       log_error("ATMOS_BOUNDARY_set_file",*) 'You need specify ATMOS_BOUNDARY_UPDATE_DT as larger than 0.0'
       call prc_abort
    endif
    update_nstep = nint( atmos_boundary_update_dt / time_dtsec )
    if ( abs(update_nstep * time_dtsec - atmos_boundary_update_dt) > 1e-10_dp ) then
       log_error("ATMOS_BOUNDARY_set_file",*) 'ATMOS_BOUNDARY_UPDATE_DT is not multiple of DT'
       call prc_abort
    end if
 
    !--- recalculate time of the run [no offset]
    run_time_startdate(:) = time_nowdate(:)
    run_time_startms      = 0.0_dp
    run_time_offset_year  = 0
 
    call calendar_date2daysec( run_time_startday,     & ! [OUT]
                               run_time_startsec,     & ! [OUT]
                               run_time_startdate(:), & ! [IN]
                               run_time_startms,      & ! [IN]
                               run_time_offset_year   ) ! [IN]
 
    run_time_nowdaysec = calendar_combine_daysec( run_time_startday, run_time_startsec )
 
    boundary_diff_daysec = run_time_nowdaysec - boundary_time_initdaysec
    boundary_timestep    = 1 + int( boundary_diff_daysec / atmos_boundary_update_dt )
    boundary_inc_offset  = mod( boundary_diff_daysec, atmos_boundary_update_dt )
    fillgaps_steps       = int( boundary_inc_offset / time_dtsec )
 
    log_info("ATMOS_BOUNDARY_set_file",*) 'BOUNDARY TIMESTEP NUMBER FOR INIT:', boundary_timestep
    log_info("ATMOS_BOUNDARY_set_file",*) 'BOUNDARY OFFSET:', boundary_inc_offset
    log_info("ATMOS_BOUNDARY_set_file",*) 'BOUNDARY FILLGAPS STEPS:', fillgaps_steps
 
 
    if ( atmos_boundary_dens_adjust ) then
       allocate( offset_time_fact(0:update_nstep) )
       total = 0.0_rp
       !$omp parallel do reduction(+:total)
       do n = 0, update_nstep
          offset_time_fact(n) = 1.0_rp - cos( 2.0_rp * pi * ( n - 1 ) / update_nstep )
          total = total + offset_time_fact(n)
       end do
       total = total / update_nstep
       !$omp parallel do
       do n = 0, update_nstep
          offset_time_fact(n) = offset_time_fact(n) / total
       end do
    end if
 
    ! read boundary data from input file
    call atmos_boundary_update_file( ref_now )
 
    if ( atmos_boundary_dens_adjust ) then
       call calc_mass( ref_now )
    end if
 
    boundary_timestep = boundary_timestep + 1
    call atmos_boundary_update_file( ref_new )
 
    if ( atmos_boundary_dens_adjust ) then
       call calc_mass( ref_new )
       call set_offset
    end if
 
    ! copy now to old
    !$omp parallel do default(none) private(i,j,k,iq) OMP_SCHEDULE_ collapse(2) &
    !$omp shared(JA,IA,KS,KE,ATMOS_BOUNDARY_ref_DENS,ref_old,ref_now,ATMOS_BOUNDARY_ref_VELX) &
    !$omp shared(ATMOS_BOUNDARY_ref_VELY,ATMOS_BOUNDARY_ref_POTT,BND_QA,ATMOS_BOUNDARY_ref_QTRC)
    do j = 1, ja
    do i = 1, ia
    do k = ks, ke
       atmos_boundary_ref_dens(k,i,j,ref_old) = atmos_boundary_ref_dens(k,i,j,ref_now)
       atmos_boundary_ref_velx(k,i,j,ref_old) = atmos_boundary_ref_velx(k,i,j,ref_now)
       atmos_boundary_ref_vely(k,i,j,ref_old) = atmos_boundary_ref_vely(k,i,j,ref_now)
       atmos_boundary_ref_pott(k,i,j,ref_old) = atmos_boundary_ref_pott(k,i,j,ref_now)
       do iq = 1, bnd_qa
          atmos_boundary_ref_qtrc(k,i,j,iq,ref_old) = atmos_boundary_ref_qtrc(k,i,j,iq,ref_now)
       end do
    end do
    end do
    end do
 
    if ( atmos_boundary_use_velz ) then
       !$omp parallel do collapse(2)
       do j = 1, ja
       do i = 1, ia
       do k = ks, ke
          atmos_boundary_ref_velz(k,i,j,:) = atmos_boundary_value_velz
       end do
       end do
       end do
    end if
 
    now_step = fillgaps_steps
 
    ! set boundary data
    call set_boundary( atmos_boundary_use_velz )
 
    return

References bnd_qa, calc_mass(), scale_calendar::calendar_combine_daysec(), scale_calendar::calendar_date2daysec(), scale_comm_cartesc_nest::comm_cartesc_nest_bnd_qa, scale_comm_cartesc_nest::comm_cartesc_nest_recvwait_issue(), scale_const::const_pi, scale_atmos_grid_cartesc_index::ke, scale_atmos_grid_cartesc_index::ks, scale_comm_cartesc_nest::online_use_velz, scale_comm_cartesc_nest::parent_dtsec, scale_prc::prc_abort(), set_boundary(), set_offset(), scale_time::time_dtsec, and scale_time::time_nowdate.

Referenced by atmos_boundary_driver_set().

Here is the call graph for this function:

Here is the caller graph for this function:

atmos_boundary_set_online()

subroutine mod_atmos_bnd_driver::atmos_boundary_set_online

Set boundary value for real case experiment [online daughter].

Definition at line 1579 of file mod_atmos_bnd_driver.F90.

    use scale_prc, only: &
       prc_abort
    use scale_const, only: &
       pi => const_pi
    use scale_time, only: &
       time_dtsec,        &
       time_nstep
    use scale_comm_cartesc_nest, only: &
       online_use_velz,          &
       parent_nstep
    implicit none
 
    ! parameters
    integer, parameter  :: handle = 2
 
    real(RP) :: total
 
    ! works
    integer  :: i, j, k, n, iq
    !---------------------------------------------------------------------------
 
    ! import data from parent domain
    boundary_timestep = 1
    log_info("ATMOS_BOUNDARY_set_online",*) 'BOUNDARY TIMESTEP NUMBER FOR INIT:', boundary_timestep
 
    call atmos_boundary_update_online_daughter( ref_now )
 
    if ( atmos_boundary_dens_adjust ) then
       call calc_mass( ref_now )
    end if
 
    boundary_timestep = boundary_timestep + 1
    log_info("ATMOS_BOUNDARY_set_online",*) 'BOUNDARY TIMESTEP NUMBER FOR INIT:', boundary_timestep
 
    call atmos_boundary_update_online_daughter( ref_new )
 
    if ( atmos_boundary_dens_adjust ) then
       call calc_mass( ref_new )
       call set_offset
    end if
 
    ! copy now to old
    !$omp parallel do collapse(2)
    do j = 1, ja
    do i = 1, ia
    do k = ks, ke
       atmos_boundary_ref_dens(k,i,j,ref_old) = atmos_boundary_ref_dens(k,i,j,ref_now)
       atmos_boundary_ref_velx(k,i,j,ref_old) = atmos_boundary_ref_velx(k,i,j,ref_now)
       atmos_boundary_ref_vely(k,i,j,ref_old) = atmos_boundary_ref_vely(k,i,j,ref_now)
       atmos_boundary_ref_pott(k,i,j,ref_old) = atmos_boundary_ref_pott(k,i,j,ref_now)
       do iq = 1, bnd_qa
          atmos_boundary_ref_qtrc(k,i,j,iq,ref_old) = atmos_boundary_ref_qtrc(k,i,j,iq,ref_now)
       end do
    end do
    end do
    end do
 
    update_nstep = nint( atmos_boundary_update_dt / time_dtsec )
    if ( update_nstep * time_dtsec /= atmos_boundary_update_dt ) then
       log_error("ATMOS_BOUNDARY_set_online",*) 'DT of the parent is not multiple of the DT'
       call prc_abort
    end if
    if ( update_nstep * parent_nstep(handle) /= time_nstep ) then
       log_error("ATMOS_BOUNDARY_set_online",*) 'DURATION must be the same as that of the parent'
       call prc_abort
    end if
 
    now_step = 0 ! should be set as zero in initialize process
 
    if ( atmos_boundary_dens_adjust ) then
       allocate( offset_time_fact(update_nstep) )
       total = 0.0_rp
       !$omp parallel do reduction(+:total)
       do n = 0, update_nstep
          offset_time_fact(n) = 1.0_rp - cos( 2.0_rp * pi * ( n - 1 ) / update_nstep )
          total = total + offset_time_fact(n)
       end do
       total = total / update_nstep
       !$omp parallel do
       do n = 0, update_nstep
          offset_time_fact(n) = offset_time_fact(n) / total
       end do
    end if
 
    ! set boundary data
    call set_boundary( online_use_velz )
 
    return

References bnd_qa, calc_mass(), scale_const::const_pi, scale_atmos_grid_cartesc_index::ia, scale_atmos_grid_cartesc_index::ja, scale_atmos_grid_cartesc_index::ke, scale_atmos_grid_cartesc_index::ks, scale_atmos_hydrometeor::n_hyd, scale_comm_cartesc_nest::online_use_velz, scale_comm_cartesc_nest::parent_nstep, scale_prc::prc_abort(), mod_atmos_phy_mp_vars::qa_mp, set_boundary(), set_offset(), scale_time::time_dtsec, and scale_time::time_nstep.

Referenced by atmos_boundary_driver_set().

Here is the call graph for this function:

Here is the caller graph for this function:

atmos_boundary_driver_finalize()

subroutine, public mod_atmos_bnd_driver::atmos_boundary_driver_finalize

Finalize boundary value.

Definition at line 1702 of file mod_atmos_bnd_driver.F90.

    use scale_comm_cartesc_nest, only: &
       comm_cartesc_nest_recvwait_issue, &
       comm_cartesc_nest_recv_cancel, &
       nestqa => comm_cartesc_nest_bnd_qa
    use scale_file_cartesc, only: &
       file_cartesc_close
    implicit none
 
    ! works
    integer :: handle
    !---------------------------------------------------------------------------
 
    if ( do_parent_process ) then !online [parent]
       handle = 1
       call comm_cartesc_nest_recvwait_issue( handle, nestqa )
    endif
 
    if ( do_daughter_process ) then !online [daughter]
       handle = 2
       call comm_cartesc_nest_recv_cancel( handle )
    endif
 
    if ( atmos_boundary_fid > 0 ) then
       call file_cartesc_close( atmos_boundary_fid )
       atmos_boundary_fid = -1
    end if
 
    return

References scale_comm_cartesc_nest::comm_cartesc_nest_bnd_qa, scale_comm_cartesc_nest::comm_cartesc_nest_recv_cancel(), scale_comm_cartesc_nest::comm_cartesc_nest_recvwait_issue(), and scale_file_cartesc::file_cartesc_close().

Referenced by mod_atmos_driver::atmos_driver_finalize().

Here is the call graph for this function:

Here is the caller graph for this function:

atmos_boundary_driver_update()

subroutine, public mod_atmos_bnd_driver::atmos_boundary_driver_update

(

logical, intent(in)

last_step

)

Update boundary value with a constant time boundary.

Definition at line 1736 of file mod_atmos_bnd_driver.F90.

    use scale_prc, only: &
       prc_abort
    use scale_comm_cartesc_nest, only: &
       online_use_velz,       &
       comm_cartesc_nest_test
    implicit none
 
    logical, intent(in) :: last_step
 
    integer :: handle
    !---------------------------------------------------------------------------
 
    if ( l_bnd ) then
 
       ! step boundary
       now_step = now_step + 1
 
       ! update referce vars
       if ( last_step ) then
          now_step = min( now_step, update_nstep-1 )
       else if ( now_step >= update_nstep ) then
          now_step          = 0
          boundary_timestep = boundary_timestep + 1
 
          call update_ref_index
 
          if ( do_daughter_process ) then !online [daughter]
             call atmos_boundary_update_online_daughter( ref_new )
          else
             call atmos_boundary_update_file( ref_new )
          end if
 
          if ( atmos_boundary_dens_adjust ) then
             call calc_mass( ref_new )
          end if
 
       end if
 
       call set_boundary( online_use_velz )
 
       if ( atmos_boundary_dens_adjust ) then
          call set_offset
       end if
 
    elseif ( do_parent_process ) then
       ! do nothing
    else
       log_error("ATMOS_BOUNDARY_update",*) '[BUG] invalid path'
       call prc_abort
    end if
 
    call history_bnd( atmos_boundary_dens, &
                      atmos_boundary_velz, &
                      atmos_boundary_velx, &
                      atmos_boundary_vely, &
                      atmos_boundary_pott, &
                      atmos_boundary_qtrc )
 
    ! To be enable to do asynchronous communicaton
    if ( do_daughter_process ) then !online [daughter]
       handle = 2
       call comm_cartesc_nest_test( handle )
    endif
 
    return

References atmos_boundary_dens, atmos_boundary_pott, atmos_boundary_qtrc, atmos_boundary_velx, atmos_boundary_vely, atmos_boundary_velz, calc_mass(), scale_comm_cartesc_nest::comm_cartesc_nest_test(), scale_comm_cartesc_nest::online_use_velz, scale_prc::prc_abort(), set_boundary(), set_offset(), and update_ref_index().

Referenced by mod_atmos_driver::atmos_driver_update().

Here is the call graph for this function:

Here is the caller graph for this function:

atmos_boundary_driver_send()

subroutine, public mod_atmos_bnd_driver::atmos_boundary_driver_send

Send data to child domain.

Definition at line 1806 of file mod_atmos_bnd_driver.F90.

    use scale_comm_cartesc_nest, only: &
       comm_cartesc_nest_test
    use mod_atmos_vars, only: &
       dens, &
       momz, &
       momx, &
       momy, &
       rhot, &
       qtrc, &
       qv,   &
       qe
    use mod_atmos_phy_mp_vars, only: &
       qs_mp, &
       qe_mp
    implicit none
 
    integer :: handle
    !---------------------------------------------------------------------------
 
    if ( do_parent_process ) then !online [parent]
       ! should be called every time step
       call atmos_boundary_update_online_parent( dens,momz,momx,momy,rhot,qtrc(:,:,:,qs_mp:qe_mp), qv, qe )
 
       ! To be enable to do asynchronous communicaton
       handle = 1
       call comm_cartesc_nest_test( handle )
    endif
 
    return

References mod_atmos_phy_mp_driver::atmos_phy_mp_driver_qhyd2qtrc(), bnd_iq, scale_comm_cartesc_nest::comm_cartesc_nest_bnd_qa, scale_comm_cartesc_nest::comm_cartesc_nest_nestdown(), scale_comm_cartesc_nest::comm_cartesc_nest_recvwait_issue(), scale_comm_cartesc_nest::comm_cartesc_nest_test(), scale_comm_cartesc_nest::daughter_ia, scale_comm_cartesc_nest::daughter_ja, scale_comm_cartesc_nest::daughter_ka, dens, scale_file_cartesc::file_cartesc_flush(), scale_atmos_grid_cartesc_index::ia, scale_atmos_grid_cartesc_index::ja, scale_atmos_grid_cartesc_index::ka, scale_atmos_grid_cartesc_index::ke, scale_atmos_grid_cartesc_index::ks, mod_atmos_vars::momx, mod_atmos_vars::momy, mod_atmos_vars::momz, scale_atmos_hydrometeor::n_hyd, scale_comm_cartesc_nest::online_boundary_diagqhyd, scale_comm_cartesc_nest::parent_ia, scale_comm_cartesc_nest::parent_ja, scale_comm_cartesc_nest::parent_ka, scale_prc::prc_abort(), scale_tracer::qa, mod_atmos_phy_mp_vars::qa_mp, mod_atmos_vars::qe, mod_atmos_phy_mp_vars::qe_mp, mod_atmos_phy_mp_vars::qs_mp, qtrc, mod_atmos_vars::qv, mod_atmos_vars::rhot, and scale_tracer::tracer_name.

Referenced by mod_atmos_driver::atmos_driver_update().

Here is the call graph for this function:

Here is the caller graph for this function:

set_boundary()

subroutine mod_atmos_bnd_driver::set_boundary

(

logical, intent(in)

use_velz

)

Definition at line 2077 of file mod_atmos_bnd_driver.F90.

    use scale_prc_cartesc, only: &
       prc_has_w,   &
       prc_has_e,   &
       prc_has_s,   &
       prc_has_n
    use mod_atmos_vars, only: &
       dens, &
       momz, &
       momx, &
       momy, &
       rhot, &
       qtrc
    implicit none
    logical, intent(in) :: use_velz
 
    real(RP) :: bnd_DENS(KA,IA,JA)        ! damping coefficient for DENS (0-1)
    real(RP) :: bnd_VELZ(KA,IA,JA)        ! damping coefficient for VELZ (0-1)
    real(RP) :: bnd_VELX(KA,IA,JA)        ! damping coefficient for VELX (0-1)
    real(RP) :: bnd_VELY(KA,IA,JA)        ! damping coefficient for VELY (0-1)
    real(RP) :: bnd_POTT(KA,IA,JA)        ! damping coefficient for POTT (0-1)
    real(RP) :: bnd_QTRC(KA,IA,JA,BND_QA) ! damping coefficient for QTRC (0-1)
 
    integer :: i, j, k, iq, iqb
 
 
    ! get boundaryal coefficients
    call get_boundary( bnd_dens(:,:,:),   & ! [OUT]
                       bnd_velz(:,:,:),   & ! [OUT]
                       bnd_velx(:,:,:),   & ! [OUT]
                       bnd_vely(:,:,:),   & ! [OUT]
                       bnd_pott(:,:,:),   & ! [OUT]
                       bnd_qtrc(:,:,:,:), & ! [OUT]
                       now_step,          & ! [IN]
                       update_nstep       ) ! [IN]
 
    ! update boundary vars
    !$omp parallel do collapse(2)
    do j = 1, ja
    do i = 1, ia
    do k = ks, ke
       atmos_boundary_dens(k,i,j) = bnd_dens(k,i,j)
       atmos_boundary_velx(k,i,j) = bnd_velx(k,i,j)
       atmos_boundary_vely(k,i,j) = bnd_vely(k,i,j)
       atmos_boundary_pott(k,i,j) = bnd_pott(k,i,j)
       do iq = 1, bnd_qa
          atmos_boundary_qtrc(k,i,j,iq) = bnd_qtrc(k,i,j,iq)
       end do
    end do
    end do
    end do
    if ( use_velz ) then
       !$omp parallel do collapse(2)
       do j = 1, ja
       do i = 1, ia
       do k = ks, ke
          atmos_boundary_velz(k,i,j) = bnd_velz(k,i,j)
       end do
       end do
       end do
    end if
 
    ! fill HALO in western region
    if ( .NOT. prc_has_w ) then
       !$omp parallel do default(none) OMP_SCHEDULE_ collapse(2) &
       !$omp shared(JA,JS,IS,KS,KE,QA,DENS,MOMX,RHOT,QTRC) &
       !$omp shared(ATMOS_BOUNDARY_DENS,ATMOS_BOUNDARY_VELX) &
       !$omp shared(ATMOS_BOUNDARY_POTT,ATMOS_BOUNDARY_QTRC) &
       !$omp shared(BND_QA,BND_IQ) &
       !$omp private(i,j,k,iq,iqb)
       do j = 1, ja
       do i = 1, is-1
       do k = ks, ke
          dens(k,i,j) = atmos_boundary_dens(k,i,j)
          momx(k,i,j) = atmos_boundary_velx(k,i,j) &
                  * ( atmos_boundary_dens(k,i,j) + atmos_boundary_dens(k,i+1,j) ) * 0.5_rp
          rhot(k,i,j) = atmos_boundary_pott(k,i,j) * atmos_boundary_dens(k,i,j)
          do iq = 1, qa
             iqb = bnd_iq(iq)
             if ( iqb > 0 ) then
                qtrc(k,i,j,iq) = atmos_boundary_qtrc(k,i,j,iqb)
             else
                qtrc(k,i,j,iq) = qtrc(k,is,j,iq)
             end if
          end do
       end do
       end do
       end do
       !$omp parallel do
       do j = 1, ja-1
       do i = 1, is-1
       do k = ks, ke
          momy(k,i,j) = atmos_boundary_vely(k,i,j) &
                  * ( atmos_boundary_dens(k,i,j) + atmos_boundary_dens(k,i,j+1) ) * 0.5_rp
       end do
       end do
       end do
       do i = 1, is-1
       do k = ks, ke
          momy(k,i,ja) = atmos_boundary_vely(k,i,ja) &
                  * atmos_boundary_dens(k,i,ja)
       end do
       end do
       if ( use_velz ) then
          !$omp parallel do
          do j = 1, ja
          do i = 1, is-1
          do k = ks, ke-1
             momz(k,i,j) = atmos_boundary_velz(k,i,j) &
                     * ( atmos_boundary_dens(k,i,j) + atmos_boundary_dens(k+1,i,j) ) * 0.5_rp
          end do
          end do
          end do
       else
          !$omp parallel do
          do j = 1, ja
          do i = 1, is-1
          do k = ks, ke-1
             momz(k,i,j) = momz(k,is,j)
          end do
          end do
          end do
       end if
    end if
 
    ! fill HALO in eastern region
    if ( .NOT. prc_has_e ) then
       !$omp parallel do default(none) OMP_SCHEDULE_ collapse(2) &
       !$omp shared(JA,IE,IA,KS,KE,QA) &
       !$omp shared(DENS,RHOT,QTRC) &
       !$omp shared(ATMOS_BOUNDARY_DENS,ATMOS_BOUNDARY_VELX) &
       !$omp shared(ATMOS_BOUNDARY_POTT,ATMOS_BOUNDARY_QTRC) &
       !$omp shared(BND_QA,BND_IQ) &
       !$omp private(i,j,k,iq,iqb)
       do j = 1, ja
       do i = ie+1, ia
       do k = ks, ke
          dens(k,i,j) = atmos_boundary_dens(k,i,j)
          rhot(k,i,j) = atmos_boundary_pott(k,i,j) * atmos_boundary_dens(k,i,j)
          do iq = 1, qa
             iqb = bnd_iq(iq)
             if ( iqb > 0 ) then
                qtrc(k,i,j,iq) = atmos_boundary_qtrc(k,i,j,iqb)
             else
                qtrc(k,i,j,iq) = qtrc(k,ie,j,iq)
             end if
          end do
       end do
       end do
       end do
       !$omp parallel do
       do j = 1, ja
       do i = ie, ia-1
       do k = ks, ke
          momx(k,i,j) = atmos_boundary_velx(k,i,j) &
                  * ( atmos_boundary_dens(k,i,j) + atmos_boundary_dens(k,i+1,j) ) * 0.5_rp
       end do
       end do
       end do
       !$omp parallel do
       do j = 1, ja
       do k = ks, ke
          momx(k,ia,j) = atmos_boundary_velx(k,ia,j) * atmos_boundary_dens(k,ia,j)
       end do
       end do
       !$omp parallel do
       do j = 1, ja-1
       do i = ie+1, ia
       do k = ks, ke
          momy(k,i,j) = atmos_boundary_vely(k,i,j) &
                  * ( atmos_boundary_dens(k,i,j) + atmos_boundary_dens(k,i,j+1) ) * 0.5_rp
       end do
       end do
       end do
       do i = ie+1, ia
       do k = ks, ke
          momy(k,i,ja) = atmos_boundary_vely(k,i,ja) &
                  * atmos_boundary_dens(k,i,ja)
       end do
       end do
       if ( use_velz ) then
          !$omp parallel do
          do j = 1, ja
          do i = ie+1, ia
          do k = ks, ke-1
             momz(k,i,j) = atmos_boundary_velz(k,i,j) &
                     * ( atmos_boundary_dens(k,i,j) + atmos_boundary_dens(k+1,i,j) ) * 0.5_rp
          end do
          end do
          end do
       else
          !$omp parallel do
          do j = 1, ja
          do i = ie+1, ia
          do k = ks, ke-1
             momz(k,i,j) = momz(k,ie,j)
          end do
          end do
          end do
       end if
    end if
 
    ! fill HALO in southern region
    if ( .NOT. prc_has_s ) then
       do j = 1, js-1
       do i = 1, ia
       do k = ks, ke
          dens(k,i,j) = atmos_boundary_dens(k,i,j)
          momy(k,i,j) = atmos_boundary_vely(k,i,j) &
                  * ( atmos_boundary_dens(k,i,j) + atmos_boundary_dens(k,i,j+1) ) * 0.5_rp
          rhot(k,i,j) = atmos_boundary_pott(k,i,j) * atmos_boundary_dens(k,i,j)
          do iq = 1, qa
             iqb = bnd_iq(iq)
             if ( iqb > 0 ) then
                qtrc(k,i,j,iq) = atmos_boundary_qtrc(k,i,j,iqb)
             else
                qtrc(k,i,j,iq) = qtrc(k,i,js,iq)
             end if
          end do
       end do
       end do
       end do
       do j = 1, js-1
       do i = 1, ia-1
       do k = ks, ke
          momx(k,i,j) = atmos_boundary_velx(k,i,j) &
                  * ( atmos_boundary_dens(k,i,j) + atmos_boundary_dens(k,i+1,j) ) * 0.5_rp
       end do
       end do
       end do
       do j = 1, js-1
       do k = ks, ke
          momx(k,ia,j) = atmos_boundary_velx(k,ia,j) &
                  * atmos_boundary_dens(k,ia,j)
       end do
       end do
       if ( use_velz ) then
          do j = 1, js-1
          do i = 1, ia
          do k = ks, ke-1
             momz(k,i,j) = atmos_boundary_velz(k,i,j) &
                     * ( atmos_boundary_dens(k,i,j) + atmos_boundary_dens(k+1,i,j) ) * 0.5_rp
          end do
          end do
          end do
       else
          do j = 1, js-1
          do i = 1, ia
          do k = ks, ke-1
             momz(k,i,j) = momz(k,i,js)
          end do
          end do
          end do
       end if
    end if
 
    ! fill HALO in northern region
    if ( .NOT. prc_has_n ) then
       do j = je+1, ja
       do i = 1, ia
       do k = ks, ke
          dens(k,i,j) = atmos_boundary_dens(k,i,j)
          rhot(k,i,j) = atmos_boundary_pott(k,i,j) * atmos_boundary_dens(k,i,j)
          do iq = 1, qa
             iqb = bnd_iq(iq)
             if ( iqb > 0 ) then
                qtrc(k,i,j,iq) = atmos_boundary_qtrc(k,i,j,iqb)
             else
                qtrc(k,i,j,iq) = qtrc(k,i,je,iq)
             end if
          end do
       end do
       end do
       end do
       do j = je, ja-1
       do i = 1, ia
       do k = ks, ke
          momy(k,i,j) = atmos_boundary_vely(k,i,j) &
                  * ( atmos_boundary_dens(k,i,j) + atmos_boundary_dens(k,i,j+1) ) * 0.5_rp
       end do
       end do
       end do
       do i = 1, ia
       do k = ks, ke
          momy(k,i,ja) = atmos_boundary_vely(k,i,ja) * atmos_boundary_dens(k,i,ja)
       end do
       end do
       do j = je+1, ja
       do i = 1, ia-1
       do k = ks, ke
          momx(k,i,j) = atmos_boundary_velx(k,i,j) &
                  * ( atmos_boundary_dens(k,i,j) + atmos_boundary_dens(k,i+1,j) ) * 0.5_rp
       end do
       end do
       end do
       do j = je+1, ja
       do k = ks, ke
          momx(k,ia,j) = atmos_boundary_velx(k,ia,j) &
                  * atmos_boundary_dens(k,ia,j)
       end do
       end do
       if ( use_velz ) then
          do j = je+1, ja
          do i = 1, ia
          do k = ks, ke-1
             momz(k,i,j) = atmos_boundary_velz(k,i,j) &
                     * ( atmos_boundary_dens(k,i,j) + atmos_boundary_dens(k+1,i,j) ) * 0.5_rp
          end do
          end do
          end do
       else
          do j = je+1, ja
          do i = 1, ia
          do k = ks, ke-1
             momz(k,i,j) = momz(k,i,je)
          end do
          end do
          end do
       end if
    end if
 
    return

References atmos_boundary_dens, atmos_boundary_pott, atmos_boundary_qtrc, atmos_boundary_velx, atmos_boundary_vely, atmos_boundary_velz, bnd_iq, bnd_qa, 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::ke, scale_atmos_grid_cartesc_index::ks, mod_atmos_vars::momx, mod_atmos_vars::momy, mod_atmos_vars::momz, scale_prc_cartesc::prc_has_e, scale_prc_cartesc::prc_has_n, scale_prc_cartesc::prc_has_s, scale_prc_cartesc::prc_has_w, scale_tracer::qa, qtrc, real, mod_atmos_vars::rhot, and scale_time::time_dtsec.

Referenced by atmos_boundary_driver_update(), atmos_boundary_set_file(), and atmos_boundary_set_online().

Here is the caller graph for this function:

update_ref_index()

subroutine mod_atmos_bnd_driver::update_ref_index

Update indices of array of boundary references.

Definition at line 2716 of file mod_atmos_bnd_driver.F90.

    implicit none
 
    ! works
    integer :: ref_tmp
    !---------------------------------------------------------------------------
 
    ref_tmp = ref_old
    ref_old = ref_now
    ref_now = ref_new
    ref_new = ref_tmp
 
    return

References bnd_iq, scale_tracer::qa, mod_atmos_phy_mp_vars::qs_mp, and scale_tracer::tracer_name.

Referenced by atmos_boundary_driver_update().

Here is the caller graph for this function:

calc_mass()

subroutine mod_atmos_bnd_driver::calc_mass

(

integer, intent(in)

ref

)

Definition at line 2768 of file mod_atmos_bnd_driver.F90.

    use scale_prc_cartesc, only: &
       prc_has_w, &
       prc_has_e, &
       prc_has_s, &
       prc_has_n
    use scale_statistics, only: &
       statistics_total
    use scale_atmos_grid_cartesc_real, only: &
       vol          => atmos_grid_cartesc_real_vol,          &
       totvol       => atmos_grid_cartesc_real_totvol,       &
       areazxv_y    => atmos_grid_cartesc_real_areazxv_y,    &
       totareazuy_x => atmos_grid_cartesc_real_totareazuy_x, &
       totareazxv_y => atmos_grid_cartesc_real_totareazxv_y
    use scale_atmos_refstate, only: &
       dens_ref => atmos_refstate_dens
    use mod_atmos_vars, only: &
       dens
    implicit none
    integer,  intent(in)  :: ref
 
    real(DP) :: masstot, masstot_current
    real(DP) :: massflx
    real(DP) :: offset_band, offset_bias
    real(DP) :: ref_tot
    real(DP) :: flx_w, flx_e, flx_s, flx_n
    real(DP) :: ref_w, ref_e, ref_s, ref_n
 
    real(RP), target :: work_x(KA,JA), work_y(KA,IA)
    real(RP), pointer :: ptr(:,:)
 
    integer :: k, i, j
 
    ! total mass
    call statistics_total( ka, ks, ke, ia, is, ie, ja, js, je, &
                           atmos_boundary_ref_dens(:,:,:,ref),      & ! (in)
                           "DENS_bnd",                              & ! (in)
                           vol(:,:,:), totvol,                      & ! (in)
                           log_suppress = .true., global = .true.,  & ! (in)
                           sum = masstot                            ) ! (out)
 
!!$    call STATISTICS_total( KA, KS, KE, IA, IS, IE, JA, JS, JE, &
!!$                           DENS(:,:,:), "DENS_bnd_update",         & ! (in)
!!$                           VOL(:,:,:), TOTVOL,                     & ! (in)
!!$                           log_suppress = .true., global = .true., & ! (in)
!!$                           sum = masstot_current                   ) ! (out)
 
 
    ! West
    if ( .NOT. prc_has_w ) then
       !$omp parallel do
       do j = js, je
       do k = ks, ke
          work_x(k,j) = atmos_boundary_ref_velx(k,is-1,j,ref) &
                      * ( atmos_boundary_ref_dens(k,is-1,j,ref) + atmos_boundary_ref_dens(k,is,j,ref) ) * 0.5_rp
       end do
       end do
       ptr => work_x
    else
       ptr => zero_x
    end if
    call statistics_total( ka, ks, ke, ja, js, je, &
                           ptr(:,:), "MFLUX_bnd_w",                & ! (in)
                           areazuy_w(:,:), totareazuy_x(is-1),     & ! (in)
                           log_suppress = .true., global = .true., & ! (in)
                           sum = flx_w                             ) ! (out)
    if ( .NOT. prc_has_w ) then
       !$omp parallel do
       do j = js, je
       do k = ks, ke
          work_x(k,j) = dens_ref(k,is,j)
       end do
       end do
    end if
    call statistics_total( ka, ks, ke, ja, js, je, &
                           ptr(:,:), "DENS_ref_w",                 & ! (in)
                           areazuy_w(:,:), totareazuy_x(is-1),     & ! (in)
                           log_suppress = .true., global = .true., & ! (in)
                           sum = ref_w                             ) ! (out)
 
    ! East
    if ( .NOT. prc_has_e ) then
       !$omp parallel do
       do j = js, je
       do k = ks, ke
          work_x(k,j) = atmos_boundary_ref_velx(k,ie,j,ref) &
                      * ( atmos_boundary_ref_dens(k,ie,j,ref) + atmos_boundary_ref_dens(k,ie+1,j,ref) ) * 0.5_rp
       end do
       end do
       ptr => work_x
    else
       ptr => zero_x
    end if
    call statistics_total( ka, ks, ke, ja, js, je, &
                           ptr(:,:), "MFLUX_bnd_e",                & ! (in)
                           areazuy_e(:,:), totareazuy_x(ie),       & ! (in)
                           log_suppress = .true., global = .true., & ! (in)
                           sum = flx_e                             ) ! (out)
    if ( .NOT. prc_has_e ) then
       !$omp parallel do
       do j = js, je
       do k = ks, ke
          work_x(k,j) = dens_ref(k,ie,j)
       end do
       end do
    end if
    call statistics_total( ka, ks, ke, ja, js, je, &
                           ptr(:,:), "DENS_ref_e",                 & ! (in)
                           areazuy_e(:,:), totareazuy_x(ie),       & ! (in)
                           log_suppress = .true., global = .true., & ! (in)
                           sum = ref_e                             ) ! (out)
 
    ! South
    if ( .NOT. prc_has_s ) then
       !$omp parallel do
       do i = is, ie
       do k = ks, ke
          work_y(k,i) = atmos_boundary_ref_vely(k,i,js-1,ref) &
                      * ( atmos_boundary_ref_dens(k,i,js-1,ref) + atmos_boundary_ref_dens(k,i,js,ref) ) * 0.5_rp
       end do
       end do
       ptr => work_y
    else
       ptr => zero_y
    end if
    call statistics_total( ka, ks, ke, ia, is, ie, &
                           ptr(:,:), "MFLUX_bnd_s",                 & ! (in)
                           areazxv_y(:,:,js-1), totareazxv_y(js-1), & ! (in)
                           log_suppress = .true., global = .true.,  & ! (in)
                           sum = flx_s                              ) ! (out)
    if ( .NOT. prc_has_s ) then
       !$omp parallel do
       do i = is, ie
       do k = ks, ke
          work_y(k,i) = dens_ref(k,i,js)
       end do
       end do
    end if
    call statistics_total( ka, ks, ke, ia, is, ie, &
                           ptr(:,:), "DENS_ref_s",                  & ! (in)
                           areazxv_y(:,:,js-1), totareazxv_y(js-1), & ! (in)
                           log_suppress = .true., global = .true.,  & ! (in)
                           sum = ref_s                              ) ! (out)
 
    ! North
    if ( .NOT. prc_has_n ) then
       !$omp parallel do
       do i = is, ie
       do k = ks, ke
          work_y(k,i) = atmos_boundary_ref_vely(k,i,je,ref) &
                      * ( atmos_boundary_ref_dens(k,i,je,ref) + atmos_boundary_ref_dens(k,i,je+1,ref) ) * 0.5_rp
       end do
       end do
       ptr => work_y
    else
       ptr => zero_y
    end if
    call statistics_total( ka, ks, ke, ia, is, ie, &
                           ptr(:,:), "MFLX_bnd_n",                 & ! (in)
                           areazxv_y(:,:,je), totareazxv_y(je),    & ! (in)
                           log_suppress = .true., global = .true., & ! (in)
                           sum = flx_n                             ) ! (out)
    if ( .NOT. prc_has_n ) then
       !$omp parallel do
       do i = is, ie
       do k = ks, ke
          work_y(k,i) = dens_ref(k,i,je)
       end do
       end do
    end if
    call statistics_total( ka, ks, ke, ia, is, ie, &
                           ptr(:,:), "DENS_ref_n",                 & ! (in)
                           areazxv_y(:,:,je), totareazxv_y(je),    & ! (in)
                           log_suppress = .true., global = .true., & ! (in)
                           sum = ref_n                             ) ! (out)
 
    massflx = flx_w - flx_e + flx_s - flx_n
 
    offset_band = ( masstot - masstot_now ) / atmos_boundary_update_dt &
                - ( massflx + massflx_now ) * 0.5_dp
!    offset_bias = ( MASSTOT_now - masstot_current ) / ATMOS_BOUNDARY_UPDATE_DT
 
    log_info("ATMOS_BOUNDARY_calc_mass",*) "Offset_band is: ", offset_band, "(", masstot, masstot_now, massflx, massflx_now, ")"
 
    ref_tot = ref_w + ref_e + ref_s + ref_n
    offset_band = offset_band / ref_tot
!    offset_bias = offset_bias / ref_tot
 
    log_info_cont(*) "          per dens  ", offset_band
 
    ! density of the reference state is used as weight
    if ( .not. prc_has_w ) then
       !$omp parallel do
       do j = js, je
       do k = ks, ke
          mflux_offset_x(k,j,1,1) = offset_band * dens_ref(k,is,j)
!          MFLUX_OFFSET_X(k,j,1,2) = offset_bias * DENS_ref(k,IS,j)
       end do
       end do
    end if
    if ( .not. prc_has_e ) then
       !$omp parallel do
       do j = js, je
       do k = ks, ke
          mflux_offset_x(k,j,2,1) = - offset_band * dens_ref(k,ie,j)
!          MFLUX_OFFSET_X(k,j,2,2) = - offset_bias * DENS_ref(k,IE,j)
       end do
       end do
    end if
    if ( .not. prc_has_s ) then
       !$omp parallel do
       do i = is, ie
       do k = ks, ke
          mflux_offset_y(k,i,1,1) = offset_band * dens_ref(k,i,js)
!          MFLUX_OFFSET_Y(k,i,1,2) = offset_bias * DENS_ref(k,i,JS)
       end do
       end do
    end if
    if ( .not. prc_has_n ) then
       !$omp parallel do
       do i = is, ie
       do k = ks, ke
          mflux_offset_y(k,i,2,1) = - offset_band * dens_ref(k,i,je)
!          MFLUX_OFFSET_Y(k,i,2,2) = - offset_bias * DENS_ref(k,i,JE)
       end do
       end do
    end if
 
    masstot_now = masstot
    massflx_now = massflx
 
    return

References scale_atmos_grid_cartesc_real::atmos_grid_cartesc_real_areazxv_y, scale_atmos_grid_cartesc_real::atmos_grid_cartesc_real_totareazuy_x, scale_atmos_grid_cartesc_real::atmos_grid_cartesc_real_totareazxv_y, scale_atmos_grid_cartesc_real::atmos_grid_cartesc_real_totvol, scale_atmos_grid_cartesc_real::atmos_grid_cartesc_real_vol, scale_atmos_refstate::atmos_refstate_dens, dens, scale_atmos_grid_cartesc_index::ie, scale_atmos_grid_cartesc_index::is, scale_atmos_grid_cartesc_index::je, scale_atmos_grid_cartesc_index::js, scale_atmos_grid_cartesc_index::ke, scale_atmos_grid_cartesc_index::ks, scale_prc_cartesc::prc_has_e, scale_prc_cartesc::prc_has_n, scale_prc_cartesc::prc_has_s, and scale_prc_cartesc::prc_has_w.

Referenced by atmos_boundary_driver_update(), atmos_boundary_set_file(), and atmos_boundary_set_online().

Here is the caller graph for this function:

set_offset()

subroutine mod_atmos_bnd_driver::set_offset

Definition at line 3003 of file mod_atmos_bnd_driver.F90.

    integer :: k, i, j, n
 
    !$omp parallel do
    do j = js, je
    do k = ks, ke
    do n = 1, 2
       atmos_boundary_mflux_offset_x(k,j,n) = mflux_offset_x(k,j,n,1) * offset_time_fact(now_step) !&
!                                            + MFLUX_OFFSET_X(k,j,n,2)
    end do
    end do
    end do
 
    !$omp parallel do
    do i = is, ie
    do k = ks, ke
    do n = 1, 2
       atmos_boundary_mflux_offset_y(k,i,n) = mflux_offset_y(k,i,n,1) * offset_time_fact(now_step) !&
!                                            + MFLUX_OFFSET_Y(k,i,n,2)
    end do
    end do
    end do
 
    return

References atmos_boundary_mflux_offset_x, atmos_boundary_mflux_offset_y, scale_atmos_grid_cartesc_index::ie, scale_atmos_grid_cartesc_index::is, scale_atmos_grid_cartesc_index::je, scale_atmos_grid_cartesc_index::js, scale_atmos_grid_cartesc_index::ke, and scale_atmos_grid_cartesc_index::ks.

Referenced by atmos_boundary_driver_update(), atmos_boundary_set_file(), and atmos_boundary_set_online().

Here is the caller graph for this function:

Variable Documentation

bnd_qa

integer, public mod_atmos_bnd_driver::bnd_qa

Definition at line 42 of file mod_atmos_bnd_driver.F90.

  integer,  public              :: BND_QA

Referenced by atmos_boundary_driver_set(), atmos_boundary_driver_setup(), atmos_boundary_set_file(), atmos_boundary_set_online(), mod_atmos_dyn_driver::atmos_dyn_driver(), and set_boundary().

of

integer allocatable public mod_atmos_bnd_driver::of

Definition at line 42 of file mod_atmos_bnd_driver.F90.

tracer

integer allocatable public mod_atmos_bnd_driver::tracer

Definition at line 42 of file mod_atmos_bnd_driver.F90.

at

integer, public mod_atmos_bnd_driver::at

Definition at line 42 of file mod_atmos_bnd_driver.F90.

boundary

integer, public mod_atmos_bnd_driver::boundary

Definition at line 42 of file mod_atmos_bnd_driver.F90.

bnd_iq

integer, dimension(:), allocatable, public mod_atmos_bnd_driver::bnd_iq

Definition at line 43 of file mod_atmos_bnd_driver.F90.

  integer,  public, allocatable :: BND_IQ(:)

Referenced by atmos_boundary_driver_send(), atmos_boundary_driver_set(), atmos_boundary_driver_setup(), mod_atmos_dyn_driver::atmos_dyn_driver(), set_boundary(), and update_ref_index().

index

integer, allocatable, public mod_atmos_bnd_driver::index

Definition at line 43 of file mod_atmos_bnd_driver.F90.

atmos_boundary_dens

real(rp), dimension(:,:,:), allocatable, public mod_atmos_bnd_driver::atmos_boundary_dens

Definition at line 45 of file mod_atmos_bnd_driver.F90.

  real(RP), public, allocatable :: ATMOS_BOUNDARY_DENS(:,:,:)

Referenced by atmos_boundary_driver_set(), atmos_boundary_driver_setup(), atmos_boundary_driver_update(), mod_atmos_dyn_driver::atmos_dyn_driver(), and set_boundary().

reference

real(rp), allocatable, public mod_atmos_bnd_driver::reference

Definition at line 45 of file mod_atmos_bnd_driver.F90.

dens

real(rp), dimension (0-1), allocatable, public mod_atmos_bnd_driver::dens

Definition at line 45 of file mod_atmos_bnd_driver.F90.

Referenced by atmos_boundary_driver_send(), atmos_boundary_driver_set(), calc_mass(), and set_boundary().

atmos_boundary_velz

real(rp), dimension(:,:,:), allocatable, public mod_atmos_bnd_driver::atmos_boundary_velz

Definition at line 46 of file mod_atmos_bnd_driver.F90.

  real(RP), public, allocatable :: ATMOS_BOUNDARY_VELZ(:,:,:)

Referenced by atmos_boundary_driver_set(), atmos_boundary_driver_setup(), atmos_boundary_driver_update(), mod_atmos_dyn_driver::atmos_dyn_driver(), and set_boundary().

velz

real(rp), dimension (0-1), allocatable, public mod_atmos_bnd_driver::velz

Definition at line 46 of file mod_atmos_bnd_driver.F90.

atmos_boundary_velx

real(rp), dimension(:,:,:), allocatable, public mod_atmos_bnd_driver::atmos_boundary_velx

Definition at line 47 of file mod_atmos_bnd_driver.F90.

  real(RP), public, allocatable :: ATMOS_BOUNDARY_VELX(:,:,:)

Referenced by atmos_boundary_driver_set(), atmos_boundary_driver_setup(), atmos_boundary_driver_update(), mod_atmos_dyn_driver::atmos_dyn_driver(), and set_boundary().

velx

real(rp), dimension (0-1), allocatable, public mod_atmos_bnd_driver::velx

Definition at line 47 of file mod_atmos_bnd_driver.F90.

atmos_boundary_vely

real(rp), dimension(:,:,:), allocatable, public mod_atmos_bnd_driver::atmos_boundary_vely

Definition at line 48 of file mod_atmos_bnd_driver.F90.

  real(RP), public, allocatable :: ATMOS_BOUNDARY_VELY(:,:,:)

Referenced by atmos_boundary_driver_set(), atmos_boundary_driver_setup(), atmos_boundary_driver_update(), mod_atmos_dyn_driver::atmos_dyn_driver(), and set_boundary().

vely

real(rp), dimension (0-1), allocatable, public mod_atmos_bnd_driver::vely

Definition at line 48 of file mod_atmos_bnd_driver.F90.

atmos_boundary_pott

real(rp), dimension(:,:,:), allocatable, public mod_atmos_bnd_driver::atmos_boundary_pott

Definition at line 49 of file mod_atmos_bnd_driver.F90.

  real(RP), public, allocatable :: ATMOS_BOUNDARY_POTT(:,:,:)

Referenced by atmos_boundary_driver_set(), atmos_boundary_driver_setup(), atmos_boundary_driver_update(), mod_atmos_dyn_driver::atmos_dyn_driver(), and set_boundary().

pott

real(rp), dimension (0-1), allocatable, public mod_atmos_bnd_driver::pott

Definition at line 49 of file mod_atmos_bnd_driver.F90.

atmos_boundary_qtrc

real(rp), dimension(:,:,:,:), allocatable, public mod_atmos_bnd_driver::atmos_boundary_qtrc

Definition at line 50 of file mod_atmos_bnd_driver.F90.

  real(RP), public, allocatable :: ATMOS_BOUNDARY_QTRC(:,:,:,:)

Referenced by atmos_boundary_driver_set(), atmos_boundary_driver_setup(), atmos_boundary_driver_update(), mod_atmos_dyn_driver::atmos_dyn_driver(), and set_boundary().

qtrc

real(rp), dimension (0-1), allocatable, public mod_atmos_bnd_driver::qtrc

Definition at line 50 of file mod_atmos_bnd_driver.F90.

Referenced by atmos_boundary_driver_send(), atmos_boundary_driver_set(), and set_boundary().

atmos_boundary_alpha_dens

real(rp), dimension(:,:,:), allocatable, public mod_atmos_bnd_driver::atmos_boundary_alpha_dens

Definition at line 52 of file mod_atmos_bnd_driver.F90.

  real(RP), public, allocatable :: ATMOS_BOUNDARY_alpha_DENS(:,:,:)

Referenced by atmos_boundary_driver_set(), atmos_boundary_driver_setup(), and mod_atmos_dyn_driver::atmos_dyn_driver().

damping

real(rp), allocatable, public mod_atmos_bnd_driver::damping

Definition at line 52 of file mod_atmos_bnd_driver.F90.

coefficient

real(rp), allocatable, public mod_atmos_bnd_driver::coefficient

Definition at line 52 of file mod_atmos_bnd_driver.F90.

for

logical allocatable public mod_atmos_bnd_driver::for

Definition at line 52 of file mod_atmos_bnd_driver.F90.

atmos_boundary_alpha_velz

real(rp), dimension(:,:,:), allocatable, public mod_atmos_bnd_driver::atmos_boundary_alpha_velz

Definition at line 53 of file mod_atmos_bnd_driver.F90.

  real(RP), public, allocatable :: ATMOS_BOUNDARY_alpha_VELZ(:,:,:)

Referenced by atmos_boundary_driver_set(), atmos_boundary_driver_setup(), and mod_atmos_dyn_driver::atmos_dyn_driver().

atmos_boundary_alpha_velx

real(rp), dimension(:,:,:), allocatable, public mod_atmos_bnd_driver::atmos_boundary_alpha_velx

Definition at line 54 of file mod_atmos_bnd_driver.F90.

  real(RP), public, allocatable :: ATMOS_BOUNDARY_alpha_VELX(:,:,:)

Referenced by atmos_boundary_driver_set(), atmos_boundary_driver_setup(), and mod_atmos_dyn_driver::atmos_dyn_driver().

atmos_boundary_alpha_vely

real(rp), dimension(:,:,:), allocatable, public mod_atmos_bnd_driver::atmos_boundary_alpha_vely

Definition at line 55 of file mod_atmos_bnd_driver.F90.

  real(RP), public, allocatable :: ATMOS_BOUNDARY_alpha_VELY(:,:,:)

Referenced by atmos_boundary_driver_set(), atmos_boundary_driver_setup(), and mod_atmos_dyn_driver::atmos_dyn_driver().

atmos_boundary_alpha_pott

real(rp), dimension(:,:,:), allocatable, public mod_atmos_bnd_driver::atmos_boundary_alpha_pott

Definition at line 56 of file mod_atmos_bnd_driver.F90.

  real(RP), public, allocatable :: ATMOS_BOUNDARY_alpha_POTT(:,:,:)

Referenced by atmos_boundary_driver_set(), atmos_boundary_driver_setup(), and mod_atmos_dyn_driver::atmos_dyn_driver().

atmos_boundary_alpha_qtrc

real(rp), dimension(:,:,:,:), allocatable, public mod_atmos_bnd_driver::atmos_boundary_alpha_qtrc

Definition at line 57 of file mod_atmos_bnd_driver.F90.

  real(RP), public, allocatable :: ATMOS_BOUNDARY_alpha_QTRC(:,:,:,:)

Referenced by atmos_boundary_driver_set(), atmos_boundary_driver_setup(), and mod_atmos_dyn_driver::atmos_dyn_driver().

atmos_boundary_mflux_offset_x

real(rp), dimension(:,:,:), allocatable, public mod_atmos_bnd_driver::atmos_boundary_mflux_offset_x

Definition at line 59 of file mod_atmos_bnd_driver.F90.

  real(RP), public, allocatable :: ATMOS_BOUNDARY_MFLUX_OFFSET_X(:,:,:)

Referenced by atmos_boundary_driver_setup(), mod_atmos_dyn_driver::atmos_dyn_driver(), and set_offset().

mass

real(rp), allocatable, public mod_atmos_bnd_driver::mass

Definition at line 59 of file mod_atmos_bnd_driver.F90.

flux

real(rp), allocatable, public mod_atmos_bnd_driver::flux

Definition at line 59 of file mod_atmos_bnd_driver.F90.

offset

real(rp), dimension (south, north), allocatable, public mod_atmos_bnd_driver::offset

Definition at line 59 of file mod_atmos_bnd_driver.F90.

atmos_boundary_mflux_offset_y

real(rp), dimension(:,:,:), allocatable, public mod_atmos_bnd_driver::atmos_boundary_mflux_offset_y

Definition at line 60 of file mod_atmos_bnd_driver.F90.

  real(RP), public, allocatable :: ATMOS_BOUNDARY_MFLUX_OFFSET_Y(:,:,:)

Referenced by atmos_boundary_driver_setup(), mod_atmos_dyn_driver::atmos_dyn_driver(), and set_offset().

atmos_boundary_smoother_fact

real(rp), public mod_atmos_bnd_driver::atmos_boundary_smoother_fact = 0.2_RP

Definition at line 62 of file mod_atmos_bnd_driver.F90.

  real(RP), public              :: ATMOS_BOUNDARY_SMOOTHER_FACT  = 0.2_rp

Referenced by atmos_boundary_driver_setup(), and mod_atmos_dyn_driver::atmos_dyn_driver().

fact

real(rp), public mod_atmos_bnd_driver::fact

Definition at line 62 of file mod_atmos_bnd_driver.F90.

smoother

real(rp), public mod_atmos_bnd_driver::smoother

Definition at line 62 of file mod_atmos_bnd_driver.F90.

to

real(rp), public mod_atmos_bnd_driver::to

Definition at line 62 of file mod_atmos_bnd_driver.F90.

atmos_boundary_update_flag

logical, public mod_atmos_bnd_driver::atmos_boundary_update_flag = .false.

Definition at line 64 of file mod_atmos_bnd_driver.F90.

  logical,  public              :: ATMOS_BOUNDARY_UPDATE_FLAG    = .false.

Referenced by atmos_boundary_driver_set(), atmos_boundary_driver_setup(), mod_atmos_driver::atmos_driver_finalize(), and mod_atmos_driver::atmos_driver_update().

switch

logical, public mod_atmos_bnd_driver::switch

Definition at line 64 of file mod_atmos_bnd_driver.F90.

real

mod_atmos_bnd_driver::real

Definition at line 64 of file mod_atmos_bnd_driver.F90.

Referenced by atmos_boundary_driver_setup(), and set_boundary().

case

logical, public mod_atmos_bnd_driver::case

Definition at line 64 of file mod_atmos_bnd_driver.F90.