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_R05.h90 in branches/dev_003_CPL/NEMO/OPA_SRC/SBC – NEMO

Context Navigation

source: branches/dev_003_CPL/NEMO/OPA_SRC/SBC/sbcrnf_ORCA_R05.h90 @ 991

Last change on this file since 991 was 991, checked in by smasson, 16 years ago
dev_003_CPL: preliminary draft (not working), see ticket #155
File size: 66.8 KB

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

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

Download in other formats:

Original Format