1 |
|
|
SUBROUTINE 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 |
|
|
|
46 |
|
|
USE PARKIND1 ,ONLY : JPIM ,JPRB |
47 |
|
|
USE YOMHOOK ,ONLY : LHOOK, DR_HOOK |
48 |
|
|
|
49 |
|
|
IMPLICIT NONE |
50 |
|
|
|
51 |
|
|
INTEGER(KIND=JPIM),INTENT(IN) :: KLEV |
52 |
|
|
REAL(KIND=JPRB) ,INTENT(IN) :: PARA |
53 |
|
|
INTEGER(KIND=JPIM),INTENT(OUT) :: KPARA |
54 |
|
|
REAL(KIND=JPRB) ,INTENT(IN) :: PSTPRE(KLEV) |
55 |
|
|
INTEGER(KIND=JPIM) :: JLEV |
56 |
|
|
REAL(KIND=JPRB) :: ZHOOK_HANDLE |
57 |
|
|
|
58 |
|
|
!* 1. SEARCH OF CORRESPONDING LEVEL |
59 |
|
|
! ----------------------------- |
60 |
|
|
IF (LHOOK) CALL DR_HOOK('SEAPRE',0,ZHOOK_HANDLE) |
61 |
|
|
KPARA=KLEV |
62 |
|
|
DO JLEV=KLEV,1,-1 |
63 |
|
|
IF(PARA <= PSTPRE(JLEV)) KPARA=JLEV |
64 |
|
|
ENDDO |
65 |
|
|
|
66 |
|
|
IF (LHOOK) CALL DR_HOOK('SEAPRE',1,ZHOOK_HANDLE) |
67 |
|
|
END SUBROUTINE SEAPRE |
68 |
|
|
|