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.

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

Context Navigation

source: trunk/NEMO/OPA_SRC/SBC/flxrnf_ORCA_R05.h90 @ 247

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

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

Download in other formats:

Original Format