1 | !!---------------------------------------------------------------------- |
---|
2 | !! *** flxrnf_ORCA_R05.h90 *** |
---|
3 | !!---------------------------------------------------------------------- |
---|
4 | !! flx_rnf : update ocean runoffs for ORCA_R05 configuration |
---|
5 | !!---------------------------------------------------------------------- |
---|
6 | !! * Modules used |
---|
7 | USE dom_oce ! ocean space and time domain |
---|
8 | USE phycst ! physical constants |
---|
9 | USE in_out_manager ! I/O manager |
---|
10 | USE daymod ! calendar |
---|
11 | |
---|
12 | IMPLICIT NONE |
---|
13 | PRIVATE |
---|
14 | |
---|
15 | !! * Routine accessibility |
---|
16 | PUBLIC flx_rnf ! routine call by step.F90 |
---|
17 | |
---|
18 | !! * Shared module variables |
---|
19 | REAL(wp), PUBLIC, DIMENSION(jpi,jpj) :: & !: |
---|
20 | runoff, & !: monthly runoff (kg/m2/s) |
---|
21 | upsadv, & !: mixed adv scheme in straits vicinity (hori.) |
---|
22 | upsrnfh !: mixed adv scheme in runoffs vicinity (hori.) |
---|
23 | REAL(wp), PUBLIC, DIMENSION(jpk) :: & !: |
---|
24 | upsrnfz !: mixed adv scheme in runoffs vicinity (vert.) |
---|
25 | INTEGER , PUBLIC :: & !: |
---|
26 | numrof = 48 , & !: logical unit for runoff data |
---|
27 | nrunoff = 0 , & !: runoff option (namelist) |
---|
28 | nrnf1, nrnf2 !: first and second record used |
---|
29 | |
---|
30 | !! * Module variable |
---|
31 | REAL(wp), DIMENSION(jpi,jpj,2) :: & !: |
---|
32 | rnfdta !: monthly runoff data array (kg/m2/s) |
---|
33 | !!---------------------------------------------------------------------- |
---|
34 | !! OPA 9.0 , LODYC-IPSL (2003) |
---|
35 | !!---------------------------------------------------------------------- |
---|
36 | |
---|
37 | CONTAINS |
---|
38 | |
---|
39 | SUBROUTINE flx_rnf( kt ) |
---|
40 | !!---------------------------------------------------------------------- |
---|
41 | !! *** ROUTINE flx_rnf *** |
---|
42 | !! |
---|
43 | !! ** Purpose : Introduce a climatological run off forcing |
---|
44 | !! |
---|
45 | !! ** Method : ORCA R05 configuration |
---|
46 | !! Initialze each mouth of river with a monthly climatology |
---|
47 | !! provided from different data. |
---|
48 | !! C a u t i o n : upward water flux, runoff is negative |
---|
49 | !! set at the last loop of the routine |
---|
50 | !! |
---|
51 | !! ** Action : |
---|
52 | !! |
---|
53 | !! References : |
---|
54 | !! J. D. Milliman and R. H. Meade, 1983 : world-wide delivery |
---|
55 | !! of river sediment to the oceans, journal of geology vol 91 |
---|
56 | !! pp 1-21. |
---|
57 | !! G. L. Russell and J. R. Miller, 1990 : global river runoff |
---|
58 | !! calculated from a global atmospheric general circulation |
---|
59 | !! model, journal of hydrology, 117(1990), pp 241-254. |
---|
60 | !! F. Van Der Leeden, Troise F. L., Todd D. K. : the water |
---|
61 | !! encyclopedia, second edition, lewis publishers. |
---|
62 | !! J. W. Weatherly, J. E. Walsh : The effects of precipitation |
---|
63 | !! and river runoff in a coupled ice-ocean model of Arctic |
---|
64 | !! Climate dynamics 1996 12:785,798 |
---|
65 | !! Jacobs et al. 1992. J. Glaciol. 38 (130) 375-387. |
---|
66 | !! |
---|
67 | !! History : |
---|
68 | !! ! 94-10 (G.Madec, M. Pontaud, M. Imbard) Original code |
---|
69 | !! ! 97-03 (G.Madec) time dependent version |
---|
70 | !! ! 98-06 (J.M. Molines) exact computation of zxy |
---|
71 | !! for months that are not 30 days |
---|
72 | !! ! 98-07 (M. Imbard) ORCA and mpp option |
---|
73 | !! ! 99-08 (J.P. Boulanger H.L.Ayina) New rivers and |
---|
74 | !! values given in m3/s |
---|
75 | !! ! 00-04 (G. Madec, K. Roberts) add antarctica ice discharge. |
---|
76 | !! ! 00-11 (R. Hordoir, E. Durand) NetCDF FORMAT |
---|
77 | !! 8.5 ! 02-09 (G. Madec) F90: Free form and module |
---|
78 | !!---------------------------------------------------------------------- |
---|
79 | !! * Modules used |
---|
80 | USE ioipsl |
---|
81 | |
---|
82 | !! * arguments |
---|
83 | INTEGER, INTENT( in ) :: kt ! ocean time step |
---|
84 | |
---|
85 | #if defined key_coupled |
---|
86 | runoff(:,:) = 0.0e0 |
---|
87 | #else |
---|
88 | !!---------------------------------------------------------------------- |
---|
89 | !! ORCA_R05 |
---|
90 | !!---------------------------------------------------------------------- |
---|
91 | !! * Local declarations |
---|
92 | INTEGER, PARAMETER :: & |
---|
93 | jpmois=12, & ! number of month in the year ! |
---|
94 | jpriv=120, & ! maximum number of rivers |
---|
95 | jpcoef=20 ! maximum number of gridpoints for mouth rivers |
---|
96 | |
---|
97 | INTEGER :: ji, jj, jl, jm, jn, jr, jc |
---|
98 | INTEGER :: i15, imois, iman |
---|
99 | INTEGER :: inbp, idbd, idmeom |
---|
100 | INTEGER :: iirnf(jpcoef,jpriv),ijrnf(jpcoef,jpriv) |
---|
101 | INTEGER :: inb(jpriv) |
---|
102 | |
---|
103 | REAL(wp) :: zconv, zxy |
---|
104 | REAL(wp) :: zrup1,zrup2,zrup3 |
---|
105 | REAL(wp) :: zrnfm(jpmois, jpriv) |
---|
106 | REAL(wp) :: zcoefr(jpi, jpj, 5) |
---|
107 | !!---------------------------------------------------------------------- |
---|
108 | !! OPA 8.5, LODYC-IPSL (2002) |
---|
109 | !!---------------------------------------------------------------------- |
---|
110 | ! ==== |
---|
111 | ! DATA |
---|
112 | ! ==== |
---|
113 | |
---|
114 | ! I NORTH AMERICA |
---|
115 | ! --------------- |
---|
116 | |
---|
117 | ! 1-Saint Lawrence (Canada) 49N30 67W00 R2 |
---|
118 | ! Old=447 km3/year=14165m3/s; CLIPPER=12430m3/s; ratio=1.14 |
---|
119 | ! --> 13975.249m3/s |
---|
120 | DATA inb(1)/3/ |
---|
121 | DATA (iirnf(jc,1),jc=1,jpcoef)/ 444, 444, 444, 17*0 / |
---|
122 | DATA (ijrnf(jc,1),jc=1,jpcoef)/ 367, 368, 369, 17*0 / |
---|
123 | DATA (zrnfm(jm,1),jm=1,jpmois)/ & |
---|
124 | 11718., 11971., 12466., 14980., 19257., 16913., & |
---|
125 | 14211., 13119., 13367., 13644., 13463., 12594./ |
---|
126 | ! 2-Hudson River (USA) 41N00 77W00 R1 |
---|
127 | ! Old=12 km3/year=380m3/s; UNESCO(64-83,p255)=408m3/s |
---|
128 | ! --> 407.416687m3/s |
---|
129 | DATA inb(2)/1/ |
---|
130 | DATA (iirnf(jc,2),jc=1,jpcoef)/ 426, 19*0 / |
---|
131 | DATA (ijrnf(jc,2),jc=1,jpcoef)/ 336, 19*0 / |
---|
132 | DATA (zrnfm(jm,2),jm=1,jpmois)/ & |
---|
133 | 364., 404., 648., 873., 593., 301., & |
---|
134 | 205., 171., 205., 290., 379., 456./ |
---|
135 | ! 3-Missisipi and Atchafalaya (USA) 29N 90W R1 |
---|
136 | ! Old=580 km3/year=18380m3/s; CLIPPER=17599m3/s |
---|
137 | ! --> 17599.334m3/s |
---|
138 | DATA inb(3)/2/ |
---|
139 | DATA (iirnf(jc,3),jc=1,jpcoef)/ 397, 398, 18*0 / |
---|
140 | DATA (ijrnf(jc,3),jc=1,jpcoef)/ 312, 313, 18*0 / |
---|
141 | DATA (zrnfm(jm,3),jm=1,jpmois)/ & |
---|
142 | 19260., 20813., 23026., 28160., 26510., 19650., & |
---|
143 | 14181., 10676., 9578., 10167., 11699., 17472./ |
---|
144 | ! 4-Brazos (USA) 29N58 95W53 R1 |
---|
145 | ! Old=7 km3/year=222m3/s; UNESCO(65-84, p298)=200m3/s |
---|
146 | ! --> 199.666656m3/s |
---|
147 | DATA inb(4)/1/ |
---|
148 | DATA (iirnf(jc,4),jc=1,jpcoef)/ 386, 19*0 / |
---|
149 | DATA (ijrnf(jc,4),jc=1,jpcoef)/ 311, 19*0 / |
---|
150 | DATA (zrnfm(jm,4),jm=1,jpmois)/ & |
---|
151 | 165., 224., 200., 266., 457., 327., & |
---|
152 | 143., 69., 102., 132., 160., 151./ |
---|
153 | ! 5-Colorado (Mexico) 31N54 114W57 R5 |
---|
154 | ! Old=20 km3/year=634m3/s |
---|
155 | ! --> 634.m3/s |
---|
156 | DATA inb(5)/2/ |
---|
157 | DATA (iirnf(jc,5),jc=1,jpcoef)/ 347, 348, 18*0 / |
---|
158 | DATA (ijrnf(jc,5),jc=1,jpcoef)/ 316, 316, 18*0 / |
---|
159 | DATA (zrnfm(jm,5),jm=1,jpmois)/12*634./ |
---|
160 | ! 6-Columbia (USA) 45N60 121W17 R1 |
---|
161 | ! Old=251 km3/year=7954m3/s; UNESCO(65-84, p312)=5368m3/s |
---|
162 | ! --> 5368.41602m3/s |
---|
163 | DATA inb(6)/2/ |
---|
164 | DATA (iirnf(jc,6),jc=1,jpcoef)/ 326, 326, 18*0 / |
---|
165 | DATA (ijrnf(jc,6),jc=1,jpcoef)/ 360, 361, 18*0 / |
---|
166 | DATA (zrnfm(jm,6),jm=1,jpmois)/ & |
---|
167 | 4752., 5077., 5633., 6075., 8422., 9822., & |
---|
168 | 6350., 4112., 3273., 3251., 3548., 4106./ |
---|
169 | ! 7-Fraser (USA) R5 |
---|
170 | ! Old=112 km3/year=3550m3/s;not found in UNESCO |
---|
171 | ! --> 3549.99976m3/s |
---|
172 | DATA inb(7)/2/ |
---|
173 | DATA (iirnf(jc,7),jc=1,jpcoef)/ 317, 317, 18*0 / |
---|
174 | DATA (ijrnf(jc,7),jc=1,jpcoef)/ 378, 379, 18*0 / |
---|
175 | DATA (zrnfm(jm,7),jm=1,jpmois)/12*3550./ |
---|
176 | ! 8-Yukon (USA) 63N00 165W00 R3 |
---|
177 | ! Old=195 km3/year=6179m3/s; UNESCO(65-77+80-84, p321)=4245; ratio=1.46 |
---|
178 | ! --> 6017.24951m3/s |
---|
179 | DATA inb(8)/2/ |
---|
180 | DATA (iirnf(jc,8),jc=1,jpcoef)/ 244, 244, 18*0 / |
---|
181 | DATA (ijrnf(jc,8),jc=1,jpcoef)/ 413, 414, 18*0 / |
---|
182 | DATA (zrnfm(jm,8),jm=1,jpmois)/ & |
---|
183 | 1459., 1213., 1065., 1100., 9067., 17675., & |
---|
184 | 10680., 11271., 8658., 5453., 2649., 1917./ |
---|
185 | ! 9-Copper (USA) 61N47 144W45 R2 |
---|
186 | ! Old=39 km3/year=1236m3/s; UNESCO(65-84, p317)=1052m3/s; ratio=1.18 |
---|
187 | ! --> 1236.25m3/s |
---|
188 | DATA inb(9)/1/ |
---|
189 | DATA (iirnf(jc,9),jc=1,jpcoef)/ 274, 19*0 / |
---|
190 | DATA (ijrnf(jc,9),jc=1,jpcoef)/ 409, 19*0 / |
---|
191 | DATA (zrnfm(jm,9),jm=1,jpmois)/ & |
---|
192 | 219., 195., 181., 214., 904., 2532., & |
---|
193 | 4026., 3627., 1632., 697., 354., 254./ |
---|
194 | ! 10-Susitna (USA) R5 |
---|
195 | ! Old=40 km3/year=1268m3/s |
---|
196 | ! --> 1268.m3/s |
---|
197 | DATA inb(10)/2/ |
---|
198 | DATA (iirnf(jc,10),jc=1,jpcoef)/ 263, 264, 18*0 / |
---|
199 | DATA (ijrnf(jc,10),jc=1,jpcoef)/ 405, 405, 18*0 / |
---|
200 | DATA (zrnfm(jm,10),jm=1,jpmois)/12*1268./ |
---|
201 | ! 11-McKenzie (Canada) 70N00 135W00 R3 |
---|
202 | ! Old=306 km3/year=9697m3/s; UNESCO(67+72+7479+82-84, p168)=8343m3/s; |
---|
203 | ! ratio=1.16 |
---|
204 | ! --> 9606.66797m3/s |
---|
205 | DATA inb(11)/3/ |
---|
206 | DATA (iirnf(jc,11),jc=1,jpcoef)/ 272, 272, 272, 17*0 / |
---|
207 | DATA (ijrnf(jc,11),jc=1,jpcoef)/ 451, 452, 453, 17*0 / |
---|
208 | DATA (zrnfm(jm,11),jm=1,jpmois)/ & |
---|
209 | 3644., 3400., 3258., 3546., 14188., 21544., & |
---|
210 | 18954., 15487., 11284., 10169., 5901., 3905./ |
---|
211 | ! 12-Kuksokwim (USA) 60N00 ~162W30 R4 |
---|
212 | ! UNESCO(65-84, p318)=1112m3/s (61N87 158E10, Crooked Creek, inland) |
---|
213 | ! --> 1111.99988m3/s |
---|
214 | DATA inb(12)/2/ |
---|
215 | DATA (iirnf(jc,12),jc=1,jpcoef)/ 244, 245, 18*0 / |
---|
216 | DATA (ijrnf(jc,12),jc=1,jpcoef)/ 413, 414, 18*0 / |
---|
217 | DATA (zrnfm(jm,12),jm=1,jpmois)/ & |
---|
218 | 335., 308., 295., 407., 2415., 2027., & |
---|
219 | 1917., 2049., 1585., 1072., 545., 389./ |
---|
220 | |
---|
221 | DATA inb(13)/0/ |
---|
222 | DATA inb(14)/0/ |
---|
223 | DATA inb(15)/0/ |
---|
224 | |
---|
225 | ! Total run-offs(I)=60974.2539m3/s |
---|
226 | |
---|
227 | ! II SOUTH AMERICA |
---|
228 | ! ----------------- |
---|
229 | ! 16-Chira (Peru) 4s54 81w08 R2 |
---|
230 | ! Old=5 km3/year=158m3/s; UNESCO (65-68, p357)=80.9; ratio=1.95 |
---|
231 | ! --> 158.000015m3/s |
---|
232 | DATA inb(16)/1/ |
---|
233 | DATA (iirnf(jc,16),jc=1,jpcoef)/ 414, 19*0 / |
---|
234 | DATA (ijrnf(jc,16),jc=1,jpcoef)/ 244, 19*0 / |
---|
235 | DATA (zrnfm(jm,16),jm=1,jpmois)/ & |
---|
236 | 70., 105., 564., 582., 230., 94., & |
---|
237 | 87., 47., 35., 32., 29., 21./ |
---|
238 | ! 17-Magdalena (Colombia) 11N06 74W51 R2 |
---|
239 | ! Old=237 km3/year=7510m3/s; UNESCO(71-73+76-79, p348)=6973m3/s; |
---|
240 | ! ratio=1.08 |
---|
241 | ! --> 7451.58301m3/s |
---|
242 | DATA inb(17)/2/ |
---|
243 | DATA (iirnf(jc,17),jc=1,jpcoef)/ 424, 424, 18*0 / |
---|
244 | DATA (ijrnf(jc,17),jc=1,jpcoef)/ 270, 271, 18*0 / |
---|
245 | DATA (zrnfm(jm,17),jm=1,jpmois)/ & |
---|
246 | 6598., 4441., 4087., 5445., 7380., 8067., & |
---|
247 | 8509., 7428., 7631., 9296., 10588., 9949./ |
---|
248 | ! 18-Orinoco (Venezuela) 8N37 62W15! R1 |
---|
249 | ! Old=1100 km3/year=34857m3/s; UNESCO(73-75, p364)=28736m3/s; |
---|
250 | ! HLA=26371m3/s |
---|
251 | ! Data from HLA: 8.1N - 63.3W 1925/89 |
---|
252 | ! --> 26371.416 |
---|
253 | DATA inb(18)/5/ |
---|
254 | DATA (iirnf(jc,18),jc=1,jpcoef)/ 451, 451, 452, 452, 452, 15*0 / |
---|
255 | DATA (ijrnf(jc,18),jc=1,jpcoef)/ 271, 272, 271, 272, 272, 15*0 / |
---|
256 | DATA (zrnfm(jm,18),jm=1,jpmois)/ & |
---|
257 | 10670., 7513., 6344., 6887., 13086., 26699., & |
---|
258 | 43938., 59957., 57435., 39951., 27410., 16567./ |
---|
259 | ! 19-Amazon (Brazil) 0S05 50W00 R1 |
---|
260 | ! Old=6300 km3/year=199634m3/s; CLIPPER=166192m3/s~UNESCO; |
---|
261 | ! HLA=184806m3/s |
---|
262 | ! UNESCO data are at Obidos which is far inland and therefore does not |
---|
263 | ! fully represent the actual discharge |
---|
264 | ! Data from HLA: 1.5N - 55.3W 1970/96 |
---|
265 | ! --> 184845.484m3/s |
---|
266 | DATA inb(19)/14/ |
---|
267 | DATA (iirnf(jc,19),jc=1,jpcoef)/ 475, 476, 476, 476, 476, & |
---|
268 | 477, 477, 477, 477, 477, & |
---|
269 | 478, 478, 478, 478, 6*0 / |
---|
270 | DATA (ijrnf(jc,19),jc=1,jpcoef)/ 254, 251, 252, 253, 254, & |
---|
271 | 250, 251, 252, 253, 254, & |
---|
272 | 250, 251, 252, 253, 6*0 / |
---|
273 | DATA (zrnfm(jm,19),jm=1,jpmois)/ & |
---|
274 | 152233., 176415., 199941., 225565., 244016., 244861., & |
---|
275 | 230711., 203128., 164531., 128486., 118214., 130045. / |
---|
276 | ! 20-Sao Francisco (Brazil) 10S30 36W24 R1 |
---|
277 | ! Old=97 km3/year=3074m3/s; UNESCO(77-79, p338)=2645m3/s |
---|
278 | ! --> 2645.3335m3/s |
---|
279 | DATA inb(20)/1/ |
---|
280 | DATA (iirnf(jc,20),jc=1,jpcoef)/ 498, 19*0 / |
---|
281 | DATA (ijrnf(jc,20),jc=1,jpcoef)/ 223, 19*0 / |
---|
282 | DATA (zrnfm(jm,20),jm=1,jpmois)/ & |
---|
283 | 3244., 4116., 5139., 4330., 2323., 1800., & |
---|
284 | 1587., 1511., 1593., 1740., 2190., 2171. / |
---|
285 | ! 21-Parana+Uruguay(Arg./Uru.) 34S12 58W18 R1 |
---|
286 | ! Old=470 km3/year=14893m3/s; UNESCO(p323/324)=22263m3/s |
---|
287 | ! --> 22203.168m3/s |
---|
288 | DATA inb(21)/5/ |
---|
289 | DATA (iirnf(jc,21),jc=1,jpcoef)/ 462, 462, 462, 462, 463, 15*0 / |
---|
290 | DATA (ijrnf(jc,21),jc=1,jpcoef)/ 172, 173, 174, 175, 175, 15*0 / |
---|
291 | DATA (zrnfm(jm,21),jm=1,jpmois)/ & |
---|
292 | 24803., 26836., 24585., 21358., 20041., 21689., & |
---|
293 | 22647., 20403., 19627., 21400., 21126., 21923./ |
---|
294 | ! 22-Negro (Argentina) 41S02 62W47 R2 |
---|
295 | ! Old=30 km3/year=951m3/s; UNESCO(65-79,p330)=765m3/s; ratio=1.24 |
---|
296 | ! --> 907.666687m3/s |
---|
297 | DATA inb(22)/2/ |
---|
298 | DATA (iirnf(jc,22),jc=1,jpcoef)/ 446, 447, 18*0 / |
---|
299 | DATA (ijrnf(jc,22),jc=1,jpcoef)/ 160, 160, 18*0 / |
---|
300 | DATA (zrnfm(jm,22),jm=1,jpmois)/ & |
---|
301 | 882., 587., 519., 482., 623., 984., & |
---|
302 | 1120., 1294., 1092., 1166., 1105., 1038./ |
---|
303 | ! 23- Tocantins : ~48W00 ~0.5S R4 |
---|
304 | ! UNESCO(76-77+80-81, p335)=12167m3/s (station far inland!) |
---|
305 | ! --> 12173.082m3/s |
---|
306 | DATA inb(23)/4/ |
---|
307 | DATA (iirnf(jc,23),jc=1,jpcoef)/ 479, 480, 481, 482, 16*0 / |
---|
308 | DATA (ijrnf(jc,23),jc=1,jpcoef)/ 249, 249, 249, 249, 16*0 / |
---|
309 | DATA (zrnfm(jm,23),jm=1,jpmois)/ & |
---|
310 | 15507., 24855., 27854., 24257., 15860., 7861., & |
---|
311 | 4440., 3082., 2728., 3417., 5511., 10705. / |
---|
312 | |
---|
313 | DATA inb(24)/0/ |
---|
314 | DATA inb(25)/0/ |
---|
315 | |
---|
316 | ! Total run-offs(II)=256755.75 |
---|
317 | ! Total run-offs(I+II)=317730 |
---|
318 | |
---|
319 | ! III EUROPE |
---|
320 | ! ---------- |
---|
321 | ! 26-Rhone (France) 43N20 4E50 R1 |
---|
322 | ! Old=49 km3/year=1553m3/s; UNESCO(65-83, p512)=1229 |
---|
323 | ! --> 1217.16663m3/s |
---|
324 | DATA inb(26)/1/ |
---|
325 | DATA (iirnf(jc,26),jc=1,jpcoef)/ 584, 19*0 / |
---|
326 | DATA (ijrnf(jc,26),jc=1,jpcoef)/ 347, 19*0 / |
---|
327 | DATA (zrnfm(jm,26),jm=1,jpmois)/ & |
---|
328 | 1385., 1848., 1568., 1311., 1420., 1191., & |
---|
329 | 937., 796., 794., 962., 1161., 1233./ |
---|
330 | ! 27-Po (Italy) 11N10 1W09 R1 |
---|
331 | ! Old=46 km3/year=1458m3/s; UNESCO(65-84, p527)=1550m3/s |
---|
332 | ! --> 1550.08337m3/s |
---|
333 | DATA inb(27)/1/ |
---|
334 | DATA (iirnf(jc,27),jc=1,jpcoef)/ 599, 19*0 / |
---|
335 | DATA (ijrnf(jc,27),jc=1,jpcoef)/ 355, 19*0 / |
---|
336 | DATA (zrnfm(jm,27),jm=1,jpmois)/ & |
---|
337 | 1304., 1605., 1856., 1549., 1958., 1760., & |
---|
338 | 1040., 1007., 1615., 1907., 1772., 1228./ |
---|
339 | ! 28-Danube (Roumania) 45N20 29E40 R1 |
---|
340 | ! Old=206 km3/year=6527m3/s; UNESCO(65-84, p544)=7011m3/s |
---|
341 | ! --> 7011.8335m3/s |
---|
342 | DATA inb(28)/2/ |
---|
343 | DATA (iirnf(jc,28),jc=1,jpcoef)/ 632, 632, 18*0 / |
---|
344 | DATA (ijrnf(jc,28),jc=1,jpcoef)/ 356, 357, 18*0 / |
---|
345 | DATA (zrnfm(jm,28),jm=1,jpmois)/ & |
---|
346 | 6625., 7352., 7986., 9279., 9482., 8618., & |
---|
347 | 7586., 6001., 5210., 4989., 4921., 6093./ |
---|
348 | |
---|
349 | DATA inb(29)/0/ |
---|
350 | DATA inb(30)/0/ |
---|
351 | |
---|
352 | ! Total run-offs(III)=9779.08301m3/s |
---|
353 | ! Total run-offs(I+II+III)=327509.094m3/s |
---|
354 | |
---|
355 | ! IV EURASIAN ARCTIC |
---|
356 | ! ------------------ |
---|
357 | ! 31-Iana (Russia) R1 |
---|
358 | ! Old=29 km3/year=919m3/s; UNESCO(65-72+74+76-84, p463)=936m3/s |
---|
359 | ! --> 0m3/s |
---|
360 | DATA inb(31)/0/ |
---|
361 | DATA (iirnf(jc,31),jc=1,jpcoef)/ 20*0 / |
---|
362 | DATA (ijrnf(jc,31),jc=1,jpcoef)/ 20*0 / |
---|
363 | DATA (zrnfm(jm,31),jm=1,jpmois)/ & |
---|
364 | 2., 1., 0., 0., 917., 3514., & |
---|
365 | 2968., 2329., 1322., 178., 36., 8./ |
---|
366 | ! 32-Ob (Russia) 66N45 69E30 R1 |
---|
367 | ! Old=385 km3/year=12200m3/s; UNESCO(65-84, p478)=12862m3/s |
---|
368 | ! --> 12861.5m3/s |
---|
369 | DATA inb(32)/8/ |
---|
370 | DATA (iirnf(jc,32),jc=1,jpcoef)/ 635, 635, 636, 636, 636,637, 637,637,12*0 / |
---|
371 | DATA (ijrnf(jc,32),jc=1,jpcoef)/ 508, 509, 509, 508, 507,509, 508,507,12*0 / |
---|
372 | DATA (zrnfm(jm,32),jm=1,jpmois)/ & |
---|
373 | 4977., 4150., 3611., 3590., 14723., 33366., & |
---|
374 | 30773., 22785., 14701., 9705., 6078., 5879./ |
---|
375 | ! 33-Yenesei (Russia) 71N50 82E40 R1 |
---|
376 | ! Old=560 km3/year=17745m3/s; UNESCO(65-84, p472)=17462m3/s |
---|
377 | ! --> 17461.666m3/s |
---|
378 | DATA inb(33)/9/ |
---|
379 | DATA (iirnf(jc,33),jc=1,jpcoef)/ 114,115,115,116, 117,117,117,118,118,11*0 / |
---|
380 | DATA (ijrnf(jc,33),jc=1,jpcoef)/ 504,504,503,503, 501,502,503,502,501,11*0 / |
---|
381 | DATA (zrnfm(jm,33),jm=1,jpmois)/ & |
---|
382 | 6659., 7007., 6895., 7089., 24183., 74965., & |
---|
383 | 25325., 15940., 16050., 12678., 6741., 6008./ |
---|
384 | ! 34-Severnaya Dvina(Russia)64N32 40E30 R5 |
---|
385 | ! Old=106 km3/year=3359m3/s; not in UNESCO data base |
---|
386 | ! --> 3359.00024m3/s |
---|
387 | DATA inb(34)/2/ |
---|
388 | DATA (iirnf(jc,34),jc=1,jpcoef)/ 620, 620, 18*0 / |
---|
389 | DATA (ijrnf(jc,34),jc=1,jpcoef)/ 448, 449, 18*0 / |
---|
390 | DATA (zrnfm(jm,34),jm=1,jpmois)/12*3359./ |
---|
391 | ! 35-Lena (Russia) 72N25 126E40 R1 |
---|
392 | ! --> 16192.333m3/s |
---|
393 | DATA inb(35)/9/ |
---|
394 | DATA (iirnf(jc,35),jc=1,jpcoef)/ 142,143,148,148,149,149,150,150,150,11*0 / |
---|
395 | DATA (ijrnf(jc,35),jc=1,jpcoef)/ 458,458,456,448,449,448,451,450,449,11*0 / |
---|
396 | DATA (zrnfm(jm,35),jm=1,jpmois)/ & |
---|
397 | 2663., 2095., 1629., 1324., 7203., 64230., & |
---|
398 | 40810., 28310., 25995., 13299., 3673., 3077./ |
---|
399 | ! 36-Kolyma (Russia) 69N30 161E00 R1 |
---|
400 | ! Old=71 km3/year=2250m3/s; UNESCO(65-84)=2138m3/s |
---|
401 | ! --> 2138.66675m3/s |
---|
402 | DATA inb(36)/2/ |
---|
403 | DATA (iirnf(jc,36),jc=1,jpcoef)/ 184, 185, 18*0 / |
---|
404 | DATA (ijrnf(jc,36),jc=1,jpcoef)/ 429, 429, 18*0 / |
---|
405 | DATA (zrnfm(jm,36),jm=1,jpmois)/ & |
---|
406 | 131., 92., 76., 65., 2576., 9143., & |
---|
407 | 4670., 3727., 3226., 1119., 630., 209./ |
---|
408 | ! 37-Indigirka (russia) 70N48 148E54 R1 |
---|
409 | ! Old=55 km3/year=1743m3/s; UNESCO(65-68+70+72-84)=1622. |
---|
410 | ! --> 1622.00012m3/s |
---|
411 | DATA inb(37)/1/ |
---|
412 | DATA (iirnf(jc,37),jc=1,jpcoef)/ 172, 19*0 / |
---|
413 | DATA (ijrnf(jc,37),jc=1,jpcoef)/ 438, 19*0 / |
---|
414 | DATA (zrnfm(jm,37),jm=1,jpmois)/ & |
---|
415 | 32., 17., 11., 7., 354., 5629., & |
---|
416 | 5210., 4617., 2823., 559., 142., 63./ |
---|
417 | |
---|
418 | DATA inb(38)/0/ |
---|
419 | DATA inb(39)/0/ |
---|
420 | DATA inb(40)/0/ |
---|
421 | |
---|
422 | ! Total run-offs(IV)=40773.6641m3/s |
---|
423 | ! Total run-offs(I+II+III+IV)=368282.75m3/s |
---|
424 | |
---|
425 | ! V ASIA |
---|
426 | ! ------ |
---|
427 | ! 41-Amur (Russia) 52N56 141E10 R2 |
---|
428 | ! Old=325 km3/year=10299m3/s; UNESCO(65-84, p451)=8792m3/s; ratio=1.17 |
---|
429 | ! --> 10299.332m3/s |
---|
430 | DATA inb(41)/2/ |
---|
431 | DATA (iirnf(jc,41),jc=1,jpcoef)/ 141, 141, 18*0 / |
---|
432 | DATA (ijrnf(jc,41),jc=1,jpcoef)/ 381, 382, 18*0 / |
---|
433 | DATA (zrnfm(jm,41),jm=1,jpmois)/ & |
---|
434 | 2108., 1543., 1252., 4088., 15464., 17336., & |
---|
435 | 16100., 20801., 21548., 15658., 5188., 2506./ |
---|
436 | ! 42-Liaohe (China) 40N40 122E09 R5 |
---|
437 | ! Old=6 km3/year=190m3/s; not found in UNESCO data base |
---|
438 | ! --> 189.999985m3/s |
---|
439 | DATA inb(42)/1/ |
---|
440 | DATA (iirnf(jc,42),jc=1,jpcoef)/ 100, 19*0 / |
---|
441 | DATA (ijrnf(jc,42),jc=1,jpcoef)/ 342, 19*0 / |
---|
442 | DATA (zrnfm(jm,42),jm=1,jpmois)/12*190./ |
---|
443 | ! 43-Daling (China) 40N56 121E43 R5 |
---|
444 | ! Old=1 km3/year=32m3/s; not found in UNESCO data base |
---|
445 | ! --> 31.9999962m3/s |
---|
446 | DATA inb(43)/2/ |
---|
447 | DATA (iirnf(jc,43),jc=1,jpcoef)/ 99, 100, 18*0 / |
---|
448 | DATA (ijrnf(jc,43),jc=1,jpcoef)/ 342, 342, 18*0 / |
---|
449 | DATA (zrnfm(jm,43),jm=1,jpmois)/12*32./ |
---|
450 | ! 44-Haiho (China) 6N17 1E23 |
---|
451 | ! Old=2 km3/year=64m3/s; not found in UNESCO data base |
---|
452 | ! --> 63.9999924m3/s |
---|
453 | DATA inb(44)/1/ |
---|
454 | DATA (iirnf(jc,44),jc=1,jpcoef)/ 121, 19*0 / |
---|
455 | DATA (ijrnf(jc,44),jc=1,jpcoef)/ 348, 19*0 / |
---|
456 | DATA (zrnfm(jm,44),jm=1,jpmois)/12*64./ |
---|
457 | ! 45-Yellow/Huangho (China) 37N32 118E19 R2 |
---|
458 | ! Old=49 km3/year=1553m3/s; UNESCO(p371)=1210m3/s; ratio=1.28 |
---|
459 | ! --> 1552.5m3/s |
---|
460 | DATA inb(45)/3/ |
---|
461 | DATA (iirnf(jc,45),jc=1,jpcoef)/ 92, 92, 92, 17*0 / |
---|
462 | DATA (ijrnf(jc,45),jc=1,jpcoef)/ 336, 337, 338, 17*0 / |
---|
463 | DATA (zrnfm(jm,45),jm=1,jpmois)/ & |
---|
464 | 697., 674., 1282., 1206., 1109., 778., & |
---|
465 | 1964., 3270., 3668., 2133., 978., 871./ |
---|
466 | ! 46-Yangtze (China) 31N48 121E10=Changjiang R2 |
---|
467 | ! Old=900 km3/year=28520m3/s; UNESCO(p372)=25032m3/s; ratio=1.14 |
---|
468 | ! --> 28519.584m3/s |
---|
469 | DATA inb(46)/2/ |
---|
470 | DATA (iirnf(jc,46),jc=1,jpcoef)/ 100, 100, 18*0 / |
---|
471 | DATA (ijrnf(jc,46),jc=1,jpcoef)/ 316, 317, 18*0 / |
---|
472 | DATA (zrnfm(jm,46),jm=1,jpmois)/ & |
---|
473 | 10994., 10548., 13254., 24097., 36972., 47197., & |
---|
474 | 53891., 42070., 37627., 30762., 22046., 12777./ |
---|
475 | ! 47-Pearl/Zhujiang (China) 23N36 106E34=Xijiang R2 |
---|
476 | ! Old=302 km3/year=9570m3/s; UNESCO(76,82, p373)=6987m3/s; ratio=1.37 |
---|
477 | ! --> 9570.75m3/s |
---|
478 | DATA inb(47)/2/ |
---|
479 | DATA (iirnf(jc,47),jc=1,jpcoef)/ 72, 73, 18*0 / |
---|
480 | DATA (ijrnf(jc,47),jc=1,jpcoef)/ 294, 294, 18*0 / |
---|
481 | DATA (zrnfm(jm,47),jm=1,jpmois)/ & |
---|
482 | 2357., 2415., 2634., 6446., 14246., 18395., & |
---|
483 | 19213., 17558., 12312., 7255., 5864., 6154./ |
---|
484 | ! 48-Red/Hungho (Vietnam) 20N17 106E34 R5 |
---|
485 | ! Old=123 km3/year=3897m3/s; not in UNESCO data base |
---|
486 | ! --> 3897.m3/s |
---|
487 | DATA inb(48)/3/ |
---|
488 | DATA (iirnf(jc,48),jc=1,jpcoef)/ 68, 68, 68, 17*0 / |
---|
489 | DATA (ijrnf(jc,48),jc=1,jpcoef)/ 288, 289, 290, 17*0 / |
---|
490 | DATA (zrnfm(jm,48),jm=1,jpmois)/12*3897./ |
---|
491 | ! 49-Mekong (Vietnam) 10N33 105E24 R5 |
---|
492 | ! Old=470 km3/year=14893m3/s; not in UNESCO data base |
---|
493 | ! --> 14893.0029m3/s |
---|
494 | DATA inb(49)/3/ |
---|
495 | DATA (iirnf(jc,49),jc=1,jpcoef)/ 69, 69, 70, 17*0 / |
---|
496 | DATA (ijrnf(jc,49),jc=1,jpcoef)/ 269, 270, 270, 17*0 / |
---|
497 | DATA (zrnfm(jm,49),jm=1,jpmois)/12*14893./ |
---|
498 | ! 50-Irrawady (Burma) 15N50 95E06 R5 |
---|
499 | ! Old=428 km3/year=13563m3/s; not in UNESCO data base |
---|
500 | ! --> 13563.m3/s |
---|
501 | DATA inb(50)/2/ |
---|
502 | DATA (iirnf(jc,50),jc=1,jpcoef)/ 44, 45, 18*0 / |
---|
503 | DATA (ijrnf(jc,50),jc=1,jpcoef)/ 283, 282, 18*0 / |
---|
504 | DATA (zrnfm(jm,50),jm=1,jpmois)/12*13563./ |
---|
505 | ! 51-Ganges+Brahmaputra (Beng.) 22N00 91E00 R1 |
---|
506 | ! Old=971 km3/year=30769m3/s; UNESCO(69-70+73-75, p367)=31760m3/s |
---|
507 | ! --> 32147.498m3/s |
---|
508 | DATA inb(51)/4/ |
---|
509 | DATA (iirnf(jc,51),jc=1,jpcoef)/ 36, 37, 38, 39, 16*0 / |
---|
510 | DATA (ijrnf(jc,51),jc=1,jpcoef)/ 295, 295, 295, 295, 16*0 / |
---|
511 | DATA (zrnfm(jm,51),jm=1,jpmois)/ & |
---|
512 | 6623., 6315., 6432., 9410., 17263., 38302., & |
---|
513 | 64688., 80338., 84802., 43387., 17888., 10322./ |
---|
514 | ! 52-Mehandi (India) R5 |
---|
515 | ! Old=67 km3/year=2123m3/s; not in UNESCO data base |
---|
516 | ! --> 0m3/s |
---|
517 | DATA inb(52)/0/ |
---|
518 | DATA (iirnf(jc,52),jc=1,jpcoef)/ 20*0 / |
---|
519 | DATA (ijrnf(jc,52),jc=1,jpcoef)/ 20*0 / |
---|
520 | DATA (zrnfm(jm,52),jm=1,jpmois)/12*2123./ |
---|
521 | ! 53-Damodar (India) R2 |
---|
522 | ! Old=10 km3/year=320m3/s; UNESCO(p386)=173m3/s; ratio=1.85 |
---|
523 | ! --> 432.916656m3/s |
---|
524 | DATA inb(53)/1/ |
---|
525 | DATA (iirnf(jc,53),jc=1,jpcoef)/ 30, 19*0 / |
---|
526 | DATA (ijrnf(jc,53),jc=1,jpcoef)/ 294, 19*0 / |
---|
527 | DATA (zrnfm(jm,53),jm=1,jpmois)/ & |
---|
528 | 51., 37., 42., 49., 85., 296., & |
---|
529 | 896., 1390., 1591., 542., 172., 44./ |
---|
530 | ! 54-Godavari (India) 17N00 81E45 R2 |
---|
531 | ! Old=84 km3/year=2662m3/s; UNESCO(p379)=1916m3/s; ratio=1.39 |
---|
532 | ! --> 3670.5m3/s |
---|
533 | DATA inb(54)/1/ |
---|
534 | DATA (iirnf(jc,54),jc=1,jpcoef)/ 19, 19*0 / |
---|
535 | DATA (ijrnf(jc,54),jc=1,jpcoef)/ 283, 19*0 / |
---|
536 | DATA (zrnfm(jm,54),jm=1,jpmois)/ & |
---|
537 | 988., 737., 280., 219., 196., 1691., & |
---|
538 | 8063., 14571., 11428., 3953., 1221., 699./ |
---|
539 | ! 55-Indus (Pakistan) 24N20 67E47 R2 |
---|
540 | ! Old=238 km3/year=7542m3/s; UNESCO (76-79, p428)=2396m3/s; ratio=3.15 |
---|
541 | ! --> 7562.m3/s |
---|
542 | DATA inb(55)/4/ |
---|
543 | DATA (iirnf(jc,55),jc=1,jpcoef)/ 711, 711, 712, 713, 16*0 / |
---|
544 | DATA (ijrnf(jc,55),jc=1,jpcoef)/ 297, 298, 297, 297, 16*0 / |
---|
545 | DATA (zrnfm(jm,55),jm=1,jpmois)/ & |
---|
546 | 988., 737., 1904., 1968., 3625., 6143., & |
---|
547 | 15969., 36656., 17173., 3215., 1251., 1115./ |
---|
548 | ! 56-Tigris and Euphrates (Iraq) 31N00 47E25 R1 |
---|
549 | ! Old=46 km3/year=1458m3/s; UNESCO(76-79, p428)=2396m3/s |
---|
550 | ! --> 2248.8335m3/s |
---|
551 | DATA inb(56)/3/ |
---|
552 | DATA (iirnf(jc,56),jc=1,jpcoef)/ 673, 673, 674, 17*0 / |
---|
553 | DATA (ijrnf(jc,56),jc=1,jpcoef)/ 312, 313, 313, 17*0 / |
---|
554 | DATA (zrnfm(jm,56),jm=1,jpmois)/ & |
---|
555 | 1872., 2127., 2962., 4944., 5036., 3078., & |
---|
556 | 1362., 900., 786., 993., 1184., 1742./ |
---|
557 | |
---|
558 | DATA inb(57)/0/ |
---|
559 | DATA inb(58)/0/ |
---|
560 | DATA inb(59)/0/ |
---|
561 | DATA inb(60)/0/ |
---|
562 | |
---|
563 | ! Total run-offs (V)=128642.922m3/s |
---|
564 | ! Total run-offs (I+II+III+IV+V)=496925.688m3/s |
---|
565 | |
---|
566 | ! VI AFRICA |
---|
567 | ! ---------- |
---|
568 | ! 61-Nile (Egypt) 30N10 31E06 R1 |
---|
569 | ! Old=30 km3/year=951m3/s; UNESCO(73-84)=1251m3/s |
---|
570 | ! --> 1251.33337m3/s |
---|
571 | DATA inb(61)/3/ |
---|
572 | DATA (iirnf(jc,61),jc=1,jpcoef)/ 636, 636, 637, 17*0 / |
---|
573 | DATA (ijrnf(jc,61),jc=1,jpcoef)/ 317, 318, 318, 17*0 / |
---|
574 | DATA (zrnfm(jm,61),jm=1,jpmois)/ & |
---|
575 | 1239., 1035., 1126., 1110., 1170., 1542., & |
---|
576 | 1742., 1559., 1191., 1103., 1075., 1124./ |
---|
577 | ! 62-Niger (Nigeria) 5N33 6E33 R1 |
---|
578 | ! Old=192 km3/year=6084m3/s; HLA=5545m3/s; |
---|
579 | ! Data from HLA: 4.5N - 6 E 1950/87 |
---|
580 | ! --> 5545.08301m3/s |
---|
581 | DATA inb(62)/3/ |
---|
582 | DATA (iirnf(jc,62),jc=1,jpcoef)/ 591, 592, 593, 17*0 / |
---|
583 | DATA (ijrnf(jc,62),jc=1,jpcoef)/ 258, 258, 258, 17*0 / |
---|
584 | DATA (zrnfm(jm,62),jm=1,jpmois)/ & |
---|
585 | 2029., 1804., 1573., 1418., 1494., 2561., & |
---|
586 | 5272., 9629., 14951., 16227., 6842., 2741./ |
---|
587 | ! 63-Zaire/Congo (Congo) 6S04 12E24 R1 |
---|
588 | ! Old=1250 km3/year=39610m3/s; UNESCO(65-83, p149)=41941m3/s; HLA |
---|
589 | ! =39720m3/s |
---|
590 | ! Data from HLA: 6S - 12 E 1902/96 |
---|
591 | ! --> 39719.5m3/s |
---|
592 | DATA inb(63)/4/ |
---|
593 | DATA (iirnf(jc,63),jc=1,jpcoef)/ 599, 599, 600, 600, 16*0 / |
---|
594 | DATA (ijrnf(jc,63),jc=1,jpcoef)/ 237, 238, 236, 237, 16*0 / |
---|
595 | DATA (zrnfm(jm,63),jm=1,jpmois)/ & |
---|
596 | 47200., 37692., 34374., 36603., 37862., 35588., & |
---|
597 | 30966., 30774., 35981., 42825., 51455., 55314./ |
---|
598 | ! 64-Orange (Sth Africa) 28S41 16E28 R5 |
---|
599 | ! Old=11 km3/year=349m3/s |
---|
600 | ! --> 349.000031m3/s |
---|
601 | DATA inb(64)/2/ |
---|
602 | DATA (iirnf(jc,64),jc=1,jpcoef)/ 609, 609, 18*0 / |
---|
603 | DATA (ijrnf(jc,64),jc=1,jpcoef)/ 188, 189, 18*0 / |
---|
604 | DATA (zrnfm(jm,64),jm=1,jpmois)/12*349./ |
---|
605 | ! 65-Zambesi (Mozambique) 18S55 36E04 R5 |
---|
606 | ! Old=223 km3/year=7066m3/s |
---|
607 | ! --> 7066.00049m3/s |
---|
608 | DATA inb(65)/3/ |
---|
609 | DATA (iirnf(jc,65),jc=1,jpcoef)/ 646, 646, 646, 17*0 / |
---|
610 | DATA (ijrnf(jc,65),jc=1,jpcoef)/ 207, 208, 209, 17*0 / |
---|
611 | DATA (zrnfm(jm,65),jm=1,jpmois)/12*7066./ |
---|
612 | ! 66-Limpopo (mozambique) 25S15 33E30 R5 |
---|
613 | ! Old=5 km3/year=158m3/s |
---|
614 | ! --> 158.m3/s |
---|
615 | DATA inb(66)/2/ |
---|
616 | DATA (iirnf(jc,66),jc=1,jpcoef)/ 642, 642, 18*0 / |
---|
617 | DATA (ijrnf(jc,66),jc=1,jpcoef)/ 196, 197, 18*0 / |
---|
618 | DATA (zrnfm(jm,66),jm=1,jpmois)/12*158./ |
---|
619 | ! 67-Rufiji (tanzania) 8S00 39E20 R5 |
---|
620 | ! Old=9 km3/year=285m3/s |
---|
621 | ! --> 285.m3/s |
---|
622 | DATA inb(67)/1/ |
---|
623 | DATA (iirnf(jc,67),jc=1,jpcoef)/ 655, 19*0 / |
---|
624 | DATA (ijrnf(jc,67),jc=1,jpcoef)/ 235, 19*0 / |
---|
625 | DATA (zrnfm(jm,67),jm=1,jpmois)/12*285./ |
---|
626 | ! Other rivers from the Gulf of Guinea |
---|
627 | ! 68-Kouilou (Congo) 4.5S-11.5E (UNESCO=p51) R1 |
---|
628 | ! Data from HLA: 4S - 11 E 1953/86 = 919m3/s |
---|
629 | ! --> 924.25m3/s |
---|
630 | DATA inb(68)/1/ |
---|
631 | DATA (iirnf(jc,68),jc=1,jpcoef)/ 599, 19*0 / |
---|
632 | DATA (ijrnf(jc,68),jc=1,jpcoef)/ 241, 19*0 / |
---|
633 | DATA (zrnfm(jm,68),jm=1,jpmois)/ & |
---|
634 | 1115., 1166., 1285., 1462., 1454., 725., & |
---|
635 | 480., 378., 327., 400., 1018., 1281./ |
---|
636 | ! 69-Nyanga (Gabon) 3S-10E (UNESCO=p66) R1 |
---|
637 | ! Data from HLA: 2.8N - 10 E 1952/82 = 264m3/s |
---|
638 | ! --> 264.166687m3/s |
---|
639 | DATA inb(69)/1/ |
---|
640 | DATA (iirnf(jc,69),jc=1,jpcoef)/ 595, 19*0 / |
---|
641 | DATA (ijrnf(jc,69),jc=1,jpcoef)/ 244, 19*0 / |
---|
642 | DATA (zrnfm(jm,69),jm=1,jpmois)/ & |
---|
643 | 126., 99., 143., 253., 316., 300., & |
---|
644 | 156., 93., 231., 557., 592., 304./ |
---|
645 | ! 70-Ogoue (Gabon) 1S-9E (UNESCO, p69) R1 |
---|
646 | ! Data from HLA: .75S - 9 E 1930/90 = 4512m3/s |
---|
647 | ! --> 4511.83301m3/s |
---|
648 | DATA inb(70)/2/ |
---|
649 | DATA (iirnf(jc,70),jc=1,jpcoef)/ 593, 593, 18*0 / |
---|
650 | DATA (ijrnf(jc,70),jc=1,jpcoef)/ 247, 248, 18*0 / |
---|
651 | DATA (zrnfm(jm,70),jm=1,jpmois)/ & |
---|
652 | 4514., 3988., 4622., 5403., 5790., 4145., & |
---|
653 | 2616., 1961., 2119., 4304., 7870., 6810./ |
---|
654 | ! 71-Sanaga (Cameroun) 3.5N-9.5E (UNESCO, p31) R1 |
---|
655 | ! Data from HLA: 3.9N - 9.6E 1944/80 = 1983m3/s |
---|
656 | ! --> 1982.58337m3/s |
---|
657 | DATA inb(71)/2/ |
---|
658 | DATA (iirnf(jc,71),jc=1,jpcoef)/ 595, 595, 18*0 / |
---|
659 | DATA (ijrnf(jc,71),jc=1,jpcoef)/ 255, 256, 18*0 / |
---|
660 | DATA (zrnfm(jm,71),jm=1,jpmois)/ & |
---|
661 | 813., 578., 517., 633., 909., 1278., & |
---|
662 | 1928., 2653., 4361., 5470., 3307., 1344./ |
---|
663 | ! 72-Nyong (Cameroun) 3.1N-9.5E (UNESCO, p36) R1 |
---|
664 | ! Data from HLA: 3.5N - 9.8E 1951/84 = 432m3/s |
---|
665 | ! --> 431.166656m3/s |
---|
666 | DATA inb(72)/1/ |
---|
667 | DATA (iirnf(jc,72),jc=1,jpcoef)/ 594, 19*0 / |
---|
668 | DATA (ijrnf(jc,72),jc=1,jpcoef)/ 257, 19*0 / |
---|
669 | DATA (zrnfm(jm,72),jm=1,jpmois)/ & |
---|
670 | 232., 128., 150., 266., 402., 459., & |
---|
671 | 353., 306., 604., 930., 878., 466./ |
---|
672 | ! 73-Ntem (Cameroun) 2.5N-9.5E (UNESCO, p36) R1 |
---|
673 | ! Data from HLA: 3N - 9.8E 1951/87 = 330m3/s |
---|
674 | ! --> 329.75m3/s |
---|
675 | DATA inb(73)/1/ |
---|
676 | DATA (iirnf(jc,73),jc=1,jpcoef)/ 594, 19*0 / |
---|
677 | DATA (ijrnf(jc,73),jc=1,jpcoef)/ 254, 19*0 / |
---|
678 | DATA (zrnfm(jm,73),jm=1,jpmois)/ & |
---|
679 | 452., 455., 492., 466., 482., 251., & |
---|
680 | 151., 108., 87., 133., 395., 485./ |
---|
681 | ! 74-Wouri (Cameroun) 4N-9.5E (UNESCO, p34) R1 |
---|
682 | ! Data from HLA: 4N - 9.5E 1951/84 = 307m3/s |
---|
683 | ! --> 306.583374m3/s |
---|
684 | DATA inb(74)/1/ |
---|
685 | DATA (iirnf(jc,74),jc=1,jpcoef)/ 594, 19*0 / |
---|
686 | DATA (ijrnf(jc,74),jc=1,jpcoef)/ 258, 19*0 / |
---|
687 | DATA (zrnfm(jm,74),jm=1,jpmois)/ & |
---|
688 | 88., 63., 68., 95., 139., 222., & |
---|
689 | 395., 633., 825., 675., 328., 148./ |
---|
690 | ! 75-Oueme (Benin) 6.5N-2.5E (UNESCO, p27) R1 |
---|
691 | ! Data from HLA: 6.2N - 2.5E 1951/84 =113m3/s |
---|
692 | ! --> 112.666664m3/s |
---|
693 | DATA inb(75)/1/ |
---|
694 | DATA (iirnf(jc,75),jc=1,jpcoef)/ 583, 19*0 / |
---|
695 | DATA (ijrnf(jc,75),jc=1,jpcoef)/ 262, 19*0 / |
---|
696 | DATA (zrnfm(jm,75),jm=1,jpmois)/ & |
---|
697 | 0., 0., 0., 0., 1., 18., & |
---|
698 | 96., 324., 531., 319., 58., 5./ |
---|
699 | ! 76-Mono(Tepetou,Benin/Togo) 6.25N-1.8E (U, p27) R1 |
---|
700 | ! Data from HLA: 6.2N - 1 E 1951/85 = 96m3/s |
---|
701 | ! --> 96.4999924m3/s |
---|
702 | DATA inb(76)/2/ |
---|
703 | DATA (iirnf(jc,76),jc=1,jpcoef)/ 578, 579, 18*0 / |
---|
704 | DATA (ijrnf(jc,76),jc=1,jpcoef)/ 262, 262, 18*0 / |
---|
705 | DATA (zrnfm(jm,76),jm=1,jpmois)/ & |
---|
706 | 2., 1., 1., 4., 8., 44., & |
---|
707 | 168., 274., 386., 215., 47., 8./ |
---|
708 | ! 77-Comoe (Ivory Coast) 5N-4W (UNESCO, p56) R1 |
---|
709 | ! Data from HLA: 5.1N - 4 W 1955/86 = 194m3/s |
---|
710 | ! --> 194.25m3/s |
---|
711 | DATA inb(77)/1/ |
---|
712 | DATA (iirnf(jc,77),jc=1,jpcoef)/ 569, 19*0 / |
---|
713 | DATA (ijrnf(jc,77),jc=1,jpcoef)/ 260, 19*0 / |
---|
714 | DATA (zrnfm(jm,77),jm=1,jpmois)/ & |
---|
715 | 14., 7., 9., 17., 29., 90., & |
---|
716 | 145., 349., 740., 688., 198., 45./ |
---|
717 | ! 78-Sassandra (Ivory Coast) 5N-6W (UNESCO, p58) R1 |
---|
718 | ! Data from HLA: 5N - 6 W 1953/86 = 449m3/s |
---|
719 | ! --> 449.333344m3/s |
---|
720 | DATA inb(78)/1/ |
---|
721 | DATA (iirnf(jc,78),jc=1,jpcoef)/ 564, 19*0 / |
---|
722 | DATA (ijrnf(jc,78),jc=1,jpcoef)/ 260, 19*0 / |
---|
723 | DATA (zrnfm(jm,78),jm=1,jpmois)/ & |
---|
724 | 144., 116., 119., 143., 156., 334., & |
---|
725 | 419., 692., 1277., 1190., 553., 249./ |
---|
726 | |
---|
727 | DATA inb(79)/0/ |
---|
728 | DATA inb(80)/0/ |
---|
729 | DATA inb(82)/0/ |
---|
730 | DATA inb(83)/0/ |
---|
731 | DATA inb(84)/0/ |
---|
732 | DATA inb(85)/0/ |
---|
733 | DATA inb(86)/0/ |
---|
734 | DATA inb(87)/0/ |
---|
735 | DATA inb(88)/0/ |
---|
736 | DATA inb(89)/0/ |
---|
737 | DATA inb(90)/0/ |
---|
738 | |
---|
739 | ! Total run-offs (VI)=63977.m3/s |
---|
740 | ! Total run-offs (I+II+III+IV+V+VI)=560902.688m3/s |
---|
741 | |
---|
742 | ! VII OCEANIA |
---|
743 | ! ------------ |
---|
744 | ! 91-Murray (Austrlia) 35S34 138E54 R5 |
---|
745 | ! Old=22 km3/year=697m3/s |
---|
746 | ! --> 696.999939m3/s |
---|
747 | DATA inb(91)/3/ |
---|
748 | DATA (iirnf(jc,91),jc=1,jpcoef)/ 129, 130, 130, 17*0 / |
---|
749 | DATA (ijrnf(jc,91),jc=1,jpcoef)/ 177, 176, 177, 17*0 / |
---|
750 | DATA (zrnfm(jm,91),jm=1,jpmois)/12*697./ |
---|
751 | ! 92-Haast (N. Zeland) 43S50 169E02 R5 |
---|
752 | ! Old=6 km3/year=190m3/s |
---|
753 | ! --> 189.999985m3/s |
---|
754 | DATA inb(92)/1/ |
---|
755 | DATA (iirnf(jc,92),jc=1,jpcoef)/ 194, 19*0 / |
---|
756 | DATA (ijrnf(jc,92),jc=1,jpcoef)/ 153, 19*0 / |
---|
757 | DATA (zrnfm(jm,92),jm=1,jpmois)/12*190./ |
---|
758 | ! 93-Fly (N. Guinea) 8S30 143E05 R5 |
---|
759 | ! Old=77 km3/year=2440m3/s |
---|
760 | ! --> 2440.m3/s |
---|
761 | DATA inb(93)/1/ |
---|
762 | DATA (iirnf(jc,93),jc=1,jpcoef)/ 143, 19*0 / |
---|
763 | DATA (ijrnf(jc,93),jc=1,jpcoef)/ 233, 19*0 / |
---|
764 | DATA (zrnfm(jm,93),jm=1,jpmois)/12*2440./ |
---|
765 | ! 94-Purari (N. Guinea) 7S25 145E05 R4 |
---|
766 | ! Old=77 km3/year=2440m3/s; UNESCO(77+83,p595)=2727m3/s |
---|
767 | ! --> 2727.33325m3/s |
---|
768 | DATA inb(94)/1/ |
---|
769 | DATA (iirnf(jc,94),jc=1,jpcoef)/ 146, 19*0 / |
---|
770 | DATA (ijrnf(jc,94),jc=1,jpcoef)/ 234, 19*0 / |
---|
771 | DATA (zrnfm(jm,94),jm=1,jpmois)/ & |
---|
772 | 1972., 3067., 2332., 2981., 2693., 2876., & |
---|
773 | 2867., 3477., 2195., 2967., 2674., 2627./ |
---|
774 | ! 95-Choshui ch''i (taiwan) 24N03 120E23 R2 |
---|
775 | ! Old=6 km3/year=190m3/s; UNESCO(65-68, p374)=120m3/s; ratio=1.58 |
---|
776 | ! --> 0.m3/s |
---|
777 | DATA inb(95)/0/ |
---|
778 | DATA (iirnf(jc,95),jc=1,jpcoef)/ 20*0 / |
---|
779 | DATA (ijrnf(jc,95),jc=1,jpcoef)/ 20*0 / |
---|
780 | DATA (zrnfm(jm,95),jm=1,jpmois)/ & |
---|
781 | 65., 63., 83., 100., 177., 535., & |
---|
782 | 332., 313., 236., 193., 101., 84./ |
---|
783 | ! 96-Kaoping (Taiwan) R5 |
---|
784 | ! Old=9 km3/year=285m3/s |
---|
785 | ! --> 0.m3/s |
---|
786 | DATA inb(96)/0/ |
---|
787 | DATA (iirnf(jc,96),jc=1,jpcoef)/ 20*0 / |
---|
788 | DATA (ijrnf(jc,96),jc=1,jpcoef)/ 20*0 / |
---|
789 | DATA (zrnfm(jm,96),jm=1,jpmois)/12*285./ |
---|
790 | ! 97-Tsengwen (Taiwan) R5 |
---|
791 | ! Old=2 km3/year=63m3/s |
---|
792 | ! --> 0.m3/s |
---|
793 | DATA inb(97)/0/ |
---|
794 | DATA (iirnf(jc,97),jc=1,jpcoef)/ 20*0 / |
---|
795 | DATA (ijrnf(jc,97),jc=1,jpcoef)/ 20*0 / |
---|
796 | DATA (zrnfm(jm,97),jm=1,jpmois)/12*63./ |
---|
797 | ! 98-Hualien (Taiwan) R5 |
---|
798 | ! Old=4 km3/year=127m3/s |
---|
799 | ! --> 0.m3/s |
---|
800 | DATA inb(98)/0/ |
---|
801 | DATA (iirnf(jc,98),jc=1,jpcoef)/ 20*0 / |
---|
802 | DATA (ijrnf(jc,98),jc=1,jpcoef)/ 20*0 / |
---|
803 | DATA (zrnfm(jm,98),jm=1,jpmois)/12*127./ |
---|
804 | ! 99-Peinan (Taiwan) 22N46 121E10 R5 |
---|
805 | ! Old=4 km3/year=127m3/s |
---|
806 | ! --> 127.000008m3/s |
---|
807 | DATA inb(99)/1/ |
---|
808 | DATA (iirnf(jc,99),jc=1,jpcoef)/ 98, 19*0 / |
---|
809 | DATA (ijrnf(jc,99),jc=1,jpcoef)/ 287, 19*0 / |
---|
810 | DATA (zrnfm(jm,99),jm=1,jpmois)/12*127./ |
---|
811 | ! 100-Hsiukuluan (Taiwan) R5 |
---|
812 | ! Old=4 km3/year=127m3/s |
---|
813 | ! --> 0.m3/s |
---|
814 | DATA inb(100)/0/ |
---|
815 | DATA (iirnf(jc,100),jc=1,jpcoef)/ 20*0 / |
---|
816 | DATA (ijrnf(jc,100),jc=1,jpcoef)/ 20*0 / |
---|
817 | DATA (zrnfm(jm,100),jm=1,jpmois)/12*127./ |
---|
818 | |
---|
819 | DATA inb(101)/0/ |
---|
820 | DATA inb(102)/0/ |
---|
821 | |
---|
822 | ! Total run-offs(VII)=6181.33301m3/s |
---|
823 | ! Total run-offs(I+II+III+IV+V+VI+VII)= 567084. m3/s |
---|
824 | ! = 1.588668056E-6 kg/m2/s |
---|
825 | ! = 17883.5586 km3/Year |
---|
826 | ! = 0.567083955 Sverdrup |
---|
827 | |
---|
828 | DATA (inb(jr),jr=103,jpriv)/18*0/ |
---|
829 | |
---|
830 | !!---------------------------------------------------------------------- |
---|
831 | !! OPA 8.5, LODYC-IPSL (2002) |
---|
832 | !!---------------------------------------------------------------------- |
---|
833 | |
---|
834 | IF( kt == nit000 ) THEN |
---|
835 | IF(lwp) WRITE(numout,*) |
---|
836 | IF(lwp) WRITE(numout,*) 'flxrnf : runoff + orca05' |
---|
837 | IF(lwp) WRITE(numout,*) '~~~~~~' |
---|
838 | IF(lwp) WRITE(numout,*) ' nrunoff = ', nrunoff |
---|
839 | IF(lwp) WRITE(numout,*) |
---|
840 | ENDIF |
---|
841 | |
---|
842 | ! 0. initialization and indexes modification for zoom version |
---|
843 | ! -------------------------------------------------------------- |
---|
844 | |
---|
845 | |
---|
846 | ! conversion factor m3/s --> kg/m2/s *e1t*e2t |
---|
847 | |
---|
848 | zconv = rauw |
---|
849 | |
---|
850 | IF( kt == nit000 ) THEN |
---|
851 | ! Set runoffs and upstream coeff to zero |
---|
852 | runoff (:,:) = 0.e0 |
---|
853 | upsrnfh(:,:) = 0.e0 |
---|
854 | upsrnfz(:) = 0.e0 |
---|
855 | upsadv (:,:) = 0.e0 |
---|
856 | zcoefr(:,:,:) = 0.e0 |
---|
857 | |
---|
858 | IF( nrunoff == 0 ) THEN |
---|
859 | IF(lwp) WRITE(numout,*) |
---|
860 | IF(lwp) WRITE(numout,*) ' NO RUNOFF in this simulation ' |
---|
861 | IF(lwp) WRITE(numout,*) |
---|
862 | ENDIF |
---|
863 | ENDIF |
---|
864 | |
---|
865 | ! 1. initialization of the runoff array |
---|
866 | ! ------------------------------------- |
---|
867 | |
---|
868 | IF( nrunoff == 1 .OR. nrunoff == 2) THEN |
---|
869 | |
---|
870 | iman = jpmois |
---|
871 | |
---|
872 | i15 = nday / 16 |
---|
873 | |
---|
874 | imois = nmonth + i15 - 1 |
---|
875 | IF( imois == 0) imois = iman |
---|
876 | IF( nleapy == 1 .AND. MOD( nyear, 4 ) == 0 ) THEN |
---|
877 | ! Number of Day Between Data |
---|
878 | idbd = nbiss(imois) |
---|
879 | ELSEIF( nleapy > 1 ) THEN |
---|
880 | idbd = nleapy |
---|
881 | ELSE |
---|
882 | idbd = nobis(imois) |
---|
883 | END IF |
---|
884 | ! Number of days between imois, 15 and the end of month |
---|
885 | idmeom = idbd - 15 |
---|
886 | |
---|
887 | ! 1.1 Calendar computation |
---|
888 | |
---|
889 | ! nrnf1 number of the first array record used in the simulation |
---|
890 | ! nrnf2 number of the last array record |
---|
891 | |
---|
892 | nrnf1 = imois |
---|
893 | nrnf2 = nrnf1 + 1 |
---|
894 | nrnf1 = MOD( nrnf1, iman ) |
---|
895 | IF( nrnf1 == 0 ) nrnf1 = iman |
---|
896 | nrnf2 = MOD( nrnf2, iman ) |
---|
897 | IF( nrnf2 == 0 ) nrnf2 = iman |
---|
898 | |
---|
899 | IF( kt == nit000 ) THEN |
---|
900 | IF(lwp) WRITE(numout,*) |
---|
901 | IF(lwp) WRITE(numout,*) 'flxrnf routine' |
---|
902 | IF(lwp) WRITE(numout,*) 'first array record used nrnf1 ',nrnf1 |
---|
903 | IF(lwp) WRITE(numout,*) 'last array record used nrnf2 ',nrnf2 |
---|
904 | ENDIF |
---|
905 | |
---|
906 | zxy = FLOAT( nday + idmeom - idbd * i15 ) / idbd |
---|
907 | |
---|
908 | DO jr = 1, jpriv-1 |
---|
909 | inbp = inb(jr) |
---|
910 | IF( inbp /= 0 ) THEN |
---|
911 | DO jl = 1, inbp |
---|
912 | DO jj = mj0(ijrnf(jl,jr)), mj1(ijrnf(jl,jr)) |
---|
913 | DO ji = mi0(iirnf(jl,jr)), mi1(iirnf(jl,jr)) |
---|
914 | runoff(ji,jj) = runoff(ji,jj) & |
---|
915 | & + ( (1.-zxy) * zrnfm(nrnf1,jr) & |
---|
916 | & + zxy * zrnfm(nrnf2,jr) ) / FLOAT(inbp) |
---|
917 | END DO |
---|
918 | END DO |
---|
919 | END DO |
---|
920 | ENDIF |
---|
921 | END DO |
---|
922 | |
---|
923 | |
---|
924 | ! 2. conversion in mm/day |
---|
925 | ! ----------------------- |
---|
926 | ! C a u t i o n : upward water flux, runoff is negative |
---|
927 | |
---|
928 | runoff(:,:) = -zconv * runoff(:,:) / ( e1t(:,:) * e2t(:,:) ) |
---|
929 | |
---|
930 | |
---|
931 | ! 3. advectve coefficient |
---|
932 | ! ----------------------- |
---|
933 | |
---|
934 | ! Creates the array coef that contains the coefficient to affect to |
---|
935 | ! the upstream scheme. advection scheme will be: |
---|
936 | ! coefr * upstream + (1- coefr) centered |
---|
937 | ! coefr must be between 0 and 1. |
---|
938 | ! 3 values are used: maximum = zrup1, between zrup2, small zrup3. |
---|
939 | ! |
---|
940 | IF( nrunoff == 2 .AND. kt == nit000 ) THEN |
---|
941 | |
---|
942 | zrup1 = 0.5 |
---|
943 | zrup2 = 0.25 |
---|
944 | zrup3 = 0.125 |
---|
945 | DO jr = 1, jpriv-1 |
---|
946 | inbp = inb(jr) |
---|
947 | IF( inbp /= 0) THEN |
---|
948 | DO jl = 1, inbp |
---|
949 | ! ji,jj |
---|
950 | DO jj = mj0(ijrnf(jl,jr)), mj1(ijrnf(jl,jr)) |
---|
951 | DO ji = mi0(iirnf(jl,jr)), mi1(iirnf(jl,jr)) |
---|
952 | DO jn = 1, 3 |
---|
953 | zcoefr(ji,jj,jn) = zrup1 |
---|
954 | END DO |
---|
955 | zcoefr(ji,jj,4) = zrup2 |
---|
956 | zcoefr(ji,jj,5) = zrup3 |
---|
957 | END DO |
---|
958 | END DO |
---|
959 | ! --> zrup2 (originally = 0.5) |
---|
960 | ! ji,jj+1 |
---|
961 | DO jj = mj0(ijrnf(jl,jr) + 1), mj1(ijrnf(jl,jr) + 1) |
---|
962 | DO ji = mi0(iirnf(jl,jr)), mi1(iirnf(jl,jr)) |
---|
963 | DO jn = 1, 3 |
---|
964 | zcoefr(ji,jj,jn) = MAX( zrup2, zcoefr(ji,jj,jn) ) |
---|
965 | END DO |
---|
966 | zcoefr(ji,jj,4) = MAX( zrup3, zcoefr(ji,jj,4) ) |
---|
967 | END DO |
---|
968 | END DO |
---|
969 | ! ji,jj-1 |
---|
970 | DO jj = mj0(ijrnf(jl,jr) - 1), mj1(ijrnf(jl,jr) - 1) |
---|
971 | DO ji = mi0(iirnf(jl,jr)), mi1(iirnf(jl,jr)) |
---|
972 | DO jn = 1, 3 |
---|
973 | zcoefr(ji,jj,jn) = MAX( zrup2, zcoefr(ji,jj,jn) ) |
---|
974 | END DO |
---|
975 | zcoefr(ji,jj,4) = MAX( zrup3, zcoefr(ji,jj,4) ) |
---|
976 | END DO |
---|
977 | END DO |
---|
978 | ! ji+1,jj-1 |
---|
979 | DO jj = mj0(ijrnf(jl,jr) - 1), mj1(ijrnf(jl,jr) - 1) |
---|
980 | DO ji = mi0(iirnf(jl,jr) + 1), mi1(iirnf(jl,jr) + 1) |
---|
981 | DO jn = 1, 3 |
---|
982 | zcoefr(ji,jj,jn) = MAX( zrup2, zcoefr(ji,jj,jn) ) |
---|
983 | END DO |
---|
984 | zcoefr(ji,jj,4) = MAX( zrup3, zcoefr(ji,jj,4) ) |
---|
985 | END DO |
---|
986 | END DO |
---|
987 | ! ji+1,jj |
---|
988 | DO jj = mj0(ijrnf(jl,jr)), mj1(ijrnf(jl,jr)) |
---|
989 | DO ji = mi0(iirnf(jl,jr) + 1), mi1(iirnf(jl,jr) + 1) |
---|
990 | DO jn = 1, 3 |
---|
991 | zcoefr(ji,jj,jn) = MAX( zrup2, zcoefr(ji,jj,jn) ) |
---|
992 | END DO |
---|
993 | zcoefr(ji,jj,4) = MAX( zrup3, zcoefr(ji,jj,4) ) |
---|
994 | END DO |
---|
995 | END DO |
---|
996 | ! ji+1,jj+1 |
---|
997 | DO jj = mj0(ijrnf(jl,jr) + 1), mj1(ijrnf(jl,jr) + 1) |
---|
998 | DO ji = mi0(iirnf(jl,jr) + 1), mi1(iirnf(jl,jr) + 1) |
---|
999 | DO jn = 1, 3 |
---|
1000 | zcoefr(ji,jj,jn) = MAX( zrup2, zcoefr(ji,jj,jn) ) |
---|
1001 | END DO |
---|
1002 | zcoefr(ji,jj,4) = MAX( zrup3, zcoefr(ji,jj,4) ) |
---|
1003 | END DO |
---|
1004 | END DO |
---|
1005 | ! ji-1,jj+1 |
---|
1006 | DO jj = mj0(ijrnf(jl,jr) + 1), mj1(ijrnf(jl,jr) + 1) |
---|
1007 | DO ji = mi0(iirnf(jl,jr) - 1), mi1(iirnf(jl,jr) - 1) |
---|
1008 | DO jn = 1, 3 |
---|
1009 | zcoefr(ji,jj,jn) = MAX( zrup2, zcoefr(ji,jj,jn) ) |
---|
1010 | END DO |
---|
1011 | zcoefr(ji,jj,4) = MAX( zrup3, zcoefr(ji,jj,4) ) |
---|
1012 | END DO |
---|
1013 | END DO |
---|
1014 | ! ji-1,jj |
---|
1015 | DO jj = mj0(ijrnf(jl,jr)), mj1(ijrnf(jl,jr)) |
---|
1016 | DO ji = mi0(iirnf(jl,jr) - 1), mi1(iirnf(jl,jr) - 1) |
---|
1017 | DO jn = 1, 3 |
---|
1018 | zcoefr(ji,jj,jn) = MAX( zrup2, zcoefr(ji,jj,jn) ) |
---|
1019 | END DO |
---|
1020 | zcoefr(ji,jj,4) = MAX( zrup3, zcoefr(ji,jj,4) ) |
---|
1021 | END DO |
---|
1022 | END DO |
---|
1023 | ! ji-1,jj-1 |
---|
1024 | DO jj = mj0(ijrnf(jl,jr) - 1), mj1(ijrnf(jl,jr) - 1) |
---|
1025 | DO ji = mi0(iirnf(jl,jr) - 1), mi1(iirnf(jl,jr) - 1) |
---|
1026 | DO jn = 1, 3 |
---|
1027 | zcoefr(ji,jj,jn) = MAX( zrup2, zcoefr(ji,jj,jn) ) |
---|
1028 | END DO |
---|
1029 | zcoefr(ji,jj,4) = MAX( zrup3, zcoefr(ji,jj,4) ) |
---|
1030 | END DO |
---|
1031 | END DO |
---|
1032 | ! --> zrup3 (originally = 0.25) |
---|
1033 | ! ji-1,jj+2 |
---|
1034 | DO jj = mj0(ijrnf(jl,jr) + 2), mj1(ijrnf(jl,jr) + 2) |
---|
1035 | DO ji = mi0(iirnf(jl,jr) - 1), mi1(iirnf(jl,jr) - 1) |
---|
1036 | DO jn = 1, 3 |
---|
1037 | zcoefr(ji,jj,jn) = MAX( zrup3, zcoefr(ji,jj,jn) ) |
---|
1038 | END DO |
---|
1039 | END DO |
---|
1040 | END DO |
---|
1041 | ! ji-1,jj-2 |
---|
1042 | DO jj = mj0(ijrnf(jl,jr) - 2), mj1(ijrnf(jl,jr) - 2) |
---|
1043 | DO ji = mi0(iirnf(jl,jr) - 1), mi1(iirnf(jl,jr) - 1) |
---|
1044 | DO jn = 1, 3 |
---|
1045 | zcoefr(ji,jj,jn) = MAX( zrup3, zcoefr(ji,jj,jn) ) |
---|
1046 | END DO |
---|
1047 | END DO |
---|
1048 | END DO |
---|
1049 | ! ji,jj-2 |
---|
1050 | DO jj = mj0(ijrnf(jl,jr) - 2), mj1(ijrnf(jl,jr) - 2) |
---|
1051 | DO ji = mi0(iirnf(jl,jr) ), mi1(iirnf(jl,jr) ) |
---|
1052 | DO jn = 1, 3 |
---|
1053 | zcoefr(ji,jj,jn) = MAX( zrup3, zcoefr(ji,jj,jn) ) |
---|
1054 | END DO |
---|
1055 | END DO |
---|
1056 | END DO |
---|
1057 | ! ji,jj+2 |
---|
1058 | DO jj = mj0(ijrnf(jl,jr) + 2), mj1(ijrnf(jl,jr) + 2) |
---|
1059 | DO ji = mi0(iirnf(jl,jr) ), mi1(iirnf(jl,jr) ) |
---|
1060 | DO jn = 1, 3 |
---|
1061 | zcoefr(ji,jj,jn) = MAX( zrup3, zcoefr(ji,jj,jn) ) |
---|
1062 | END DO |
---|
1063 | END DO |
---|
1064 | END DO |
---|
1065 | ! ji+1,jj-2 |
---|
1066 | DO jj = mj0(ijrnf(jl,jr) - 2), mj1(ijrnf(jl,jr) - 2) |
---|
1067 | DO ji = mi0(iirnf(jl,jr) + 1), mi1(iirnf(jl,jr) + 1) |
---|
1068 | DO jn = 1, 3 |
---|
1069 | zcoefr(ji,jj,jn) = MAX( zrup3, zcoefr(ji,jj,jn) ) |
---|
1070 | END DO |
---|
1071 | END DO |
---|
1072 | END DO |
---|
1073 | ! ji+1,jj+2 |
---|
1074 | DO jj = mj0(ijrnf(jl,jr) + 2), mj1(ijrnf(jl,jr) + 2) |
---|
1075 | DO ji = mi0(iirnf(jl,jr) + 1), mi1(iirnf(jl,jr) + 1) |
---|
1076 | DO jn = 1, 3 |
---|
1077 | zcoefr(ji,jj,jn) = MAX( zrup3, zcoefr(ji,jj,jn) ) |
---|
1078 | END DO |
---|
1079 | END DO |
---|
1080 | END DO |
---|
1081 | ! ji+2,jj-2 |
---|
1082 | DO jj = mj0(ijrnf(jl,jr) - 2), mj1(ijrnf(jl,jr) - 2) |
---|
1083 | DO ji = mi0(iirnf(jl,jr) + 2), mi1(iirnf(jl,jr) + 2) |
---|
1084 | DO jn = 1, 3 |
---|
1085 | zcoefr(ji,jj,jn) = MAX( zrup3, zcoefr(ji,jj,jn) ) |
---|
1086 | END DO |
---|
1087 | END DO |
---|
1088 | END DO |
---|
1089 | ! ji+2,jj+2 |
---|
1090 | DO jj = mj0(ijrnf(jl,jr) + 2), mj1(ijrnf(jl,jr) + 2) |
---|
1091 | DO ji = mi0(iirnf(jl,jr) + 1), mi1(iirnf(jl,jr) + 1) |
---|
1092 | DO jn = 1, 3 |
---|
1093 | zcoefr(ji,jj,jn) = MAX( zrup3, zcoefr(ji,jj,jn) ) |
---|
1094 | END DO |
---|
1095 | END DO |
---|
1096 | END DO |
---|
1097 | ! ji+2,jj-1 |
---|
1098 | DO jj = mj0(ijrnf(jl,jr) - 1), mj1(ijrnf(jl,jr) - 1) |
---|
1099 | DO ji = mi0(iirnf(jl,jr) + 2), mi1(iirnf(jl,jr) + 2) |
---|
1100 | DO jn = 1, 3 |
---|
1101 | zcoefr(ji,jj,jn) = MAX( zrup3, zcoefr(ji,jj,jn) ) |
---|
1102 | END DO |
---|
1103 | END DO |
---|
1104 | END DO |
---|
1105 | ! ji+2,jj+1 |
---|
1106 | DO jj = mj0(ijrnf(jl,jr) + 1), mj1(ijrnf(jl,jr) + 1) |
---|
1107 | DO ji = mi0(iirnf(jl,jr) + 2), mi1(iirnf(jl,jr) + 2) |
---|
1108 | DO jn = 1, 3 |
---|
1109 | zcoefr(ji,jj,jn) = MAX( zrup3, zcoefr(ji,jj,jn) ) |
---|
1110 | END DO |
---|
1111 | END DO |
---|
1112 | END DO |
---|
1113 | ! ji+2,jj |
---|
1114 | DO jj = mj0(ijrnf(jl,jr) ), mj1(ijrnf(jl,jr) ) |
---|
1115 | DO ji = mi0(iirnf(jl,jr) + 2), mi1(iirnf(jl,jr) + 2) |
---|
1116 | DO jn = 1, 3 |
---|
1117 | zcoefr(ji,jj,jn) = MAX( zrup3, zcoefr(ji,jj,jn) ) |
---|
1118 | END DO |
---|
1119 | END DO |
---|
1120 | END DO |
---|
1121 | ! ji-2,jj |
---|
1122 | DO jj = mj0(ijrnf(jl,jr) ), mj1(ijrnf(jl,jr) ) |
---|
1123 | DO ji = mi0(iirnf(jl,jr) - 2), mi1(iirnf(jl,jr) - 2) |
---|
1124 | DO jn = 1, 3 |
---|
1125 | zcoefr(ji,jj,jn) = MAX( zrup3, zcoefr(ji,jj,jn) ) |
---|
1126 | END DO |
---|
1127 | END DO |
---|
1128 | END DO |
---|
1129 | ! ji-2,jj-2 |
---|
1130 | DO jj = mj0(ijrnf(jl,jr) - 2), mj1(ijrnf(jl,jr) - 2) |
---|
1131 | DO ji = mi0(iirnf(jl,jr) - 2), mi1(iirnf(jl,jr) - 2) |
---|
1132 | DO jn = 1, 3 |
---|
1133 | zcoefr(ji,jj,jn) = MAX( zrup3, zcoefr(ji,jj,jn) ) |
---|
1134 | END DO |
---|
1135 | END DO |
---|
1136 | END DO |
---|
1137 | ! ji-2,jj+2 |
---|
1138 | DO jj = mj0(ijrnf(jl,jr) + 2), mj1(ijrnf(jl,jr) + 2) |
---|
1139 | DO ji = mi0(iirnf(jl,jr) - 2), mi1(iirnf(jl,jr) - 2) |
---|
1140 | DO jn = 1, 3 |
---|
1141 | zcoefr(ji,jj,jn) = MAX( zrup3, zcoefr(ji,jj,jn) ) |
---|
1142 | END DO |
---|
1143 | END DO |
---|
1144 | END DO |
---|
1145 | ! ji-2,jj-1 |
---|
1146 | DO jj = mj0(ijrnf(jl,jr) - 1), mj1(ijrnf(jl,jr) - 1) |
---|
1147 | DO ji = mi0(iirnf(jl,jr) - 2), mi1(iirnf(jl,jr) - 2) |
---|
1148 | DO jn = 1, 3 |
---|
1149 | zcoefr(ji,jj,jn) = MAX( zrup3, zcoefr(ji,jj,jn) ) |
---|
1150 | END DO |
---|
1151 | END DO |
---|
1152 | END DO |
---|
1153 | ! ji-2,jj+1 |
---|
1154 | DO jj = mj0(ijrnf(jl,jr) + 1), mj1(ijrnf(jl,jr) + 1) |
---|
1155 | DO ji = mi0(iirnf(jl,jr) - 2), mi1(iirnf(jl,jr) - 2) |
---|
1156 | DO jn = 1, 3 |
---|
1157 | zcoefr(ji,jj,jn) = MAX( zrup3, zcoefr(ji,jj,jn) ) |
---|
1158 | END DO |
---|
1159 | END DO |
---|
1160 | END DO |
---|
1161 | END DO |
---|
1162 | ENDIF |
---|
1163 | END DO |
---|
1164 | upsrnfh(:,:) = zcoefr(:,:,1) |
---|
1165 | upsrnfz(:) = 0.e0 |
---|
1166 | upsrnfz(1) = zrup1 |
---|
1167 | upsrnfz(2) = zrup1 |
---|
1168 | upsrnfz(3) = zrup1 |
---|
1169 | upsrnfz(4) = zrup2 |
---|
1170 | upsrnfz(5) = zrup3 |
---|
1171 | ENDIF |
---|
1172 | ENDIF |
---|
1173 | |
---|
1174 | #endif |
---|
1175 | |
---|
1176 | END SUBROUTINE flx_rnf |
---|