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