1 | include(dom.inc)
2 |
3 | SUBROUTINE FINDI(LICO,LICO2,LIH2O,ICO,ICO2,IH2O,WVNB,DWVNB)
4 |
5 | !*****************************************************************
6 | ! *
7 | ! THIS SUBROUTINE SEARCHS THE PARAMETER INDEXES *
8 | ! CORRESPONDING THE WAVE NUMBER 'WVNB' *
9 | ! *
10 | !*****************************************************************
11 |
12 | IMPLICIT NONE
13 |
14 | include 'dom_constants.h'
15 |
16 | DOM_INT :: ICO, ICO2, IH2O
17 | DOM_REAL :: WVNB, DWVNB
18 |
19 | LOGICAL LICO,LICO2,LIH2O
20 |
21 | !-----------------------------!
22 | ! LOOKING FOR THE CO2 INDEXES !
23 | !-----------------------------!
24 |
25 | ICO2=-9
26 |
27 | IF (WVNB.ge.450.0.and.WVNB.le.1200) THEN
28 | ICO2=INT((WVNB-450.)/25.)+1
29 | ENDIF
30 |
31 | IF (WVNB.ge.1950.0.and.WVNB.le.2450.) THEN
32 | ICO2=INT((WVNB-1950.)/25.)+32
33 | ENDIF
34 |
35 | IF(WVNB.ge.3300.0.and.WVNB.le.3800.) THEN
36 | ICO2=INT((WVNB-3300.)/25.)+53
37 | ENDIF
38 |
39 | IF(WVNB.ge.4700.0.and.WVNB.le.5250.) THEN
40 | ICO2=INT((WVNB-4700.)/25.)+74
41 | ENDIF
42 |
43 | !-----------------------------!
44 | ! LOOKING FOR THE H2O INDEXES !
45 | !-----------------------------!
46 |
47 | IH2O=-9
48 | IF(WVNB.ge.150.and.WVNB.le.9300.) THEN
49 | IH2O=INT((WVNB-150.)/25.)+1
50 | ENDIF
51 |
52 | !----------------------------!
53 | ! LOOKING FOR THE CO INDEXES !
54 | !----------------------------!
55 |
56 | ICO=-9
57 |
58 | IF (WVNB.ge.1750.0.and.WVNB.le.2335.) THEN
59 | ICO=INT((WVNB-1750.)/25.)+1
60 | ENDIF
61 |
62 | IF (WVNB.ge.3775.0.and.WVNB.le.4350.) THEN
63 | ICO=INT((WVNB-3775.)/25.)+25
64 | ENDIF
65 |
66 | LICO = ICO.GT.0
67 | LICO2 = ICO2.GT.0
68 | LIH2O = IH2O.GT.0
69 |
70 | END SUBROUTINE FINDI
findi.F could be called by: