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_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...

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_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_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

Detailed Description

module ATMOSPHERE / Physics Aerosol Microphysics

Description

Container for mod_atmos_phy_ae

Author

Team SCALE

History

• 2014-05-04 (H.Yashiro) [new]

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_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_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

No history output

Function/Subroutine Documentation

atmos_phy_ae_vars_setup()

subroutine, public mod_atmos_phy_ae_vars::atmos_phy_ae_vars_setup

(

)

Setup.

Definition at line 86 of file mod_atmos_phy_ae_vars.f90.

References atmos_phy_ae_ccn, atmos_phy_ae_ccn_t, atmos_phy_ae_emit, atmos_phy_ae_restart_in_basename, atmos_phy_ae_restart_in_postfix_timelabel, 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_grid_index::ia, scale_stdio::io_fid_conf, scale_stdio::io_fid_log, scale_stdio::io_fid_nml, scale_stdio::io_l, scale_stdio::io_nml, scale_grid_index::ja, scale_grid_index::ka, scale_process::prc_mpistop(), scale_atmos_phy_ae::qa_ae, scale_atmos_phy_ae::qe_ae, and scale_atmos_phy_ae::qs_ae.

Referenced by mod_atmos_vars::atmos_vars_setup().

    use scale_process, only: &
       prc_mpistop
    use scale_const, only: &
       undef => const_undef
    use scale_atmos_phy_ae, only: &
       qa_ae, &
       qs_ae, &
       qe_ae
    implicit none

    namelist / param_atmos_phy_ae_vars / &
       atmos_phy_ae_restart_in_basename,           &
       atmos_phy_ae_restart_in_postfix_timelabel,  &
       atmos_phy_ae_restart_output,                &
       atmos_phy_ae_restart_out_basename,          &
       atmos_phy_ae_restart_out_postfix_timelabel, &
       atmos_phy_ae_restart_out_title,             &
       atmos_phy_ae_restart_out_dtype

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

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

    allocate( atmos_phy_ae_rhoq_t(ka,ia,ja,qs_ae:qe_ae) )
    atmos_phy_ae_rhoq_t(:,:,:, :) = undef

    allocate( atmos_phy_ae_ccn(ka,ia,ja) )
    atmos_phy_ae_ccn(:,:,:) = undef

    allocate( atmos_phy_ae_ccn_t(ka,ia,ja) )
    atmos_phy_ae_ccn_t(:,:,:) = undef

    allocate( atmos_phy_ae_emit(ka,ia,ja,qa_ae) )
    atmos_phy_ae_emit(:,:,:,:) = 0.0_rp

    !--- read namelist
    rewind(io_fid_conf)
    read(io_fid_conf,nml=param_atmos_phy_ae_vars,iostat=ierr)
    if( ierr < 0 ) then !--- missing
       if( io_l ) write(io_fid_log,*) '*** Not found namelist. Default used.'
    elseif( ierr > 0 ) then !--- fatal error
       write(*,*) 'xxx Not appropriate names in namelist PARAM_ATMOS_PHY_AE_VARS. Check!'
       call prc_mpistop
    endif
    if( io_nml ) write(io_fid_nml,nml=param_atmos_phy_ae_vars)

    if( io_l ) write(io_fid_log,*)
    if( io_l ) write(io_fid_log,*) '*** [ATMOS_PHY_AE] prognostic/diagnostic variables'
    if( io_l ) write(io_fid_log,'(1x,A,A24,A,A48,A,A12,A)') &
               '***       |', 'VARNAME                 ','|', &
               'DESCRIPTION                                     ', '[', 'UNIT        ', ']'
    do iv = 1, vmax
       if( io_l ) write(io_fid_log,'(1x,A,I3,A,A24,A,A48,A,A12,A)') &
                  '*** NO.',iv,'|',var_name(iv),'|',var_desc(iv),'[',var_unit(iv),']'
    enddo

    if( io_l ) write(io_fid_log,*)
    if ( atmos_phy_ae_restart_in_basename /= '' ) then
       if( io_l ) write(io_fid_log,*) '*** Restart input?  : YES, file = ', trim(atmos_phy_ae_restart_in_basename)
       if( io_l ) write(io_fid_log,*) '*** Add timelabel?  : ', atmos_phy_ae_restart_in_postfix_timelabel
    else
       if( io_l ) write(io_fid_log,*) '*** Restart input?  : NO'
    endif
    if (       atmos_phy_ae_restart_output             &
         .AND. atmos_phy_ae_restart_out_basename /= '' ) then
       if( io_l ) write(io_fid_log,*) '*** Restart output? : YES, file = ', trim(atmos_phy_ae_restart_out_basename)
       if( io_l ) write(io_fid_log,*) '*** Add timelabel?  : ', atmos_phy_ae_restart_out_postfix_timelabel
    else
       if( io_l ) write(io_fid_log,*) '*** Restart output? : NO'
       atmos_phy_ae_restart_output = .false.
    endif

    return

Here is the call graph for this function:

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 167 of file mod_atmos_phy_ae_vars.f90.

References atmos_phy_ae_ccn, scale_grid_index::ie, scale_grid_index::is, scale_grid_index::je, scale_grid_index::js, scale_grid_index::ka, scale_grid_index::ke, and scale_grid_index::ks.

Referenced by atmos_phy_ae_vars_restart_read(), and atmos_phy_ae_vars_restart_write().

    use scale_comm, only: &
       comm_vars8, &
       comm_wait
    implicit none

    integer :: i, j
    !---------------------------------------------------------------------------

    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

    call comm_vars8( atmos_phy_ae_ccn(:,:,:), 1 )
    call comm_wait ( atmos_phy_ae_ccn(:,:,:), 1 )

    return

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 191 of file mod_atmos_phy_ae_vars.f90.

References atmos_phy_ae_restart_in_basename, atmos_phy_ae_restart_in_postfix_timelabel, scale_fileio::fileio_open(), scale_stdio::io_fid_log, scale_stdio::io_l, and scale_time::time_gettimelabel().

Referenced by mod_atmos_vars::atmos_vars_restart_open().

    use scale_time, only: &
       time_gettimelabel
    use scale_fileio, only: &
       fileio_open
    implicit none

    character(len=19)     :: timelabel
    character(len=H_LONG) :: basename
    !---------------------------------------------------------------------------

    if( io_l ) write(io_fid_log,*)
    if( io_l ) write(io_fid_log,*) '*** 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

       if( io_l ) write(io_fid_log,*) '*** basename: ', trim(basename)

       call fileio_open( restart_fid, basename )
    else
       if( io_l ) write(io_fid_log,*) '*** restart file for ATMOS_PHY_AE is not specified.'
    endif

    return

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 226 of file mod_atmos_phy_ae_vars.f90.

References atmos_phy_ae_ccn, atmos_phy_ae_vars_fillhalo(), scale_fileio::fileio_flush(), scale_grid_index::ia, scale_stdio::io_aggregate, scale_stdio::io_fid_log, scale_stdio::io_l, scale_grid_index::ja, scale_grid_index::ka, scale_grid_index::ke, scale_grid_index::ks, and scale_rm_statistics::statistics_checktotal.

Referenced by mod_atmos_vars::atmos_vars_restart_read().

    use scale_rm_statistics, only: &
       statistics_checktotal, &
       stat_total
    use scale_fileio, only: &
       fileio_read, &
       fileio_flush
    implicit none

    real(RP) :: total
    integer  :: i, j
    !---------------------------------------------------------------------------

    if ( restart_fid /= -1 ) then
       if( io_l ) write(io_fid_log,*)
       if( io_l ) write(io_fid_log,*) '*** Read from restart file (ATMOS_PHY_AE) ***'

       call fileio_read( atmos_phy_ae_ccn(:,:,:),                & ! [OUT]
                         restart_fid, var_name(1), 'ZXY', step=1 ) ! [IN]

       if ( io_aggregate ) then
          call fileio_flush( restart_fid ) ! X/Y halos have been read from file

          ! fill K halos
          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
       else
          call atmos_phy_ae_vars_fillhalo
       end if

       if ( statistics_checktotal ) then
          call stat_total( total, atmos_phy_ae_ccn(:,:,:), var_name(1) )
       end if
    else
       if( io_l ) write(io_fid_log,*) '*** invlaid restart file ID for ATMOS_PHY_AE.'
    endif

    return

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 272 of file mod_atmos_phy_ae_vars.f90.

References 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_fileio::fileio_create(), scale_stdio::io_fid_log, scale_stdio::io_l, and scale_time::time_gettimelabel().

Referenced by mod_atmos_vars::atmos_vars_restart_create().

    use scale_time, only: &
       time_gettimelabel
    use scale_fileio, only: &
       fileio_create
    implicit none

    character(len=19)     :: timelabel
    character(len=H_LONG) :: basename
    !---------------------------------------------------------------------------

    if ( atmos_phy_ae_restart_out_basename /= '' ) then

       if( io_l ) write(io_fid_log,*)
       if( io_l ) write(io_fid_log,*) '*** 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

       if( io_l ) write(io_fid_log,*) '*** basename: ', trim(basename)

       call fileio_create( restart_fid,                                                             & ! [OUT]
                           basename, atmos_phy_ae_restart_out_title, atmos_phy_ae_restart_out_dtype ) ! [IN]

    endif

    return

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 307 of file mod_atmos_phy_ae_vars.f90.

References scale_fileio::fileio_enddef().

Referenced by mod_atmos_vars::atmos_vars_restart_enddef().

    use scale_fileio, only: &
       fileio_enddef
    implicit none

    if ( restart_fid /= -1 ) then
       call fileio_enddef( restart_fid ) ! [IN]
    endif

    return

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 321 of file mod_atmos_phy_ae_vars.f90.

References scale_fileio::fileio_close(), scale_stdio::io_fid_log, and scale_stdio::io_l.

Referenced by mod_atmos_vars::atmos_vars_restart_close().

    use scale_fileio, only: &
       fileio_close
    implicit none
    !---------------------------------------------------------------------------

    if ( restart_fid /= -1 ) then
       if( io_l ) write(io_fid_log,*)
       if( io_l ) write(io_fid_log,*) '*** Close restart file (ATMOS_PHY_AE) ***'

       call fileio_close( restart_fid ) ! [IN]

       restart_fid = -1
    endif

    return

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 341 of file mod_atmos_phy_ae_vars.f90.

References atmos_phy_ae_restart_out_dtype, and scale_fileio::fileio_def_var().

Referenced by mod_atmos_vars::atmos_vars_restart_def_var().

    use scale_fileio, only: &
       fileio_def_var
    implicit none
    !---------------------------------------------------------------------------

    if ( restart_fid /= -1 ) then
       call fileio_def_var( restart_fid, var_id(1), var_name(1), var_desc(1), var_unit(1), &
                            'ZXY', atmos_phy_ae_restart_out_dtype  ) ! [IN]
    endif

    return

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 357 of file mod_atmos_phy_ae_vars.f90.

References atmos_phy_ae_ccn, atmos_phy_ae_vars_fillhalo(), and scale_rm_statistics::statistics_checktotal.

Referenced by mod_atmos_vars::atmos_vars_restart_write().

    use scale_rm_statistics, only: &
       statistics_checktotal, &
       stat_total
    use scale_fileio, only: &
       fileio_write_var
    implicit none

    real(RP) :: total
    !---------------------------------------------------------------------------

    if ( restart_fid /= -1 ) then

       call atmos_phy_ae_vars_fillhalo

       if ( statistics_checktotal ) then
          call stat_total( total, atmos_phy_ae_ccn(:,:,:), var_name(1) )
       end if

       call fileio_write_var( restart_fid, var_id(1), atmos_phy_ae_ccn(:,:,:), &
                              var_name(1), 'ZXY' ) ! [IN]

    endif

    return

Here is the call graph for this function:

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 46 of file mod_atmos_phy_ae_vars.f90.

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

  logical,               public :: ATMOS_PHY_AE_RESTART_OUTPUT                = .false.                

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 48 of file mod_atmos_phy_ae_vars.f90.

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

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

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 49 of file mod_atmos_phy_ae_vars.f90.

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

  logical,                public :: ATMOS_PHY_AE_RESTART_IN_POSTFIX_TIMELABEL  = .false.                

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 50 of file mod_atmos_phy_ae_vars.f90.

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

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

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 51 of file mod_atmos_phy_ae_vars.f90.

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

  logical,                public :: ATMOS_PHY_AE_RESTART_OUT_POSTFIX_TIMELABEL = .true.                 

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 52 of file mod_atmos_phy_ae_vars.f90.

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

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

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 53 of file mod_atmos_phy_ae_vars.f90.

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().

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

atmos_phy_ae_rhoq_t

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

Definition at line 55 of file mod_atmos_phy_ae_vars.f90.

Referenced by mod_atmos_phy_ae_driver::atmos_phy_ae_driver(), and atmos_phy_ae_vars_setup().

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

atmos_phy_ae_ccn

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

Definition at line 57 of file mod_atmos_phy_ae_vars.f90.

Referenced by mod_atmos_phy_ae_driver::atmos_phy_ae_driver(), atmos_phy_ae_vars_fillhalo(), atmos_phy_ae_vars_restart_read(), atmos_phy_ae_vars_restart_write(), and atmos_phy_ae_vars_setup().

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

atmos_phy_ae_ccn_t

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

Definition at line 58 of file mod_atmos_phy_ae_vars.f90.

Referenced by mod_atmos_phy_ae_driver::atmos_phy_ae_driver(), atmos_phy_ae_vars_setup(), and mod_atmos_phy_mp_driver::atmos_phy_mp_driver().

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

atmos_phy_ae_emit

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

Definition at line 59 of file mod_atmos_phy_ae_vars.f90.

Referenced by mod_atmos_phy_ae_driver::atmos_phy_ae_driver(), and atmos_phy_ae_vars_setup().

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