Convert between buoyancy fluxes and others (Source File: convert_fluxes.F90)


INTERFACE:

   subroutine  convert_fluxes(nlev,g,cp,rho_0,heat,p_e,rad,T,S,           &
                             tFlux,sFlux,btFlux,bsFlux,tRad,bRad)
DESCRIPTION:

This subroutine computes the buoyancy fluxes that are due to

  1. the surface heat flux,
  2. the surface salinity flux caused by the value of P-E (precipitation-evaporation),
  3. and the short wave radiative flux.
Additionally, it outputs the temperature flux (tFlux) corresponding to the surface heat flux, the salinity flux (sFlux) corresponding to the value P-E, and the profile of the temperature flux (tRad) corresponding to the profile of the radiative heat flux.

This function is only called when the KPP turbulence model is used. When you call the KPP routines from another model outside GOTM, you are on your own in computing the fluxes required by the KPP model, because they have to be consistent with the equation of state used in your model.


USES:

   use eqstate
   IMPLICIT NONE
INPUT PARAMETERS:
   integer,  intent(in)                :: nlev
   REALTYPE, intent(in)                :: g,cp,rho_0
   REALTYPE, intent(in)                :: heat,p_e
   REALTYPE, intent(in)                :: rad(0:nlev)
   REALTYPE, intent(in)                :: T(0:nlev)
   REALTYPE, intent(in)                :: S(0:nlev)
OUTPUT PARAMETERS:
   REALTYPE, intent(out)              ::  tFlux,sFlux
   REALTYPE, intent(out)               :: btFlux,bsFlux
   REALTYPE, intent(out)               :: tRad(0:nlev)
   REALTYPE, intent(out)               :: bRad(0:nlev)
REVISION HISTORY:
    Original author(s): Lars Umlauf
    $Log: convert_fluxes.F90,v $
    Revision 1.2  2005-08-11 12:34:32  lars
    corrected indentation for Protex
    Revision 1.1  2005/06/27 10:54:33  kbk
    new files needed

Karsten Bolding 2012-01-24