source: LMDZ6/trunk/libf/phylmd/rrtm/srtm_cmbgb28.F90 @ 5450

Last change on this file since 5450 was 1990, checked in by Laurent Fairhead, 11 years ago

Corrections à la version r1989 pour permettre la compilation avec RRTM
Inclusion de la licence CeCILL_V2 pour RRTM


Changes to revision r1989 to enable RRTM code compilation
RRTM part put under CeCILL_V2 licence

  • Property copyright set to
    Name of program: LMDZ
    Creation date: 1984
    Version: LMDZ5
    License: CeCILL version 2
    Holder: Laboratoire de m\'et\'eorologie dynamique, CNRS, UMR 8539
    See the license file in the root directory
File size: 1.6 KB
Line 
1SUBROUTINE SRTM_CMBGB28
2
3!     BAND 28:  38000-50000 cm-1 (low - O3,O2; high - O3,O2)
4!-----------------------------------------------------------------------
5
6USE PARKIND1  ,ONLY : JPIM , JPRB
7USE YOMHOOK   ,ONLY : LHOOK, DR_HOOK
8
9USE YOESRTWN , ONLY : NGC, NGS, NGN, RWGT
10USE YOESRTA28, ONLY : KA, KB, SFLUXREF, &
11                    & KAC, KBC, SFLUXREFC
12
13IMPLICIT NONE
14
15! Local variables
16INTEGER(KIND=JPIM) :: JN, JT, JP, IGC, IPR, IPRSM
17REAL(KIND=JPRB)    :: ZSUMK, ZSUMF
18
19REAL(KIND=JPRB) :: ZHOOK_HANDLE
20!     ------------------------------------------------------------------
21IF (LHOOK) CALL DR_HOOK('SRTM_CMBGB28',0,ZHOOK_HANDLE)
22
23DO JN = 1,9
24  DO JT = 1,5
25    DO JP = 1,13
26      IPRSM = 0
27      DO IGC = 1,NGC(13)
28        ZSUMK = 0.
29        DO IPR = 1, NGN(NGS(12)+IGC)
30          IPRSM = IPRSM + 1
31          ZSUMK = ZSUMK + KA(JN,JT,JP,IPRSM)*RWGT(IPRSM+192)
32        ENDDO
33        KAC(JN,JT,JP,IGC) = ZSUMK
34      ENDDO
35    ENDDO
36  ENDDO
37ENDDO
38
39DO JN = 1,5
40  DO JT = 1,5
41    DO JP = 13,59
42      IPRSM = 0
43      DO IGC = 1,NGC(13)
44        ZSUMK = 0.
45        DO IPR = 1, NGN(NGS(12)+IGC)
46          IPRSM = IPRSM + 1
47          ZSUMK = ZSUMK + KB(JN,JT,JP,IPRSM)*RWGT(IPRSM+192)
48        ENDDO
49        KBC(JN,JT,JP,IGC) = ZSUMK
50      ENDDO
51    ENDDO
52  ENDDO
53ENDDO
54
55DO JP = 1,5
56  IPRSM = 0
57  DO IGC = 1,NGC(13)
58    ZSUMF = 0.
59    DO IPR = 1, NGN(NGS(12)+IGC)
60      IPRSM = IPRSM + 1
61      ZSUMF = ZSUMF + SFLUXREF(IPRSM,JP)
62    ENDDO
63    SFLUXREFC(IGC,JP) = ZSUMF
64  ENDDO
65ENDDO
66
67!     -----------------------------------------------------------------
68IF (LHOOK) CALL DR_HOOK('SRTM_CMBGB28',1,ZHOOK_HANDLE)
69END SUBROUTINE SRTM_CMBGB28
70
Note: See TracBrowser for help on using the repository browser.