source: LMDZ5/branches/IPSLCM5A2.1_ISO/libf/phyiso/rrtm/rrtm_cmbgb11.F90

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

Add modification for isotopes

  • Property svn:executable set to *
File size: 2.0 KB
Line 
1!***************************************************************************
2SUBROUTINE RRTM_CMBGB11
3!***************************************************************************
4
5!     BAND 11:  1480-1800 cm-1 (low - H2O; high - H2O)
6!***************************************************************************
7
8! Parameters
9USE PARKIND1  ,ONLY : JPIM     ,JPRB
10USE YOMHOOK   ,ONLY : LHOOK,   DR_HOOK
11
12USE YOERRTO11, ONLY : KAO     ,KBO     ,SELFREFO    ,FRACREFAO ,FRACREFBO
13USE YOERRTA11, ONLY : KA      ,KB      ,SELFREF     ,FRACREFA  ,FRACREFB
14USE YOERRTRWT, ONLY : FREFA    ,FREFB    ,RWGT
15USE YOERRTFTR, ONLY : NGC      ,NGS      ,NGN     
16
17IMPLICIT NONE
18
19INTEGER(KIND=JPIM) :: IGC, IPR, IPRSM, JP, JT
20
21REAL(KIND=JPRB) :: Z_SUMF1, Z_SUMF2, Z_SUMK
22REAL(KIND=JPRB) :: ZHOOK_HANDLE
23
24IF (LHOOK) CALL DR_HOOK('RRTM_CMBGB11',0,ZHOOK_HANDLE)
25DO JT = 1,5
26  DO JP = 1,13
27    IPRSM = 0
28    DO IGC = 1,NGC(11)
29      Z_SUMK = 0.0_JPRB
30      DO IPR = 1, NGN(NGS(10)+IGC)
31        IPRSM = IPRSM + 1
32
33        Z_SUMK = Z_SUMK + KAO(JT,JP,IPRSM)*RWGT(IPRSM+160)
34      ENDDO
35
36      KA(JT,JP,IGC) = Z_SUMK
37    ENDDO
38  ENDDO
39ENDDO
40DO JT = 1,5
41  DO JP = 13,59
42    IPRSM = 0
43    DO IGC = 1,NGC(11)
44      Z_SUMK = 0.0_JPRB
45      DO IPR = 1, NGN(NGS(10)+IGC)
46        IPRSM = IPRSM + 1
47
48        Z_SUMK = Z_SUMK + KBO(JT,JP,IPRSM)*RWGT(IPRSM+160)
49      ENDDO
50
51      KB(JT,JP,IGC) = Z_SUMK
52    ENDDO
53  ENDDO
54ENDDO
55
56DO JT = 1,10
57  IPRSM = 0
58  DO IGC = 1,NGC(11)
59    Z_SUMK = 0.0_JPRB
60    DO IPR = 1, NGN(NGS(10)+IGC)
61      IPRSM = IPRSM + 1
62
63      Z_SUMK = Z_SUMK + SELFREFO(JT,IPRSM)*RWGT(IPRSM+160)
64    ENDDO
65
66    SELFREF(JT,IGC) = Z_SUMK
67  ENDDO
68ENDDO
69
70IPRSM = 0
71DO IGC = 1,NGC(11)
72  Z_SUMF1= 0.0_JPRB
73  Z_SUMF2= 0.0_JPRB
74  DO IPR = 1, NGN(NGS(10)+IGC)
75    IPRSM = IPRSM + 1
76
77    Z_SUMF1= Z_SUMF1+ FRACREFAO(IPRSM)
78    Z_SUMF2= Z_SUMF2+ FRACREFBO(IPRSM)
79  ENDDO
80
81  FRACREFA(IGC) = Z_SUMF1
82  FRACREFB(IGC) = Z_SUMF2
83ENDDO
84
85DO IGC = 1,NGC(11)
86
87  FREFA(NGS(10)+IGC,1) = FRACREFA(IGC)
88  FREFB(NGS(10)+IGC,1) = FRACREFB(IGC)
89ENDDO
90
91IF (LHOOK) CALL DR_HOOK('RRTM_CMBGB11',1,ZHOOK_HANDLE)
92END SUBROUTINE RRTM_CMBGB11
Note: See TracBrowser for help on using the repository browser.