source: LMDZ5/branches/IPSLCM5A2.1_ISO/libf/phyiso/rrtm/seapre.F90 @ 5308

Last change on this file since 5308 was 3331, checked in by acozic, 6 years ago

Add modification for isotopes

  • Property svn:executable set to *
File size: 1.7 KB
Line 
1SUBROUTINE SEAPRE(PARA,KPARA,PSTPRE,KLEV)
2
3!**** *SEAPRE*   - Search corresponding level to input pressure
4
5!     Purpose.
6!     --------
7!           Serach corresponding level of the model to an input pressure
8!           through standard atmosphere.
9
10!**   Interface.
11!     ----------
12!        *CALL* *SEAPRE(PARA,KPARA,PSTPRE,KLEV)
13
14!        Explicit arguments :
15!        --------------------
16!        PARA  : Pressure                           (input)
17!        KPARA  : Level                              (output)
18!        PSTPRE : Standard atmosphere                (input)
19!        KLEV   : Number of level of the model       (input)
20
21!        Implicit arguments :
22!        --------------------
23!        none
24
25!     Method.
26!     -------
27!        See documentation
28
29!     Externals.
30!     ----------
31
32!     Reference.
33!     ----------
34!        Documentation ARPEGE
35
36!     Author.
37!     -------
38!        A. Lasserre-Bigorry
39
40!     Modifications.
41!     --------------
42!        Original : 91-06-10
43!        M.Hamrud      01-Oct-2003 CY28 Cleaning
44!     ------------------------------------------------------------------
45
46USE PARKIND1  ,ONLY : JPIM     ,JPRB
47USE YOMHOOK   ,ONLY : LHOOK,   DR_HOOK
48
49IMPLICIT NONE
50
51INTEGER(KIND=JPIM),INTENT(IN)    :: KLEV
52REAL(KIND=JPRB)   ,INTENT(IN)    :: PARA
53INTEGER(KIND=JPIM),INTENT(OUT)   :: KPARA
54REAL(KIND=JPRB)   ,INTENT(IN)    :: PSTPRE(KLEV)
55INTEGER(KIND=JPIM) :: JLEV
56REAL(KIND=JPRB) :: ZHOOK_HANDLE
57
58!*      1.  SEARCH OF CORRESPONDING LEVEL
59!           -----------------------------
60IF (LHOOK) CALL DR_HOOK('SEAPRE',0,ZHOOK_HANDLE)
61KPARA=KLEV
62DO JLEV=KLEV,1,-1
63  IF(PARA <= PSTPRE(JLEV)) KPARA=JLEV
64ENDDO
65
66IF (LHOOK) CALL DR_HOOK('SEAPRE',1,ZHOOK_HANDLE)
67END SUBROUTINE SEAPRE
68
Note: See TracBrowser for help on using the repository browser.