scale_tracer.F90

Go to the documentation of this file.

!-------------------------------------------------------------------------------
!-------------------------------------------------------------------------------
! Warning: This file was generated from common/scale_tracer.F90.erb.
!          Do not edit this file.
!-------------------------------------------------------------------------------
module scale_tracer
  !-----------------------------------------------------------------------------
  !
  !++ used modules
  !
  use scale_precision
  use scale_stdio
  !-----------------------------------------------------------------------------
  implicit none
  private
  !-----------------------------------------------------------------------------
  !
  !++ Public procedure
  !
  !-----------------------------------------------------------------------------
  !
  !++ Public parameters & variables
  !

#ifdef TRACER
#define EXTM(name) inc_tracer_ ## name ## .h
#define NAME(name) EXTM(name)
#define XSTR(s) # s
#define STR(s) XSTR(s)
  include str(name(tracer))
#else
  integer, public :: qa

  integer, public :: qa_mp
  integer, public :: i_qv
  integer, public :: i_qc
  integer, public :: i_qr
  integer, public :: i_qi
  integer, public :: i_qs
  integer, public :: i_qg
  integer, public :: i_nc
  integer, public :: i_nr
  integer, public :: i_ni
  integer, public :: i_ns
  integer, public :: i_ng

  integer, public :: qqa ! mass tracer (water)
  integer, public :: qqs ! start index for mass tracer
  integer, public :: qqe ! end   index for mass tracer

  integer, public :: qws ! start index for water tracer
  integer, public :: qwe ! end   index for water tracer
  integer, public :: qis ! start index for ice tracer
  integer, public :: qie ! end   index for ice tracer

  integer, public :: qa_ae   ! number of tracer for aerosol
  integer, public :: ae_ctg  ! category number of aerosol
  integer, public :: gas_ctg ! category number of gas
  integer, public :: n_atr   ! number of attribute
  integer, public, allocatable :: nkap(:)
  integer, public, allocatable :: nsiz(:)
  integer, public :: ic_mix
  integer, public :: ic_sea
  integer, public :: ic_dus
  integer, public :: ig_h2so4
  integer, public :: ig_cgas
  integer, public :: qaes ! start index for aerosol tracer
  integer, public :: qaee ! end   index for aerosol tracer

  character(len=H_SHORT), public, allocatable :: aq_name(:)
  character(len=H_SHORT), public, allocatable :: aq_mp_name(:)
  character(len=H_SHORT), public, allocatable :: aq_ae_name(:)
  character(len=H_MID)  , public, allocatable :: aq_desc(:)
  character(len=H_MID)  , public, allocatable :: aq_mp_desc(:)
  character(len=H_MID)  , public, allocatable :: aq_ae_desc(:)
  character(len=H_SHORT), public, allocatable :: aq_unit(:)
  character(len=H_SHORT), public, allocatable :: aq_mp_unit(:)
  character(len=H_SHORT), public, allocatable :: aq_ae_unit(:)

  !-----------------------------------------------------------------------------
  !
  !++ tracer index & relationship (MP_dry+AE_dummy+RD_mstrnx)
  !
  !-----------------------------------------------------------------------------

  integer, public :: mp_qa ! number of hydrometeor tracer
  integer, public, allocatable :: i_mp2all(:)

  integer, public, allocatable :: i_mp2rd(:)

  integer, public :: ae_qa ! number of aerosol tracer
  integer, public :: i_ae_dummy
  integer, public, allocatable :: i_ae2all(:)

  integer, public, allocatable :: i_ae2rd(:)

  character(len=H_SHORT), public :: tracer_type = 'DRY'
  character(len=H_SHORT), public :: aetracer_type = 'NONE'
#endif
  !-----------------------------------------------------------------------------
  !
  !++ Private procedure
  !
  public :: tracer_setup

  !-----------------------------------------------------------------------------
  !
  !++ Private parameters & variables
  !
  !-----------------------------------------------------------------------------
contains
  !-----------------------------------------------------------------------------
  subroutine tracer_setup
#ifdef TRACER
#define EXTM2(name) scale_tracer_ ## name
#define _MODNAME(name) EXTM2(name)
#define EXTM3(name) tracer_ ## name ## _setup
#define _SETUP(name) EXTM3(name)
  use _modname(tracer), only: &
    _setup(tracer)
#else
  use scale_process, only: &
    prc_mpistop
  use scale_tracer_dry, only: &
    qa_mp_dry => qa_mp, &
    i_qv_dry => i_qv, &
    i_qc_dry => i_qc, &
    i_qr_dry => i_qr, &
    i_qi_dry => i_qi, &
    i_qs_dry => i_qs, &
    i_qg_dry => i_qg, &
    i_nc_dry => i_nc, &
    i_nr_dry => i_nr, &
    i_ni_dry => i_ni, &
    i_ns_dry => i_ns, &
    i_ng_dry => i_ng, &
    qqa_dry => qqa, &
    qqs_dry => qqs, &
    qqe_dry => qqe, &
    qws_dry => qws, &
    qwe_dry => qwe, &
    qis_dry => qis, &
    qie_dry => qie, &
    mp_qa_dry => mp_qa, &
    aq_mp_name_dry => aq_mp_name, &
    aq_mp_desc_dry => aq_mp_desc, &
    aq_mp_unit_dry => aq_mp_unit, &
    i_mp2all_dry => i_mp2all, &
    i_mp2rd_dry => i_mp2rd, &
    tracer_dry_setup
  use scale_tracer_kessler, only: &
    qa_mp_kessler => qa_mp, &
    i_qv_kessler => i_qv, &
    i_qc_kessler => i_qc, &
    i_qr_kessler => i_qr, &
    i_qi_kessler => i_qi, &
    i_qs_kessler => i_qs, &
    i_qg_kessler => i_qg, &
    i_nc_kessler => i_nc, &
    i_nr_kessler => i_nr, &
    i_ni_kessler => i_ni, &
    i_ns_kessler => i_ns, &
    i_ng_kessler => i_ng, &
    qqa_kessler => qqa, &
    qqs_kessler => qqs, &
    qqe_kessler => qqe, &
    qws_kessler => qws, &
    qwe_kessler => qwe, &
    qis_kessler => qis, &
    qie_kessler => qie, &
    mp_qa_kessler => mp_qa, &
    aq_mp_name_kessler => aq_mp_name, &
    aq_mp_desc_kessler => aq_mp_desc, &
    aq_mp_unit_kessler => aq_mp_unit, &
    i_mp2all_kessler => i_mp2all, &
    i_mp2rd_kessler => i_mp2rd, &
    tracer_kessler_setup
  use scale_tracer_tomita08, only: &
    qa_mp_tomita08 => qa_mp, &
    i_qv_tomita08 => i_qv, &
    i_qc_tomita08 => i_qc, &
    i_qr_tomita08 => i_qr, &
    i_qi_tomita08 => i_qi, &
    i_qs_tomita08 => i_qs, &
    i_qg_tomita08 => i_qg, &
    i_nc_tomita08 => i_nc, &
    i_nr_tomita08 => i_nr, &
    i_ni_tomita08 => i_ni, &
    i_ns_tomita08 => i_ns, &
    i_ng_tomita08 => i_ng, &
    qqa_tomita08 => qqa, &
    qqs_tomita08 => qqs, &
    qqe_tomita08 => qqe, &
    qws_tomita08 => qws, &
    qwe_tomita08 => qwe, &
    qis_tomita08 => qis, &
    qie_tomita08 => qie, &
    mp_qa_tomita08 => mp_qa, &
    aq_mp_name_tomita08 => aq_mp_name, &
    aq_mp_desc_tomita08 => aq_mp_desc, &
    aq_mp_unit_tomita08 => aq_mp_unit, &
    i_mp2all_tomita08 => i_mp2all, &
    i_mp2rd_tomita08 => i_mp2rd, &
    tracer_tomita08_setup
  use scale_tracer_sn14, only: &
    qa_mp_sn14 => qa_mp, &
    i_qv_sn14 => i_qv, &
    i_qc_sn14 => i_qc, &
    i_qr_sn14 => i_qr, &
    i_qi_sn14 => i_qi, &
    i_qs_sn14 => i_qs, &
    i_qg_sn14 => i_qg, &
    i_nc_sn14 => i_nc, &
    i_nr_sn14 => i_nr, &
    i_ni_sn14 => i_ni, &
    i_ns_sn14 => i_ns, &
    i_ng_sn14 => i_ng, &
    qqa_sn14 => qqa, &
    qqs_sn14 => qqs, &
    qqe_sn14 => qqe, &
    qws_sn14 => qws, &
    qwe_sn14 => qwe, &
    qis_sn14 => qis, &
    qie_sn14 => qie, &
    mp_qa_sn14 => mp_qa, &
    aq_mp_name_sn14 => aq_mp_name, &
    aq_mp_desc_sn14 => aq_mp_desc, &
    aq_mp_unit_sn14 => aq_mp_unit, &
    i_mp2all_sn14 => i_mp2all, &
    i_mp2rd_sn14 => i_mp2rd, &
    tracer_sn14_setup
  use scale_tracer_suzuki10, only: &
    qa_mp_suzuki10 => qa_mp, &
    i_qv_suzuki10 => i_qv, &
    i_qc_suzuki10 => i_qc, &
    i_qr_suzuki10 => i_qr, &
    i_qi_suzuki10 => i_qi, &
    i_qs_suzuki10 => i_qs, &
    i_qg_suzuki10 => i_qg, &
    i_nc_suzuki10 => i_nc, &
    i_nr_suzuki10 => i_nr, &
    i_ni_suzuki10 => i_ni, &
    i_ns_suzuki10 => i_ns, &
    i_ng_suzuki10 => i_ng, &
    qqa_suzuki10 => qqa, &
    qqs_suzuki10 => qqs, &
    qqe_suzuki10 => qqe, &
    qws_suzuki10 => qws, &
    qwe_suzuki10 => qwe, &
    qis_suzuki10 => qis, &
    qie_suzuki10 => qie, &
    mp_qa_suzuki10 => mp_qa, &
    aq_mp_name_suzuki10 => aq_mp_name, &
    aq_mp_desc_suzuki10 => aq_mp_desc, &
    aq_mp_unit_suzuki10 => aq_mp_unit, &
    i_mp2all_suzuki10 => i_mp2all, &
    i_mp2rd_suzuki10 => i_mp2rd, &
    tracer_suzuki10_setup
  use scale_tracer_sdm, only: &
    qa_mp_sdm => qa_mp, &
    i_qv_sdm => i_qv, &
    i_qc_sdm => i_qc, &
    i_qr_sdm => i_qr, &
    i_qi_sdm => i_qi, &
    i_qs_sdm => i_qs, &
    i_qg_sdm => i_qg, &
    i_nc_sdm => i_nc, &
    i_nr_sdm => i_nr, &
    i_ni_sdm => i_ni, &
    i_ns_sdm => i_ns, &
    i_ng_sdm => i_ng, &
    qqa_sdm => qqa, &
    qqs_sdm => qqs, &
    qqe_sdm => qqe, &
    qws_sdm => qws, &
    qwe_sdm => qwe, &
    qis_sdm => qis, &
    qie_sdm => qie, &
    mp_qa_sdm => mp_qa, &
    aq_mp_name_sdm => aq_mp_name, &
    aq_mp_desc_sdm => aq_mp_desc, &
    aq_mp_unit_sdm => aq_mp_unit, &
    i_mp2all_sdm => i_mp2all, &
    i_mp2rd_sdm => i_mp2rd, &
    tracer_sdm_setup
  use scale_aetracer_none, only: &
    qa_ae_none => qa_ae, &
    ae_ctg_none => ae_ctg, &
    n_atr_none => n_atr, &
    gas_ctg_none => gas_ctg, &
    nsiz_none => nsiz, &
    nkap_none => nkap, &
    ic_mix_none => ic_mix, &
    ic_sea_none => ic_sea, &
    ic_dus_none => ic_dus, &
    ig_h2so4_none => ig_h2so4, &
    ig_cgas_none => ig_cgas, &
    ae_qa_none => ae_qa, &
    i_ae_dummy_none => i_ae_dummy, &
    aq_ae_name_none => aq_ae_name, &
    aq_ae_desc_none => aq_ae_desc, &
    aq_ae_unit_none => aq_ae_unit, &
    i_ae2all_none => i_ae2all, &
    i_ae2rd_none => i_ae2rd, &
    aetracer_none_setup
  use scale_aetracer_kajino13, only: &
    qa_ae_kajino13 => qa_ae, &
    ae_ctg_kajino13 => ae_ctg, &
    n_atr_kajino13 => n_atr, &
    gas_ctg_kajino13 => gas_ctg, &
    nsiz_kajino13 => nsiz, &
    nkap_kajino13 => nkap, &
    ic_mix_kajino13 => ic_mix, &
    ic_sea_kajino13 => ic_sea, &
    ic_dus_kajino13 => ic_dus, &
    ig_h2so4_kajino13 => ig_h2so4, &
    ig_cgas_kajino13 => ig_cgas, &
    ae_qa_kajino13 => ae_qa, &
    i_ae_dummy_kajino13 => i_ae_dummy, &
    aq_ae_name_kajino13 => aq_ae_name, &
    aq_ae_desc_kajino13 => aq_ae_desc, &
    aq_ae_unit_kajino13 => aq_ae_unit, &
    i_ae2all_kajino13 => i_ae2all, &
    i_ae2rd_kajino13 => i_ae2rd, &
    aetracer_kajino13_setup
#endif
    implicit none
#ifndef TRACER
    namelist / param_tracer / &
         tracer_type, &
         aetracer_type
#endif

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

    if( io_l ) write(io_fid_log,*)
    if( io_l ) write(io_fid_log,*) '+++ Module[TRACER]/Categ[COMMON]'

#ifdef TRACER
    call _setup(tracer)
#else
    !--- read namelist
    rewind(io_fid_conf)
    read(io_fid_conf,nml=param_tracer,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_TRACER. Check!'
       call prc_mpistop
    endif
    if( io_l ) write(io_fid_log,nml=param_tracer)

    select case (tracer_type)
    case ("DRY")
       call tracer_dry_setup
       qa_mp = qa_mp_dry
       i_qv = i_qv_dry
       i_qc = i_qc_dry
       i_qr = i_qr_dry
       i_qi = i_qi_dry
       i_qs = i_qs_dry
       i_qg = i_qg_dry
       i_nc = i_nc_dry
       i_nr = i_nr_dry
       i_ni = i_ni_dry
       i_ns = i_ns_dry
       i_ng = i_ng_dry
       qqa = qqa_dry
       qqs = qqs_dry
       qqe = qqe_dry
       qws = qws_dry
       qwe = qwe_dry
       qis = qis_dry
       qie = qie_dry
       mp_qa = mp_qa_dry
       allocate( aq_mp_name(qa_mp) )
       allocate( aq_mp_desc(qa_mp) )
       allocate( aq_mp_unit(qa_mp) )
       allocate( i_mp2all(mp_qa) )
       allocate( i_mp2rd(mp_qa) )
       aq_mp_name = aq_mp_name_dry
       aq_mp_desc = aq_mp_desc_dry
       aq_mp_unit = aq_mp_unit_dry
       i_mp2all = i_mp2all_dry
       i_mp2rd = i_mp2rd_dry
    case ("KESSLER")
       call tracer_kessler_setup
       qa_mp = qa_mp_kessler
       i_qv = i_qv_kessler
       i_qc = i_qc_kessler
       i_qr = i_qr_kessler
       i_qi = i_qi_kessler
       i_qs = i_qs_kessler
       i_qg = i_qg_kessler
       i_nc = i_nc_kessler
       i_nr = i_nr_kessler
       i_ni = i_ni_kessler
       i_ns = i_ns_kessler
       i_ng = i_ng_kessler
       qqa = qqa_kessler
       qqs = qqs_kessler
       qqe = qqe_kessler
       qws = qws_kessler
       qwe = qwe_kessler
       qis = qis_kessler
       qie = qie_kessler
       mp_qa = mp_qa_kessler
       allocate( aq_mp_name(qa_mp) )
       allocate( aq_mp_desc(qa_mp) )
       allocate( aq_mp_unit(qa_mp) )
       allocate( i_mp2all(mp_qa) )
       allocate( i_mp2rd(mp_qa) )
       aq_mp_name = aq_mp_name_kessler
       aq_mp_desc = aq_mp_desc_kessler
       aq_mp_unit = aq_mp_unit_kessler
       i_mp2all = i_mp2all_kessler
       i_mp2rd = i_mp2rd_kessler
    case ("TOMITA08")
       call tracer_tomita08_setup
       qa_mp = qa_mp_tomita08
       i_qv = i_qv_tomita08
       i_qc = i_qc_tomita08
       i_qr = i_qr_tomita08
       i_qi = i_qi_tomita08
       i_qs = i_qs_tomita08
       i_qg = i_qg_tomita08
       i_nc = i_nc_tomita08
       i_nr = i_nr_tomita08
       i_ni = i_ni_tomita08
       i_ns = i_ns_tomita08
       i_ng = i_ng_tomita08
       qqa = qqa_tomita08
       qqs = qqs_tomita08
       qqe = qqe_tomita08
       qws = qws_tomita08
       qwe = qwe_tomita08
       qis = qis_tomita08
       qie = qie_tomita08
       mp_qa = mp_qa_tomita08
       allocate( aq_mp_name(qa_mp) )
       allocate( aq_mp_desc(qa_mp) )
       allocate( aq_mp_unit(qa_mp) )
       allocate( i_mp2all(mp_qa) )
       allocate( i_mp2rd(mp_qa) )
       aq_mp_name = aq_mp_name_tomita08
       aq_mp_desc = aq_mp_desc_tomita08
       aq_mp_unit = aq_mp_unit_tomita08
       i_mp2all = i_mp2all_tomita08
       i_mp2rd = i_mp2rd_tomita08
    case ("SN14")
       call tracer_sn14_setup
       qa_mp = qa_mp_sn14
       i_qv = i_qv_sn14
       i_qc = i_qc_sn14
       i_qr = i_qr_sn14
       i_qi = i_qi_sn14
       i_qs = i_qs_sn14
       i_qg = i_qg_sn14
       i_nc = i_nc_sn14
       i_nr = i_nr_sn14
       i_ni = i_ni_sn14
       i_ns = i_ns_sn14
       i_ng = i_ng_sn14
       qqa = qqa_sn14
       qqs = qqs_sn14
       qqe = qqe_sn14
       qws = qws_sn14
       qwe = qwe_sn14
       qis = qis_sn14
       qie = qie_sn14
       mp_qa = mp_qa_sn14
       allocate( aq_mp_name(qa_mp) )
       allocate( aq_mp_desc(qa_mp) )
       allocate( aq_mp_unit(qa_mp) )
       allocate( i_mp2all(mp_qa) )
       allocate( i_mp2rd(mp_qa) )
       aq_mp_name = aq_mp_name_sn14
       aq_mp_desc = aq_mp_desc_sn14
       aq_mp_unit = aq_mp_unit_sn14
       i_mp2all = i_mp2all_sn14
       i_mp2rd = i_mp2rd_sn14
    case ("SUZUKI10")
       call tracer_suzuki10_setup
       qa_mp = qa_mp_suzuki10
       i_qv = i_qv_suzuki10
       i_qc = i_qc_suzuki10
       i_qr = i_qr_suzuki10
       i_qi = i_qi_suzuki10
       i_qs = i_qs_suzuki10
       i_qg = i_qg_suzuki10
       i_nc = i_nc_suzuki10
       i_nr = i_nr_suzuki10
       i_ni = i_ni_suzuki10
       i_ns = i_ns_suzuki10
       i_ng = i_ng_suzuki10
       qqa = qqa_suzuki10
       qqs = qqs_suzuki10
       qqe = qqe_suzuki10
       qws = qws_suzuki10
       qwe = qwe_suzuki10
       qis = qis_suzuki10
       qie = qie_suzuki10
       mp_qa = mp_qa_suzuki10
       allocate( aq_mp_name(qa_mp) )
       allocate( aq_mp_desc(qa_mp) )
       allocate( aq_mp_unit(qa_mp) )
       allocate( i_mp2all(mp_qa) )
       allocate( i_mp2rd(mp_qa) )
       aq_mp_name = aq_mp_name_suzuki10
       aq_mp_desc = aq_mp_desc_suzuki10
       aq_mp_unit = aq_mp_unit_suzuki10
       i_mp2all = i_mp2all_suzuki10
       i_mp2rd = i_mp2rd_suzuki10
    case ("SDM")
       call tracer_sdm_setup
       qa_mp = qa_mp_sdm
       i_qv = i_qv_sdm
       i_qc = i_qc_sdm
       i_qr = i_qr_sdm
       i_qi = i_qi_sdm
       i_qs = i_qs_sdm
       i_qg = i_qg_sdm
       i_nc = i_nc_sdm
       i_nr = i_nr_sdm
       i_ni = i_ni_sdm
       i_ns = i_ns_sdm
       i_ng = i_ng_sdm
       qqa = qqa_sdm
       qqs = qqs_sdm
       qqe = qqe_sdm
       qws = qws_sdm
       qwe = qwe_sdm
       qis = qis_sdm
       qie = qie_sdm
       mp_qa = mp_qa_sdm
       allocate( aq_mp_name(qa_mp) )
       allocate( aq_mp_desc(qa_mp) )
       allocate( aq_mp_unit(qa_mp) )
       allocate( i_mp2all(mp_qa) )
       allocate( i_mp2rd(mp_qa) )
       aq_mp_name = aq_mp_name_sdm
       aq_mp_desc = aq_mp_desc_sdm
       aq_mp_unit = aq_mp_unit_sdm
       i_mp2all = i_mp2all_sdm
       i_mp2rd = i_mp2rd_sdm
     case default
        write(*,*) 'xxx Unsupported TRACER_TYPE (', trim(tracer_type), '). Check!'
        call prc_mpistop
     end select
     select case (aetracer_type)
     case ("NONE")
       call aetracer_none_setup
       qa_ae = qa_ae_none
       ae_ctg = ae_ctg_none
       n_atr = n_atr_none
       gas_ctg = gas_ctg_none
       allocate( nsiz(ae_ctg) )
       allocate( nkap(ae_ctg) )
       nsiz = nsiz_none
       nkap = nkap_none
       ic_mix = ic_mix_none
       ic_sea = ic_sea_none
       ic_dus = ic_dus_none
       ig_h2so4 = ig_h2so4_none
       ig_cgas = ig_cgas_none
       ae_qa = ae_qa_none
       i_ae_dummy = i_ae_dummy_none
       allocate( aq_ae_name(qa_ae) )
       allocate( aq_ae_desc(qa_ae) )
       allocate( aq_ae_unit(qa_ae) )
       allocate( i_ae2all(ae_qa) )
       allocate( i_ae2rd(ae_qa) )
       aq_ae_name = aq_ae_name_none
       aq_ae_desc = aq_ae_desc_none
       aq_ae_unit = aq_ae_unit_none
       i_ae2all = i_ae2all_none
       i_ae2rd = i_ae2rd_none
     case ("KAJINO13")
       call aetracer_kajino13_setup
       qa_ae = qa_ae_kajino13
       ae_ctg = ae_ctg_kajino13
       n_atr = n_atr_kajino13
       gas_ctg = gas_ctg_kajino13
       allocate( nsiz(ae_ctg) )
       allocate( nkap(ae_ctg) )
       nsiz = nsiz_kajino13
       nkap = nkap_kajino13
       ic_mix = ic_mix_kajino13
       ic_sea = ic_sea_kajino13
       ic_dus = ic_dus_kajino13
       ig_h2so4 = ig_h2so4_kajino13
       ig_cgas = ig_cgas_kajino13
       ae_qa = ae_qa_kajino13
       i_ae_dummy = i_ae_dummy_kajino13
       allocate( aq_ae_name(qa_ae) )
       allocate( aq_ae_desc(qa_ae) )
       allocate( aq_ae_unit(qa_ae) )
       allocate( i_ae2all(ae_qa) )
       allocate( i_ae2rd(ae_qa) )
       aq_ae_name = aq_ae_name_kajino13
       aq_ae_desc = aq_ae_desc_kajino13
       aq_ae_unit = aq_ae_unit_kajino13
       i_ae2all = i_ae2all_kajino13
       i_ae2rd = i_ae2rd_kajino13
     case default
        write(*,*) 'xxx Unsupported AETRACER_TYPE (', trim(aetracer_type), '). Check!'
        call prc_mpistop
     end select
     qaes = qa_mp+1
     qaee = qa_mp+qa_ae
     qa = qa_mp + qa_ae
     allocate( aq_name(qa) )
     allocate( aq_desc(qa) )
     allocate( aq_unit(qa) )
     aq_name(1:qa_mp) = aq_mp_name(1:qa_mp)
     aq_desc(1:qa_mp) = aq_mp_desc(1:qa_mp)
     aq_unit(1:qa_mp) = aq_mp_unit(1:qa_mp)
     aq_name(qa_mp+1:qa_mp+qa_ae) = aq_ae_name(1:qa_ae)
     aq_desc(qa_mp+1:qa_mp+qa_ae) = aq_ae_desc(1:qa_ae)
     aq_unit(qa_mp+1:qa_mp+qa_ae) = aq_ae_unit(1:qa_ae)
#endif

    return
  end subroutine tracer_setup

end module scale_tracer