source: LMDZ5/branches/IPSLCM5A2.1_ISO/libf/phyiso/rrtm/rrtm_cmbgb10.F90 @ 3331

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

Add modification for isotopes

  • Property svn:executable set to *
File size: 1.8 KB
Line 
1!***************************************************************************
2SUBROUTINE RRTM_CMBGB10
3!***************************************************************************
4
5!     BAND 10:  1390-1480 cm-1 (low - H2O; high - H2O)
6!***************************************************************************
7
8! Parameters
9USE PARKIND1  ,ONLY : JPIM     ,JPRB
10USE YOMHOOK   ,ONLY : LHOOK,   DR_HOOK
11
12USE YOERRTO10, ONLY : KAO     ,KBO      ,FRACREFAO   ,FRACREFBO
13USE YOERRTA10, ONLY : KA      ,KB       ,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_CMBGB10',0,ZHOOK_HANDLE)
25DO JT = 1,5
26  DO JP = 1,13
27    IPRSM = 0
28    DO IGC = 1,NGC(10)
29      Z_SUMK = 0.0_JPRB
30      DO IPR = 1, NGN(NGS(9)+IGC)
31        IPRSM = IPRSM + 1
32
33        Z_SUMK = Z_SUMK + KAO(JT,JP,IPRSM)*RWGT(IPRSM+144)
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(10)
44      Z_SUMK = 0.0_JPRB
45      DO IPR = 1, NGN(NGS(9)+IGC)
46        IPRSM = IPRSM + 1
47
48        Z_SUMK = Z_SUMK + KBO(JT,JP,IPRSM)*RWGT(IPRSM+144)
49      ENDDO
50
51      KB(JT,JP,IGC) = Z_SUMK
52    ENDDO
53  ENDDO
54ENDDO
55
56IPRSM = 0
57DO IGC = 1,NGC(10)
58  Z_SUMF1= 0.0_JPRB
59  Z_SUMF2= 0.0_JPRB
60  DO IPR = 1, NGN(NGS(9)+IGC)
61    IPRSM = IPRSM + 1
62
63    Z_SUMF1= Z_SUMF1+ FRACREFAO(IPRSM)
64    Z_SUMF2= Z_SUMF2+ FRACREFBO(IPRSM)
65  ENDDO
66
67  FRACREFA(IGC) = Z_SUMF1
68  FRACREFB(IGC) = Z_SUMF2
69ENDDO
70
71DO IGC = 1,NGC(10)
72
73  FREFA(NGS(9)+IGC,1) = FRACREFA(IGC)
74  FREFB(NGS(9)+IGC,1) = FRACREFB(IGC)
75ENDDO
76
77IF (LHOOK) CALL DR_HOOK('RRTM_CMBGB10',1,ZHOOK_HANDLE)
78END SUBROUTINE RRTM_CMBGB10
Note: See TracBrowser for help on using the repository browser.