lw.F90

Go to the documentation of this file.

SUBROUTINE lw &
 & ( kidia, kfdia , klon  , klev  , kmode,&
 & pcco2, pcldld, pcldlu,&
 & pdp  , pdt0  , pemis , pemiw,&
 & ppmb , pqof  , ptl,&
 & paer , ptave , pview , pwv,&
 & pemit, pflux , pfluc &
 & )  

!**** *LW*   - ORGANIZES THE LONGWAVE CALCULATIONS

!     PURPOSE.
!     --------
!           COMPUTES LONGWAVE FLUXES 

!**   INTERFACE.
!     ----------

!        *LW* IS CALLED FROM *RADLSW*

!        EXPLICIT ARGUMENTS :
!        --------------------
! PAER   : (KLON,6,KLEV)     ; OPTICAL THICKNESS OF THE AEROSOLS
! PCCO2  :                   ; CONCENTRATION IN CO2 (PA/PA)
! PCLDLD : (KLON,KLEV)       ; DOWNWARD EFFECTIVE CLOUD FRACTION
! PCLDLU : (KLON,KLEV)       ; UPWARD EFFECTIVE CLOUD FRACTION
! PDP    : (KLON,KLEV)       ; LAYER PRESSURE THICKNESS
! PDT0   : (KLON)            ; SURFACE TEMPERATURE DISCONTINUITY  
! PEMIS  : (KLON)            ; SURFACE LW EMISSIVITY
! PEMIW  : (KLON)            ; SURFACE LW WINDOW EMISSIVITY
! PPMB   : (KLON,KLEV+1)     ; HALF LEVEL PRESSURE
! PQOF   : (KLON,KLEV)       ; CONCENTRATION IN OZONE (PA/PA)
! PTAVE  : (KLON,KLEV)       ; TEMPERATURE
! PTL    : (KLON,KLEV+1)     ; HALF LEVEL TEMPERATURE
! PVIEW  : (KLON)            ; COSECANT OF VIEWING ANGLE
! PWV    : (KLON,KLEV)       ; SPECIFIC HUMIDITY  (PA/PA)
!     ==== OUTPUTS ===
! PEMIT(KLON)                ; SURFACE TOTAL LW EMISSIVITY
! PFLUX(KLON,2,KLEV+1)       ; RADIATIVE FLUXES :
!                     1  ==>  UPWARD   FLUX TOTAL
!                     2  ==>  DOWNWARD FLUX TOTAL
! PFLUC(KLON,2,KLEV+1)       ; RADIATIVE FLUXES CLEAR SKY:
!                     1  ==>  UPWARD   FLUX TOTAL
!                     2  ==>  DOWNWARD FLUX TOTAL

!        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
!        99-05-25   JJMorcrette    Revised aerosols
!        M.Hamrud      01-Oct-2003 CY28 Cleaning

!-----------------------------------------------------------------------

USE parkind1  ,ONLY : jpim     ,jprb
USE yomhook   ,ONLY : lhook,   dr_hook

USE yoelw    , ONLY : nua
IMPLICIT NONE

INTEGER(KIND=JPIM),INTENT(IN)    :: KLON 
INTEGER(KIND=JPIM),INTENT(IN)    :: KLEV 
INTEGER(KIND=JPIM),INTENT(IN)    :: KIDIA 
INTEGER(KIND=JPIM),INTENT(IN)    :: KFDIA 
INTEGER(KIND=JPIM),INTENT(IN)    :: KMODE 
REAL(KIND=JPRB)   ,INTENT(IN)    :: PCCO2 
REAL(KIND=JPRB)   ,INTENT(IN)    :: PCLDLD(klon,klev) 
REAL(KIND=JPRB)   ,INTENT(IN)    :: PCLDLU(klon,klev) 
REAL(KIND=JPRB)   ,INTENT(IN)    :: PDP(klon,klev) 
REAL(KIND=JPRB)   ,INTENT(IN)    :: PDT0(klon) 
REAL(KIND=JPRB)   ,INTENT(IN)    :: PEMIS(klon) 
REAL(KIND=JPRB)   ,INTENT(IN)    :: PEMIW(klon) 
REAL(KIND=JPRB)   ,INTENT(IN)    :: PPMB(klon,klev+1) 
REAL(KIND=JPRB)   ,INTENT(IN)    :: PQOF(klon,klev) 
REAL(KIND=JPRB)   ,INTENT(IN)    :: PTL(klon,klev+1) 
REAL(KIND=JPRB)   ,INTENT(IN)    :: PAER(klon,6,klev) 
REAL(KIND=JPRB)   ,INTENT(IN)    :: PTAVE(klon,klev) 
REAL(KIND=JPRB)   ,INTENT(IN)    :: PVIEW(klon) 
REAL(KIND=JPRB)   ,INTENT(IN)    :: PWV(klon,klev) 
REAL(KIND=JPRB)   ,INTENT(INOUT) :: PEMIT(klon) 
REAL(KIND=JPRB)   ,INTENT(OUT)   :: PFLUX(klon,2,klev+1) 
REAL(KIND=JPRB)   ,INTENT(INOUT) :: PFLUC(klon,2,klev+1) 
!-----------------------------------------------------------------------

!*       0.1   ARGUMENTS
!              ---------

!-------------------------------------------------------------------------

!              ------------
REAL(KIND=JPRB) :: ZABCU(klon,nua,3*klev+1)&
 & ,  ZBINT(KLON,KLEV+1)        , ZBSUI(KLON)&
 & ,  ZCNTRB(KLON,KLEV+1,KLEV+1)  

REAL(KIND=JPRB) :: ZHOOK_HANDLE

#include "lwbv.intfb.h"
#include "lwc.intfb.h"
#include "lwu.intfb.h"

!     ------------------------------------------------------------------

!*         1.    INITIALIZATION
!                --------------

!100  CONTINUE

!     ------------------------------------------------------------------

!*         1.1   COMPUTES ABSORBER AMOUNTS
!                -------------------------

IF (lhook) CALL dr_hook('LW',0,zhook_handle)
print *,'       LW: Avant LWU'
CALL lwu &
 & (  kidia, kfdia, klon, klev,&
 & paer , pcco2, pdp , ppmb, pqof , ptave, pview, pwv,&
 & zabcu &
 & )  

!     ------------------------------------------------------------------

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

print *,'       LW: Avant LWBV'
CALL lwbv &
 & ( kidia, kfdia, klon , klev  , kmode,&
 & pdt0 , pemis, pemiw, ptl   , ptave,&
 & pemit, pfluc,&
 & zabcu, zbint, zbsui, zcntrb &
 & )  

!     ------------------------------------------------------------------

!*         4.    INTRODUCES THE EFFECTS OF CLOUDS
!                --------------------------------

print *,'       LW: Avant LWC'
CALL lwc &
 & ( kidia , kfdia, klon  , klev,&
 & zbint , zbsui, pcldld, pcldlu,&
 & zcntrb, pemit, pfluc,&
 & pflux    &
 & )  

!     ------------------------------------------------------------------

IF (lhook) CALL dr_hook('LW',1,zhook_handle)
END SUBROUTINE lw