source: LMDZ6/trunk/libf/phylmd/ecrad.v1.5.1/yoesw.F90 @ 5322

Last change on this file since 5322 was 3908, checked in by idelkadi, 3 years ago

Online implementation of the radiative transfer code ECRAD in the LMDZ model.

  • Inclusion of the ecrad directory containing the sources of the ECRAD code
    • interface routine : radiation_scheme.F90
  • Adaptation of compilation scripts :
    • compilation under CPP key CPP_ECRAD
    • compilation with option "-rad ecard" or "-ecard true"
    • The "-rad old/rtm/ecran" build option will need to replace the "-rrtm true" and "-ecrad true" options in the future.
  • Runing LMDZ simulations with ecrad, you need :
    • logical key iflag_rrtm = 2 in physiq.def
    • namelist_ecrad (DefLists?)
    • the directory "data" containing the configuration files is temporarily placed in ../libfphylmd/ecrad/
  • Compilation and execution are tested in the 1D case. The repository under svn would allow to continue the implementation work: tests, verification of the results, ...
File size: 7.4 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) :: REFFIA
88
89REAL(KIND=JPRB) :: RTAUA(6,6)
90REAL(KIND=JPRB) :: RPIZA(6,6)
91REAL(KIND=JPRB) :: RCGA(6,6)
92REAL(KIND=JPRB) :: RAER(6,6)
93
94INTEGER(KIND=JPIM) :: NMPSRTM(14), NTYPS
95
96REAL(KIND=JPRB) :: RADJUST
97
98!        * E.C.M.W.F. PHYSICS PACKAGE *
99
100!     J.-J. MORCRETTE       E.C.M.W.F.      89/07/14
101
102!  NAME     TYPE     PURPOSE
103!  ----  :  ----   : ---------------------------------------------------
104!  APAD  :  REAL     PADE APPROXIMANTS NUMERATOR
105!  BPAD  :  REAL     PADE APPROXIMANTS DENOMINATOR
106!  D     :  REAL     TRANSMISSION LIMIT FOR INFINITE ABSORBER AMOUNT
107!  RRAY  :  REAL     RAYLEIGH SCATTERING COEFFICIENTS
108!  RSUN  :  REAL     SOLAR FRACTION IN SPECTRAL INTERVALS
109!  RPDH1 :  1 + EXPONENT PRESSURE DEPENDENCE H2O
110!  RPDU1 :  1 + EXPONENT PRESSURE DEPENDENCE UNIFORMLY MIXED GASES
111!  RPNH  :  REFERENCE PRESSURE FACTOR FOR H2O
112!  RPNU  :  REFERENCE PRESSURE FACTOR FOR UNIFORMLY MIXED GASES
113!  RSWCE :  E-TYPE, H2O CONTINUUM ABSORPTION COEFFICIENT
114!  RSWCP :  P-TYPE, H2O CONTINUUM ABSORPTION COEFFICIENT
115!  RTDH2O:  EXPONENT TEMPERATURE DEPENDENCE H2O
116!  RTDUMG:  EXPONENT TEMPERATURE DEPENDENCE UNIFORMLY MIXED GASES
117!  RTH2O :  REFERENCE TEMPERATURE H2O
118!  RTUMG :  REFERENCE TEMPERATURE UNIFORMLY MIXED GASES
119!     -----------------------------------------------------------------
120
121!        * E.C.M.W.F. PHYSICS PACKAGE *
122
123!     J.-J. MORCRETTE       E.C.M.W.F.      89/07/14
124
125!  NAME     TYPE     PURPOSE
126!  ----  :  ----   : ---------------------------------------------------
127!*    FOUQUART (1987) WATER CLOUD OPTICAL PROPERTIES
128
129! RYFWCA :  REAL   : C1 IN OPTICAL THICKNESS FORMULA
130! RYFWCB :  REAL   : C2 IN OPTICAL THICKNESS FORMULA
131! RYFWCC :  REAL   : SINGLE SCATTERING ALBEDO PARAMETER
132! RYFWCD :  REAL   : SINGLE SCATTERING ALBEDO PARAMETER
133! RYFWCE :  REAL   : SINGLE SCATTERING ALBEDO PARAMETER
134! RYFWCF :  REAL   : ASSYMETRY FACTOR
135
136!*    SLINGO (1989) WATER CLOUD OPTICAL PROPERTIES
137
138! RASWCA :  REAL   : C1 IN OPTICAL THICKNESS FORMULA
139! RASWCB :  REAL   : C2 IN OPTICAL THICKNESS FORMULA
140! RASWCC :  REAL   : SINGLE SCATTERING ALBEDO PARAMETER
141! RASWCD :  REAL   : SINGLE SCATTERING ALBEDO PARAMETER
142! RASWCE :  REAL   : SINGLE SCATTERING ALBEDO PARAMETER
143! RASWCF :  REAL   : ASSYMETRY FACTOR
144
145!*   LINDNER,LI (2000) WATER CLOUD OPTICAL PROPERTIES (RRTM)
146
147! RLILIA : REAL    : MASS ABSORPTION COEFFICIENTS (POLYNOMIAL DEVELOPM)
148! RLILIB : REAL    : 1-SSA COEFFICIENTS  (POLYNOMIAL DEVELOPM)
149
150!*    ICE CLOUD OPTICAL PROPERTIES DERIVED FROM EBERT-CURRY (1992)
151
152! REBCUA :  REAL   : C1 IN OPTICAL THICKNESS FORMULA
153! REBCUB :  REAL   : C2 IN OPTICAL THICKNESS FORMULA
154! REBCUC :  REAL   : 1-C3  IN SINGLE SCATTERING ALBEDO FORMULA
155! REBCUD :  REAL   : C4 IN SINGLE SCATTERING ALBEDO FORMULA
156! REBCUE :  REAL   : C5 IN ASSYMETRY FACTOR FORMULA
157! REBCUF :  REAL   : C6 IN ASSYMETRY FACTOR FORMULA
158! REBCUG :  REAL   : C7 IN MASS ABSORPTION COEFFICIENT FORMULA
159! REBCUH :  REAL   : C8 IN MASS ABSORPTION COEFFICIENT FORMULA
160! REBCUI :  REAL   : C7 IN MASS ABSORPTION COEFFICIENT SPECTRAL FORMULA
161! REBCUJ :  REAL   : C8 IN MASS ABSORPTION COEFFICIENT SPECTRAL FORMULA
162
163!*    ICE CLOUD OPTICAL PROPERTIES DERIVED FROM SUN-SHINE (1995)
164
165! RSHSUE :  REAL   : E IN SINGLE SCATTERING ALBEDO FORMULA
166! RSHSUF :  REAL   : F IN SINGLE SCATTERING ALBEDO FORMULA
167! RSHSUH :  REAL   : H IN ASSYMETRY FACTOR FORMULA
168! RSHSUK :  REAL   : K IN ASSYMETRY FACTOR FORMULA
169! RSHSUA :  REAL   : ALPHA IN SSA CORRECTION FACTOR FORMULA
170! RSHSUG :  REAL   : GAMMA IN ASSYMETRY CORRECTION FACTOR FORMULA
171! RSHSUFA:  REAL   : COEFFICIENTS IN TEMPERATURE CORRECTION FACTOR
172
173! REFFIA :  REAL   : C9  IN EFFECTIVE RADIUS FORMULA
174
175!*    ICE CLOUD OPTICAL PROPERTIES DERIVED FROM FU-LIOU (1993)
176
177! RFULIO :  REAL   : COEFFICIENTS IN EXPRESSION FOR LW EXTINCTION COEFF.
178! RFLAA  :  REAL   : COEFFICIENTS IN EXPRESSION FOR SW EXTINCTION COEFF.
179! RFLBB  :  REAL   : COEFFICIENTS IN EXPRESSION FOR SW SINGLE SCATT.ALB.
180! RFLCC  :  REAL   : COEFFICIENTS IN EXPRESSION FOR SW ASSYMETRY FACTOR
181! RFLDD  :  REAL   : COEFFICIENTS IN EXPRESSION FOR SW ASSYMETRY FACTOR
182
183!*    ICE CLOUD OPTICAL PROPERTIES DERIVED FROM FU (1996) & FU ET AL. (1998)
184
185! RFUETA :  REAL   : COEFFICIENTS IN EXPRESSION FOR LW EXTINCTION COEFF.
186! RFUAA  :  REAL   : COEFFICIENTS IN EXPRESSION FOR SW EXTINCTION COEFF.
187! RFUBB  :  REAL   : COEFFICIENTS IN EXPRESSION FOR SW SINGLE SCATT.ALB.
188! RFUCC  :  REAL   : COEFFICIENTS IN EXPRESSION FOR SW ASSYMETRY FACTOR
189
190!     -----------------------------------------------------------------
191
192!        * E.C.M.W.F. PHYSICS PACKAGE *
193
194!     J.-J. MORCRETTE       E.C.M.W.F.      89/07/14
195
196!  NAME     TYPE     PURPOSE
197!  ----  :  ----   : -------
198!  RTAUA :  REAL     S.W. NORMALIZED OPTICAL THICKNESS AT 0.55 MICRON
199!  RPIZA :  REAL     S.W. SINGLE SCATTERING ALBEDO
200!  RCGA  :  REAL     S.W. ASSYMETRY FACTOR
201!  RAER  :  REAL     L.W. ABSORPTION COEFFICIENTS
202!     -----------------------------------------------------------------
203
204!        * E.C.M.W.F. PHYSICS PACKAGE *
205
206!     J.-J. MORCRETTE       E.C.M.W.F.      89/07/14
207
208!  NAME     TYPE     PURPOSE
209!  ----  :  ----   : -------
210!RTWEIGHT:  REAL     S.W. INTEGRATED WEIGHT
211! NMPSRTM: INTEGER  : Indices for mapping SW[1:6] albedo into SRTM[1:14] 
212!     -----------------------------------------------------------------
213END MODULE YOESW
Note: See TracBrowser for help on using the repository browser.