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 64 of file scale_tracer.F90.

References scale_stdio::io_fid_log, scale_stdio::io_l, scale_process::prc_mpistop(), 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(), scale_atmos_phy_ae_kajino13::atmos_phy_ae_kajino13_config(), scale_atmos_phy_mp_sn14::atmos_phy_mp_sn14_config(), scale_atmos_phy_mp_suzuki10::atmos_phy_mp_suzuki10_config(), scale_atmos_phy_tb_d1980::atmos_phy_tb_d1980_config(), scale_atmos_phy_tb_mynn::atmos_phy_tb_mynn_config(), scale_atmos_phy_tb_smg::atmos_phy_tb_smg_config(), and mod_user::user_config().

     use scale_process, only: &
       prc_mpistop
     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
        write(*,*) 'xxx total number of tracer must be less or equal to ', qa_max
        call prc_mpistop
     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
 
     if( io_l ) write(io_fid_log,*)
     do n = 1, nq
 
        name_trim = trim(name(n))
 
        if( io_l ) write(io_fid_log,'(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