source: dynamico_lmdz/aquaplanet/LMDZ5/libf/phymar/susw.F90 @ 4031

Last change on this file since 4031 was 3990, checked in by millour, 8 years ago

An updated LMDZ5 (based on LMDZ rev 2786) to play with.
EM
File size: 4.7 KB
RevLine 
[3990]1SUBROUTINE SUSW
2
3!**** *SUSW*   - INITIALIZE COMMON YOESW
4
5!     PURPOSE.
6!     --------
7!           INITIALIZE YOESW, THE COMMON THAT CONTAINS COEFFICIENTS
8!           NEEDED TO RUN THE SHORTWAVE RADIATION SUBROUTINES
9
10!**   INTERFACE.
11!     ----------
12!        *CALL* *SUSW
13
14!        EXPLICIT ARGUMENTS :
15!        --------------------
16!        NONE
17
18!        IMPLICIT ARGUMENTS :
19!        --------------------
20!        COMMON YOESW
21
22!     METHOD.
23!     -------
24!        SEE DOCUMENTATION
25
26!     EXTERNALS.
27!     ----------
28
29!     REFERENCE.
30!     ----------
31!        ECMWF RESEARCH DEPARTMENT DOCUMENTATION OF THE IFS
32
33!     AUTHOR.
34!     -------
35!        JEAN-JACQUES MORCRETTE *ECMWF*
36
37!     MODIFICATIONS.
38!     --------------
39!        ORIGINAL : 88-12-15
40
41!     ------------------------------------------------------------------
42
43#include "tsmbkind.h"
44
45USE YOMCST   , ONLY : RG
46USE YOESW    , ONLY : APAD     ,BPAD     ,RRAY     ,RSUN     ,&
47            &RPDH1    ,RPDU1    ,RPNH     ,RPNU     ,RSWCE    ,&
48            &RSWCP    ,RTDH2O   ,RTDUMG   ,RTH2O    ,RTUMG    ,D
49
50!      ----------------------------------------------------------------
51
52!*       1.    SET VALUES.
53!              -----------
54
55
56
57IMPLICIT NONE
58REAL_B   :: ZPDH2O=0.8_JPRB
59REAL_B   :: ZPDUMG=0.75_JPRB
60REAL_B   :: ZPRH2O=30000._JPRB
61REAL_B   :: ZPRUMG=30000._JPRB
62
63!     LOCAL REAL SCALARS
64REAL_B :: ZH2O, ZUMG
65
66RPDH1=ZPDH2O+_ONE_
67RPDU1=ZPDUMG+_ONE_
68ZH2O=_ONE_/( 10._JPRB* RG * RPDH1 )
69ZUMG=_ONE_/( 10._JPRB* RG * RPDU1 )
70RPNU = ZUMG/(ZPRUMG**ZPDUMG)
71RPNH = ZH2O/(ZPRH2O**ZPDH2O)
72
73!      RSWCP=0.002*RSWCE
74
75
76RTDH2O = 0.40_JPRB
77RTDUMG = 0.375_JPRB
78RTH2O = 240._JPRB
79RTUMG = 240._JPRB
80RSWCE = 0._JPRB
81RSWCP = 0._JPRB
82
83RSUN(1) = 0.441676_JPRB
84D(1,:) = (/ 0.00_JPRB, 0.00_JPRB, 0.00_JPRB /)
85!* DERIVED FROM HITRAN APRIL 1991
86!       H2O:  Pref=300 hPa, Tref=240K, Pdep=0.8
87!       O3 :  unchanged
88
89APAD(1, 1, :) = (/&
90 &0.912418292E+05_JPRB,&
91 &0.723613782E+05_JPRB,&
92 &0.596037057E+04_JPRB,&
93 &0.000000000E-00_JPRB,&
94 &0.000000000E-00_JPRB,&
95 &0.000000000E-00_JPRB,&
96 &0.000000000E-00_JPRB/)
97APAD(1, 2, :) = (/&
98 &0.000000000E-00_JPRB,&
99 &0.000000000E-00_JPRB,&
100 &0.000000000E-00_JPRB,&
101 &0.000000000E-00_JPRB,&
102 &0.000000000E-00_JPRB,&
103 &0.000000000E-00_JPRB,&
104 &0.000000000E-00_JPRB/)
105APAD(1, 3, :) = (/&
106 &0.925887084E-04_JPRB,&
107 &0.129353723E-01_JPRB,&
108 &0.800821928E+00_JPRB,&
109 &0.242715973E+02_JPRB,&
110 &0.878331486E+02_JPRB,&
111 &0.191559725E+02_JPRB,&
112 &0.000000000E+00_JPRB/)
113
114BPAD(1, 1, :) = (/&
115 &0.912418292E+05_JPRB,&
116 &0.724555318E+05_JPRB,&
117 &0.602593328E+04_JPRB,&
118 &0.100000000E+01_JPRB,&
119 &0.000000000E-00_JPRB,&
120 &0.000000000E-00_JPRB,&
121 &0.000000000E-00_JPRB/)
122BPAD(1, 2, :) = (/&
123 &0.000000000E-00_JPRB,&
124 &0.000000000E-00_JPRB,&
125 &0.000000000E-00_JPRB,&
126 &0.000000000E-00_JPRB,&
127 &0.000000000E-00_JPRB,&
128 &0.000000000E-00_JPRB,&
129 &0.000000000E-00_JPRB/)
130BPAD(1, 3, :) = (/&
131 &0.925887084E-04_JPRB,&
132 &0.131812683E-01_JPRB,&
133 &0.812706117E+00_JPRB,&
134 &0.249863591E+02_JPRB,&
135 &0.931071925E+02_JPRB,&
136 &0.252233437E+02_JPRB,&
137 &0.100000000E+01_JPRB/)
138
139RRAY(1,1:6) = (/&
140 &.428937E-01_JPRB, .890743E+00_JPRB,-.288555E+01_JPRB,&
141 &.522744E+01_JPRB,-.469173E+01_JPRB, .161645E+01_JPRB/)
142
143RSUN(2) = 0.558324_JPRB
144
145D(2,1:3) = (/ 0.000000000_JPRB, 0.000000000_JPRB, 0.800000000_JPRB /)
146
147!* INTERVAL 2:  0.68 - 4.00 MICRONS
148!* DERIVED FROM HITRAN APRIL 1991
149!       H2O:  Pref=300 hPa, Tref=240K, Pdep=0.80
150!       UMG:  Pref=300 hPa, Tref=240K, Pdep=0.75 (CO2+O2+CH4+N2O+CO)
151!       O3 :  unchanged
152
153APAD(2, 1, :) = (/&
154 &0.376655383E-08_JPRB,&
155 &0.978576773E-04_JPRB,&
156 &0.387714006E+00_JPRB,&
157 &0.118461660E+03_JPRB,&
158 &0.119079797E+04_JPRB,&
159 &0.293353397E+03_JPRB,&
160 &0.000000000E+00_JPRB/)
161APAD(2, 2, :) = (/&
162 &0.739646016E-08_JPRB,&
163 &0.131849595E-03_JPRB,&
164 &0.437772681E+00_JPRB,&
165 &0.151345118E+03_JPRB,&
166 &0.233628890E+04_JPRB,&
167 &0.797219934E+03_JPRB,&
168 &0.000000000E+00_JPRB/)
169APAD(2, 3, :) = (/&
170 &0.410177786E+03_JPRB,&
171 &0.672595424E+02_JPRB,&
172 &0.000000000E-00_JPRB,&
173 &0.000000000E-00_JPRB,&
174 &0.000000000E-00_JPRB,&
175 &0.000000000E-00_JPRB,&
176 &0.000000000E+00_JPRB/)
177
178BPAD(2, 1, :) = (/&
179 &0.376655383E-08_JPRB,&
180 &0.979023421E-04_JPRB,&
181 &0.388611139E+00_JPRB,&
182 &0.120291383E+03_JPRB,&
183 &0.130531005E+04_JPRB,&
184 &0.415049409E+03_JPRB,&
185 &0.100000000E+01_JPRB/)
186BPAD(2, 2, :) = (/&
187 &0.739646016E-08_JPRB,&
188 &0.131861712E-03_JPRB,&
189 &0.437949001E+00_JPRB,&
190 &0.151692730E+03_JPRB,&
191 &0.237071130E+04_JPRB,&
192 &0.867914360E+03_JPRB,&
193 &0.100000000E+01_JPRB/)
194BPAD(2, 3, :) = (/&
195 &0.410177786E+03_JPRB,&
196 &0.731185438E+02_JPRB,&
197 &0.100000000E+01_JPRB,&
198 &0.000000000E+00_JPRB,&
199 &0.000000000E+00_JPRB,&
200 &0.000000000E+00_JPRB,&
201 &0.000000000E+00_JPRB/)
202
203RRAY(2,1:6) = (/&
204 &.697200E-02_JPRB, .173297E-01_JPRB,-.850903E-01_JPRB,&
205 &.248261E+00_JPRB,-.302031E+00_JPRB, .129662E+00_JPRB/)
206RETURN
207END SUBROUTINE SUSW
Note: See TracBrowser for help on using the repository browser.