mod_atmos_phy_ae_vars Module Reference

module ATMOSPHERE / Physics Aerosol Microphysics More...

Functions/Subroutines
subroutine, public	atmos_phy_ae_vars_setup
	Setup. More...
subroutine, public	atmos_phy_ae_vars_finalize
	Finalize. More...
subroutine, public	atmos_phy_ae_vars_fillhalo
	HALO Communication. More...
subroutine, public	atmos_phy_ae_vars_restart_open
	Open restart file for read. More...
subroutine, public	atmos_phy_ae_vars_restart_read
	Read restart. More...
subroutine, public	atmos_phy_ae_vars_restart_create
	Create restart file. More...
subroutine, public	atmos_phy_ae_vars_restart_enddef
	Exit netCDF define mode. More...
subroutine, public	atmos_phy_ae_vars_restart_close
	Close restart file. More...
subroutine, public	atmos_phy_ae_vars_restart_def_var
	Write restart. More...
subroutine, public	atmos_phy_ae_vars_restart_write
	Write restart. More...
subroutine, public	atmos_phy_ae_vars_history (QTRC, RH)
subroutine, public	atmos_phy_ae_vars_get_diagnostic (QTRC, RH, Re, Qe)
subroutine, public	atmos_phy_ae_vars_reset_diagnostics

Variables
logical, public	atmos_phy_ae_restart_output = .false.
	output restart file? More...
character(len=h_long), public	atmos_phy_ae_restart_in_basename = ''
	Basename of the input file. More...
logical, public	atmos_phy_ae_restart_in_aggregate
	Switch to use aggregate file. More...
logical, public	atmos_phy_ae_restart_in_postfix_timelabel = .false.
	Add timelabel to the basename of input file? More...
character(len=h_long), public	atmos_phy_ae_restart_out_basename = ''
	Basename of the output file. More...
logical, public	atmos_phy_ae_restart_out_aggregate
	Switch to use aggregate file. More...
logical, public	atmos_phy_ae_restart_out_postfix_timelabel = .true.
	Add timelabel to the basename of output file? More...
character(len=h_mid), public	atmos_phy_ae_restart_out_title = 'ATMOS_PHY_AE restart'
	title of the output file More...
character(len=h_short), public	atmos_phy_ae_restart_out_dtype = 'DEFAULT'
	REAL4 or REAL8. More...
real(rp), dimension(:,:,:,:), allocatable, public	atmos_phy_ae_rhoq_t
real(rp), dimension(:,:,:), allocatable, public	atmos_phy_ae_ccn
real(rp), dimension(:,:,:), allocatable, public	atmos_phy_ae_ccn_t
real(rp), dimension(:,:,:,:), allocatable, public	atmos_phy_ae_emit
integer, public	qa_ae = 0
integer, public	qs_ae = -1
integer, public	qe_ae = -2

Detailed Description

module ATMOSPHERE / Physics Aerosol Microphysics

Description

Container for mod_atmos_phy_ae

Author

Team SCALE

NAMELIST

• PARAM_ATMOS_PHY_AE_VARS

  name	type	default value	comment
  ATMOS_PHY_AE_RESTART_IN_BASENAME	character(len=H_LONG)	''	Basename of the input file
  ATMOS_PHY_AE_RESTART_IN_AGGREGATE	logical		Switch to use aggregate file
  ATMOS_PHY_AE_RESTART_IN_POSTFIX_TIMELABEL	logical	.false.	Add timelabel to the basename of input file?
  ATMOS_PHY_AE_RESTART_OUTPUT	logical	.false.	output restart file?
  ATMOS_PHY_AE_RESTART_OUT_BASENAME	character(len=H_LONG)	''	Basename of the output file
  ATMOS_PHY_AE_RESTART_OUT_AGGREGATE	logical		Switch to use aggregate file
  ATMOS_PHY_AE_RESTART_OUT_POSTFIX_TIMELABEL	logical	.true.	Add timelabel to the basename of output file?
  ATMOS_PHY_AE_RESTART_OUT_TITLE	character(len=H_MID)	'ATMOS_PHY_AE restart'	title of the output file
  ATMOS_PHY_AE_RESTART_OUT_DTYPE	character(len=H_SHORT)	'DEFAULT'	REAL4 or REAL8

  
  

History Output

name	description	unit	variable
Qe_{AE_NAME}	mass mixing ratio of {AE_NAME}; {AE_NAME} is AD, ASO, AVA, AS, AR, ASS, AU, AT, AOC.	kg/kg	{'Qe_'//trim}
Re_{AE_NAME}	effective radius of {AE_NAME}; {AE_NAME} is AD, ASO, AVA, AS, AR, ASS, AU, AT, AOC.	cm	{'Re_'//trim}

Function/Subroutine Documentation

atmos_phy_ae_vars_setup()

subroutine, public mod_atmos_phy_ae_vars::atmos_phy_ae_vars_setup

Setup.

Definition at line 110 of file mod_atmos_phy_ae_vars.F90.

    use scale_prc, only: &
       prc_abort
    use scale_const, only: &
       undef => const_undef
    use scale_atmos_aerosol, only: &
       n_ae, &
       ae_name, &
       ae_desc
    use scale_file_history, only: &
       file_history_reg
    implicit none
 
    namelist / param_atmos_phy_ae_vars / &
       atmos_phy_ae_restart_in_basename,           &
       atmos_phy_ae_restart_in_aggregate,          &
       atmos_phy_ae_restart_in_postfix_timelabel,  &
       atmos_phy_ae_restart_output,                &
       atmos_phy_ae_restart_out_basename,          &
       atmos_phy_ae_restart_out_aggregate,         &
       atmos_phy_ae_restart_out_postfix_timelabel, &
       atmos_phy_ae_restart_out_title,             &
       atmos_phy_ae_restart_out_dtype
 
    integer :: ierr
    integer :: iv
    !---------------------------------------------------------------------------
 
    log_newline
    log_info("ATMOS_PHY_AE_vars_setup",*) 'Setup'
 
    allocate( atmos_phy_ae_rhoq_t(ka,ia,ja,qs_ae:qe_ae) )
    atmos_phy_ae_rhoq_t(:,:,:, :) = undef
    !$acc enter data create(ATMOS_PHY_AE_RHOQ_t)
 
    allocate( atmos_phy_ae_ccn(ka,ia,ja) )
    atmos_phy_ae_ccn(:,:,:) = undef
    !$acc enter data create(ATMOS_PHY_AE_CCN)
 
    allocate( atmos_phy_ae_ccn_t(ka,ia,ja) )
    atmos_phy_ae_ccn_t(:,:,:) = undef
    !$acc enter data create(ATMOS_PHY_AE_CCN_t)
 
    allocate( atmos_phy_ae_emit(ka,ia,ja,qs_ae:qe_ae) )
    atmos_phy_ae_emit(:,:,:,:) = 0.0_rp
    !$acc enter data copyin(ATMOS_PHY_AE_EMIT)
 
    !--- read namelist
    rewind(io_fid_conf)
    read(io_fid_conf,nml=param_atmos_phy_ae_vars,iostat=ierr)
    if( ierr < 0 ) then !--- missing
       log_info("ATMOS_PHY_AE_vars_setup",*) 'Not found namelist. Default used.'
    elseif( ierr > 0 ) then !--- fatal error
       log_error("ATMOS_PHY_AE_vars_setup",*) 'Not appropriate names in namelist PARAM_ATMOS_PHY_AE_VARS. Check!'
       call prc_abort
    endif
    log_nml(param_atmos_phy_ae_vars)
 
    log_newline
    log_info("ATMOS_PHY_AE_vars_setup",*) '[ATMOS_PHY_AE] prognostic/diagnostic variables'
    log_info_cont('(1x,A,A24,A,A48,A,A12,A)') &
               '      |', 'VARNAME                 ','|', &
               'DESCRIPTION                                     ', '[', 'UNIT        ', ']'
    do iv = 1, vmax
       log_info_cont('(1x,A,I3,A,A24,A,A48,A,A12,A)') &
                  'NO.',iv,'|',var_name(iv),'|',var_desc(iv),'[',var_unit(iv),']'
    enddo
 
    log_newline
    if ( atmos_phy_ae_restart_in_basename /= '' ) then
       log_info("ATMOS_PHY_AE_vars_setup",*) 'Restart input?  : YES, file = ', trim(atmos_phy_ae_restart_in_basename)
       log_info("ATMOS_PHY_AE_vars_setup",*) 'Add timelabel?  : ', atmos_phy_ae_restart_in_postfix_timelabel
    else
       log_info("ATMOS_PHY_AE_vars_setup",*) 'Restart input?  : NO'
    endif
    if (       atmos_phy_ae_restart_output             &
         .AND. atmos_phy_ae_restart_out_basename /= '' ) then
       log_info("ATMOS_PHY_AE_vars_setup",*) 'Restart output? : YES, file = ', trim(atmos_phy_ae_restart_out_basename)
       log_info("ATMOS_PHY_AE_vars_setup",*) 'Add timelabel?  : ', atmos_phy_ae_restart_out_postfix_timelabel
    else
       log_info("ATMOS_PHY_AE_vars_setup",*) 'Restart output? : NO'
       atmos_phy_ae_restart_output = .false.
    endif
 
    ! diagnostices
    allocate( atmos_phy_ae_re(ka,ia,ja,n_ae) )
    allocate( atmos_phy_ae_qe(ka,ia,ja,n_ae) )
!OCL XFILL
    atmos_phy_ae_re(:,:,:,:) = undef
!OCL XFILL
    atmos_phy_ae_qe(:,:,:,:) = undef
    diag_re = .false.
    diag_qe = .false.
    !$acc enter data create(ATMOS_PHY_AE_Re,ATMOS_PHY_AE_Qe)
 
    ! history
    hist_re = .false.
    allocate( hist_re_id(n_ae) )
    do iv = 1, n_ae
       call file_history_reg( 'Re_'//trim(ae_name(iv)), 'effective radius of '//trim(ae_desc(iv)), 'cm', hist_re_id(iv), fill_halo=.true., dim_type='ZXY' )
       if ( hist_re_id(iv) > 0 ) hist_re = .true.
    end do
 
    hist_qe = .false.
    allocate( hist_qe_id(n_ae) )
    do iv = 1, n_ae
       call file_history_reg( 'Qe_'//trim(ae_name(iv)), 'mass mixing ratio of '//trim(ae_desc(iv)), 'kg/kg', hist_qe_id(iv), fill_halo=.true., dim_type='ZXY' )
       if ( hist_qe_id(iv) > 0 ) hist_qe = .true.
    end do
 
    return

References scale_atmos_aerosol::ae_desc, scale_atmos_aerosol::ae_name, atmos_phy_ae_ccn, atmos_phy_ae_ccn_t, atmos_phy_ae_emit, atmos_phy_ae_restart_in_aggregate, atmos_phy_ae_restart_in_basename, atmos_phy_ae_restart_in_postfix_timelabel, atmos_phy_ae_restart_out_aggregate, atmos_phy_ae_restart_out_basename, atmos_phy_ae_restart_out_dtype, atmos_phy_ae_restart_out_postfix_timelabel, atmos_phy_ae_restart_out_title, atmos_phy_ae_restart_output, atmos_phy_ae_rhoq_t, scale_const::const_undef, scale_file_history::file_history_reg(), scale_atmos_grid_cartesc_index::ia, scale_io::io_fid_conf, scale_atmos_grid_cartesc_index::ja, scale_atmos_grid_cartesc_index::ka, scale_atmos_aerosol::n_ae, scale_prc::prc_abort(), qe_ae, and qs_ae.

Referenced by mod_atmos_vars::atmos_vars_setup().

Here is the call graph for this function:

Here is the caller graph for this function:

atmos_phy_ae_vars_finalize()

subroutine, public mod_atmos_phy_ae_vars::atmos_phy_ae_vars_finalize

Finalize.

Definition at line 225 of file mod_atmos_phy_ae_vars.F90.

    implicit none
    !---------------------------------------------------------------------------
 
    log_newline
    log_info("ATMOS_PHY_AE_vars_finalize",*) 'Finalize'
 
    !$acc exit data delete(ATMOS_PHY_AE_RHOQ_t)
    deallocate( atmos_phy_ae_rhoq_t )
 
    !$acc exit data delete(ATMOS_PHY_AE_CCN)
    deallocate( atmos_phy_ae_ccn )
 
    !$acc exit data delete(ATMOS_PHY_AE_CCN_t)
    deallocate( atmos_phy_ae_ccn_t )
 
    !$acc exit data delete(ATMOS_PHY_AE_EMIT)
    deallocate( atmos_phy_ae_emit )
 
    ! diagnostices
    !$acc exit data delete(ATMOS_PHY_AE_Re,ATMOS_PHY_AE_Qe)
    deallocate( atmos_phy_ae_re )
    deallocate( atmos_phy_ae_qe )
 
    ! history
    deallocate( hist_re_id )
 
    deallocate( hist_qe_id )
 
    return

References atmos_phy_ae_ccn, atmos_phy_ae_ccn_t, atmos_phy_ae_emit, and atmos_phy_ae_rhoq_t.

Referenced by mod_atmos_vars::atmos_vars_finalize().

Here is the caller graph for this function:

atmos_phy_ae_vars_fillhalo()

subroutine, public mod_atmos_phy_ae_vars::atmos_phy_ae_vars_fillhalo

HALO Communication.

Definition at line 259 of file mod_atmos_phy_ae_vars.F90.

    use scale_comm_cartesc, only: &
       comm_vars8, &
       comm_wait
    implicit none
 
    integer :: i, j
    !---------------------------------------------------------------------------
 
    !$omp parallel do
    !$acc kernels
    !$acc loop independent collapse(2)
    do j  = js, je
    do i  = is, ie
       atmos_phy_ae_ccn(   1:ks-1,i,j) = atmos_phy_ae_ccn(ks,i,j)
       atmos_phy_ae_ccn(ke+1:ka,  i,j) = atmos_phy_ae_ccn(ke,i,j)
    enddo
    enddo
    !$acc end kernels
 
    call comm_vars8( atmos_phy_ae_ccn(:,:,:), 1 )
    call comm_wait ( atmos_phy_ae_ccn(:,:,:), 1 )
 
    return

References atmos_phy_ae_ccn, 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::ka, scale_atmos_grid_cartesc_index::ke, and scale_atmos_grid_cartesc_index::ks.

Referenced by atmos_phy_ae_vars_restart_read(), and atmos_phy_ae_vars_restart_write().

Here is the caller graph for this function:

atmos_phy_ae_vars_restart_open()

subroutine, public mod_atmos_phy_ae_vars::atmos_phy_ae_vars_restart_open

Open restart file for read.

Definition at line 287 of file mod_atmos_phy_ae_vars.F90.

    use scale_time, only: &
       time_gettimelabel
    use scale_file_cartesc, only: &
       file_cartesc_open
    implicit none
 
    character(len=19)     :: timelabel
    character(len=H_LONG) :: basename
    !---------------------------------------------------------------------------
 
    log_newline
    log_info("ATMOS_PHY_AE_vars_restart_open",*) 'Open restart file (ATMOS_PHY_AE) '
 
    if ( atmos_phy_ae_restart_in_basename /= '' ) then
 
       if ( atmos_phy_ae_restart_in_postfix_timelabel ) then
          call time_gettimelabel( timelabel )
          basename = trim(atmos_phy_ae_restart_in_basename)//'_'//trim(timelabel)
       else
          basename = trim(atmos_phy_ae_restart_in_basename)
       endif
 
       log_info("ATMOS_PHY_AE_vars_restart_open",*) 'basename: ', trim(basename)
 
       call file_cartesc_open( basename, restart_fid, aggregate=atmos_phy_ae_restart_in_aggregate )
    else
       log_info("ATMOS_PHY_AE_vars_restart_open",*) 'restart file for ATMOS_PHY_AE is not specified.'
    endif
 
    return

References atmos_phy_ae_restart_in_aggregate, atmos_phy_ae_restart_in_basename, atmos_phy_ae_restart_in_postfix_timelabel, scale_file_cartesc::file_cartesc_open(), and scale_time::time_gettimelabel().

Referenced by mod_atmos_vars::atmos_vars_restart_open().

Here is the call graph for this function:

Here is the caller graph for this function:

atmos_phy_ae_vars_restart_read()

subroutine, public mod_atmos_phy_ae_vars::atmos_phy_ae_vars_restart_read

Read restart.

Definition at line 322 of file mod_atmos_phy_ae_vars.F90.

    use scale_file, only: &
       file_get_aggregate
    use scale_file_cartesc, only: &
       file_cartesc_read, &
       file_cartesc_flush
    implicit none
 
    integer  :: i, j
    !---------------------------------------------------------------------------
 
    if ( restart_fid /= -1 ) then
       log_newline
       log_info("ATMOS_PHY_AE_vars_restart_read",*) 'Read from restart file (ATMOS_PHY_AE) '
 
       call file_cartesc_read( restart_fid, var_name(1), 'ZXY', & ! [IN]
                               atmos_phy_ae_ccn(:,:,:)          ) ! [OUT]
 
       if ( file_get_aggregate(restart_fid) ) then
          call file_cartesc_flush( restart_fid ) ! X/Y halos have been read from file
          !$acc update device(ATMOS_PHY_AE_CCN)
 
          ! fill K halos
          !$omp parallel do
          !$acc kernels
          !$acc loop independent collapse(2)
          do j  = 1, ja
          do i  = 1, ia
             atmos_phy_ae_ccn(   1:ks-1,i,j) = atmos_phy_ae_ccn(ks,i,j)
             atmos_phy_ae_ccn(ke+1:ka,  i,j) = atmos_phy_ae_ccn(ke,i,j)
          enddo
          enddo
          !$acc end kernels
       else
          call atmos_phy_ae_vars_fillhalo
       end if
 
       call atmos_phy_ae_vars_check
 
    else
       log_info("ATMOS_PHY_AE_vars_restart_read",*) 'invlaid restart file ID for ATMOS_PHY_AE.'
    endif
 
    return

References atmos_phy_ae_ccn, atmos_phy_ae_vars_fillhalo(), scale_file_cartesc::file_cartesc_flush(), scale_file::file_get_aggregate(), scale_atmos_grid_cartesc_index::ia, scale_atmos_grid_cartesc_index::ja, scale_atmos_grid_cartesc_index::ka, scale_atmos_grid_cartesc_index::ke, and scale_atmos_grid_cartesc_index::ks.

Referenced by mod_atmos_vars::atmos_vars_restart_read().

Here is the call graph for this function:

Here is the caller graph for this function:

atmos_phy_ae_vars_restart_create()

subroutine, public mod_atmos_phy_ae_vars::atmos_phy_ae_vars_restart_create

Create restart file.

Definition at line 370 of file mod_atmos_phy_ae_vars.F90.

    use scale_time, only: &
       time_gettimelabel
    use scale_file_cartesc, only: &
       file_cartesc_create
    implicit none
 
    character(len=19)     :: timelabel
    character(len=H_LONG) :: basename
    !---------------------------------------------------------------------------
 
    if ( atmos_phy_ae_restart_out_basename /= '' ) then
 
       log_newline
       log_info("ATMOS_PHY_AE_vars_restart_create",*) 'Create restart file (ATMOS_PHY_AE) '
 
       if ( atmos_phy_ae_restart_out_postfix_timelabel ) then
          call time_gettimelabel( timelabel )
          basename = trim(atmos_phy_ae_restart_out_basename)//'_'//trim(timelabel)
       else
          basename = trim(atmos_phy_ae_restart_out_basename)
       endif
 
       log_info("ATMOS_PHY_AE_vars_restart_create",*) 'basename: ', trim(basename)
 
       call file_cartesc_create( &
            basename, atmos_phy_ae_restart_out_title, atmos_phy_ae_restart_out_dtype, & ! [IN]
            restart_fid,                                                              & ! [OUT]
            aggregate=atmos_phy_ae_restart_out_aggregate                              ) ! [IN]
 
    endif
 
    return

References atmos_phy_ae_restart_out_aggregate, atmos_phy_ae_restart_out_basename, atmos_phy_ae_restart_out_dtype, atmos_phy_ae_restart_out_postfix_timelabel, atmos_phy_ae_restart_out_title, scale_file_cartesc::file_cartesc_create(), and scale_time::time_gettimelabel().

Referenced by mod_atmos_vars::atmos_vars_restart_create().

Here is the call graph for this function:

Here is the caller graph for this function:

atmos_phy_ae_vars_restart_enddef()

subroutine, public mod_atmos_phy_ae_vars::atmos_phy_ae_vars_restart_enddef

Exit netCDF define mode.

Definition at line 407 of file mod_atmos_phy_ae_vars.F90.

    use scale_file_cartesc, only: &
       file_cartesc_enddef
    implicit none
 
    if ( restart_fid /= -1 ) then
       call file_cartesc_enddef( restart_fid ) ! [IN]
    endif
 
    return

References scale_file_cartesc::file_cartesc_enddef().

Referenced by mod_atmos_vars::atmos_vars_restart_enddef().

Here is the call graph for this function:

Here is the caller graph for this function:

atmos_phy_ae_vars_restart_close()

subroutine, public mod_atmos_phy_ae_vars::atmos_phy_ae_vars_restart_close

Close restart file.

Definition at line 421 of file mod_atmos_phy_ae_vars.F90.

    use scale_file_cartesc, only: &
       file_cartesc_close
    implicit none
    !---------------------------------------------------------------------------
 
    if ( restart_fid /= -1 ) then
       log_newline
       log_info("ATMOS_PHY_AE_vars_restart_close",*) 'Close restart file (ATMOS_PHY_AE) '
 
       call file_cartesc_close( restart_fid ) ! [IN]
 
       restart_fid = -1
    endif
 
    return

References scale_file_cartesc::file_cartesc_close().

Referenced by mod_atmos_vars::atmos_vars_restart_close().

Here is the call graph for this function:

Here is the caller graph for this function:

atmos_phy_ae_vars_restart_def_var()

subroutine, public mod_atmos_phy_ae_vars::atmos_phy_ae_vars_restart_def_var

Write restart.

Definition at line 441 of file mod_atmos_phy_ae_vars.F90.

    use scale_file_cartesc, only: &
       file_cartesc_def_var
    implicit none
    !---------------------------------------------------------------------------
 
    if ( restart_fid /= -1 ) then
       call file_cartesc_def_var( restart_fid, var_name(1), var_desc(1), var_unit(1), 'ZXY', atmos_phy_ae_restart_out_dtype, &
                                  var_id(1) )
    endif
 
    return

References atmos_phy_ae_restart_out_dtype, and scale_file_cartesc::file_cartesc_def_var().

Referenced by mod_atmos_vars::atmos_vars_restart_def_var().

Here is the call graph for this function:

Here is the caller graph for this function:

atmos_phy_ae_vars_restart_write()

subroutine, public mod_atmos_phy_ae_vars::atmos_phy_ae_vars_restart_write

Write restart.

Definition at line 457 of file mod_atmos_phy_ae_vars.F90.

    use scale_file_cartesc, only: &
       file_cartesc_write_var
    implicit none
    !---------------------------------------------------------------------------
 
    if ( restart_fid /= -1 ) then
 
       call atmos_phy_ae_vars_fillhalo
 
       call atmos_phy_ae_vars_check
 
       call file_cartesc_write_var( restart_fid, var_id(1), atmos_phy_ae_ccn(:,:,:), var_name(1), 'ZXY' ) ! [IN]
 
    endif
 
    return

References atmos_phy_ae_ccn, and atmos_phy_ae_vars_fillhalo().

Referenced by mod_atmos_vars::atmos_vars_restart_write().

Here is the call graph for this function:

Here is the caller graph for this function:

atmos_phy_ae_vars_history()

subroutine, public mod_atmos_phy_ae_vars::atmos_phy_ae_vars_history	(	real(rp), dimension(ka,ia,ja,qa), intent(in)	QTRC,
		real(rp), dimension(ka,ia,ja), intent(in)	RH
	)

Definition at line 477 of file mod_atmos_phy_ae_vars.F90.

    use scale_tracer, only: &
       qa
    use scale_atmos_aerosol, only: &
       n_ae
    use scale_file_history, only: &
       file_history_put
    real(RP), intent(in) :: QTRC(KA,IA,JA,QA)
    real(RP), intent(in) :: RH(KA,IA,JA)
 
    real(RP) :: WORK(KA,IA,JA,N_AE)
    integer  :: iv
 
    if ( hist_re ) then
       call atmos_phy_ae_vars_get_diagnostic( &
            qtrc(:,:,:,:), rh(:,:,:), & ! [IN]
            re=work(:,:,:,:)          ) ! [OUT]
       do iv = 1, n_ae
          if ( hist_re_id(iv) > 0 ) &
               call file_history_put( hist_re_id(iv), work(:,:,:,iv) )
       end do
    end if
 
    if ( hist_qe ) then
       call atmos_phy_ae_vars_get_diagnostic( &
            qtrc(:,:,:,:), rh(:,:,:), & ! [IN]
            qe=work(:,:,:,:)          ) ! [OUT]
       do iv = 1, n_ae
          if ( hist_qe_id(iv) > 0 ) &
               call file_history_put( hist_qe_id(iv), work(:,:,:,iv) )
       end do
    end if
 
    return

References atmos_phy_ae_vars_get_diagnostic(), scale_atmos_aerosol::n_ae, and scale_tracer::qa.

Referenced by mod_atmos_vars::atmos_vars_history().

Here is the call graph for this function:

Here is the caller graph for this function:

atmos_phy_ae_vars_get_diagnostic()

subroutine, public mod_atmos_phy_ae_vars::atmos_phy_ae_vars_get_diagnostic	(	real(rp), dimension(ka,ia,ja,qa), intent(in)	QTRC,
		real(rp), dimension(ka,ia,ja), intent(in)	RH,
		real(rp), dimension(ka,ia,ja,n_ae), intent(out), optional	Re,
		real(rp), dimension(ka,ia,ja,n_ae), intent(out), optional	Qe
	)

Definition at line 515 of file mod_atmos_phy_ae_vars.F90.

    use scale_time, only: &
       time_nowdaysec
    use scale_tracer, only: &
       qa
    use scale_atmos_aerosol, only: &
       n_ae
    use scale_atmos_phy_ae_kajino13, only: &
       atmos_phy_ae_kajino13_effective_radius
    use scale_atmos_phy_ae_offline, only: &
       atmos_phy_ae_offline_effective_radius, &
       atmos_phy_ae_offline_qtrc2qaero
    use mod_atmos_admin, only: &
       atmos_phy_ae_type
 
    real(RP), intent(in) :: QTRC(KA,IA,JA,QA)
    real(RP), intent(in) :: RH(KA,IA,JA)
    real(RP), intent(out), optional :: Re(KA,IA,JA,N_AE)
    real(RP), intent(out), optional :: Qe(KA,IA,JA,N_AE)
 
    if ( present(re) ) then
       !$acc data copyout(Re)
       if ( .not. diag_re ) then
 
          select case ( atmos_phy_ae_type )
          case ( 'KAJINO13' )
             !$acc update host(QTRC(:,:,:,QS_AE:QE_AE),RH)
             call atmos_phy_ae_kajino13_effective_radius( &
                     ka, ia, ja, qa_ae, &
                     qtrc(:,:,:,qs_ae:qe_ae), rh(:,:,:), & ! [IN]
                     atmos_phy_ae_re(:,:,:,:)            ) ! [OUT]
             !$acc update device(ATMOS_PHY_AE_Re)
          case ( 'OFFLINE' )
             !$acc update host(RH)
             call atmos_phy_ae_offline_effective_radius( &
                     ka, ia, ja,              &
                     rh(:,:,:),               & ! [IN]
                     atmos_phy_ae_re(:,:,:,:) ) ! [OUT]
             !$acc update device(ATMOS_PHY_AE_Re)
          case default
             !$acc kernels
             atmos_phy_ae_re(:,:,:,:) = 0.0_rp
             !$acc end kernels
          end select
 
          diag_re = .true.
       end if
!OCL XFILL
       !$acc kernels
       re(:,:,:,:) = atmos_phy_ae_re(:,:,:,:)
       !$acc end kernels
       !$acc end data
    end if
 
    if ( present(qe) ) then
       !$acc data copyout(Qe)
       if ( .not. diag_qe ) then
 
          select case ( atmos_phy_ae_type )
          case ( 'OFFLINE' )
             call atmos_phy_ae_offline_qtrc2qaero( &
                     ka, ia, ja,              &
                     time_nowdaysec,          &
                     atmos_phy_ae_qe(:,:,:,:) ) ! [OUT]
             !$acc update device(ATMOS_PHY_AE_Qe)
          case default
             !$acc kernels
             atmos_phy_ae_qe(:,:,:,:) = 0.0_rp
             !$acc end kernels
          end select
 
          diag_qe = .true.
       end if
!OCL XFILL
       !$acc kernels
       qe(:,:,:,:) = atmos_phy_ae_qe(:,:,:,:)
       !$acc end kernels
       !$acc end data
    end if
 
    return

References scale_atmos_phy_ae_kajino13::atmos_phy_ae_kajino13_effective_radius(), scale_atmos_phy_ae_offline::atmos_phy_ae_offline_effective_radius(), scale_atmos_phy_ae_offline::atmos_phy_ae_offline_qtrc2qaero(), mod_atmos_admin::atmos_phy_ae_type, scale_atmos_grid_cartesc_index::ia, scale_atmos_grid_cartesc_index::ja, scale_atmos_grid_cartesc_index::ka, scale_atmos_aerosol::n_ae, scale_tracer::qa, qa_ae, qe_ae, qs_ae, and scale_time::time_nowdaysec.

Referenced by atmos_phy_ae_vars_history(), and mod_atmos_phy_rd_driver::atmos_phy_rd_driver_calc_tendency().

Here is the call graph for this function:

Here is the caller graph for this function:

atmos_phy_ae_vars_reset_diagnostics()

subroutine, public mod_atmos_phy_ae_vars::atmos_phy_ae_vars_reset_diagnostics

Definition at line 598 of file mod_atmos_phy_ae_vars.F90.

    diag_re      = .false.
    diag_qe      = .false.
 
    return

References scale_atmos_grid_cartesc_real::atmos_grid_cartesc_real_totvol, scale_atmos_grid_cartesc_real::atmos_grid_cartesc_real_vol, atmos_phy_ae_ccn, scale_atmos_grid_cartesc_index::ia, scale_atmos_grid_cartesc_index::ie, scale_atmos_grid_cartesc_index::is, scale_atmos_grid_cartesc_index::ja, scale_atmos_grid_cartesc_index::je, scale_atmos_grid_cartesc_index::js, scale_atmos_grid_cartesc_index::ka, scale_atmos_grid_cartesc_index::ke, scale_atmos_grid_cartesc_index::ks, and scale_statistics::statistics_checktotal.

Referenced by mod_atmos_vars::atmos_vars_calc_diagnostics().

Here is the caller graph for this function:

Variable Documentation

atmos_phy_ae_restart_output

logical, public mod_atmos_phy_ae_vars::atmos_phy_ae_restart_output = .false.

output restart file?

Definition at line 50 of file mod_atmos_phy_ae_vars.F90.

  logical,                public :: ATMOS_PHY_AE_RESTART_OUTPUT                = .false.                

Referenced by mod_admin_restart::admin_restart_setup(), and atmos_phy_ae_vars_setup().

atmos_phy_ae_restart_in_basename

character(len=h_long), public mod_atmos_phy_ae_vars::atmos_phy_ae_restart_in_basename = ''

Basename of the input file.

Definition at line 52 of file mod_atmos_phy_ae_vars.F90.

  character(len=H_LONG),  public :: ATMOS_PHY_AE_RESTART_IN_BASENAME           = ''

Referenced by mod_admin_restart::admin_restart_setup(), atmos_phy_ae_vars_restart_open(), and atmos_phy_ae_vars_setup().

atmos_phy_ae_restart_in_aggregate

logical, public mod_atmos_phy_ae_vars::atmos_phy_ae_restart_in_aggregate

Switch to use aggregate file.

Definition at line 53 of file mod_atmos_phy_ae_vars.F90.

  logical,                public :: ATMOS_PHY_AE_RESTART_IN_AGGREGATE

Referenced by mod_admin_restart::admin_restart_setup(), atmos_phy_ae_vars_restart_open(), and atmos_phy_ae_vars_setup().

atmos_phy_ae_restart_in_postfix_timelabel

logical, public mod_atmos_phy_ae_vars::atmos_phy_ae_restart_in_postfix_timelabel = .false.

Add timelabel to the basename of input file?

Definition at line 54 of file mod_atmos_phy_ae_vars.F90.

  logical,                public :: ATMOS_PHY_AE_RESTART_IN_POSTFIX_TIMELABEL  = .false.                

Referenced by mod_admin_restart::admin_restart_setup(), atmos_phy_ae_vars_restart_open(), and atmos_phy_ae_vars_setup().

atmos_phy_ae_restart_out_basename

character(len=h_long), public mod_atmos_phy_ae_vars::atmos_phy_ae_restart_out_basename = ''

Basename of the output file.

Definition at line 55 of file mod_atmos_phy_ae_vars.F90.

  character(len=H_LONG),  public :: ATMOS_PHY_AE_RESTART_OUT_BASENAME          = ''

Referenced by mod_admin_restart::admin_restart_setup(), atmos_phy_ae_vars_restart_create(), and atmos_phy_ae_vars_setup().

atmos_phy_ae_restart_out_aggregate

logical, public mod_atmos_phy_ae_vars::atmos_phy_ae_restart_out_aggregate

Switch to use aggregate file.

Definition at line 56 of file mod_atmos_phy_ae_vars.F90.

  logical,                public :: ATMOS_PHY_AE_RESTART_OUT_AGGREGATE

Referenced by mod_admin_restart::admin_restart_setup(), atmos_phy_ae_vars_restart_create(), and atmos_phy_ae_vars_setup().

atmos_phy_ae_restart_out_postfix_timelabel

logical, public mod_atmos_phy_ae_vars::atmos_phy_ae_restart_out_postfix_timelabel = .true.

Add timelabel to the basename of output file?

Definition at line 57 of file mod_atmos_phy_ae_vars.F90.

  logical,                public :: ATMOS_PHY_AE_RESTART_OUT_POSTFIX_TIMELABEL = .true.                 

Referenced by mod_admin_restart::admin_restart_setup(), atmos_phy_ae_vars_restart_create(), and atmos_phy_ae_vars_setup().

atmos_phy_ae_restart_out_title

character(len=h_mid), public mod_atmos_phy_ae_vars::atmos_phy_ae_restart_out_title = 'ATMOS_PHY_AE restart'

title of the output file

Definition at line 58 of file mod_atmos_phy_ae_vars.F90.

  character(len=H_MID),   public :: ATMOS_PHY_AE_RESTART_OUT_TITLE             = 'ATMOS_PHY_AE restart'

Referenced by mod_admin_restart::admin_restart_setup(), atmos_phy_ae_vars_restart_create(), and atmos_phy_ae_vars_setup().

atmos_phy_ae_restart_out_dtype

character(len=h_short), public mod_atmos_phy_ae_vars::atmos_phy_ae_restart_out_dtype = 'DEFAULT'

REAL4 or REAL8.

Definition at line 59 of file mod_atmos_phy_ae_vars.F90.

  character(len=H_SHORT), public :: ATMOS_PHY_AE_RESTART_OUT_DTYPE             = 'DEFAULT'

Referenced by mod_admin_restart::admin_restart_setup(), atmos_phy_ae_vars_restart_create(), atmos_phy_ae_vars_restart_def_var(), and atmos_phy_ae_vars_setup().

atmos_phy_ae_rhoq_t

real(rp), dimension(:,:,:,:), allocatable, public mod_atmos_phy_ae_vars::atmos_phy_ae_rhoq_t

Definition at line 61 of file mod_atmos_phy_ae_vars.F90.

  real(RP), public, allocatable :: ATMOS_PHY_AE_RHOQ_t(:,:,:,:)                               ! tendency rho*QTRC [kg/kg/s]

Referenced by mod_atmos_phy_ae_driver::atmos_phy_ae_driver_calc_tendency(), atmos_phy_ae_vars_finalize(), and atmos_phy_ae_vars_setup().

atmos_phy_ae_ccn

real(rp), dimension(:,:,:), allocatable, public mod_atmos_phy_ae_vars::atmos_phy_ae_ccn

Definition at line 63 of file mod_atmos_phy_ae_vars.F90.

  real(RP), public, allocatable :: ATMOS_PHY_AE_CCN(:,:,:)                                    ! cloud condensation nuclei [/m3]

Referenced by mod_atmos_phy_ae_driver::atmos_phy_ae_driver_calc_tendency(), atmos_phy_ae_vars_fillhalo(), atmos_phy_ae_vars_finalize(), atmos_phy_ae_vars_reset_diagnostics(), atmos_phy_ae_vars_restart_read(), atmos_phy_ae_vars_restart_write(), and atmos_phy_ae_vars_setup().

atmos_phy_ae_ccn_t

real(rp), dimension(:,:,:), allocatable, public mod_atmos_phy_ae_vars::atmos_phy_ae_ccn_t

Definition at line 64 of file mod_atmos_phy_ae_vars.F90.

  real(RP), public, allocatable :: ATMOS_PHY_AE_CCN_t(:,:,:)                                  ! tendency CCN [/m3/s]

Referenced by mod_atmos_phy_ae_driver::atmos_phy_ae_driver_calc_tendency(), atmos_phy_ae_vars_finalize(), atmos_phy_ae_vars_setup(), and mod_atmos_phy_mp_driver::atmos_phy_mp_driver_calc_tendency().

atmos_phy_ae_emit

real(rp), dimension(:,:,:,:), allocatable, public mod_atmos_phy_ae_vars::atmos_phy_ae_emit

Definition at line 65 of file mod_atmos_phy_ae_vars.F90.

  real(RP), public, allocatable :: ATMOS_PHY_AE_EMIT(:,:,:,:)                                 ! emission of aerosol and gas

Referenced by mod_atmos_phy_ae_driver::atmos_phy_ae_driver_calc_tendency(), atmos_phy_ae_vars_finalize(), and atmos_phy_ae_vars_setup().

qa_ae

integer, public mod_atmos_phy_ae_vars::qa_ae = 0

Definition at line 67 of file mod_atmos_phy_ae_vars.F90.

  integer, public :: QA_AE = 0

Referenced by mod_atmos_phy_ae_driver::atmos_phy_ae_driver_adjustment(), mod_atmos_phy_ae_driver::atmos_phy_ae_driver_calc_tendency(), mod_atmos_phy_ae_driver::atmos_phy_ae_driver_tracer_setup(), atmos_phy_ae_vars_get_diagnostic(), and mod_mkinit::rect_setup().

qs_ae

integer, public mod_atmos_phy_ae_vars::qs_ae = -1

Definition at line 68 of file mod_atmos_phy_ae_vars.F90.

  integer, public :: QS_AE = -1

Referenced by mod_atmos_phy_ae_driver::atmos_phy_ae_driver_adjustment(), mod_atmos_phy_ae_driver::atmos_phy_ae_driver_calc_tendency(), mod_atmos_phy_ae_driver::atmos_phy_ae_driver_tracer_setup(), atmos_phy_ae_vars_get_diagnostic(), atmos_phy_ae_vars_setup(), and mod_mkinit::rect_setup().

qe_ae

integer, public mod_atmos_phy_ae_vars::qe_ae = -2

Definition at line 69 of file mod_atmos_phy_ae_vars.F90.

  integer, public :: QE_AE = -2

Referenced by mod_atmos_phy_ae_driver::atmos_phy_ae_driver_adjustment(), mod_atmos_phy_ae_driver::atmos_phy_ae_driver_calc_tendency(), mod_atmos_phy_ae_driver::atmos_phy_ae_driver_tracer_setup(), atmos_phy_ae_vars_get_diagnostic(), atmos_phy_ae_vars_setup(), and mod_mkinit::rect_setup().