module TRACER More...

Functions/Subroutines
subroutine, public	tracer_regist (QS, NQ, NAME, DESC, UNIT, CV, CP, R, ADVC, MASS)
	Regist tracer. More...

Variables
integer, public	qa = 0

real(rp), dimension(qa_max), public	tracer_cp

real(rp), dimension(qa_max), public	tracer_cv

real(rp), dimension(qa_max), public	tracer_r

real(rp), dimension(qa_max), public	tracer_mass

logical, dimension(qa_max), public	tracer_advc

character(len=h_short), dimension(qa_max), public	tracer_name

character(len=h_mid), dimension(qa_max), public	tracer_desc

character(len=h_short), dimension(qa_max), public	tracer_unit

Detailed Description

module TRACER

Description: Tracer module

Author: Team SCALE

History

2013-12-04 (S.Nishizawa) [new]
2016-08-02 (S.Nishizawa) [mod] add register

Function/Subroutine Documentation

◆ tracer_regist()

subroutine, public scale_tracer::tracer_regist	(	integer, intent(out)	QS,
		integer, intent(in)	NQ,
		character(len=*), dimension(nq), intent(in)	NAME,
		character(len=*), dimension(nq), intent(in)	DESC,
		character(len=*), dimension(nq), intent(in)	UNIT,
		real(rp), dimension (nq), intent(in), optional	CV,
		real(rp), dimension (nq), intent(in), optional	CP,
		real(rp), dimension (nq), intent(in), optional	R,
		logical, dimension(nq), intent(in), optional	ADVC,
		logical, dimension(nq), intent(in), optional	MASS
	)

Regist tracer.

Parameters

[in]	advc	if .true., the tracer is advected in the dynamical process. (default is .true.)
[in]	mass	if .true., the tracer has mass. (default is .false.)

Definition at line 66 of file scale_tracer.F90.

References scale_prc::prc_abort(), qa, tracer_advc, tracer_cp, tracer_cv, tracer_desc, tracer_mass, tracer_name, tracer_r, and tracer_unit.

Referenced by scale_atmos_hydrometeor::atmos_hydrometeor_regist(), mod_atmos_phy_ae_driver::atmos_phy_ae_driver_tracer_setup(), mod_atmos_phy_bl_driver::atmos_phy_bl_driver_tracer_setup(), mod_atmos_phy_ch_driver::atmos_phy_ch_driver_tracer_setup(), mod_atmos_phy_mp_driver::atmos_phy_mp_driver_tracer_setup(), scale_atmos_phy_tb_d1980::atmos_phy_tb_d1980_config(), and mod_user::user_tracer_setup().

     use scale_prc, only: &
       prc_abort
     implicit none
 
     integer,          intent(out)          :: qs
     integer,          intent(in)           :: nq
     character(len=*), intent(in)           :: name(nq)
     character(len=*), intent(in)           :: desc(nq)
     character(len=*), intent(in)           :: unit(nq)
     real(RP),         intent(in), optional :: cv  (nq)
     real(RP),         intent(in), optional :: cp  (nq)
     real(RP),         intent(in), optional :: r   (nq)
     logical,          intent(in), optional :: advc(nq)
     logical,          intent(in), optional :: mass(nq)
 
     real(RP) :: cv_  (nq)
     real(RP) :: cp_  (nq)
     real(RP) :: r_   (nq)
     logical  :: advc_(nq)
     logical  :: mass_(nq)
 
     character(len=24) :: name_trim
 
     integer  :: n
     !---------------------------------------------------------------------------
 
     if ( qa + nq > qa_max ) then
        log_error("TRACER_regist",*) 'total number of tracer must be less or equal to ', qa_max
        call prc_abort
     end if
 
     if ( present(cv) ) then
        cv_(:) = cv(:)
     else
        cv_(:) = 0.0_rp
     end if
 
     if ( present(cp) ) then
        cp_(:) = cp(:)
     else
        cp_(:) = 0.0_rp
     end if
 
     if ( present(r) ) then
        r_(:) = r(:)
     else
        r_(:) = 0.0_rp
     end if
 
     if ( present(advc) ) then
        advc_(:) = advc(:)
     else
        advc_(:) = .true.
     end if
 
     if ( present(mass) ) then
        mass_(:) = mass(:)
     else
        mass_(:) = .false.
     end if
 
     log_newline
     do n = 1, nq
 
        name_trim = trim(name(n))
 
        log_info("TRACER_regist",'(1x,A,I3,A,A,A,F6.1,A,F6.1,A,L1,A,L1)') &
                                       '] Register tracer : No.', qa+n,      &
                                                     ', NAME = ', name_trim, &
                                                       ', CV = ', cv_(n),  &
                                                       ', CP = ', cp_(n),  &
                                                     ', ADVC = ', advc_(n),  &
                                                     ', MASS = ', mass_(n)
 
        tracer_name(qa+n) = name(n)
        tracer_desc(qa+n) = desc(n)
        tracer_unit(qa+n) = unit(n)
        tracer_cv(qa+n) = cv_(n)
        tracer_cp(qa+n) = cp_(n)
        tracer_r(qa+n) = r_(n)
        tracer_advc(qa+n) = advc_(n)
 
        if ( mass_(n) ) then
           tracer_mass(qa+n) = 1.0_rp
        else
           tracer_mass(qa+n) = 0.0_rp
        end if
     end do
 
     qs = qa + 1
     qa = qa + nq
 
     return