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 trunk/NEMO/OPA_SRC/SBC – NEMO

source: trunk/NEMO/OPA_SRC/SBC/sbcrnf_ORCA_R05.h90 @ 912

Last change on this file since 912 was 888, checked in by ctlod, 16 years ago

merge dev_001_SBC branche with the trunk to include the New Surface Module package, see ticket: #113

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