presf.f90

      !+ computes fulllevel pressures.
      !+ $Id: presf.f90,v 1.5 1999/01/26 09:31:06 m214003 Exp $
      
      SUBROUTINE presf(pf,kdimp,ph,klen)
   5: 
        ! Description:
        !
        ! Compute full-level pressures from half-level values.
        !
  10:   ! Method:
        !
        ! Full-level pressures are defined as the arithmetic
        ! average of the two adjoining half-level pressures.
        !
  15:   ! *presf* is called from *physc*. 
        ! Parameters are:
        !    *pf*        *computed full-level pressures.
        !    *ph*        *half-level pressures.
        !    *kdimp*     *first dimension of 2-d arrays *pf* and *ph*
  20:   !    *klen*      *number of points for which computation is
        !                 performed.
        !
        ! Authors:
        !
  25:   ! M. Jarraud, ECMWF, May 1982, original source
        ! L. Kornblueh, MPI, May 1998, f90 rewrite
        ! U. Schulzweida, MPI, May 1998, f90 rewrite
        ! 
        ! for more details see file AUTHORS
  30:   !
      
        USE mo_control, ONLY: nlev
      
        IMPLICIT NONE
  35: 
        !  Scalar arguments 
        INTEGER ,INTENT(in) :: kdimp, klen
      
        !  Array arguments 
  40:   REAL    ,INTENT(in) :: ph(kdimp,*)
        REAL    ,INTENT(out):: pf(kdimp,*)
      
        !  Local scalars: 
        INTEGER :: jlev, jlon
  45: 
      
        !  Executable statements 
      
      !-- 1. Compute full-level pressure values
  50: 
        DO jlev = 1, nlev
          DO jlon = 1, klen
            pf(jlon,jlev) = (ph(jlon,jlev)+ph(jlon,jlev+1))*.5
          END DO
  55:   END DO
      
        RETURN
      END SUBROUTINE presf


Info Section
uses: mo_control
back to top
ECHAM 4 vf90 (C) 1998 Max-Planck-Institut für Meteorologie, Hamburg
Wed Nov 24 01:25:21 CST 1999

HTML derived from FORTRAN source by f2html.pl v0.3 (C) 1997,98 Beroud Jean-Marc.