Convert between buoyancy fluxes and others

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