source: LMDZ5/branches/IPSLCM5A2.1_ISO/libf/phyiso/rrtm/yoesw.F90 @ 5308

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

Add modification for isotopes

  • Property svn:executable set to *
File size: 10.1 KB
Line 
1MODULE YOESW
2
3USE PARKIND1  ,ONLY : JPIM     ,JPRB
4
5IMPLICIT NONE
6
7SAVE
8
9!     ------------------------------------------------------------------
10!*    ** *YOESW* - COEFFICIENTS FOR SHORTWAVE RADIATION TRANSFER
11!     ------------------------------------------------------------------
12
13REAL(KIND=JPRB) :: APAD(6,3,7)
14REAL(KIND=JPRB) :: BPAD(6,3,7)
15REAL(KIND=JPRB) :: RRAY(6,6)
16REAL(KIND=JPRB), ALLOCATABLE :: RSUN(:)
17!$OMP THREADPRIVATE(RSUN)
18REAL(KIND=JPRB) :: RPDH1
19REAL(KIND=JPRB) :: RPDU1
20REAL(KIND=JPRB) :: RPNH
21REAL(KIND=JPRB) :: RPNU
22REAL(KIND=JPRB) :: RSWCE(6)
23REAL(KIND=JPRB) :: RSWCP(6)
24REAL(KIND=JPRB) :: RTDH2O
25REAL(KIND=JPRB) :: RTDUMG
26REAL(KIND=JPRB) :: RTH2O
27REAL(KIND=JPRB) :: RTUMG
28REAL(KIND=JPRB) :: D(6,3)
29REAL(KIND=JPRB) :: REXPO3(6,2,7)
30INTEGER(KIND=JPIM) :: NEXPO3(6)
31
32REAL(KIND=JPRB) :: RYFWCA(6)
33REAL(KIND=JPRB) :: RYFWCB(6)
34REAL(KIND=JPRB) :: RYFWCC(6)
35REAL(KIND=JPRB) :: RYFWCD(6)
36REAL(KIND=JPRB) :: RYFWCE(6)
37REAL(KIND=JPRB) :: RYFWCF(6)
38
39REAL(KIND=JPRB) :: REBCUA(6)
40REAL(KIND=JPRB) :: REBCUB(6)
41REAL(KIND=JPRB) :: REBCUC(6)
42REAL(KIND=JPRB) :: REBCUD(6)
43REAL(KIND=JPRB) :: REBCUE(6)
44REAL(KIND=JPRB) :: REBCUF(6)
45REAL(KIND=JPRB) :: REBCUG(16)
46REAL(KIND=JPRB) :: REBCUH(16)
47REAL(KIND=JPRB) :: REBCUI(6)
48REAL(KIND=JPRB) :: REBCUJ(6)
49
50REAL(KIND=JPRB) :: RASWCA(6)
51REAL(KIND=JPRB) :: RASWCB(6)
52REAL(KIND=JPRB) :: RASWCC(6)
53REAL(KIND=JPRB) :: RASWCD(6)
54REAL(KIND=JPRB) :: RASWCE(6)
55REAL(KIND=JPRB) :: RASWCF(6)
56
57REAL(KIND=JPRB) :: RFUETA(16,3),RFUETB(16,4), RFUETC(16,4)
58REAL(KIND=JPRB) :: RFULIO(16,3)
59REAL(KIND=JPRB) :: RHSAVI(16,3)
60REAL(KIND=JPRB) :: RLILIA(16,5),RLILIB(16,4)
61
62REAL(KIND=JPRB) :: RFLAA0(6)
63REAL(KIND=JPRB) :: RFLAA1(6)
64REAL(KIND=JPRB) :: RFLBB0(6)
65REAL(KIND=JPRB) :: RFLBB1(6)
66REAL(KIND=JPRB) :: RFLBB2(6)
67REAL(KIND=JPRB) :: RFLBB3(6)
68REAL(KIND=JPRB) :: RFLCC0(6)
69REAL(KIND=JPRB) :: RFLCC1(6)
70REAL(KIND=JPRB) :: RFLCC2(6)
71REAL(KIND=JPRB) :: RFLCC3(6)
72
73REAL(KIND=JPRB) :: RFUAA0(6)
74REAL(KIND=JPRB) :: RFUAA1(6)
75REAL(KIND=JPRB) :: RFUBB0(6)
76REAL(KIND=JPRB) :: RFUBB1(6)
77REAL(KIND=JPRB) :: RFUBB2(6)
78REAL(KIND=JPRB) :: RFUBB3(6)
79REAL(KIND=JPRB) :: RFUCC0(6)
80REAL(KIND=JPRB) :: RFUCC1(6)
81REAL(KIND=JPRB) :: RFUCC2(6)
82REAL(KIND=JPRB) :: RFUCC3(6)
83REAL(KIND=JPRB) :: RFLDD0(6)
84REAL(KIND=JPRB) :: RFLDD1(6)
85REAL(KIND=JPRB) :: RFLDD2(6)
86REAL(KIND=JPRB) :: RFLDD3(6)
87
88REAL(KIND=JPRB) :: RSUSHE(6)
89REAL(KIND=JPRB) :: RSUSHF(6)
90REAL(KIND=JPRB) :: RSUSHH(6)
91REAL(KIND=JPRB) :: RSUSHK(6)
92REAL(KIND=JPRB) :: RSUSHA(6)
93REAL(KIND=JPRB) :: RSUSHG(6)
94REAL(KIND=JPRB) :: RSUSHFA(4)
95REAL(KIND=JPRB) :: RSUSHC
96REAL(KIND=JPRB) :: RSUSHD
97
98REAL(KIND=JPRB) :: REFFIA
99REAL(KIND=JPRB) :: REFFIB
100REAL(KIND=JPRB) :: RTIW
101REAL(KIND=JPRB) :: RRIW
102REAL(KIND=JPRB) :: RROMA(6)
103REAL(KIND=JPRB) :: RROMB(6)
104REAL(KIND=JPRB) :: RRASY(6)
105
106REAL(KIND=JPRB) :: RHSRA(6)
107REAL(KIND=JPRB) :: RHSRB(6)
108REAL(KIND=JPRB) :: RHSRC(6)
109REAL(KIND=JPRB) :: RHSRD(6)
110REAL(KIND=JPRB) :: RHSRE(6)
111REAL(KIND=JPRB) :: RHSRF(6)
112REAL(KIND=JPRB) :: RHSRTA
113REAL(KIND=JPRB) :: RHSRTB
114
115REAL(KIND=JPRB) :: RTAUA(6,6)
116REAL(KIND=JPRB) :: RPIZA(6,6)
117REAL(KIND=JPRB) :: RCGA(6,6)
118REAL(KIND=JPRB) :: RAER(6,6)
119
120REAL(KIND=JPRB) :: RTWEIGHT(8)
121REAL(KIND=JPRB) :: RWEIGS(6)
122REAL(KIND=JPRB) :: RWEIGV(6)
123
124REAL(KIND=JPRB) :: RWEIGHT(6,8)
125INTEGER(KIND=JPIM) :: NMPSRTM(14), NTYPS
126
127REAL(KIND=JPRB) :: RADJUST
128
129!        * E.C.M.W.F. PHYSICS PACKAGE *
130
131!     J.-J. MORCRETTE       E.C.M.W.F.      89/07/14
132
133!  NAME     TYPE     PURPOSE
134!  ----  :  ----   : ---------------------------------------------------
135!  APAD  :  REAL     PADE APPROXIMANTS NUMERATOR
136!  BPAD  :  REAL     PADE APPROXIMANTS DENOMINATOR
137!  D     :  REAL     TRANSMISSION LIMIT FOR INFINITE ABSORBER AMOUNT
138!  RRAY  :  REAL     RAYLEIGH SCATTERING COEFFICIENTS
139!  RSUN  :  REAL     SOLAR FRACTION IN SPECTRAL INTERVALS
140!  RPDH1 :  1 + EXPONENT PRESSURE DEPENDENCE H2O
141!  RPDU1 :  1 + EXPONENT PRESSURE DEPENDENCE UNIFORMLY MIXED GASES
142!  RPNH  :  REFERENCE PRESSURE FACTOR FOR H2O
143!  RPNU  :  REFERENCE PRESSURE FACTOR FOR UNIFORMLY MIXED GASES
144!  RSWCE :  E-TYPE, H2O CONTINUUM ABSORPTION COEFFICIENT
145!  RSWCP :  P-TYPE, H2O CONTINUUM ABSORPTION COEFFICIENT
146!  RTDH2O:  EXPONENT TEMPERATURE DEPENDENCE H2O
147!  RTDUMG:  EXPONENT TEMPERATURE DEPENDENCE UNIFORMLY MIXED GASES
148!  RTH2O :  REFERENCE TEMPERATURE H2O
149!  RTUMG :  REFERENCE TEMPERATURE UNIFORMLY MIXED GASES
150!     -----------------------------------------------------------------
151
152!        * E.C.M.W.F. PHYSICS PACKAGE *
153
154!     J.-J. MORCRETTE       E.C.M.W.F.      89/07/14
155
156!  NAME     TYPE     PURPOSE
157!  ----  :  ----   : ---------------------------------------------------
158!*    FOUQUART (1987) WATER CLOUD OPTICAL PROPERTIES
159
160! RYFWCA :  REAL   : C1 IN OPTICAL THICKNESS FORMULA
161! RYFWCB :  REAL   : C2 IN OPTICAL THICKNESS FORMULA
162! RYFWCC :  REAL   : SINGLE SCATTERING ALBEDO PARAMETER
163! RYFWCD :  REAL   : SINGLE SCATTERING ALBEDO PARAMETER
164! RYFWCE :  REAL   : SINGLE SCATTERING ALBEDO PARAMETER
165! RYFWCF :  REAL   : ASSYMETRY FACTOR
166
167!*    SLINGO (1989) WATER CLOUD OPTICAL PROPERTIES
168
169! RASWCA :  REAL   : C1 IN OPTICAL THICKNESS FORMULA
170! RASWCB :  REAL   : C2 IN OPTICAL THICKNESS FORMULA
171! RASWCC :  REAL   : SINGLE SCATTERING ALBEDO PARAMETER
172! RASWCD :  REAL   : SINGLE SCATTERING ALBEDO PARAMETER
173! RASWCE :  REAL   : SINGLE SCATTERING ALBEDO PARAMETER
174! RASWCF :  REAL   : ASSYMETRY FACTOR
175
176!*   SAVIJARVI (1998) WATER CLOUD OPTICAL PROPERTIES (RRTM)
177
178! RHSAVI : REAL    : MASS ABSORPTION COEFFICIENTS (POLYNOMIAL DEVELOPM)
179
180!*   LINDNER,LI (2000) WATER CLOUD OPTICAL PROPERTIES (RRTM)
181
182! RLILIA : REAL    : MASS ABSORPTION COEFFICIENTS (POLYNOMIAL DEVELOPM)
183! RLILIB : REAL    : 1-SSA COEFFICIENTS  (POLYNOMIAL DEVELOPM)
184
185!*    ICE CLOUD OPTICAL PROPERTIES DERIVED FROM EBERT-CURRY (1992)
186
187! REBCUA :  REAL   : C1 IN OPTICAL THICKNESS FORMULA
188! REBCUB :  REAL   : C2 IN OPTICAL THICKNESS FORMULA
189! REBCUC :  REAL   : 1-C3  IN SINGLE SCATTERING ALBEDO FORMULA
190! REBCUD :  REAL   : C4 IN SINGLE SCATTERING ALBEDO FORMULA
191! REBCUE :  REAL   : C5 IN ASSYMETRY FACTOR FORMULA
192! REBCUF :  REAL   : C6 IN ASSYMETRY FACTOR FORMULA
193! REBCUG :  REAL   : C7 IN MASS ABSORPTION COEFFICIENT FORMULA
194! REBCUH :  REAL   : C8 IN MASS ABSORPTION COEFFICIENT FORMULA
195! REBCUI :  REAL   : C7 IN MASS ABSORPTION COEFFICIENT SPECTRAL FORMULA
196! REBCUJ :  REAL   : C8 IN MASS ABSORPTION COEFFICIENT SPECTRAL FORMULA
197
198!*    ICE CLOUD OPTICAL PROPERTIES DERIVED FROM SUN-SHINE (1995)
199
200! RSHSUE :  REAL   : E IN SINGLE SCATTERING ALBEDO FORMULA
201! RSHSUF :  REAL   : F IN SINGLE SCATTERING ALBEDO FORMULA
202! RSHSUH :  REAL   : H IN ASSYMETRY FACTOR FORMULA
203! RSHSUK :  REAL   : K IN ASSYMETRY FACTOR FORMULA
204! RSHSUA :  REAL   : ALPHA IN SSA CORRECTION FACTOR FORMULA
205! RSHSUG :  REAL   : GAMMA IN ASSYMETRY CORRECTION FACTOR FORMULA
206! RSHSUFA:  REAL   : COEFFICIENTS IN TEMPERATURE CORRECTION FACTOR
207
208! REFFIA :  REAL   : C9  IN EFFECTIVE RADIUS FORMULA
209! REFFIB :  REAL   : C10 IN EFFECTIVE RADIUS FORMULA
210
211!*    ICE CLOUD OPTICAL PROPERTIES DERIVED FROM FU-LIOU (1993)
212
213! RFULIO :  REAL   : COEFFICIENTS IN EXPRESSION FOR LW EXTINCTION COEFF.
214! RFLAA  :  REAL   : COEFFICIENTS IN EXPRESSION FOR SW EXTINCTION COEFF.
215! RFLBB  :  REAL   : COEFFICIENTS IN EXPRESSION FOR SW SINGLE SCATT.ALB.
216! RFLCC  :  REAL   : COEFFICIENTS IN EXPRESSION FOR SW ASSYMETRY FACTOR
217! RFLDD  :  REAL   : COEFFICIENTS IN EXPRESSION FOR SW ASSYMETRY FACTOR
218
219!*    ICE CLOUD OPTICAL PROPERTIES DERIVED FROM FU (1996) & FU ET AL. (1998)
220
221! RFUETA :  REAL   : COEFFICIENTS IN EXPRESSION FOR LW EXTINCTION COEFF.
222! RFUAA  :  REAL   : COEFFICIENTS IN EXPRESSION FOR SW EXTINCTION COEFF.
223! RFUBB  :  REAL   : COEFFICIENTS IN EXPRESSION FOR SW SINGLE SCATT.ALB.
224! RFUCC  :  REAL   : COEFFICIENTS IN EXPRESSION FOR SW ASSYMETRY FACTOR
225
226!*    TRANSITION BETWEEN LIQUID AND SOLID WATER
227
228! RTIW   :  REAL   : TEMPERATURE THRESHOLD
229! RRIW   :  REAL   : TRANSITION RANGE
230
231!*    RAIN OPTICAL PROPERTIES FROM SAVIJARVI (1996)
232
233! RROMA  :  REAL   : COEFFICIENTS FOR SINGLE SCATTERING ALBEDO
234! RROMB  :  REAL   : COEFFICIENTS FOR SINGLE SCATTERING ALBEDO
235! RRASY  :  REAL   : COEFFICIENTS FOR ASSYMETRY FACTOR
236! RHSRA  :  REAL   : COEFFICIENTS FOR OPTICAL THICKNESS
237! RHSRB  :  REAL   : COEFFICIENTS FOR OPTICAL THICKNESS
238! RHSRC  :  REAL   : COEFFICIENTS FOR SINGLE SCATTERING ALBEDO
239! RHSRD  :  REAL   : COEFFICIENTS FOR SINGLE SCATTERING ALBEDO
240! RHSRE  :  REAL   : COEFFICIENTS FOR ASSYMETRY FACTOR
241! RHSRF  :  REAL   : COEFFICIENTS FOR ASSYMETRY FACTOR
242! RHSRTA :  REAL   : COEFFICIENTS FOR OPTICAL THICKNESS
243! RHSRTB :  REAL   : COEFFICIENTS FOR OPTICAL THICKNESS
244!     -----------------------------------------------------------------
245
246!        * E.C.M.W.F. PHYSICS PACKAGE *
247
248!     J.-J. MORCRETTE       E.C.M.W.F.      89/07/14
249
250!  NAME     TYPE     PURPOSE
251!  ----  :  ----   : -------
252!  RTAUA :  REAL     S.W. NORMALIZED OPTICAL THICKNESS AT 0.55 MICRON
253!  RPIZA :  REAL     S.W. SINGLE SCATTERING ALBEDO
254!  RCGA  :  REAL     S.W. ASSYMETRY FACTOR
255!  RAER  :  REAL     L.W. ABSORPTION COEFFICIENTS
256!     -----------------------------------------------------------------
257
258!        * E.C.M.W.F. PHYSICS PACKAGE *
259
260!     J.-J. MORCRETTE       E.C.M.W.F.      89/07/14
261
262!  NAME     TYPE     PURPOSE
263!  ----  :  ----   : -------
264!RTWEIGHT:  REAL     S.W. INTEGRATED WEIGHT
265! NMPSRTM: INTEGER  : Indices for mapping SW[1:6] albedo into SRTM[1:14] 
266!     -----------------------------------------------------------------
267!$OMP THREADPRIVATE(apad,bpad,d,nexpo3,nmpsrtm,ntyps,radjust,raer,raswca)
268!$OMP THREADPRIVATE(raswcb,raswcc,raswcd,raswce,raswcf,rcga,rebcua,rebcub)
269!$OMP THREADPRIVATE(rebcuc,rebcud,rebcue,rebcuf,rebcug,rebcuh,rebcui,rebcuj)
270!$OMP THREADPRIVATE(reffia,reffib,rexpo3,rflaa0,rflaa1,rflbb0,rflbb1,rflbb2)
271!$OMP THREADPRIVATE(rflbb3,rflcc0,rflcc1,rflcc2,rflcc3,rfldd0,rfldd1,rfldd2)
272!$OMP THREADPRIVATE(rfldd3,rfuaa0,rfuaa1,rfubb0,rfubb1,rfubb2,rfubb3,rfucc0)
273!$OMP THREADPRIVATE(rfucc1,rfucc2,rfucc3,rfueta,rfuetb,rfuetc,rfulio,rhsavi)
274!$OMP THREADPRIVATE(rhsra,rhsrb,rhsrc,rhsrd,rhsre,rhsrf,rhsrta,rhsrtb,rlilia)
275!$OMP THREADPRIVATE(rlilib,rpdh1,rpdu1,rpiza,rpnh,rpnu,rrasy,rray,rriw,rroma)
276!$OMP THREADPRIVATE(rromb,rsusha,rsushc,rsushd,rsushe,rsushf,rsushfa,rsushg)
277!$OMP THREADPRIVATE(rsushh,rsushk,rswce,rswcp,rtaua,rtdh2o,rtdumg,rth2o,rtiw)
278!$OMP THREADPRIVATE(rtumg,rtweight,rweight,rweigs,rweigv,ryfwca,ryfwcb,ryfwcc,ryfwcd,ryfwce,ryfwcf)
279END MODULE YOESW
Note: See TracBrowser for help on using the repository browser.