New URL for NEMO forge! http://forge.nemo-ocean.eu

Since March 2022 along with NEMO 4.2 release, the code development moved to a self-hosted GitLab.
This present forge is now archived and remained online for history.

sbcrnf_ORCA_R2.h90 in trunk/NEMO/OPA_SRC/SBC – NEMO

Context Navigation

source: trunk/NEMO/OPA_SRC/SBC/sbcrnf_ORCA_R2.h90 @ 702

Last change on this file since 702 was 702, checked in by smasson, 17 years ago
add first set of new surface module, see ticket:3
Property svn:executable set to ``*
File size: 58.8 KB

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

Note: See TracBrowser for help on using the repository browser.

Download in other formats:

Original Format