source: LMDZ5/branches/IPSLCM5A2.1_ISO/libf/phyiso/rrtm/srtm_cmbgb28.F90 @ 4249

Last change on this file since 4249 was 3331, checked in by acozic, 7 years ago

Add modification for isotopes

  • Property svn:executable set to *
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.