scale_tracer_suzuki10.F90

Go to the documentation of this file.

!-------------------------------------------------------------------------------
!-------------------------------------------------------------------------------
module scale_tracer_suzuki10
  !-----------------------------------------------------------------------------
  !
  !++ used modules
  !
  use scale_precision
  use scale_stdio
  !-----------------------------------------------------------------------------
  implicit none
  private
  !-----------------------------------------------------------------------------
  !
  !++ Public procedure
  !
  public :: tracer_suzuki10_setup

  include "inc_tracer_suzuki10.h"

  !-----------------------------------------------------------------------------
contains
  !-----------------------------------------------------------------------------
  subroutine tracer_suzuki10_setup
    use scale_process, only: &
      prc_mpistop
    implicit none

    namelist / param_bin / &
       nbin, &
       nccn, &
       iceflg, &
       kphase

    if( io_l ) write(io_fid_log,*)
    if( io_l ) write(io_fid_log,*) '+++ READ BIN NUMBER'

    rewind(io_fid_conf)
    read(io_fid_conf,nml=param_bin,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_BIN, Check!'
     call prc_mpistop
    end if

    if( io_l ) write(io_fid_log,nml=param_bin)

    if( iceflg == 0 ) then
     nspc = 1
    elseif( iceflg == 1 ) then
     nspc = 7
    else
     write(*,*) "ICEFLG should be 0(warm rain) or 1(mixed rain) check!!"
     call prc_mpistop
    endif

    nccn1 = max(nccn,1)

    !-- setup QA_MP ...
    qa_mp   = i_qv+nbin*nspc+nccn
    qqa     = i_qv+nbin*nspc
    qqs     = i_qv
    qqe     = i_qv+nbin*nspc
    qws     =  i_qv+1
    qwe     =  i_qv+nbin
    if( iceflg == 0 ) then
      qis =  0
      qie = -1
    elseif( iceflg == 1 ) then
      qis =  (i_qv+nbin+1)*iceflg
      qie =  (i_qv+nbin*nspc)*iceflg
    endif

    allocate( aq_mp_name(qa_mp) )
    allocate( aq_mp_desc(qa_mp) )
    allocate( aq_mp_unit(qa_mp) )


    !-----------------------------------------------------------------------------
    !
    !++ calculate each category and aerosol
    !
    !-----------------------------------------------------------------------------
    do n = 1, qa_mp
     write(aq_mp_unit(n),'(a)')  'kg/kg'
    enddo

    write(aq_mp_name(i_qv),'(a)') 'QV'
    do m = 1, nspc
    do n = 1, nbin
     write(aq_mp_name(i_qv+nbin*(m-1)+n),'(a,i0)') trim(namspc(m)), n
    enddo
    enddo

    do n = 1, nccn
     write(aq_mp_name(i_qv+nbin*nspc+n),'(a,i0)') trim(namspc(8)), n
    enddo

    write(aq_mp_desc(i_qv),'(a)')  'Water Vapor mixing ratio'
    do m = 1, nspc
    do n = 1, nbin
     write(aq_mp_desc(i_qv+nbin*(m-1)+n),'(a,i0)') trim(lnamspc(m)), n
    enddo
    enddo

    do n = 1, nccn
     write(aq_mp_desc(i_qv+nbin*nspc+n),'(a,i0)') trim(lnamspc(8)), n
    enddo

    return
  end subroutine tracer_suzuki10_setup

end module scale_tracer_suzuki10