gppre.F90

Go to the documentation of this file.

SUBROUTINE gppre(KPROMA,KSTART,KPROF,KFLEV,PVAH,PVBH,PRESH,PRESF)

!**** *GPPRE* - Computes half and full level pressure

!     Purpose.
!     --------
!           Computes pressures at half and full model levels.

!**   Interface.
!     ----------
!        *CALL* *GPPRE(KPROMA,KSTART,KPROF,KFLEV,PVAH,PVBH,PRESH,PRESF)

!        Explicit arguments :
!        --------------------
!                              KPROMA :  dimensioning
!                              KSTART :  start of work
!                              KPROF  :  depth of work
!                              KFLEV     : vert. dimensioning
!                              PRESH(KPROMA,0:KFLEV) - HALF LEVEL PRESSURE
!                              PRESF(KPROMA,KFLEV)   - FULL LEVEL PRESSURE
!                              PVAH(KFLEV),PVBH(KFLEV)- vertical coordinate
!        Implicit arguments :  NONE.
!        --------------------

!     Method.
!     -------
!        See documentation

!     Externals.  None.
!     ----------

!     Reference.
!     ----------
!        ECMWF Research Department documentation of the IFS
!     Half level P:     PHk = Ak + Bk * Ps

!                                PHk*ln(PHk) - PHk-1*ln(PHk-1)
!     Full level P: ln(PFk) = [ ------------------------------- - 1. ]
!                                        PHk - PHk-1

!     Author.
!     -------
!        Mats Hamrud and Philippe Courtier  *ECMWF*

!     Modifications.
!     --------------
!        Original : 88-02-04
!     Erik Andersson  920326: Altered computation of full level pressure
!     Erik Andersson  930225: Use GPPREH/F.
!        M.Hamrud      01-Oct-2003 CY28 Cleaning
!     ------------------------------------------------------------------

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

USE yomgem   , ONLY : vc       ,vdelb

IMPLICIT NONE

INTEGER(KIND=JPIM),INTENT(IN)    :: KPROMA 
INTEGER(KIND=JPIM),INTENT(IN)    :: KFLEV 
INTEGER(KIND=JPIM),INTENT(IN)    :: KSTART 
INTEGER(KIND=JPIM),INTENT(IN)    :: KPROF 
REAL(KIND=JPRB)   ,INTENT(IN)    :: PVAH(0:kflev) 
REAL(KIND=JPRB)   ,INTENT(IN)    :: PVBH(0:kflev) 
REAL(KIND=JPRB)   ,INTENT(INOUT) :: PRESH(kproma,0:kflev) 
REAL(KIND=JPRB)   ,INTENT(OUT)   :: PRESF(kproma,kflev) 
REAL(KIND=JPRB) ::    ZLNPR(kproma,kflev) ,  ZALPH (kproma,kflev)
REAL(KIND=JPRB) ::    ZDELP(kproma,kflev) ,  ZRDELP(kproma,kflev)
REAL(KIND=JPRB) ::    ZRTGR(kproma,kflev) ,  ZRPRES(kproma,kflev)
REAL(KIND=JPRB) ::    ZRPP (kproma,kflev)
REAL(KIND=JPRB) :: ZHOOK_HANDLE

#include "gppref.intfb.h"
#include "gppreh.intfb.h"
#include "gpxyb.intfb.h"

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

!*       1.    COMPUTES HALF AND FULL LEVEL PRESSURES
!              --------------------------------------

IF (lhook) CALL dr_hook('GPPRE',0,zhook_handle)
CALL gppreh(kproma,kstart,kprof,kflev,pvah,pvbh,presh)
CALL gpxyb(kproma,kstart,kprof,kflev,vdelb,vc,presh,zdelp,&
 & zrdelp,zlnpr,zalph,zrtgr,zrpres,zrpp)  
CALL gppref(kproma,kstart,kprof,kflev,pvah,pvbh,zalph,presh,presf)

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

IF (lhook) CALL dr_hook('GPPRE',1,zhook_handle)
END SUBROUTINE gppre