zm_convr_tendency_diagnostics.F90

module zm_convr_tendency_diagnostics
   use ccpp_kinds, only:  kind_phys

   implicit none
   private
   save

   public :: zm_convr_tendency_diagnostics_init ! init routine
   public :: zm_convr_tendency_diagnostics_run  ! main routine

CONTAINS

 !> \section arg_table_zm_convr_tendency_diagnostics_init  Argument Table
 !! \htmlinclude zm_convr_tendency_diagnostics_init.html
 subroutine zm_convr_tendency_diagnostics_init(errmsg, errflg)
    use cam_history,         only: history_add_field
    use cam_history_support, only: horiz_only

    character(len=512), intent(out) :: errmsg
    integer,            intent(out) :: errflg

    ! Local variables:

    errmsg = ''
    errflg = 0

    call history_add_field ('ZMDT', 'T tendency - Zhang-McFarlane moist convection', 'lev', 'avg', 'K s-1')
    call history_add_field ('ZMDQ', 'Q tendency - Zhang-McFarlane moist convection', 'lev', 'avg', 'kg kg-1 s-1')

   end subroutine zm_convr_tendency_diagnostics_init

   !> \section arg_table_zm_convr_tendency_diagnostics_run  Argument Table
   !! \htmlinclude zm_convr_tendency_diagnostics_run.html
   subroutine zm_convr_tendency_diagnostics_run(ncol, pver, cpair, heat, qtnd, errmsg, errflg)

      use cam_history, only: history_out_field

      !------------------------------------------------
      !   Input / output parameters
      !------------------------------------------------
      integer, intent(in) :: ncol
      integer, intent(in) :: pver

      real(kind_phys), intent(in) :: cpair
      real(kind_phys), intent(in) :: heat(:,:)
      real(kind_phys), intent(in) :: qtnd(:,:)

      ! CCPP error handling variables
      character(len=512), intent(out) :: errmsg
      integer,            intent(out) :: errflg


      real(kind_phys) :: ftem(ncol,pver)

      errmsg = ''
      errflg = 0

      ftem(:,:) = 0._kind_phys

      ftem(:ncol,:pver) = heat(:ncol,:pver)/cpair
      call history_out_field('ZMDT    ', ftem)
      call history_out_field('ZMDQ    ', qtnd)

   end subroutine zm_convr_tendency_diagnostics_run

   !=======================================================================

end module zm_convr_tendency_diagnostics