Compute various diagnostic/integrated variables

INTERFACE:

    subroutine do_diagnostics(n,nlev,BuoyMeth,dt,u_taus,u_taub,I_0,heat)

DESCRIPTION:

This subroutine calculates the following diagnostic/integrated variables.

USES:

    use airsea,       only: sst
    use meanflow,     only: gravity,rho_0,cp
    use meanflow,     only: h,u,v,s,t,NN,SS,buoy,rad
    use turbulence,   only: kappa
    use turbulence,   only: tke
    use observations, only: tprof,b_obs_sbf
    use eqstate,      only: eqstate1
    IMPLICIT NONE

INPUT PARAMETERS:

    integer, intent(in)                 :: n,nlev,BuoyMeth
    REALTYPE, intent(in)                :: dt
    REALTYPE, intent(in)                :: u_taus,u_taub
    REALTYPE, intent(in)                :: I_0,heat

REVISION HISTORY:

    Original author(s): Karsten Bolding & Hans Burchard
    See output module