source: LMDZ5/trunk/libf/phymar/rrtm_cmbgb3.F90 @ 5304

Last change on this file since 5304 was 2089, checked in by Laurent Fairhead, 10 years ago

Inclusion de la physique de MAR

Integration of MAR physics
File size: 4.4 KB
RevLine 
[2089]1!***************************************************************************
2SUBROUTINE RRTM_CMBGB3
3!***************************************************************************
4
5!     BAND 3:  500-630 cm-1 (low - H2O,CO2; high - H2O,CO2)
6!***************************************************************************
7
8! Parameters
9#include "tsmbkind.h"
10
11USE PARRRTM  , ONLY : JPBAND   ,JPG      ,JPXSEC   ,JPGPT
12
13USE YOERRTO3 , ONLY : KAO     ,KBO     ,SELFREFO   ,FRACREFAO  ,&
14           &FRACREFBO  ,FORREFO    ,ABSN2OAO   ,ABSN2OBO
15USE YOERRTA3 , ONLY : KA      ,KB      ,SELFREF    ,FRACREFA   ,&
16             &        FRACREFB,FORREF  ,ABSN2OA    ,ABSN2OB    ,&
17             &        ABSA    ,ABSB    ,NG3
18USE YOERRTRWT, ONLY : FREFA    ,FREFB    ,FREFADF  ,FREFBDF   ,RWGT
19USE YOERRTFTR, ONLY : NGC      ,NGS      ,NGN      ,NGB       ,NGM     , WT
20
21IMPLICIT NONE
22
23!     LOCAL INTEGER SCALARS
24INTEGER_M :: IGC, IPR, IPRSM, JN, JP, JT
25INTEGER_M :: MEQ, NEQ                    ! To force equivalence, HG, 13-DEC-2003
26
27!     LOCAL REAL SCALARS
28REAL_B :: SUMF, SUMK, SUMK1, SUMK2, SUMK3
29
30
31DO JN = 1,10
32  DO JT = 1,5
33    DO JP = 1,13
34      IPRSM = 0
35      DO IGC = 1,NGC(3)
36        SUMK = _ZERO_
37        DO IPR = 1, NGN(NGS(2)+IGC)
38          IPRSM = IPRSM + 1
39
40          SUMK = SUMK + KAO(JN,JT,JP,IPRSM)*RWGT(IPRSM+32)
41        ENDDO
42
43        KA(JN,JT,JP,IGC) = SUMK
44      ENDDO
45    ENDDO
46  ENDDO
47ENDDO
48DO JN = 1,5
49  DO JT = 1,5
50    DO JP = 13,59
51      IPRSM = 0
52      DO IGC = 1,NGC(3)
53        SUMK = _ZERO_
54        DO IPR = 1, NGN(NGS(2)+IGC)
55          IPRSM = IPRSM + 1
56
57          SUMK = SUMK + KBO(JN,JT,JP,IPRSM)*RWGT(IPRSM+32)
58        ENDDO
59
60        KB(JN,JT,JP,IGC) = SUMK
61      ENDDO
62    ENDDO
63  ENDDO
64ENDDO
65
66DO JT = 1,10
67  IPRSM = 0
68  DO IGC = 1,NGC(3)
69    SUMK = _ZERO_
70    SUMF = _ZERO_
71    DO IPR = 1, NGN(NGS(2)+IGC)
72      IPRSM = IPRSM + 1
73
74
75      SUMK = SUMK + SELFREFO(JT,IPRSM)*RWGT(IPRSM+32)
76      SUMF = SUMF + FRACREFAO(IPRSM,JT)
77    ENDDO
78
79
80    SELFREF(JT,IGC) = SUMK
81    FRACREFA(IGC,JT) = SUMF
82  ENDDO
83ENDDO
84
85DO JP = 1,5
86  IPRSM = 0
87  DO IGC = 1,NGC(3)
88    SUMF = _ZERO_
89    DO IPR = 1, NGN(NGS(2)+IGC)
90      IPRSM = IPRSM + 1
91
92      SUMF = SUMF + FRACREFBO(IPRSM,JP)
93    ENDDO
94
95    FRACREFB(IGC,JP) = SUMF
96  ENDDO
97ENDDO
98
99IPRSM = 0
100DO IGC = 1,NGC(3)
101  SUMK1= _ZERO_
102  SUMK2= _ZERO_
103  SUMK3= _ZERO_
104  DO IPR = 1, NGN(NGS(2)+IGC)
105    IPRSM = IPRSM + 1
106
107
108
109    SUMK1= SUMK1+ FORREFO(IPRSM)*RWGT(IPRSM+32)
110    SUMK2= SUMK2+ ABSN2OAO(IPRSM)*RWGT(IPRSM+32)
111    SUMK3= SUMK3+ ABSN2OBO(IPRSM)*RWGT(IPRSM+32)
112  ENDDO
113
114
115
116  FORREF(IGC) = SUMK1
117  ABSN2OA(IGC) = SUMK2
118  ABSN2OB(IGC) = SUMK3
119ENDDO
120
121DO JP = 1,10
122  DO IGC = 1,NGC(3)
123
124    FREFA(NGS(2)+IGC,JP) = FRACREFA(IGC,JP)
125  ENDDO
126ENDDO
127DO JP = 1,9
128  DO IGC = 1,NGC(3)
129
130
131    FREFADF(NGS(2)+IGC,JP) = FRACREFA(IGC,JP+1) -FRACREFA(IGC,JP)
132  ENDDO
133ENDDO
134DO JP = 1,5
135  DO IGC = 1,NGC(3)
136
137    FREFB(NGS(2)+IGC,JP) = FRACREFB(IGC,JP)
138  ENDDO
139ENDDO
140DO JP = 1,4
141  DO IGC = 1,NGC(3)
142
143
144    FREFBDF(NGS(2)+IGC,JP) = FRACREFB(IGC,JP+1) -FRACREFB(IGC,JP)
145  ENDDO
146ENDDO
147
148
149! +--Force the equivalence: BEGIN (HG, 13-DEC-2003)
150! +  ============================
151
152! +--ABSA
153! +  ^^^^
154         JN  = 0
155         JT  = 1
156         JP  = 1
157         IGC = 1
158      DO NEQ=1,NG3
159      DO MEQ=1,650
160             JN =  JN  + 1
161      IF   ( JN == 10  + 1)                                         THEN
162             JN =  1
163             JT =  JT  + 1
164       IF  ( JT == 5   + 1 )                                        THEN
165             JT =  1
166             JP =  JP  + 1
167        IF ( JP == 13  + 1 )                                        THEN
168             JP =  1
169             IGC=  IGC + 1
170        END IF
171       END IF
172      END IF
173             ABSA(MEQ,NEQ) = KA(JN,JT,JP,IGC)
174      ENDDO
175      ENDDO
176
177! +--ABSB
178! +  ^^^^
179         JN  = 0
180         JT  = 1
181         JP  = 13
182         IGC = 1
183      DO NEQ=1,NG3
184      DO MEQ=1,1175
185             JN =  JN  + 1
186      IF   ( JN == 5   + 1)                                         THEN
187             JN =  1
188             JT =  JT  + 1
189       IF  ( JT == 5   + 1 )                                        THEN
190             JT =  1
191             JP =  JP  + 1
192        IF ( JP == 59  + 1 )                                        THEN
193             JP =  13
194             IGC=  IGC + 1
195        END IF
196       END IF
197      END IF
198             ABSB(MEQ,NEQ) = KB(JN,JT,JP,IGC)
199      ENDDO
200      ENDDO
201
202! +--Force the equivalence: END   (HG, 13-DEC-2003)
203! +  ==========================
204
205
206RETURN
207END SUBROUTINE RRTM_CMBGB3
Note: See TracBrowser for help on using the repository browser.