Add T, U, V tendency due to dynamics diagnostics

[jimmielin]

DTCORE, DQCORE, UTEND_CORE, VTEND_CORE diagnostics should be implemented after the energy fixer (check_energy_fix_run)

Defined half in check_energy.F90

    call addfld('DTCORE', (/ 'lev' /), 'A', 'K/s' , 'T tendency due to dynamical core')
    call addfld('DQCORE', (/ 'lev' /), 'A', 'kg/kg/s' , 'Water vapor tendency due to dynamical core')

... the other half in physpkg.F90

    call addfld('UTEND_CORE', (/ 'lev' /), 'A', 'm/s2' , 'Zonal wind tendency due to dynamical core')
    call addfld('VTEND_CORE', (/ 'lev' /), 'A', 'm/s2' , 'Meridional wind tendency due to dynamical core')
    call register_vector_field('UTEND_CORE','VTEND_CORE')

It is computed in physpkg.F90

    ! T, U, V tendency due to dynamics
    if( nstep > dyn_time_lvls-1 ) then
       dtcore(:ncol,:pver) = (state%t(:ncol,:pver) - dtcore(:ncol,:pver))/ztodt
       dqcore(:ncol,:pver) = (state%q(:ncol,:pver,ixq) - dqcore(:ncol,:pver))/ztodt
       ducore(:ncol,:pver) = (state%u(:ncol,:pver) - ducore(:ncol,:pver))/ztodt
       dvcore(:ncol,:pver) = (state%v(:ncol,:pver) - dvcore(:ncol,:pver))/ztodt
       call outfld( 'DTCORE', dtcore, pcols, lchnk )
       call outfld( 'DQCORE', dqcore, pcols, lchnk )
       call outfld( 'UTEND_CORE', ducore, pcols, lchnk )
       call outfld( 'VTEND_CORE', dvcore, pcols, lchnk )
    end if

"adjusted" (after dycore_energy_consistency_adjust_run) state%T is stored in the pbuf for use in the next step:

    ! store T, U, and V in buffer for use in computing dynamics T-tendency in next timestep
    do k = 1,pver
       dtcore(:ncol,k) = state%t(:ncol,k)
       dqcore(:ncol,k) = state%q(:ncol,k,ixq)
       ducore(:ncol,k) = state%u(:ncol,k)
       dvcore(:ncol,k) = state%v(:ncol,k)
    end do