source: LMDZ5/branches/LMDZ5_SPLA/libf/phylmd/rrtm/yoesw.F90 @ 4934

Last change on this file since 4934 was 2010, checked in by Laurent Fairhead, 11 years ago

Modifications pour OpenMP

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