scale_tracer_suzuki10 Module Reference

module TRACER / suzuki10 More...

Functions/Subroutines
subroutine, public	tracer_suzuki10_setup

Detailed Description

module TRACER / suzuki10

Description

Tracer suzuki10 module

Author

Team SCALE

History

• 2013-12-21 (Y.Sato) [new] imported from inc_tracer_suzuki10.f90

NAMELIST

• PARAM_BIN

  name	type	default value	comment
  NBIN	integer	33
  NCCN	integer	0
  ICEFLG	integer	1
  KPHASE	integer	0

  
  

History Output

No history output

Function/Subroutine Documentation

tracer_suzuki10_setup()

subroutine, public scale_tracer_suzuki10::tracer_suzuki10_setup

(

)

Definition at line 36 of file scale_tracer_suzuki10.F90.

References scale_stdio::io_fid_conf, scale_stdio::io_fid_log, scale_stdio::io_l, and scale_process::prc_mpistop().

Referenced by scale_tracer::tracer_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

Here is the call graph for this function:

Here is the caller graph for this function: