olw.F90

Go to the documentation of this file.

SUBROUTINE olw &
     & ( kidia, kfdia , klon  , klev &
     & , pcco2, pcldld, pcldlu &
     & , pdp  , pdt0  , pemis  &
     & , paph , pqof  , pth &
     & , paer , pt    , pview , pwv &
     & , pcolr, pcolc , pflux, pfluc &
     & )
!
!**** *LW*   - ORGANIZES THE LONGWAVE CALCULATIONS
!
!     PURPOSE.
!     --------
!           DEPENDING ON KMODE, COMPUTES LONGWAVE FLUXES AND/OR
!           RADIANCES
!
!**   INTERFACE.
!     ----------
!
!        *LW* IS CALLED FROM *RADLSW*
!
!        EXPLICIT ARGUMENTS :
!        --------------------
! PAER   : (KLON,KLEV,6)     ; OPTICAL THICKNESS OF THE AEROSOLS
! PCCO2  :                   ; CONCENTRATION IN CO2 (PA/PA)
! PCLDLD : (KLON,KLEV)       ; DOWNWARD EFFECTIVE FRACTIONAL COVER
! PCLDLU : (KLON,KLEV)       ; UPWARD EFFECTIVE FRACTIONAL COVER
! PDP    : (KLON,KLEV)       ; LAYER PRESSURE THICKNESS
! PDT0   : (KLON)            ; SURFACE TEMPERATURE DISCONTINUITY  
! PEMIS  : (KLON)            ; SURFACE EMISSIVITY
! PAPH   : (KLON,KLEV+1)     ; HALF LEVEL PRESSURE
! PQOF   : (KLON,KLEV)       ; CONCENTRATION IN OZONE (PA/PA)
! PT     : (KLON,KLEV)       ; TEMPERATURE
! PTH    : (KLON,KLEV+1)     ; HALF LEVEL TEMPERATURE
! PVIEW  : (KLON)            ; COSECANT OF VIEWING ANGLE
! PWV    : (KLON,KLEV)       ; SPECIFIC HUMIDITY  (PA/PA)
!     ==== OUTPUTS ===
!  IF KMODE = 0, 1, 2
! PFLUX(KLON,2,KLEV)         ; RADIATIVE FLUXES :
!                     1  ==>  UPWARD   FLUX TOTAL
!                     2  ==>  DOWNWARD FLUX TOTAL
! PFLUC(KLON,2,KLEV)         ; RADIATIVE FLUXES CLEAR SKY:
!                     1  ==>  UPWARD   FLUX TOTAL
!                     2  ==>  DOWNWARD FLUX TOTAL
! PCOLR(KLON,KLEV)           ; LONG-WAVE TENDENCY
! PCOLC(KLON,KLEV)           ; LONG-WAVE TENDENCY CLEAR SKY
!
!        IMPLICIT ARGUMENTS :   NONE
!        --------------------
!
!     METHOD.
!     -------
!
!          1. COMPUTES THE PRESSURE AND TEMPERATURE WEIGHTED AMOUNTS OF
!     ABSORBERS.
!          2. COMPUTES THE PLANCK FUNCTIONS ON THE INTERFACES AND THE
!     GRADIENT OF PLANCK FUNCTIONS IN THE LAYERS.
!          3. PERFORMS THE VERTICAL INTEGRATION DISTINGUISHING THE CON-
!     TRIBUTIONS OF THE ADJACENT AND DISTANT LAYERS AND THOSE FROM THE
!     BOUNDARIES.
!          4. COMPUTES THE CLEAR-SKY DOWNWARD AND UPWARD EMISSIVITIES.
!          5. INTRODUCES THE EFFECTS OF THE CLOUDS ON THE FLUXES.
!
!     EXTERNALS.
!     ----------
!
!          *LWU*, *LWBV*, *LWC*
!
!     REFERENCE.
!     ----------
!
!        SEE RADIATION'S PART OF THE MODEL'S DOCUMENTATION AND
!        ECMWF RESEARCH DEPARTMENT DOCUMENTATION OF THE IFS
!
!     AUTHOR.
!     -------
!        JEAN-JACQUES MORCRETTE  *ECMWF*
!
!     MODIFICATIONS.
!     --------------
!        ORIGINAL : 89-07-14
!-----------------------------------------------------------------------

#include "tsmbkind.h"

USE yoeolw   , ONLY : nua
USE yoerdu   , ONLY : rcday
USE yoedbug  , ONLY : ldebug


IMPLICIT NONE

!     DUMMY INTEGER SCALARS
integer_m :: kfdia
integer_m :: kidia
integer_m :: klev
integer_m :: klon

!     DUMMY REAL SCALARS
real_b :: pcco2, zdfnet

!-----------------------------------------------------------------------
!
!*       0.1   ARGUMENTS
!              ---------
!
real_b :: pcldld(klon,klev)      , pcldlu(klon,klev) &
     &  ,  pdp(klon,klev)   , pdt0(klon) &
     &  ,  pemis(klon)      , paph(klon,klev+1) &
     &  ,  pqof(klon,klev)  , pth(klon,klev+1) &
     &  ,  paer(klon,klev,6), pt(klon,klev) &
     &  ,  pview(klon)      , pwv(klon,klev)
!
real_b :: pcolr(klon,klev)       , pcolc(klon,klev) &
     &  ,  pflux(klon,2,klev+1), pfluc(klon,2,klev+1)
!
!-------------------------------------------------------------------------
!
!*       0.2   LOCAL ARRAYS
!              ------------
real_b :: zabcu(klon,nua,3*klev+1) &
     &  ,  zbint(klon,klev+1) &
     &  ,  zbsui(klon) &
     &  ,  zcntrb(klon,klev+1,klev+1) &
     &  ,  zfdn(klon,klev+1) &
     &  ,  zfup(klon,klev+1)
     
!     LOCAL INTEGER SCALARS
integer_m :: jk, jkl, jl, ilim, iua, klevt, jk1, jk2

     
!
!     ------------------------------------------------------------------
!
!*         1.    INITIALIZATION
!                --------------
!
!     ------------------------------------------------------------------
!
!*         1.1   COMPUTES ABSORBER AMOUNTS
!                -------------------------
!

CALL olwu ( kidia, kfdia, klon, klev &
     &  ,  paer, pcco2, pdp, paph, pqof, pt, pview, pwv &
     &  ,  zabcu   )

!
!     ------------------------------------------------------------------
!
!*         2.    COMPUTES PLANCK FUNCTIONS
!                -------------------------
!                PERFORMS THE VERTICAL INTEGRATION
!                ---------------------------------

DO jk1=1,klev+1
  DO jk2=1,klev+1
    DO jl=kidia,kfdia
      zcntrb(jl,jk1,jk2)=0.
    END DO
  END DO
END DO      
!
CALL olwbv( kidia,kfdia,klon,klev &
     &          , pdp,pdt0,pemis,pth &
     &          , pt &
     &          , pcolc,pfluc &
     &          , zabcu,zbint,zbsui,zcntrb,zfdn,zfup)
     
!     ------------------------------------------------------------------
!
!*         4.    INTRODUCES THE EFFECTS OF CLOUDS
!                --------------------------------
!
CALL olwc ( kidia,kfdia,klon,klev &
     &  , zbint,zbsui,pcldld,pcldlu,zcntrb,pemis,zfdn,zfup &
     &  , pflux                                                )
!
DO jkl = 1 , klev
  jk = klev+1 - jkl

  DO jl = kidia,kfdia
    zdfnet = pflux(jl,1,jk+1) + pflux(jl,2,jk+1) &
     &        -pflux(jl,1,jk  ) - pflux(jl,2,jk  )
    pcolr(jl,jk) = rcday * zdfnet / pdp(jl,jkl)
    
    zdfnet = pfluc(jl,1,jk+1) + pfluc(jl,2,jk+1) &
    &              -pfluc(jl,1,jk  ) - pfluc(jl,2,jk  )
    pcolc(jl,jk) = rcday * zdfnet / pdp(jl,jkl)
  END DO
  
END DO
!
!     ------------------------------------------------------------------
!
RETURN
END SUBROUTINE olw