Index: trunk/NEMO/OPA_SRC/SBC/sbcrnf_ORCA_R05.h90
===================================================================
--- trunk/NEMO/OPA_SRC/SBC/sbcrnf_ORCA_R05.h90	(revision 1122)
+++ 	(revision )
@@ -1,1537 +1,0 @@
-   !!----------------------------------------------------------------------
-   !!                  *** flxrnf_ORCA_R05.h90  ***
-   !!----------------------------------------------------------------------
-   !!   flx_rnf      : update ocean runoffs for ORCA_R05 configuration
-   !!----------------------------------------------------------------------
-   !! * Modules used
-   USE dom_oce         ! ocean space and time domain
-   USE phycst          ! physical constants
-   USE in_out_manager  ! I/O manager
-   USE daymod          ! calendar
-
-   IMPLICIT NONE
-   PRIVATE
-
-   !! * Routine accessibility
-   PUBLIC flx_rnf         ! routine call by step.F90
-
-   !! * Shared module variables
-   REAL(wp), PUBLIC, DIMENSION(jpi,jpj) ::   &  !:
-      runoff,            &  !: monthly runoff (kg/m2/s)
-      upsadv,            &  !: mixed adv scheme in straits vicinity (hori.)
-      upsrnfh               !: mixed adv scheme in runoffs vicinity (hori.)
-   REAL(wp), PUBLIC, DIMENSION(jpk) ::   &  !:
-      upsrnfz               !: mixed adv scheme in runoffs vicinity (vert.)
-   INTEGER , PUBLIC ::   &  !:
-      nrunoff =  0          !: runoff option (namelist)
-
-   !! * Module variable
-   REAL(wp), DIMENSION(jpi,jpj,2) ::   &  !:
-      rnfdta                !: monthly runoff data array (kg/m2/s)
-   INTEGER    ::         &  !:
-      nrnf1, nrnf2          !: first and second record used
-   !!----------------------------------------------------------------------
-   !!   OPA 9.0 , LOCEAN-IPSL (2005) 
-   !! $ Id: $ 
-   !! Software governed by the CeCILL licence (modipsl/doc/NEMO_CeCILL.txt)
-   !!----------------------------------------------------------------------
-
-CONTAINS
-
-   SUBROUTINE flx_rnf( kt )
-      !!----------------------------------------------------------------------
-      !!                  ***  ROUTINE flx_rnf  ***
-      !!       
-      !! ** Purpose :   Introduce a climatological run off forcing
-      !!
-      !! ** Method  :   ORCA R05 configuration
-      !!      Initialze each mouth of river with a monthly climatology 
-      !!      provided from different data.
-      !!     C a u t i o n : upward water flux, runoff is negative
-      !!                     set at the last loop of the routine
-      !!
-      !! ** Action :
-      !!
-      !! References : 
-      !!       J. D. Milliman and R. H. Meade, 1983 : world-wide delivery
-      !!          of river sediment to the oceans, journal of geology vol 91
-      !!          pp 1-21.
-      !!       G. L. Russell and J. R. Miller, 1990 : global river runoff
-      !!          calculated from a global atmospheric general circulation
-      !!          model, journal of hydrology, 117(1990), pp 241-254.
-      !!       F. Van Der Leeden, Troise F. L., Todd D. K. : the water
-      !!          encyclopedia, second edition, lewis publishers.
-      !!       J. W. Weatherly, J. E. Walsh : The effects of precipitation
-      !!          and river runoff in a coupled ice-ocean model of Arctic
-      !!          Climate dynamics 1996 12:785,798
-      !!       Jacobs et al. 1992. J. Glaciol. 38 (130) 375-387.
-      !!
-      !! History :
-      !!        !  94-10  (G.Madec, M. Pontaud, M. Imbard)  Original code
-      !!        !  97-03  (G.Madec)  time dependent version
-      !!        !  98-06  (J.M. Molines)  exact computation of zxy 
-      !!                         for months that are not 30 days
-      !!        !  98-07  (M. Imbard)  ORCA and mpp option
-      !!        !  99-08  (J.P. Boulanger H.L.Ayina)  New rivers and 
-      !!                         values given in m3/s 
-      !!        !  00-04  (G. Madec, K. Roberts) add antarctica ice discharge.
-      !!        !  00-11  (R. Hordoir, E. Durand)  NetCDF FORMAT
-      !!   8.5  !  02-09  (G. Madec)  F90: Free form and module
-      !!----------------------------------------------------------------------
-      !! * Modules used
-
-      !! * arguments
-      INTEGER, INTENT( in  ) ::   kt ! ocean time step
-
-#if defined key_coupled
-      runoff(:,:) = 0.0e0
-#else
-      !!----------------------------------------------------------------------
-      !!  ORCA_R05
-      !!----------------------------------------------------------------------
-      !! * Local declarations
-      INTEGER, PARAMETER ::   &
-         jpmois=12,   &  ! number of month in the year !
-         jpriv=200,   &  ! maximum number of rivers
-         jpcoef=20       ! maximum number of gridpoints for mouth rivers
-
-      INTEGER  ::   ji, jj, jl, jm, jn, jr, jc
-      INTEGER  ::   i15, imois, iman
-      INTEGER  ::   inbp, idbd, idmeom
-      INTEGER  ::   iirnf(jpcoef,jpriv),ijrnf(jpcoef,jpriv)
-      INTEGER  ::   inb(jpriv)
-
-      REAL(wp) ::   zconv, zxy
-      REAL(wp) ::   zrup1,zrup2,zrup3
-      REAL(wp) ::   zrnfm(jpmois, jpriv)
-      REAL(wp) ::   zcoefr(jpi, jpj, 5)
-      !!----------------------------------------------------------------------
-      !!  OPA 8.5, LODYC-IPSL (2002)
-      !!----------------------------------------------------------------------
-      ! ====
-      ! DATA
-      ! ====
-
-      ! I NORTH AMERICA
-      ! ---------------
-
-      ! 1-Saint Lawrence (Canada) 49N30 67W00               R2
-      ! Old=447 km3/year=14165m3/s; CLIPPER=12430m3/s; ratio=1.14
-      ! --> 13975.249m3/s
-      DATA inb(1)/3/
-      DATA (iirnf(jc,1),jc=1,jpcoef)/ 444, 444, 444, 17*0 /
-      DATA (ijrnf(jc,1),jc=1,jpcoef)/ 367, 368, 369, 17*0 /
-      DATA (zrnfm(jm,1),jm=1,jpmois)/   &
-          11718.,  11971.,  12466.,  14980.,  19257.,  16913.,   &
-          14211.,  13119.,  13367.,  13644.,  13463.,  12594./
-      ! 2-Hudson River (USA) 41N00 77W00                    R1
-      ! Old=12 km3/year=380m3/s; UNESCO(64-83,p255)=408m3/s
-      ! --> 407.416687m3/s
-      DATA inb(2)/1/
-      DATA (iirnf(jc,2),jc=1,jpcoef)/ 426, 19*0 /
-      DATA (ijrnf(jc,2),jc=1,jpcoef)/ 336, 19*0 /
-      DATA (zrnfm(jm,2),jm=1,jpmois)/   &
-          364.,    404.,    648.,    873.,    593.,    301.,   &
-          205.,    171.,    205.,    290.,    379.,    456./
-      ! 3-Missisipi and Atchafalaya (USA) 29N 90W           R1 
-      ! Old=580 km3/year=18380m3/s; CLIPPER=17599m3/s
-      ! --> 17599.334m3/s
-      DATA inb(3)/2/
-      DATA (iirnf(jc,3),jc=1,jpcoef)/ 397, 398, 18*0 /
-      DATA (ijrnf(jc,3),jc=1,jpcoef)/ 312, 313, 18*0 /
-      DATA (zrnfm(jm,3),jm=1,jpmois)/   &
-          19260.,  20813.,  23026.,  28160.,  26510.,  19650.,   &
-          14181.,  10676.,   9578.,  10167.,  11699.,  17472./
-      ! 4-Brazos (USA) 29N58 95W53                          R1
-      ! Old=7 km3/year=222m3/s; UNESCO(65-84, p298)=200m3/s
-      ! --> 199.666656m3/s
-      DATA inb(4)/1/
-      DATA (iirnf(jc,4),jc=1,jpcoef)/ 386, 19*0 /
-      DATA (ijrnf(jc,4),jc=1,jpcoef)/ 311, 19*0 /
-      DATA (zrnfm(jm,4),jm=1,jpmois)/   &
-          165.,    224.,    200.,    266.,    457.,    327.,   &
-          143.,     69.,    102.,    132.,    160.,    151./
-      ! 5-Colorado (Mexico) 31N54 114W57                    R5
-      ! Old=20 km3/year=634m3/s
-      ! --> 634.m3/s
-      DATA inb(5)/2/
-      DATA (iirnf(jc,5),jc=1,jpcoef)/ 347, 348, 18*0 /
-      DATA (ijrnf(jc,5),jc=1,jpcoef)/ 316, 316, 18*0 /
-      DATA (zrnfm(jm,5),jm=1,jpmois)/12*634./ 
-      ! 6-Columbia (USA) 45N60 121W17                       R1
-      ! Old=251 km3/year=7954m3/s; UNESCO(65-84, p312)=5368m3/s
-      ! --> 5368.41602m3/s
-      DATA inb(6)/2/
-      DATA (iirnf(jc,6),jc=1,jpcoef)/ 326, 326, 18*0 /
-      DATA (ijrnf(jc,6),jc=1,jpcoef)/ 360, 361, 18*0 /
-      DATA (zrnfm(jm,6),jm=1,jpmois)/   &
-          4752.,    5077.,    5633.,    6075.,    8422.,    9822.,   &
-          6350.,    4112.,    3273.,    3251.,    3548.,    4106./ 
-      ! 7-Fraser (USA)                                      R5
-      ! Old=112 km3/year=3550m3/s;not found in UNESCO
-      ! --> 3549.99976m3/s
-      DATA inb(7)/2/
-      DATA (iirnf(jc,7),jc=1,jpcoef)/ 317, 317, 18*0 /
-      DATA (ijrnf(jc,7),jc=1,jpcoef)/ 378, 379, 18*0 /
-      DATA (zrnfm(jm,7),jm=1,jpmois)/12*3550./ 
-      ! 8-Yukon (USA) 63N00 165W00                          R3
-      ! Old=195 km3/year=6179m3/s; UNESCO(65-77+80-84, p321)=4245; ratio=1.46
-      ! --> 6017.24951m3/s
-      DATA inb(8)/2/
-      DATA (iirnf(jc,8),jc=1,jpcoef)/ 244, 244, 18*0 /
-      DATA (ijrnf(jc,8),jc=1,jpcoef)/ 413, 414, 18*0 /
-      DATA (zrnfm(jm,8),jm=1,jpmois)/   &
-           1459.,    1213.,    1065.,    1100.,    9067.,   17675.,   &
-          10680.,   11271.,    8658.,    5453.,    2649.,    1917./ 
-      ! 9-Copper (USA) 61N47 144W45                         R2
-      ! Old=39 km3/year=1236m3/s; UNESCO(65-84, p317)=1052m3/s; ratio=1.18
-      ! --> 1236.25m3/s
-      DATA inb(9)/1/
-      DATA (iirnf(jc,9),jc=1,jpcoef)/ 274, 19*0 /
-      DATA (ijrnf(jc,9),jc=1,jpcoef)/ 409, 19*0 /
-      DATA (zrnfm(jm,9),jm=1,jpmois)/   &
-           219.,    195.,    181.,    214.,    904.,   2532.,   &
-          4026.,   3627.,   1632.,    697.,    354.,    254./ 
-      ! 10-Susitna (USA)                                    R5
-      ! Old=40 km3/year=1268m3/s
-      ! --> 1268.m3/s
-      DATA inb(10)/2/
-      DATA (iirnf(jc,10),jc=1,jpcoef)/ 263, 264, 18*0 /
-      DATA (ijrnf(jc,10),jc=1,jpcoef)/ 405, 405, 18*0 /
-      DATA (zrnfm(jm,10),jm=1,jpmois)/12*1268./ 
-      ! 11-McKenzie (Canada) 70N00 135W00                   R3
-      ! Old=306 km3/year=9697m3/s; UNESCO(67+72+7479+82-84, p168)=8343m3/s;
-      ! ratio=1.16
-      ! --> 9606.66797m3/s
-      DATA inb(11)/3/
-      DATA (iirnf(jc,11),jc=1,jpcoef)/ 272, 272, 272, 17*0 /
-      DATA (ijrnf(jc,11),jc=1,jpcoef)/ 451, 452, 453, 17*0 /
-      DATA (zrnfm(jm,11),jm=1,jpmois)/   &
-           3644.,    3400.,    3258.,    3546.,   14188.,   21544.,   &
-          18954.,   15487.,   11284.,   10169.,    5901.,    3905./
-      ! 12-Kuksokwim (USA) 60N00 ~162W30                    R4
-      ! UNESCO(65-84, p318)=1112m3/s (61N87 158E10, Crooked Creek, inland)
-      ! --> 1111.99988m3/s
-      DATA inb(12)/2/
-      DATA (iirnf(jc,12),jc=1,jpcoef)/ 244, 244, 18*0 /
-      DATA (ijrnf(jc,12),jc=1,jpcoef)/ 413, 414, 18*0 /
-      DATA (zrnfm(jm,12),jm=1,jpmois)/   &
-          335.,     308.,     295.,     407.,    2415.,    2027.,   &
-          1917.,   2049.,    1585.,    1072.,     545.,     389./     
-
-      DATA inb(13)/0/
-      DATA inb(14)/0/
-      DATA inb(15)/0/
-
-      ! Total run-offs(I)=60974.2539m3/s
-
-      ! II SOUTH AMERICA
-      ! -----------------
-      ! 16-Chira (Peru) 4s54 81w08                          R2
-      ! Old=5 km3/year=158m3/s; UNESCO (65-68, p357)=80.9; ratio=1.95
-      ! --> 158.000015m3/s
-      DATA inb(16)/1/
-      DATA (iirnf(jc,16),jc=1,jpcoef)/ 414, 19*0 /
-      DATA (ijrnf(jc,16),jc=1,jpcoef)/ 244, 19*0 /
-      DATA (zrnfm(jm,16),jm=1,jpmois)/   &
-          70.,     105.,     564.,     582.,     230.,      94.,   &
-          87.,      47.,      35.,      32.,      29.,      21./
-      ! 17-Magdalena (Colombia) 11N06 74W51                 R2
-      ! Old=237 km3/year=7510m3/s; UNESCO(71-73+76-79, p348)=6973m3/s;
-      ! ratio=1.08
-      ! --> 7451.58301m3/s
-      DATA inb(17)/2/
-      DATA (iirnf(jc,17),jc=1,jpcoef)/ 424, 424, 18*0 /
-      DATA (ijrnf(jc,17),jc=1,jpcoef)/ 270, 271, 18*0 /
-      DATA (zrnfm(jm,17),jm=1,jpmois)/   &
-          6598.,    4441.,    4087.,    5445.,    7380.,    8067.,   &
-          8509.,    7428.,    7631.,    9296.,   10588.,    9949./
-      ! 18-Orinoco (Venezuela) 8N37 62W15!                  R1
-      ! Old=1100 km3/year=34857m3/s; UNESCO(73-75, p364)=28736m3/s;
-      ! HLA=26371m3/s
-      ! Data from HLA: 8.1N - 63.3W   1925/89
-      ! --> 26371.416
-      DATA inb(18)/5/
-      DATA (iirnf(jc,18),jc=1,jpcoef)/ 451, 451, 452, 452, 452, 15*0 /
-      DATA (ijrnf(jc,18),jc=1,jpcoef)/ 271, 272, 271, 272, 272, 15*0 /
-      DATA (zrnfm(jm,18),jm=1,jpmois)/   &
-          10670.,    7513.,    6344.,    6887.,   13086.,   26699.,   &
-          43938.,   59957.,   57435.,   39951.,   27410.,   16567./ 
-      ! 19-Amazon (Brazil) 0S05 50W00                       R1
-      ! Old=6300 km3/year=199634m3/s; CLIPPER=166192m3/s~UNESCO;
-      ! HLA=184806m3/s
-      ! UNESCO data are at Obidos which is far inland and therefore does not
-      ! fully represent the actual discharge
-      ! Data from HLA: 1.5N - 55.3W   1970/96
-      ! --> 184845.484m3/s
-      DATA inb(19)/14/
-      DATA (iirnf(jc,19),jc=1,jpcoef)/ 475, 476, 476, 476, 476,  &
-                                       477, 477, 477, 477, 477,  &
-                                       478, 478, 478, 478, 6*0 /
-      DATA (ijrnf(jc,19),jc=1,jpcoef)/ 254, 251, 252, 253, 254,  &
-                                       250, 251, 252, 253, 254,  &
-                                       250, 251, 252, 253, 6*0 /
-      DATA (zrnfm(jm,19),jm=1,jpmois)/   &
-          152233., 176415., 199941., 225565., 244016., 244861.,   &
-          230711., 203128., 164531., 128486., 118214., 130045. /
-      ! 20-Sao Francisco (Brazil) 10S30 36W24               R1
-      ! Old=97 km3/year=3074m3/s; UNESCO(77-79, p338)=2645m3/s
-      ! --> 2645.3335m3/s
-      DATA inb(20)/1/
-      DATA (iirnf(jc,20),jc=1,jpcoef)/ 498, 19*0 /
-      DATA (ijrnf(jc,20),jc=1,jpcoef)/ 223, 19*0 /
-      DATA (zrnfm(jm,20),jm=1,jpmois)/   &
-          3244.,   4116.,   5139.,   4330.,   2323.,   1800.,     &
-          1587.,   1511.,   1593.,   1740.,   2190.,   2171. /
-      ! 21-Parana+Uruguay(Arg./Uru.) 34S12 58W18            R1 
-      ! Old=470 km3/year=14893m3/s; UNESCO(p323/324)=22263m3/s
-      ! --> 22203.168m3/s 
-      DATA inb(21)/5/
-      DATA (iirnf(jc,21),jc=1,jpcoef)/ 462, 462, 462, 462, 463, 15*0 /
-      DATA (ijrnf(jc,21),jc=1,jpcoef)/ 172, 173, 174, 175, 175, 15*0 /
-      DATA (zrnfm(jm,21),jm=1,jpmois)/   &
-          24803.,   26836.,   24585.,   21358.,   20041.,   21689.,   &
-          22647.,   20403.,   19627.,   21400.,   21126.,   21923./
-      ! 22-Negro (Argentina)  41S02  62W47                  R2
-      ! Old=30 km3/year=951m3/s; UNESCO(65-79,p330)=765m3/s; ratio=1.24
-      ! --> 907.666687m3/s
-      DATA inb(22)/2/
-      DATA (iirnf(jc,22),jc=1,jpcoef)/ 446, 447, 18*0 /
-      DATA (ijrnf(jc,22),jc=1,jpcoef)/ 160, 160, 18*0 /
-      DATA (zrnfm(jm,22),jm=1,jpmois)/   &
-          882.,     587.,     519.,     482.,     623.,     984.,   &
-          1120.,   1294.,    1092.,    1166.,    1105.,    1038./
-      ! 23- Tocantins :  ~48W00 ~0.5S                       R4
-      ! UNESCO(76-77+80-81, p335)=12167m3/s (station far inland!)
-      ! --> 12173.082m3/s
-      DATA inb(23)/4/
-      DATA (iirnf(jc,23),jc=1,jpcoef)/ 479, 480, 481, 482, 16*0 /
-      DATA (ijrnf(jc,23),jc=1,jpcoef)/ 249, 249, 249, 249, 16*0 /
-      DATA (zrnfm(jm,23),jm=1,jpmois)/   &
-          15507.,  24855.,  27854.,  24257.,  15860.,   7861.,   &
-          4440.,   3082.,   2728.,   3417.,   5511.,  10705. /
-
-      DATA inb(24)/0/
-      DATA inb(25)/0/
-
-      ! Total run-offs(II)=256755.75
-      ! Total run-offs(I+II)=317730
-
-      ! III EUROPE
-      ! ----------
-      ! 26-Rhone (France) 43N20 4E50                        R1 
-      ! Old=49 km3/year=1553m3/s; UNESCO(65-83, p512)=1229
-      ! --> 1217.16663m3/s
-      DATA inb(26)/1/
-      DATA (iirnf(jc,26),jc=1,jpcoef)/ 584, 19*0 /
-      DATA (ijrnf(jc,26),jc=1,jpcoef)/ 347, 19*0 /
-      DATA (zrnfm(jm,26),jm=1,jpmois)/   &
-          1385.,    1848.,    1568.,    1311.,    1420.,    1191.,   &
-           937.,     796.,     794.,     962.,    1161.,    1233./
-      ! 27-Po (Italy) 11N10 1W09                            R1 
-      ! Old=46 km3/year=1458m3/s; UNESCO(65-84, p527)=1550m3/s
-      ! --> 1550.08337m3/s
-      DATA inb(27)/1/
-      DATA (iirnf(jc,27),jc=1,jpcoef)/ 599, 19*0 /
-      DATA (ijrnf(jc,27),jc=1,jpcoef)/ 355, 19*0 /
-      DATA (zrnfm(jm,27),jm=1,jpmois)/   &
-          1304.,    1605.,    1856.,    1549.,    1958.,    1760.,   &
-          1040.,    1007.,    1615.,    1907.,    1772.,    1228./
-      ! 28-Danube (Roumania) 45N20 29E40                    R1 
-      ! Old=206 km3/year=6527m3/s; UNESCO(65-84, p544)=7011m3/s
-      ! --> 7011.8335m3/s
-      DATA inb(28)/2/
-      DATA (iirnf(jc,28),jc=1,jpcoef)/ 632, 632, 18*0 /
-      DATA (ijrnf(jc,28),jc=1,jpcoef)/ 356, 357, 18*0 /
-      DATA (zrnfm(jm,28),jm=1,jpmois)/   &
-          6625.,    7352.,    7986.,    9279.,    9482.,    8618.,   &
-          7586.,    6001.,    5210.,    4989.,    4921.,    6093./
-
-      DATA inb(29)/0/
-      DATA inb(30)/0/
-
-      ! Total run-offs(III)=9779.08301m3/s
-      ! Total run-offs(I+II+III)=327509.094m3/s
-      
-      ! IV EURASIAN ARCTIC
-      ! ------------------
-      ! 31-Iana (Russia)                                    R1
-      ! Old=29 km3/year=919m3/s; UNESCO(65-72+74+76-84, p463)=936m3/s
-      ! --> 0m3/s
-      DATA inb(31)/0/
-      DATA (iirnf(jc,31),jc=1,jpcoef)/ 20*0 /
-      DATA (ijrnf(jc,31),jc=1,jpcoef)/ 20*0 /
-      DATA (zrnfm(jm,31),jm=1,jpmois)/   &
-             2.,       1.,       0.,       0.,     917.,    3514.,   &
-          2968.,    2329.,    1322.,     178.,      36.,       8./ 
-      ! 32-Ob (Russia) 66N45 69E30                          R1 
-      ! Old=385 km3/year=12200m3/s; UNESCO(65-84, p478)=12862m3/s
-      ! --> 12861.5m3/s
-      DATA inb(32)/8/
-      DATA (iirnf(jc,32),jc=1,jpcoef)/ 635, 635, 636, 636, 636,637, 637,637,12*0 /
-      DATA (ijrnf(jc,32),jc=1,jpcoef)/ 508, 509, 509, 508, 507,509, 508,507,12*0 /
-     DATA (zrnfm(jm,32),jm=1,jpmois)/12*0./
-!CT bug     DATA (zrnfm(jm,32),jm=1,jpmois)/   &
-!CT bug          4977.,    4150.,    3611.,    3590.,   14723.,   33366.,   &
-!CT bug         30773.,   22785.,   14701.,    9705.,    6078.,    5879./
-      ! 33-Yenesei (Russia) 71N50 82E40                     R1 
-      ! Old=560 km3/year=17745m3/s; UNESCO(65-84, p472)=17462m3/s
-      ! --> 17461.666m3/s
-      DATA inb(33)/9/
-      DATA (iirnf(jc,33),jc=1,jpcoef)/ 114,115,115,116, 117,117,117,118,118,11*0 /
-      DATA (ijrnf(jc,33),jc=1,jpcoef)/ 504,504,503,503, 501,502,503,502,501,11*0 /
-      DATA (zrnfm(jm,33),jm=1,jpmois)/   &
-           6659.,    7007.,    6895.,    7089.,   24183.,   74965.,   &
-          25325.,   15940.,   16050.,   12678.,    6741.,    6008./
-      ! 34-Severnaya Dvina(Russia)64N32 40E30               R5 
-      ! Old=106 km3/year=3359m3/s; not in UNESCO data base
-      ! --> 3359.00024m3/s
-      DATA inb(34)/2/
-      DATA (iirnf(jc,34),jc=1,jpcoef)/ 620, 620, 18*0 /
-      DATA (ijrnf(jc,34),jc=1,jpcoef)/ 448, 449, 18*0 /
-      DATA (zrnfm(jm,34),jm=1,jpmois)/12*3359./
-      ! 35-Lena (Russia) 72N25 126E40                       R1 
-      ! --> 16192.333m3/s
-      DATA inb(35)/9/
-      DATA (iirnf(jc,35),jc=1,jpcoef)/ 142,143,148,148,149,149,150,150,150,11*0 /
-      DATA (ijrnf(jc,35),jc=1,jpcoef)/ 458,458,456,448,449,448,451,450,449,11*0 /
-      DATA (zrnfm(jm,35),jm=1,jpmois)/   &
-           2663.,    2095.,    1629.,    1324.,    7203.,   64230.,   &
-          40810.,   28310.,   25995.,   13299.,    3673.,    3077./
-      ! 36-Kolyma (Russia) 69N30 161E00                     R1
-      ! Old=71 km3/year=2250m3/s; UNESCO(65-84)=2138m3/s
-      ! --> 2138.66675m3/s
-      DATA inb(36)/2/
-      DATA (iirnf(jc,36),jc=1,jpcoef)/ 184, 185, 18*0 /
-      DATA (ijrnf(jc,36),jc=1,jpcoef)/ 429, 429, 18*0 /
-      DATA (zrnfm(jm,36),jm=1,jpmois)/   &
-           131.,      92.,      76.,      65.,    2576.,    9143.,   &
-          4670.,    3727.,    3226.,    1119.,     630.,     209./
-      ! 37-Indigirka (russia) 70N48 148E54                  R1 
-      ! Old=55 km3/year=1743m3/s; UNESCO(65-68+70+72-84)=1622.
-      ! --> 1622.00012m3/s
-      DATA inb(37)/1/
-      DATA (iirnf(jc,37),jc=1,jpcoef)/ 172, 19*0 /
-      DATA (ijrnf(jc,37),jc=1,jpcoef)/ 438, 19*0 /
-      DATA (zrnfm(jm,37),jm=1,jpmois)/   &
-            32.,      17.,      11.,       7.,     354.,    5629.,   &
-          5210.,    4617.,    2823.,     559.,     142.,      63./
-
-      DATA inb(38)/0/
-      DATA inb(39)/0/
-      DATA inb(40)/0/
-
-      ! Total run-offs(IV)=40773.6641m3/s
-      ! Total run-offs(I+II+III+IV)=368282.75m3/s
-
-      ! V ASIA
-      ! ------
-      ! 41-Amur (Russia) 52N56 141E10                       R2
-      ! Old=325 km3/year=10299m3/s; UNESCO(65-84, p451)=8792m3/s; ratio=1.17
-      ! --> 10299.332m3/s
-      DATA inb(41)/2/
-      DATA (iirnf(jc,41),jc=1,jpcoef)/ 141, 141, 18*0 /
-      DATA (ijrnf(jc,41),jc=1,jpcoef)/ 381, 382, 18*0 /
-      DATA (zrnfm(jm,41),jm=1,jpmois)/   &
-           2108.,    1543.,    1252.,    4088.,   15464.,   17336.,   &
-          16100.,   20801.,   21548.,   15658.,    5188.,    2506./
-      ! 42-Liaohe (China) 40N40 122E09                      R5
-      ! Old=6 km3/year=190m3/s; not found in UNESCO data base
-      ! --> 189.999985m3/s
-      DATA inb(42)/1/
-      DATA (iirnf(jc,42),jc=1,jpcoef)/ 100, 19*0 /
-      DATA (ijrnf(jc,42),jc=1,jpcoef)/ 342, 19*0 /
-      DATA (zrnfm(jm,42),jm=1,jpmois)/12*190./
-      ! 43-Daling (China) 40N56 121E43                      R5
-      ! Old=1 km3/year=32m3/s; not found in UNESCO data base
-      ! --> 31.9999962m3/s
-      DATA inb(43)/2/
-      DATA (iirnf(jc,43),jc=1,jpcoef)/  99, 100, 18*0 /
-      DATA (ijrnf(jc,43),jc=1,jpcoef)/ 342, 342, 18*0 /
-      DATA (zrnfm(jm,43),jm=1,jpmois)/12*32./
-      ! 44-Haiho (China) 6N17   1E23
-      ! Old=2 km3/year=64m3/s; not found in UNESCO data base
-      ! --> 63.9999924m3/s
-      DATA inb(44)/1/
-      DATA (iirnf(jc,44),jc=1,jpcoef)/ 121, 19*0 /
-      DATA (ijrnf(jc,44),jc=1,jpcoef)/ 348, 19*0 /
-      DATA (zrnfm(jm,44),jm=1,jpmois)/12*64./
-      ! 45-Yellow/Huangho (China) 37N32 118E19              R2 
-      ! Old=49 km3/year=1553m3/s; UNESCO(p371)=1210m3/s; ratio=1.28
-      ! --> 1552.5m3/s
-      DATA inb(45)/3/
-      DATA (iirnf(jc,45),jc=1,jpcoef)/  92,  92,  92, 17*0 /
-      DATA (ijrnf(jc,45),jc=1,jpcoef)/ 336, 337, 338, 17*0 /
-      DATA (zrnfm(jm,45),jm=1,jpmois)/   &
-           697.,     674.,    1282.,    1206.,    1109.,     778.,   &
-          1964.,    3270.,    3668.,    2133.,     978.,     871./
-      ! 46-Yangtze (China) 31N48 121E10=Changjiang          R2
-      ! Old=900 km3/year=28520m3/s; UNESCO(p372)=25032m3/s; ratio=1.14
-      ! --> 28519.584m3/s
-      DATA inb(46)/2/
-      DATA (iirnf(jc,46),jc=1,jpcoef)/ 100, 100, 18*0 /
-      DATA (ijrnf(jc,46),jc=1,jpcoef)/ 316, 317, 18*0 /
-      DATA (zrnfm(jm,46),jm=1,jpmois)/   &
-          10994.,   10548.,   13254.,   24097.,   36972.,   47197.,   &
-          53891.,   42070.,   37627.,   30762.,   22046.,   12777./
-      ! 47-Pearl/Zhujiang (China) 23N36 106E34=Xijiang      R2 
-      ! Old=302 km3/year=9570m3/s; UNESCO(76,82, p373)=6987m3/s; ratio=1.37
-      ! --> 9570.75m3/s
-      DATA inb(47)/2/
-      DATA (iirnf(jc,47),jc=1,jpcoef)/  72,  73, 18*0 /
-      DATA (ijrnf(jc,47),jc=1,jpcoef)/ 294, 294, 18*0 /
-      DATA (zrnfm(jm,47),jm=1,jpmois)/   &
-           2357.,    2415.,    2634.,    6446.,   14246.,   18395.,   &
-          19213.,   17558.,   12312.,    7255.,    5864.,    6154./
-      ! 48-Red/Hungho (Vietnam) 20N17 106E34                R5 
-      ! Old=123 km3/year=3897m3/s; not in UNESCO data base
-      ! --> 3897.m3/s
-      DATA inb(48)/3/
-      DATA (iirnf(jc,48),jc=1,jpcoef)/  68,  68,  68, 17*0 /
-      DATA (ijrnf(jc,48),jc=1,jpcoef)/ 288, 289, 290, 17*0 /
-      DATA (zrnfm(jm,48),jm=1,jpmois)/12*3897./
-      ! 49-Mekong (Vietnam) 10N33 105E24                    R5
-      ! Old=470 km3/year=14893m3/s; not in UNESCO data base
-      ! --> 14893.0029m3/s
-      DATA inb(49)/3/
-      DATA (iirnf(jc,49),jc=1,jpcoef)/  69,  69,  70, 17*0 /
-      DATA (ijrnf(jc,49),jc=1,jpcoef)/ 269, 270, 270, 17*0 /
-      DATA (zrnfm(jm,49),jm=1,jpmois)/12*14893./
-
-      DATA inb(50)/0/
-      DATA inb(51)/0/
-      DATA inb(52)/0/
-
-      ! 53-Damodar (India)                                  R2
-      ! Old=10 km3/year=320m3/s; UNESCO(p386)=173m3/s; ratio=1.85
-      ! --> 432.916656m3/s
-      DATA inb(53)/1/
-      DATA (iirnf(jc,53),jc=1,jpcoef)/  30, 19*0 /
-      DATA (ijrnf(jc,53),jc=1,jpcoef)/ 294, 19*0 /
-      DATA (zrnfm(jm,53),jm=1,jpmois)/   &
-           51.,      37.,      42.,      49.,      85.,     296.,   &
-          896.,    1390.,    1591.,     542.,     172.,      44./   
-
-      DATA inb(54)/0/
-      DATA inb(55)/0/
-      DATA inb(56)/0/
-      DATA inb(57)/0/
-      DATA inb(58)/0/
-      DATA inb(59)/0/
-      DATA inb(60)/0/
-
-      ! Total run-offs (V)=128642.922m3/s
-      ! Total run-offs (I+II+III+IV+V)=496925.688m3/s
-
-      ! VI AFRICA
-      ! ----------
-      ! 61-Nile (Egypt) 30N10  31E06                        R1
-      ! Old=30 km3/year=951m3/s; UNESCO(73-84)=1251m3/s
-      ! --> 1251.33337m3/s
-      DATA inb(61)/3/
-      DATA (iirnf(jc,61),jc=1,jpcoef)/ 636, 636, 637, 17*0 /
-      DATA (ijrnf(jc,61),jc=1,jpcoef)/ 317, 318, 318, 17*0 /
-      DATA (zrnfm(jm,61),jm=1,jpmois)/   &
-          1239.,    1035.,    1126.,    1110.,    1170.,    1542.,   &
-          1742.,    1559.,    1191.,    1103.,    1075.,    1124./   
-      ! 62-Niger (Nigeria) 5N33   6E33                      R1
-      ! Old=192 km3/year=6084m3/s; HLA=5545m3/s;
-      ! Data from HLA: 4.5N -  6  E   1950/87
-      ! --> 5545.08301m3/s
-      DATA inb(62)/3/
-      DATA (iirnf(jc,62),jc=1,jpcoef)/ 591, 592, 593, 17*0 /
-      DATA (ijrnf(jc,62),jc=1,jpcoef)/ 258, 258, 258, 17*0 /
-      DATA (zrnfm(jm,62),jm=1,jpmois)/   &
-          2029.,    1804.,    1573.,    1418.,    1494.,    2561.,   &
-          5272.,    9629.,   14951.,   16227.,    6842.,    2741./   
-      ! 63-Zaire/Congo (Congo) 6S04  12E24                  R1      
-      ! Old=1250 km3/year=39610m3/s; UNESCO(65-83, p149)=41941m3/s; HLA
-      ! =39720m3/s
-      ! Data from HLA: 6S - 12  E   1902/96
-      ! --> 39719.5m3/s
-      DATA inb(63)/4/
-      DATA (iirnf(jc,63),jc=1,jpcoef)/ 599, 599, 600, 600, 16*0 /
-      DATA (ijrnf(jc,63),jc=1,jpcoef)/ 237, 238, 236, 237, 16*0 /
-      DATA (zrnfm(jm,63),jm=1,jpmois)/   &
-          47200.,   37692.,   34374.,   36603.,   37862.,   35588.,   &
-          30966.,   30774.,   35981.,   42825.,   51455.,   55314./  
-      ! 64-Orange (Sth Africa) 28S41 16E28                  R5
-      ! Old=11 km3/year=349m3/s
-      ! --> 349.000031m3/s
-      DATA inb(64)/2/
-      DATA (iirnf(jc,64),jc=1,jpcoef)/ 609, 609, 18*0 /
-      DATA (ijrnf(jc,64),jc=1,jpcoef)/ 188, 189, 18*0 /
-      DATA (zrnfm(jm,64),jm=1,jpmois)/12*349./
-      ! 65-Zambesi (Mozambique) 18S55 36E04                 R5
-      ! Old=223 km3/year=7066m3/s
-      ! --> 7066.00049m3/s
-      DATA inb(65)/3/
-      DATA (iirnf(jc,65),jc=1,jpcoef)/ 646, 646, 646, 17*0 /
-      DATA (ijrnf(jc,65),jc=1,jpcoef)/ 207, 208, 209, 17*0 /
-      DATA (zrnfm(jm,65),jm=1,jpmois)/12*7066./
-      ! 66-Limpopo (mozambique) 25S15 33E30                 R5
-      ! Old=5 km3/year=158m3/s
-      ! --> 158.m3/s
-      DATA inb(66)/2/
-      DATA (iirnf(jc,66),jc=1,jpcoef)/ 642, 642, 18*0 /
-      DATA (ijrnf(jc,66),jc=1,jpcoef)/ 196, 197, 18*0 /
-      DATA (zrnfm(jm,66),jm=1,jpmois)/12*158./
-      ! 67-Rufiji (tanzania) 8S00 39E20                     R5
-      ! Old=9 km3/year=285m3/s
-      ! --> 285.m3/s
-      DATA inb(67)/1/
-      DATA (iirnf(jc,67),jc=1,jpcoef)/ 655, 19*0 /
-      DATA (ijrnf(jc,67),jc=1,jpcoef)/ 235, 19*0 /
-      DATA (zrnfm(jm,67),jm=1,jpmois)/12*285./
-      ! Other rivers from the Gulf of Guinea
-      ! 68-Kouilou (Congo)  4.5S-11.5E (UNESCO=p51)         R1
-      ! Data from HLA:  4S - 11  E   1953/86 = 919m3/s
-      ! --> 924.25m3/s
-      DATA inb(68)/1/
-      DATA (iirnf(jc,68),jc=1,jpcoef)/ 599, 19*0 /
-      DATA (ijrnf(jc,68),jc=1,jpcoef)/ 240, 19*0 /
-      DATA (zrnfm(jm,68),jm=1,jpmois)/   &
-          1115.,    1166.,    1285.,    1462.,    1454.,     725.,   &
-           480.,     378.,     327.,     400.,    1018.,    1281./
-      ! 69-Nyanga (Gabon)  3S-10E (UNESCO=p66)              R1
-      ! Data from HLA: 2.8N - 10  E   1952/82 = 264m3/s
-      ! --> 264.166687m3/s
-      DATA inb(69)/1/
-      DATA (iirnf(jc,69),jc=1,jpcoef)/ 595, 19*0 /
-      DATA (ijrnf(jc,69),jc=1,jpcoef)/ 244, 19*0 /
-      DATA (zrnfm(jm,69),jm=1,jpmois)/   &
-          126.,      99.,     143.,     253.,     316.,     300.,   &
-          156.,      93.,     231.,     557.,     592.,     304./
-      ! 70-Ogoue (Gabon) 1S-9E (UNESCO, p69)                R1
-      ! Data from HLA: .75S -  9  E   1930/90 = 4512m3/s
-      ! --> 4511.83301m3/s
-      DATA inb(70)/2/
-      DATA (iirnf(jc,70),jc=1,jpcoef)/ 593, 593, 18*0 /
-      DATA (ijrnf(jc,70),jc=1,jpcoef)/ 247, 248, 18*0 /
-      DATA (zrnfm(jm,70),jm=1,jpmois)/   &
-          4514.,    3988.,    4622.,    5403.,    5790.,    4145.,   &
-          2616.,    1961.,    2119.,    4304.,    7870.,    6810./
-      ! 71-Sanaga (Cameroun) 3.5N-9.5E (UNESCO, p31)        R1
-      ! Data from HLA: 3.9N -  9.6E   1944/80 = 1983m3/s
-      ! --> 1982.58337m3/s
-      DATA inb(71)/2/
-      DATA (iirnf(jc,71),jc=1,jpcoef)/ 595, 595, 18*0 /
-      DATA (ijrnf(jc,71),jc=1,jpcoef)/ 255, 256, 18*0 /
-      DATA (zrnfm(jm,71),jm=1,jpmois)/   &
-           813.,     578.,     517.,     633.,     909.,    1278.,   &
-          1928.,    2653.,    4361.,    5470.,    3307.,    1344./
-      ! 72-Nyong (Cameroun) 3.1N-9.5E (UNESCO, p36)         R1
-      ! Data from HLA: 3.5N -  9.8E   1951/84 = 432m3/s
-      ! --> 431.166656m3/s      
-      DATA inb(72)/1/
-      DATA (iirnf(jc,72),jc=1,jpcoef)/ 594, 19*0 /
-      DATA (ijrnf(jc,72),jc=1,jpcoef)/ 257, 19*0 /
-      DATA (zrnfm(jm,72),jm=1,jpmois)/   &
-          232.,     128.,     150.,     266.,     402.,     459.,   &
-          353.,     306.,     604.,     930.,     878.,     466./
-      ! 73-Ntem (Cameroun) 2.5N-9.5E (UNESCO, p36)          R1
-      ! Data from HLA: 3N -  9.8E   1951/87 = 330m3/s      
-      ! --> 329.75m3/s
-      DATA inb(73)/1/
-      DATA (iirnf(jc,73),jc=1,jpcoef)/ 594, 19*0 /
-      DATA (ijrnf(jc,73),jc=1,jpcoef)/ 254, 19*0 /
-      DATA (zrnfm(jm,73),jm=1,jpmois)/   &
-          452.,     455.,     492.,     466.,     482.,     251.,   &
-          151.,     108.,      87.,     133.,     395.,     485./
-      ! 74-Wouri (Cameroun) 4N-9.5E (UNESCO, p34)           R1
-      ! Data from HLA: 4N -  9.5E   1951/84  = 307m3/s
-      ! --> 306.583374m3/s
-      DATA inb(74)/1/
-      DATA (iirnf(jc,74),jc=1,jpcoef)/ 594, 19*0 /
-      DATA (ijrnf(jc,74),jc=1,jpcoef)/ 257, 19*0 /
-      DATA (zrnfm(jm,74),jm=1,jpmois)/   &
-           88.,      63.,      68.,      95.,     139.,     222.,   &
-          395.,     633.,     825.,     675.,     328.,     148./
-      ! 75-Oueme (Benin) 6.5N-2.5E (UNESCO, p27)            R1
-      ! Data from HLA:  6.2N -  2.5E   1951/84 =113m3/s
-      ! --> 112.666664m3/s
-      DATA inb(75)/1/
-      DATA (iirnf(jc,75),jc=1,jpcoef)/ 583, 19*0 /
-      DATA (ijrnf(jc,75),jc=1,jpcoef)/ 262, 19*0 /
-      DATA (zrnfm(jm,75),jm=1,jpmois)/   &
-           0.,       0.,       0.,       0.,       1.,      18.,   &
-          96.,     324.,     531.,     319.,      58.,       5./
-      ! 76-Mono(Tepetou,Benin/Togo) 6.25N-1.8E (U, p27)     R1
-      ! Data from HLA: 6.2N -  1  E   1951/85  = 96m3/s       
-      ! --> 96.4999924m3/s
-      DATA inb(76)/2/
-      DATA (iirnf(jc,76),jc=1,jpcoef)/ 578, 579, 18*0 /
-      DATA (ijrnf(jc,76),jc=1,jpcoef)/ 262, 262, 18*0 /
-      DATA (zrnfm(jm,76),jm=1,jpmois)/   &
-            2.,       1.,       1.,       4.,       8.,      44.,   &
-          168.,     274.,     386.,     215.,      47.,       8./
-      ! 77-Comoe (Ivory Coast) 5N-4W (UNESCO, p56)          R1
-      ! Data from HLA: 5.1N -  4  W   1955/86 = 194m3/s       
-      ! --> 194.25m3/s
-      DATA inb(77)/1/
-      DATA (iirnf(jc,77),jc=1,jpcoef)/ 569, 19*0 /
-      DATA (ijrnf(jc,77),jc=1,jpcoef)/ 260, 19*0 /
-      DATA (zrnfm(jm,77),jm=1,jpmois)/   &
-           14.,       7.,       9.,      17.,      29.,      90.,   &
-          145.,     349.,     740.,     688.,     198.,      45./
-      ! 78-Sassandra (Ivory Coast) 5N-6W (UNESCO, p58)      R1
-      ! Data from HLA: 5N -  6  W   1953/86 = 449m3/s
-      ! --> 449.333344m3/s
-      DATA inb(78)/1/
-      DATA (iirnf(jc,78),jc=1,jpcoef)/ 564, 19*0 /
-      DATA (ijrnf(jc,78),jc=1,jpcoef)/ 260, 19*0 /
-      DATA (zrnfm(jm,78),jm=1,jpmois)/   &
-          144.,     116.,     119.,     143.,     156.,     334.,   &
-          419.,     692.,    1277.,    1190.,     553.,     249./
-      
-      DATA inb(79)/0/
-      DATA inb(80)/0/
-      DATA inb(81)/0/
-      DATA inb(82)/0/
-      DATA inb(83)/0/
-      DATA inb(84)/0/
-      DATA inb(85)/0/
-      DATA inb(86)/0/
-      DATA inb(87)/0/
-      DATA inb(88)/0/
-      DATA inb(89)/0/
-      DATA inb(90)/0/
-
-      ! Total run-offs (VI)=63977.m3/s
-      ! Total run-offs (I+II+III+IV+V+VI)=560902.688m3/s
-
-      ! VII OCEANIA
-      ! ------------
-      ! 91-Murray (Austrlia) 35S34 138E54                   R5
-      ! Old=22 km3/year=697m3/s
-      ! --> 696.999939m3/s
-      DATA inb(91)/3/
-      DATA (iirnf(jc,91),jc=1,jpcoef)/ 129, 130, 130, 17*0 /
-      DATA (ijrnf(jc,91),jc=1,jpcoef)/ 177, 176, 177, 17*0 /
-      DATA (zrnfm(jm,91),jm=1,jpmois)/12*697./
-      ! 92-Haast (N. Zeland) 43S50 169E02                   R5
-      ! Old=6 km3/year=190m3/s
-      ! --> 189.999985m3/s
-      DATA inb(92)/1/
-      DATA (iirnf(jc,92),jc=1,jpcoef)/ 194, 19*0 /
-      DATA (ijrnf(jc,92),jc=1,jpcoef)/ 153, 19*0 /
-      DATA (zrnfm(jm,92),jm=1,jpmois)/12*190./
-      ! 93-Fly (N. Guinea) 8S30 143E05                      R5
-      ! Old=77 km3/year=2440m3/s
-      ! --> 2440.m3/s
-      DATA inb(93)/1/
-      DATA (iirnf(jc,93),jc=1,jpcoef)/ 143, 19*0 /
-      DATA (ijrnf(jc,93),jc=1,jpcoef)/ 233, 19*0 /
-      DATA (zrnfm(jm,93),jm=1,jpmois)/12*2440./
-      ! 94-Purari (N. Guinea) 7S25 145E05                   R4
-      ! Old=77 km3/year=2440m3/s; UNESCO(77+83,p595)=2727m3/s
-      ! --> 2727.33325m3/s
-      DATA inb(94)/1/
-      DATA (iirnf(jc,94),jc=1,jpcoef)/ 146, 19*0 /
-      DATA (ijrnf(jc,94),jc=1,jpcoef)/ 234, 19*0 /
-      DATA (zrnfm(jm,94),jm=1,jpmois)/   &
-          1972.,    3067.,    2332.,    2981.,    2693.,    2876.,   &
-          2867.,    3477.,    2195.,    2967.,    2674.,    2627./
-      ! 95-Choshui ch''i (taiwan) 24N03 120E23              R2
-      ! Old=6 km3/year=190m3/s; UNESCO(65-68, p374)=120m3/s; ratio=1.58 
-      ! --> 0.m3/s
-      DATA inb(95)/0/
-      DATA (iirnf(jc,95),jc=1,jpcoef)/  20*0 /
-      DATA (ijrnf(jc,95),jc=1,jpcoef)/  20*0 /
-      DATA (zrnfm(jm,95),jm=1,jpmois)/   &
-           65.,      63.,      83.,     100.,     177.,     535.,   &
-          332.,     313.,     236.,     193.,     101.,      84./
-      ! 96-Kaoping (Taiwan)                                 R5  
-      ! Old=9 km3/year=285m3/s
-      ! --> 0.m3/s
-      DATA inb(96)/0/
-      DATA (iirnf(jc,96),jc=1,jpcoef)/ 20*0 /
-      DATA (ijrnf(jc,96),jc=1,jpcoef)/ 20*0 /
-      DATA (zrnfm(jm,96),jm=1,jpmois)/12*285./
-      ! 97-Tsengwen (Taiwan)                                R5
-      ! Old=2 km3/year=63m3/s
-      ! --> 0.m3/s
-      DATA inb(97)/0/
-      DATA (iirnf(jc,97),jc=1,jpcoef)/ 20*0 /
-      DATA (ijrnf(jc,97),jc=1,jpcoef)/ 20*0 /
-      DATA (zrnfm(jm,97),jm=1,jpmois)/12*63./
-      ! 98-Hualien (Taiwan)                                 R5
-      ! Old=4 km3/year=127m3/s
-      ! --> 0.m3/s
-      DATA inb(98)/0/
-      DATA (iirnf(jc,98),jc=1,jpcoef)/ 20*0 /
-      DATA (ijrnf(jc,98),jc=1,jpcoef)/ 20*0 /
-      DATA (zrnfm(jm,98),jm=1,jpmois)/12*127./
-      ! 99-Peinan (Taiwan) 22N46 121E10                     R5
-      ! Old=4 km3/year=127m3/s
-      ! --> 127.000008m3/s
-      DATA inb(99)/1/
-      DATA (iirnf(jc,99),jc=1,jpcoef)/  98, 19*0 /
-      DATA (ijrnf(jc,99),jc=1,jpcoef)/ 288, 19*0 /
-      DATA (zrnfm(jm,99),jm=1,jpmois)/12*127./
-      ! 100-Hsiukuluan (Taiwan)                             R5
-      ! Old=4 km3/year=127m3/s
-      ! --> 0.m3/s
-      DATA inb(100)/0/
-      DATA (iirnf(jc,100),jc=1,jpcoef)/ 20*0 /
-      DATA (ijrnf(jc,100),jc=1,jpcoef)/ 20*0 /
-      DATA (zrnfm(jm,100),jm=1,jpmois)/12*127./
-
-
-
-
-
-
-
-       ! I ASIA
-       ! ------
-
-       ! 1-Indus (Pakistan)
-       ! --> 3949m3/s
-       DATA inb(101)/4/
-       DATA (iirnf(jc,101),jc=1,jpcoef)/ 711, 711, 712, 713, 16*0 /
-       DATA (ijrnf(jc,101),jc=1,jpcoef)/ 297, 298, 297, 297, 16*0 /
-       DATA (zrnfm(jm,101),jm=1,jpmois)/  &
-     &      1691.,    2121.,    2083.,    2940.,    3738.,    4198.,  &
-     &      5334.,    8247.,    7833.,    4908.,    2681.,    1614./
-       ! 2-Sabarmati+Mahi+Narmada (India)
-       ! --> 3383m3/s
-       DATA inb(102)/1/
-!CT bug point terre       DATA (iirnf(jc,102),jc=1,jpcoef)/ 720, 19*0 /
-!CT bug point terre       DATA (ijrnf(jc,102),jc=1,jpcoef)/ 292, 19*0 /
-       DATA (iirnf(jc,102),jc=1,jpcoef)/ 719, 19*0 /
-       DATA (ijrnf(jc,102),jc=1,jpcoef)/ 291, 19*0 /
-       DATA (zrnfm(jm,102),jm=1,jpmois)/  &
-     &      1343.,     110.,       0.,       0.,      13.,    3051.,  &
-     &      9453.,   11655.,    7619.,    3547.,    2081.,    1727./
-       ! 3-
-       ! --> 182m3/s
-       DATA inb(103)/1/
-!CT bug point terre       DATA (iirnf(jc,103),jc=1,jpcoef)/ 721, 19*0 /
-!CT bug point terre       DATA (ijrnf(jc,103),jc=1,jpcoef)/ 286, 19*0 /
-       DATA (iirnf(jc,103),jc=1,jpcoef)/ 720, 19*0 /
-       DATA (ijrnf(jc,103),jc=1,jpcoef)/ 286, 19*0 /
-       DATA (zrnfm(jm,103),jm=1,jpmois)/  &
-     &      0.,       0.,       0.,       0.,     225.,     669.,  &
-     &      526.,     327.,     187.,     100.,      62.,      92./
-       ! 4-Mandovi + Zuari + Kalinadi
-       ! --> 347m3/s
-       DATA inb(104)/1/
-       DATA (iirnf(jc,104),jc=1,jpcoef)/ 2, 19*0 /
-       DATA (ijrnf(jc,104),jc=1,jpcoef)/ 282, 19*0 /
-       DATA (zrnfm(jm,104),jm=1,jpmois)/  &
-     &      7.,       0.,       0.,       0.,     428.,    1339.,  &
-     &      966.,     601.,     350.,     190.,     117.,     166./
-       ! 5-                    R2
-       ! --> 190m3/s
-       DATA inb(105)/1/
-       DATA (iirnf(jc,105),jc=1,jpcoef)/ 4, 19*0 /
-       DATA (ijrnf(jc,105),jc=1,jpcoef)/ 279, 19*0 /
-       DATA (zrnfm(jm,105),jm=1,jpmois)/  &
-     &      0.,       0.,       0.,       0.,     363.,     710.,   &
-     &      484.,     315.,     172.,      95.,      59.,      87./
-       ! 6-
-       ! --> 1101m3/s
-       DATA inb(106)/1/
-       DATA (iirnf(jc,106),jc=1,jpcoef)/ 4, 19*0 /
-       DATA (ijrnf(jc,106),jc=1,jpcoef)/ 277, 19*0 /
-       DATA (zrnfm(jm,106),jm=1,jpmois)/   &
-     &      79.,       0.,       0.,       0.,    1457.,    3639.,   &
-     &      3084.,    2031.,    1282.,     696.,     414.,     535. /
-       ! 7-
-       ! --> 948m3/s
-       DATA inb(107)/1/
-       DATA (iirnf(jc,107),jc=1,jpcoef)/ 5, 19*0 /
-       DATA (ijrnf(jc,107),jc=1,jpcoef)/ 275, 19*0 /
-       DATA (zrnfm(jm,107),jm=1,jpmois)/  &
-     &      72.,       0.,       0.,       0.,    1151.,    3038.,     & 
-     &      2583.,    1603.,    1216.,     749.,     426.,     539./
-
-       ! 8-
-       ! --> 551m3/s
-       DATA inb(108)/1/
-       DATA (iirnf(jc,108),jc=1,jpcoef)/ 7, 19*0 /
-       DATA (ijrnf(jc,108),jc=1,jpcoef)/ 271, 19*0 /
-       DATA (zrnfm(jm,108),jm=1,jpmois)/  &
-     &      45.,       0.,       2.,       7.,     680.,    1621.,  &
-     &      1234.,     744.,     784.,     696.,     359.,     442. /
-       ! 9-
-       ! --> 103m3/s
-       DATA inb(109)/1/
-       DATA (iirnf(jc,109),jc=1,jpcoef)/ 8, 19*0 /
-       DATA (ijrnf(jc,109),jc=1,jpcoef)/ 269, 19*0 /
-       DATA (zrnfm(jm,109),jm=1,jpmois)/  &
-     &      11.,       0.,       0.,       0.,     160.,     250.,   &
-     &      200.,     145.,     168.,     146.,      74.,      89. /
-       ! 10-
-       ! --> 99m3/s
-       DATA inb(110)/1/
-       DATA (iirnf(jc,110),jc=1,jpcoef)/12, 19*0 /
-       DATA (ijrnf(jc,110),jc=1,jpcoef)/267, 19*0 /
-       DATA (zrnfm(jm,110),jm=1,jpmois)/  &
-     &      50.,       1.,       1.,       0.,      25.,      69.,  &
-     &      51.,      32.,     176.,     329.,     233.,     222./
-       ! 11- Kaveri (India)
-       ! --> 173m3/s
-       DATA inb(111)/1/
-       DATA (iirnf(jc,111),jc=1,jpcoef)/16, 19*0 /
-       DATA (ijrnf(jc,111),jc=1,jpcoef)/271, 19*0 /
-       DATA (zrnfm(jm,111),jm=1,jpmois)/  &
-     &       97.,       2.,       1.,       2.,      19.,      67.,   &
-     &      108.,     134.,     218.,     515.,     498.,     416. /
-       ! 12-
-       ! --> 116m3/s
-       DATA inb(112)/1/
-       DATA (iirnf(jc,112),jc=1,jpcoef)/16, 19*0 /
-       DATA (ijrnf(jc,112),jc=1,jpcoef)/279, 19*0 /
-       DATA (zrnfm(jm,112),jm=1,jpmois)/  &
-     &       43.,       0.,       0.,       0.,       0.,       0.,   &
-     &        0.,       0.,      99.,     568.,     343.,     342. /
-       ! 13- Krishna
-       ! --> 2864m3/s
-       DATA inb(113)/1/
-       DATA (iirnf(jc,113),jc=1,jpcoef)/ 16, 19*0 /
-       DATA (ijrnf(jc,113),jc=1,jpcoef)/ 281, 19*0 /
-       DATA (zrnfm(jm,113),jm=1,jpmois)/   &
-     &       1368.,     174.,       0.,       0.,     166.,    2362.,  &
-     &       9501.,    8639.,    6053.,    3173.,    1763.,    1175. /
-       ! 14- (SriLanka)
-       ! --> 371m3/s
-       DATA inb(114)/1/
-       DATA (iirnf(jc,114),jc=1,jpcoef)/ 16, 19*0 /
-       DATA (ijrnf(jc,114),jc=1,jpcoef)/ 261, 19*0 /
-       DATA (zrnfm(jm,114),jm=1,jpmois)/  &
-     &      171.,     144.,     292.,     469.,     414.,     309.,  &
-     &      269.,     244.,     526.,     560.,     476.,     588.  /
-       ! 15- (SriLanka)
-       ! --> 305m3/s
-       DATA inb(115)/1/
-       DATA (iirnf(jc,115),jc=1,jpcoef)/18, 19*0 /
-       DATA (ijrnf(jc,115),jc=1,jpcoef)/262, 19*0 /
-       DATA (zrnfm(jm,115),jm=1,jpmois)/  &
-     &      116.,     111.,     225.,     379.,     366.,     300.,  &
-     &      256.,     239.,     385.,     437.,     388.,     467. /
-       ! 16- (SriLanka)
-       ! --> 722m3/s
-       DATA inb(116)/1/
-       DATA (iirnf(jc,116),jc=1,jpcoef)/20, 19*0 /
-       DATA (ijrnf(jc,116),jc=1,jpcoef)/265, 19*0 /
-       DATA (zrnfm(jm,116),jm=1,jpmois)/  &
-     &      711.,     347.,     477.,     556.,     430.,     302.,  &
-     &      228.,     200.,     629.,    1156.,    1579.,    2055.   /
-       ! 17- (SriLanka)
-       ! --> 188m3/s
-       DATA inb(117)/1/
-       DATA (iirnf(jc,117),jc=1,jpcoef)/19, 19*0 /
-       DATA (ijrnf(jc,117),jc=1,jpcoef)/267, 19*0 /
-       DATA (zrnfm(jm,117),jm=1,jpmois)/  &
-     &      189.,      50.,      27.,      17.,      10.,       6.,  &
-     &       4.,       2.,       1.,     236.,     829.,     894. /
-       ! 18- Godavari (India) 17N00  81E45
-       ! --> 2709m3/s
-       DATA inb(118)/1/
-       DATA (iirnf(jc,118),jc=1,jpcoef)/ 19, 19*0 /
-       DATA (ijrnf(jc,118),jc=1,jpcoef)/283, 19*0 /
-       DATA (zrnfm(jm,118),jm=1,jpmois)/  &
-     &      1151.,      57.,       0.,       0.,      38.,    1490.,  &
-     &      6714.,    9065.,    7208.,    3294.,    1874.,    1629. /
-       ! 19-
-       ! --> 116m3/s
-       DATA inb(119)/1/
-       DATA (iirnf(jc,119),jc=1,jpcoef)/23, 19*0 /
-       DATA (ijrnf(jc,119),jc=1,jpcoef)/286, 19*0 /
-       DATA (zrnfm(jm,119),jm=1,jpmois)/  &
-     &         17.,       0.,       0.,       0.,       0.,     103.,  &
-     &        259.,     341.,     330.,     155.,      87.,     110. /
-       ! 20- Mahanadi (India)
-       ! --> 2390m3/s
-       DATA inb(120)/1/
-       DATA (iirnf(jc,120),jc=1,jpcoef)/28, 19*0 /
-       DATA (ijrnf(jc,120),jc=1,jpcoef)/290, 19*0 /
-       DATA (zrnfm(jm,120),jm=1,jpmois)/  &
-     &       809.,       5.,       0.,       0.,      60.,    1401.,  &
-     &      6652.,    8828.,    5410.,    2561.,    1503.,    1452. /
-       ! 21-
-       ! --> 247m3/s
-       DATA inb(121)/1/
-       DATA (iirnf(jc,121),jc=1,jpcoef)/30, 19*0 /
-       DATA (ijrnf(jc,121),jc=1,jpcoef)/294, 19*0 /
-       DATA (zrnfm(jm,121),jm=1,jpmois)/  &
-     &         15.,       0.,       0.,       0.,      73.,     367.,   &
-     &        744.,     722.,     491.,     227.,     137.,     188. /
-       ! 22-
-       ! --> 167m3/s
-       DATA inb(122)/1/
-       DATA (iirnf(jc,122),jc=1,jpcoef)/31, 19*0 /
-       DATA (ijrnf(jc,122),jc=1,jpcoef)/294, 19*0 /
-       DATA (zrnfm(jm,122),jm=1,jpmois)/  &
-     &         16.,       0.,       0.,       0.,      95.,     337.,  &
-     &        488.,     453.,     285.,     143.,      84.,     109. /
-       ! 23-
-       ! --> 206m3/s
-       DATA inb(123)/1/
-       DATA (iirnf(jc,123),jc=1,jpcoef)/32, 19*0 /
-       DATA (ijrnf(jc,123),jc=1,jpcoef)/294, 19*0 /
-       DATA (zrnfm(jm,123),jm=1,jpmois)/  &
-     &         12.,       0.,       0.,       0.,      35.,     400.,  &
-     &        652.,     616.,     347.,     172.,     105.,     144. /
-       ! 24- Ganga + Brahmaputra (Bangladesh) 22N00  91E00
-       ! --> 42436m3/s
-       DATA inb(124)/9/
-       DATA (iirnf(jc,124),jc=1,jpcoef)/ 37, 38, 39, 32, 33, 34, 35  &
-     &      , 36, 37, 11*0 /
-       DATA (ijrnf(jc,124),jc=1,jpcoef)/ 295,295,295,294,294,294,294  &
-     &      , 294, 294, 11*0 /
-       DATA (zrnfm(jm,124),jm=1,jpmois)/  &
-     &        14293.,   4223.,   3004.,  11743.,  32210.,  56899.,  &
-     &        82263.,  95069.,  77997.,  44194.,  23724.,  17272. /
-
-       DATA inb(125)/0/
-      DATA inb(126)/0/
-
-       ! 27- Kuladan
-       ! --> 1541m3/s
-       DATA inb(127)/1/
-       DATA (iirnf(jc,127),jc=1,jpcoef)/42, 19*0 /
-       DATA (ijrnf(jc,127),jc=1,jpcoef)/290, 19*0 /
-       DATA (zrnfm(jm,127),jm=1,jpmois)/  &
-     &        140.,       0.,       0.,      27.,    2280.,    4133.,  &
-     &       4306.,    3049.,    1993.,    1115.,     646.,     812. /
-       ! 28-
-       ! --> 618m3/s
-       DATA inb(128)/1/
-       DATA (iirnf(jc,128),jc=1,jpcoef)/44, 19*0 /
-       DATA (ijrnf(jc,128),jc=1,jpcoef)/287, 19*0 /
-       DATA (zrnfm(jm,128),jm=1,jpmois)/  &
-     &        0.,       0.,       0.,      73.,    1139.,    1755.,  &
-     &      1725.,    1139.,     677.,     362.,     221.,     326. /
-       ! 29-
-       ! --> 158m3/s
-       DATA inb(129)/1/
-       DATA (iirnf(jc,129),jc=1,jpcoef)/44, 19*0 /
-       DATA (ijrnf(jc,129),jc=1,jpcoef)/284, 19*0 /
-       DATA (zrnfm(jm,129),jm=1,jpmois)/  &
-     &        0.,       0.,       0.,       0.,     248.,     419.,  &
-     &      445.,     334.,     201.,     103.,      63.,      93. /
-       ! 30-
-       ! --> 136m3/s
-       DATA inb(130)/1/
-       DATA (iirnf(jc,130),jc=1,jpcoef)/32, 19*0 /
-       DATA (ijrnf(jc,130),jc=1,jpcoef)/277, 19*0 /
-       DATA (zrnfm(jm,130),jm=1,jpmois)/  &
-     &         0.,       0.,       0.,       0.,     223.,     261.,  &
-     &       269.,     315.,     223.,     147.,      86.,     116. /
-       ! 31-
-       ! --> 142m3/s
-       DATA inb(131)/1/
-       DATA (iirnf(jc,131),jc=1,jpcoef)/42, 19*0 /
-       DATA (ijrnf(jc,131),jc=1,jpcoef)/275, 19*0 /
-       DATA (zrnfm(jm,131),jm=1,jpmois)/  &
-     &         0.,       0.,       0.,       0.,     241.,     260.,  &
-     &       270.,     325.,     231.,     158.,      96.,     126.  /
-       ! 32- Irrawady (Myanmar) 15N50  95E06
-       ! --> 16751m3/s
-       DATA inb(132)/3/
-       DATA (iirnf(jc,132),jc=1,jpcoef)/  45,  46, 47, 17*0 /
-       DATA (ijrnf(jc,132),jc=1,jpcoef)/ 282, 281, 282, 17*0 /
-       DATA (zrnfm(jm,132),jm=1,jpmois)/  &
-     &       7174.,     2281.,     80.,    687.,    8133.,   23980.,  &
-     &       38452.,   41442.,  34497.,  23150.,   12529.,    8629. /
-
-       DATA inb(133)/0/
-       DATA inb(134)/0/
-       DATA inb(135)/0/
-       DATA inb(136)/0/
-
-       ! 37- Tenasserim (Myanmar)
-       ! --> 1369m3/s
-       DATA inb(137)/1/
-       DATA (iirnf(jc,137),jc=1,jpcoef)/51, 19*0 /
-       DATA (ijrnf(jc,137),jc=1,jpcoef)/279, 19*0 /
-       DATA (zrnfm(jm,137),jm=1,jpmois)/  &
-     &        65.,       0.,       0.,     540.,    2127.,    3260.,  &
-     &      3540.,    2966.,    1747.,     900.,     544.,     741. /
-       ! 38-
-       ! --> 413m3/s
-       DATA inb(138)/1/
-       DATA (iirnf(jc,138),jc=1,jpcoef)/52, 19*0 /
-       DATA (ijrnf(jc,138),jc=1,jpcoef)/277, 19*0 /
-       DATA (zrnfm(jm,138),jm=1,jpmois)/  &
-     &          0.,       0.,       0.,     202.,     773.,    1000.,  &
-     &       1032.,     831.,     492.,     248.,     153.,     226. /
-       ! 39-
-       ! --> 810m3/s
-       DATA inb(139)/1/
-       DATA (iirnf(jc,139),jc=1,jpcoef)/52, 19*0 /
-       DATA (ijrnf(jc,139),jc=1,jpcoef)/275, 19*0 /
-       DATA (zrnfm(jm,139),jm=1,jpmois)/  &
-     &         59.,       0.,       0.,      22.,    1104.,    1857.,  &
-     &       2069.,    1858.,    1252.,     637.,     373.,     493. /
-       ! 40-
-       ! --> 896m3/s
-       DATA inb(140)/1/
-       DATA (iirnf(jc,140),jc=1,jpcoef)/52, 19*0 /
-       DATA (ijrnf(jc,140),jc=1,jpcoef)/273, 19*0 /
-       DATA (zrnfm(jm,140),jm=1,jpmois)/  &
-     &         94.,       0.,       0.,     196.,    1145.,    1716.,  &
-     &       1936.,    1899.,    1519.,    1000.,     562.,     688. /
-       ! 41-
-       ! --> 559m3/s
-       DATA inb(141)/1/
-       DATA (iirnf(jc,141),jc=1,jpcoef)/52, 19*0 /
-       DATA (ijrnf(jc,141),jc=1,jpcoef)/271, 19*0 /
-       DATA (zrnfm(jm,141),jm=1,jpmois)/  &
-     &       38.,       0.,       0.,     227.,     805.,    1051.,  &
-     &      1105.,    1236.,     941.,     575.,     321.,     416. /
-       ! 42-
-       ! --> 1070m3/s
-       DATA inb(142)/1/
-       DATA (iirnf(jc,142),jc=1,jpcoef)/53, 19*0 /
-       DATA (ijrnf(jc,142),jc=1,jpcoef)/265, 19*0 /
-       DATA (zrnfm(jm,142),jm=1,jpmois)/  &
-     &       356.,      18.,      14.,     160.,     745.,    1123.,  &
-     &      1347.,    1732.,    2076.,    2057.,    1574.,    1639. /
-       ! 43-
-       ! --> 513m3/s
-       DATA inb(143)/1/
-       DATA (iirnf(jc,143),jc=1,jpcoef)/54, 19*0 /
-       DATA (ijrnf(jc,143),jc=1,jpcoef)/263, 19*0 /
-       DATA (zrnfm(jm,143),jm=1,jpmois)/  &
-     &       231.,      12.,     121.,     236.,     216.,     239.,  &
-     &       351.,     517.,    1008.,    1324.,     975.,     929.  /
-       ! 44-
-       ! --> 881m3/s
-       DATA inb(144)/1/
-       DATA (iirnf(jc,144),jc=1,jpcoef)/56, 19*0 /
-       DATA (ijrnf(jc,144),jc=1,jpcoef)/259, 19*0 /
-       DATA (zrnfm(jm,144),jm=1,jpmois)/  &
-     &      558.,     632.,     980.,     902.,     575.,     421.,  &
-     &      432.,     731.,    1276.,    1405.,    1188.,    1480. /
-       ! 45-
-       ! --> 905m3/s
-       DATA inb(145)/1/
-       DATA (iirnf(jc,145),jc=1,jpcoef)/52, 19*0 /
-       DATA (ijrnf(jc,145),jc=1,jpcoef)/259, 19*0 /
-       DATA (zrnfm(jm,145),jm=1,jpmois)/  &
-     &      604.,     552.,     781.,     853.,     627.,     488.,  &
-     &      503.,     781.,    1245.,    1356.,    1331.,    1747. /
-       ! 46-
-       ! --> 136m3/s
-       DATA inb(146)/1/
-       DATA (iirnf(jc,146),jc=1,jpcoef)/50, 19*0 /
-       DATA (ijrnf(jc,146),jc=1,jpcoef)/261, 19*0 /
-       DATA (zrnfm(jm,146),jm=1,jpmois)/  &
-     &      103.,     111.,     126.,     112.,      66.,      40.,  &
-     &       35.,      51.,     126.,     227.,     279.,     362.  /
-       ! 47-
-       ! --> 222m3/s
-       DATA inb(147)/1/
-       DATA (iirnf(jc,147),jc=1,jpcoef)/47, 19*0 /
-       DATA (ijrnf(jc,147),jc=1,jpcoef)/258, 19*0 /
-       DATA (zrnfm(jm,147),jm=1,jpmois)/  &
-     &      176.,     160.,     209.,     197.,     130.,      99.,  &
-     &       90.,     151.,     240.,     342.,     373.,     502./
-       ! 48-
-       ! --> 326m3/s
-       DATA inb(148)/1/
-       DATA (iirnf(jc,148),jc=1,jpcoef)/49, 19*0 /
-       DATA (ijrnf(jc,148),jc=1,jpcoef)/257, 19*0 /
-       DATA (zrnfm(jm,148),jm=1,jpmois)/  &
-     &       226.,     278.,     363.,     301.,     181.,     130.,  &
-     &       131.,     194.,     407.,     505.,     519.,     689. /
-       ! 49-
-       ! --> 203m3/s
-       DATA inb(149)/1/
-       DATA (iirnf(jc,149),jc=1,jpcoef)/50, 19*0 /
-       DATA (ijrnf(jc,149),jc=1,jpcoef)/256, 19*0 /
-       DATA (zrnfm(jm,149),jm=1,jpmois)/  &
-     &      120.,     218.,     233.,     189.,     134.,     128.,  &
-     &      115.,     178.,     244.,     272.,     247.,     360. /
-       ! 50-
-       ! --> 153m3/s
-       DATA inb(150)/1/
-       DATA (iirnf(jc,150),jc=1,jpcoef)/50, 19*0 /
-       DATA (ijrnf(jc,150),jc=1,jpcoef)/255, 19*0 /
-       DATA (zrnfm(jm,150),jm=1,jpmois)/  &
-     &         97.,     108.,     191.,     141.,      97.,      82.,  &
-     &        107.,     118.,     188.,     230.,     208.,     278. /
-       ! 51-
-       ! --> 196m3/s
-       DATA inb(151)/1/
-       DATA (iirnf(jc,151),jc=1,jpcoef)/53, 19*0 /
-       DATA (ijrnf(jc,151),jc=1,jpcoef)/253, 19*0 /
-       DATA (zrnfm(jm,151),jm=1,jpmois)/  &
-     &        74.,     125.,     165.,     123.,     117.,     121.,  &
-     &       160.,     196.,     272.,     322.,     296.,     389. /
-       ! 52-
-       ! --> 166m3/s
-       DATA inb(152)/1/
-       DATA (iirnf(jc,152),jc=1,jpcoef)/53, 19*0 /
-       DATA (ijrnf(jc,152),jc=1,jpcoef)/252, 19*0 /
-       DATA (zrnfm(jm,152),jm=1,jpmois)/  &
-     &      136.,     140.,     192.,     169.,     108.,      80.,  &
-     &      75.,     145.,     200.,     219.,     221.,     307. /
-       ! 53-Tigris+Euphrates (Irak) 31N00  47E25
-       ! --> 4762m3/s
-       DATA inb(153)/3/
-       DATA (iirnf(jc,153),jc=1,jpcoef)/ 673, 673, 674, 17*0 /
-       DATA (ijrnf(jc,153),jc=1,jpcoef)/ 312, 313, 313, 17*0 /
-       DATA (zrnfm(jm,153),jm=1,jpmois)/  &
-     &       6056.,    7229.,    8377.,    8505.,    7966.,    6217.,  &
-     &       3396.,    1898.,    1147.,     869.,    1714.,    3779. /
-
-       DATA (inb(jr),jr=154,200)/47*0/
-
-
-      !  Total run-offs(VII)=6181.33301m3/s
-      !  Total run-offs(I+II+III+IV+V+VI+VII)= 567084. m3/s
-      !                                      = 1.588668056E-6 kg/m2/s
-      !                                      = 17883.5586 km3/Year
-      !                                      = 0.567083955 Sverdrup
-
-      !!----------------------------------------------------------------------
-      !!  OPA 8.5, LODYC-IPSL (2002)
-      !!----------------------------------------------------------------------
-
-      IF( kt == nit000 ) THEN
-         IF(lwp) WRITE(numout,*)
-         IF(lwp) WRITE(numout,*) 'flxrnf : runoff + orca05'
-         IF(lwp) WRITE(numout,*) '~~~~~~'
-         IF(lwp) WRITE(numout,*) '   nrunoff = ', nrunoff
-         IF(lwp) WRITE(numout,*)
-      ENDIF
-
-      ! 0. initialization and  indexes  modification for zoom version
-      ! --------------------------------------------------------------
-      
-      
-      ! conversion factor m3/s --> kg/m2/s *e1t*e2t
-
-      zconv = rauw
-
-      IF( kt == nit000 ) THEN 
-         ! Set runoffs and upstream coeff to zero
-         runoff (:,:) = 0.e0                                                                                  
-         upsrnfh(:,:) = 0.e0
-         upsrnfz(:)   = 0.e0
-         upsadv (:,:) = 0.e0
-         zcoefr(:,:,:) = 0.e0
-
-         IF( nrunoff == 0 ) THEN
-            IF(lwp) WRITE(numout,*)
-            IF(lwp) WRITE(numout,*) ' NO RUNOFF in this simulation '
-            IF(lwp) WRITE(numout,*)
-         ENDIF
-      ENDIF
-
-      ! 1. initialization of the runoff array
-      ! -------------------------------------
-
-      IF( nrunoff == 1 .OR. nrunoff == 2) THEN
-
-         iman  = jpmois
- !!! better but change the results      i15 = INT( 2*FLOAT( nday ) / ( FLOAT( nobis(nmonth) ) + 0.5 ) )
-         i15 = nday / 16
-         imois = nmonth + i15 - 1
-         IF( imois == 0) imois = iman
-         IF( nn_leap_year == 1 .AND. MOD( nyear, 4 ) == 0 ) THEN
-            ! Number of Day Between Data
-            idbd = nbiss(imois)
-         ELSEIF( nn_leap_year > 1 ) THEN
-            idbd = nn_leap_year
-         ELSE
-            idbd = nobis(imois)
-         END IF
-         ! Number of days between imois, 15 and the end of month
-         idmeom = idbd - 15
-         
-         ! 1.1 Calendar computation
-         
-         ! nrnf1 number of the first array record used in the simulation
-         ! nrnf2 number of the last  array record
-         
-         nrnf1 = imois
-         nrnf2 = nrnf1 + 1
-         nrnf1 = MOD( nrnf1, iman )
-         IF( nrnf1 == 0 ) nrnf1 = iman
-         nrnf2 = MOD( nrnf2, iman )
-         IF( nrnf2 == 0 ) nrnf2 = iman
-         
-         IF( kt == nit000 ) THEN
-            IF(lwp) WRITE(numout,*)
-            IF(lwp) WRITE(numout,*) 'flxrnf routine'
-            IF(lwp) WRITE(numout,*) 'first array record used nrnf1 ',nrnf1
-            IF(lwp) WRITE(numout,*) 'last  array record used nrnf2 ',nrnf2
-         ENDIF
-
-         zxy = FLOAT( nday + idmeom - idbd * i15 ) / idbd
-
-         DO jr = 1, jpriv-1
-            inbp = inb(jr)
-            IF( inbp /= 0 ) THEN
-               DO jl = 1, inbp
-                  DO jj = mj0(ijrnf(jl,jr)), mj1(ijrnf(jl,jr))
-                     DO ji = mi0(iirnf(jl,jr)), mi1(iirnf(jl,jr))
-                        runoff(ji,jj) = runoff(ji,jj)   & 
-                           &          + ( (1.-zxy) * zrnfm(nrnf1,jr)   &
-                           &          +       zxy  * zrnfm(nrnf2,jr)  ) / FLOAT(inbp)
-                     END DO
-                  END DO
-               END DO
-            ENDIF
-         END DO
-
-
-         ! 2. conversion in mm/day
-         ! -----------------------
-         !     C a u t i o n : upward water flux, runoff is negative
-
-         runoff(:,:) = -zconv * runoff(:,:) / ( e1t(:,:) * e2t(:,:) )
-         
-
-         ! 3. advectve coefficient
-         ! -----------------------
-
-         !  Creates the array coef that contains the coefficient to affect to
-         !  the upstream scheme. advection scheme will be:
-         !  coefr * upstream + (1- coefr) centered
-         !  coefr must be between 0 and 1.
-         !  3 values are used: maximum = zrup1, between zrup2, small zrup3.
-         ! 
-         IF( nrunoff == 2 .AND. kt == nit000 ) THEN
-
-            zrup1 = 0.5
-            zrup2 = 0.25
-            zrup3 = 0.125
-            DO jr = 1, jpriv-1
-               inbp = inb(jr)
-               IF( inbp /= 0) THEN
-                  DO jl = 1, inbp
-                     ! ji,jj
-                     DO jj = mj0(ijrnf(jl,jr)), mj1(ijrnf(jl,jr))
-                        DO ji = mi0(iirnf(jl,jr)), mi1(iirnf(jl,jr))
-                           DO jn = 1, 3
-                              zcoefr(ji,jj,jn) = zrup1
-                           END DO
-                           zcoefr(ji,jj,4) = zrup2
-                           zcoefr(ji,jj,5) = zrup3
-                        END DO
-                     END DO
-                     ! --> zrup2 (originally = 0.5)
-                     ! ji,jj+1
-                     DO jj = mj0(ijrnf(jl,jr) + 1), mj1(ijrnf(jl,jr) + 1)
-                        DO ji = mi0(iirnf(jl,jr)), mi1(iirnf(jl,jr))
-                           DO jn = 1, 3
-                              zcoefr(ji,jj,jn) = MAX( zrup2, zcoefr(ji,jj,jn) )
-                           END DO
-                           zcoefr(ji,jj,4) = MAX( zrup3, zcoefr(ji,jj,4) )
-                        END DO
-                     END DO
-                     ! ji,jj-1
-                     DO jj = mj0(ijrnf(jl,jr) - 1), mj1(ijrnf(jl,jr) - 1)
-                        DO ji = mi0(iirnf(jl,jr)), mi1(iirnf(jl,jr))
-                           DO jn = 1, 3
-                              zcoefr(ji,jj,jn) = MAX( zrup2, zcoefr(ji,jj,jn) )
-                           END DO
-                           zcoefr(ji,jj,4) = MAX( zrup3, zcoefr(ji,jj,4) )
-                        END DO
-                     END DO
-                     ! ji+1,jj-1
-                     DO jj = mj0(ijrnf(jl,jr) - 1), mj1(ijrnf(jl,jr) - 1)
-                        DO ji = mi0(iirnf(jl,jr) + 1), mi1(iirnf(jl,jr) + 1)
-                           DO jn = 1, 3
-                              zcoefr(ji,jj,jn) = MAX( zrup2, zcoefr(ji,jj,jn) )
-                           END DO
-                           zcoefr(ji,jj,4) = MAX( zrup3, zcoefr(ji,jj,4) )
-                        END DO
-                     END DO
-                     ! ji+1,jj
-                     DO jj = mj0(ijrnf(jl,jr)), mj1(ijrnf(jl,jr))
-                        DO ji = mi0(iirnf(jl,jr) + 1), mi1(iirnf(jl,jr) + 1)
-                           DO jn = 1, 3
-                              zcoefr(ji,jj,jn) = MAX( zrup2, zcoefr(ji,jj,jn) )
-                           END DO
-                           zcoefr(ji,jj,4) = MAX( zrup3, zcoefr(ji,jj,4) )
-                        END DO
-                     END DO
-                     ! ji+1,jj+1
-                     DO jj = mj0(ijrnf(jl,jr) + 1), mj1(ijrnf(jl,jr) + 1)
-                        DO ji = mi0(iirnf(jl,jr) + 1), mi1(iirnf(jl,jr) + 1)
-                           DO jn = 1, 3
-                              zcoefr(ji,jj,jn) = MAX( zrup2, zcoefr(ji,jj,jn) )
-                           END DO
-                           zcoefr(ji,jj,4) = MAX( zrup3, zcoefr(ji,jj,4) )
-                        END DO
-                     END DO
-                     ! ji-1,jj+1
-                     DO jj = mj0(ijrnf(jl,jr) + 1), mj1(ijrnf(jl,jr) + 1)
-                        DO ji = mi0(iirnf(jl,jr) - 1), mi1(iirnf(jl,jr) - 1)
-                           DO jn = 1, 3
-                              zcoefr(ji,jj,jn) = MAX( zrup2, zcoefr(ji,jj,jn) )
-                           END DO
-                           zcoefr(ji,jj,4) = MAX( zrup3, zcoefr(ji,jj,4) )
-                        END DO
-                     END DO
-                     ! ji-1,jj
-                     DO jj = mj0(ijrnf(jl,jr)), mj1(ijrnf(jl,jr))
-                        DO ji = mi0(iirnf(jl,jr) - 1), mi1(iirnf(jl,jr) - 1)
-                           DO jn = 1, 3
-                              zcoefr(ji,jj,jn) = MAX( zrup2, zcoefr(ji,jj,jn) )
-                           END DO
-                           zcoefr(ji,jj,4) = MAX( zrup3, zcoefr(ji,jj,4) )
-                        END DO
-                     END DO
-                     ! ji-1,jj-1
-                     DO jj = mj0(ijrnf(jl,jr) - 1), mj1(ijrnf(jl,jr) - 1)
-                        DO ji = mi0(iirnf(jl,jr) - 1), mi1(iirnf(jl,jr) - 1)
-                           DO jn = 1, 3
-                              zcoefr(ji,jj,jn) = MAX( zrup2, zcoefr(ji,jj,jn) )
-                           END DO
-                           zcoefr(ji,jj,4) = MAX( zrup3, zcoefr(ji,jj,4) )
-                        END DO
-                     END DO
-                     ! --> zrup3 (originally = 0.25)
-                     ! ji-1,jj+2
-                     DO jj = mj0(ijrnf(jl,jr) + 2), mj1(ijrnf(jl,jr) + 2)
-                        DO ji = mi0(iirnf(jl,jr) - 1), mi1(iirnf(jl,jr) - 1)
-                           DO jn = 1, 3
-                              zcoefr(ji,jj,jn) = MAX( zrup3, zcoefr(ji,jj,jn) )
-                           END DO
-                        END DO
-                     END DO
-                     ! ji-1,jj-2
-                     DO jj = mj0(ijrnf(jl,jr) - 2), mj1(ijrnf(jl,jr) - 2)
-                        DO ji = mi0(iirnf(jl,jr) - 1), mi1(iirnf(jl,jr) - 1)
-                           DO jn = 1, 3 
-                              zcoefr(ji,jj,jn) = MAX( zrup3, zcoefr(ji,jj,jn) )
-                           END DO
-                        END DO
-                     END DO
-                     ! ji,jj-2
-                     DO jj = mj0(ijrnf(jl,jr) - 2), mj1(ijrnf(jl,jr) - 2)
-                        DO ji = mi0(iirnf(jl,jr)    ), mi1(iirnf(jl,jr)    )
-                           DO jn = 1, 3 
-                              zcoefr(ji,jj,jn) = MAX( zrup3, zcoefr(ji,jj,jn) )
-                           END DO
-                        END DO
-                     END DO
-                     ! ji,jj+2
-                     DO jj = mj0(ijrnf(jl,jr) + 2), mj1(ijrnf(jl,jr) + 2)
-                        DO ji = mi0(iirnf(jl,jr)    ), mi1(iirnf(jl,jr)    )
-                           DO jn = 1, 3
-                              zcoefr(ji,jj,jn) = MAX( zrup3, zcoefr(ji,jj,jn) )
-                           END DO
-                        END DO
-                     END DO
-                     ! ji+1,jj-2
-                     DO jj = mj0(ijrnf(jl,jr) - 2), mj1(ijrnf(jl,jr) - 2)
-                        DO ji = mi0(iirnf(jl,jr) + 1), mi1(iirnf(jl,jr) + 1)
-                           DO jn = 1, 3 
-                              zcoefr(ji,jj,jn) = MAX( zrup3, zcoefr(ji,jj,jn) )
-                           END DO
-                        END DO
-                     END DO
-                     ! ji+1,jj+2
-                     DO jj = mj0(ijrnf(jl,jr) + 2), mj1(ijrnf(jl,jr) + 2)
-                        DO ji = mi0(iirnf(jl,jr) + 1), mi1(iirnf(jl,jr) + 1)
-                           DO jn = 1, 3 
-                              zcoefr(ji,jj,jn) = MAX( zrup3, zcoefr(ji,jj,jn) )
-                           END DO
-                        END DO
-                     END DO
-                     ! ji+2,jj-2
-                     DO jj = mj0(ijrnf(jl,jr) - 2), mj1(ijrnf(jl,jr) - 2)
-                        DO ji = mi0(iirnf(jl,jr) + 2), mi1(iirnf(jl,jr) + 2)
-                           DO jn = 1, 3 
-                              zcoefr(ji,jj,jn) = MAX( zrup3, zcoefr(ji,jj,jn) )
-                           END DO
-                        END DO
-                     END DO
-                     ! ji+2,jj+2
-                     DO jj = mj0(ijrnf(jl,jr) + 2), mj1(ijrnf(jl,jr) + 2)
-                        DO ji = mi0(iirnf(jl,jr) + 2), mi1(iirnf(jl,jr) + 2)
-                           DO jn = 1, 3  
-                              zcoefr(ji,jj,jn) = MAX( zrup3, zcoefr(ji,jj,jn) )
-                           END DO
-                        END DO
-                     END DO
-                     ! ji+2,jj-1
-                     DO jj = mj0(ijrnf(jl,jr) - 1), mj1(ijrnf(jl,jr) - 1)
-                        DO ji = mi0(iirnf(jl,jr) + 2), mi1(iirnf(jl,jr) + 2)
-                           DO jn = 1, 3
-                              zcoefr(ji,jj,jn) = MAX( zrup3, zcoefr(ji,jj,jn) )
-                           END DO
-                        END DO
-                     END DO
-                     ! ji+2,jj+1
-                     DO jj = mj0(ijrnf(jl,jr) + 1), mj1(ijrnf(jl,jr) + 1)
-                        DO ji = mi0(iirnf(jl,jr) + 2), mi1(iirnf(jl,jr) + 2)
-                           DO jn = 1, 3 
-                              zcoefr(ji,jj,jn) = MAX( zrup3, zcoefr(ji,jj,jn) )
-                           END DO
-                        END DO
-                     END DO
-                     ! ji+2,jj
-                     DO jj = mj0(ijrnf(jl,jr)    ), mj1(ijrnf(jl,jr)    )
-                        DO ji = mi0(iirnf(jl,jr) + 2), mi1(iirnf(jl,jr) + 2)
-                           DO jn = 1, 3
-                              zcoefr(ji,jj,jn) = MAX( zrup3, zcoefr(ji,jj,jn) )
-                           END DO
-                        END DO
-                     END DO
-                     ! ji-2,jj
-                     DO jj = mj0(ijrnf(jl,jr)    ), mj1(ijrnf(jl,jr)    )
-                        DO ji = mi0(iirnf(jl,jr) - 2), mi1(iirnf(jl,jr) - 2)
-                           DO jn = 1, 3        
-                              zcoefr(ji,jj,jn) = MAX( zrup3, zcoefr(ji,jj,jn) )
-                           END DO
-                        END DO
-                     END DO
-                     ! ji-2,jj-2
-                     DO jj = mj0(ijrnf(jl,jr) - 2), mj1(ijrnf(jl,jr) - 2)
-                        DO ji = mi0(iirnf(jl,jr) - 2), mi1(iirnf(jl,jr) - 2)
-                           DO jn = 1, 3 
-                              zcoefr(ji,jj,jn) = MAX( zrup3, zcoefr(ji,jj,jn) )
-                           END DO
-                        END DO
-                     END DO
-                     ! ji-2,jj+2
-                     DO jj = mj0(ijrnf(jl,jr) + 2), mj1(ijrnf(jl,jr) + 2)
-                        DO ji = mi0(iirnf(jl,jr) - 2), mi1(iirnf(jl,jr) - 2)
-                           DO jn = 1, 3
-                              zcoefr(ji,jj,jn) = MAX( zrup3, zcoefr(ji,jj,jn) )
-                           END DO
-                        END DO
-                     END DO
-                     ! ji-2,jj-1
-                     DO jj = mj0(ijrnf(jl,jr) - 1), mj1(ijrnf(jl,jr) - 1)
-                        DO ji = mi0(iirnf(jl,jr) - 2), mi1(iirnf(jl,jr) - 2)
-                           DO jn = 1, 3 
-                              zcoefr(ji,jj,jn) = MAX( zrup3, zcoefr(ji,jj,jn) )
-                           END DO
-                        END DO
-                     END DO
-                     ! ji-2,jj+1
-                     DO jj = mj0(ijrnf(jl,jr) + 1), mj1(ijrnf(jl,jr) + 1)
-                        DO ji = mi0(iirnf(jl,jr) - 2), mi1(iirnf(jl,jr) - 2)
-                           DO jn = 1, 3 
-                              zcoefr(ji,jj,jn) = MAX( zrup3, zcoefr(ji,jj,jn) )
-                           END DO
-                        END DO
-                     END DO
-                  END DO
-               ENDIF
-            END DO
-            upsrnfh(:,:) = zcoefr(:,:,1)
-            upsrnfz(:)   = 0.e0
-            upsrnfz(1)   = zrup1
-            upsrnfz(2)   = zrup1
-            upsrnfz(3)   = zrup1
-            upsrnfz(4)   = zrup2
-            upsrnfz(5)   = zrup3
-         ENDIF
-      ENDIF
-    
-#endif
-
-   END SUBROUTINE flx_rnf
Index: trunk/NEMO/OPA_SRC/SBC/sbcrnf_ORCA_R2.h90
===================================================================
--- trunk/NEMO/OPA_SRC/SBC/sbcrnf_ORCA_R2.h90	(revision 1122)
+++ 	(revision )
@@ -1,1340 +1,0 @@
-   !!----------------------------------------------------------------------
-   !!                    ***  flxrnf_ORCA_R2.h90  ***
-   !!----------------------------------------------------------------------
-
-   !!----------------------------------------------------------------------
-   !!                  OLD ORCA R2 runoff specification
-   !!                      kept here for memory
-   !!                         ( 20/04/2000 )
-   !!----------------------------------------------------------------------
-   !! History :         !  94-10  (G.Madec, M. Pontaud, M. Imbard)  Original code
-   !!                   !  97-03  (G.Madec)  time dependent version
-   !!                   !  98-06  (J.M. Molines)  exact computation of zxy 
-   !!                                      for months that are not 30 days
-   !!                   !  98-07  (M. Imbard)  ORCA and mpp option
-   !!                   !  99-08  (J.P. Boulanger H.L.Ayina)  New rivers and 
-   !!                                                   values given in m3/s 
-   !!                   !  00-04  (G. Madec, K. Roberts) add antarctica ice discharge.
-   !!----------------------------------------------------------------------
-   !! References : 
-   !!       J. D. Milliman and R. H. Meade, 1983 : world-wide delivery
-   !!          of river sediment to the oceans, journal of geology vol 91
-   !!          pp 1-21.
-   !!       G. L. Russell and J. R. Miller, 1990 : global river runoff
-   !!          calculated from a global atmospheric general circulation
-   !!          model, journal of hydrology, 117(1990), pp 241-254.
-   !!       F. Van Der Leeden, Troise F. L., Todd D. K. : the water
-   !!          encyclopedia, second edition, lewis publishers.
-   !!       J. W. Weatherly, J. E. Walsh : The effects of precipitation
-   !!          and river runoff in a coupled ice-ocean model of Arctic
-   !!          Climate dynamics 1996 12:785,798
-   !!       Jacobs et al. 1992. J. Glaciol. 38 (130) 375-387.
-   !!----------------------------------------------------------------------
-   !!  OPA 9.0 , LOCEAN-IPSL (2005) 
-   !! $ Id: $
-   !! Software governed by the CeCILL licence (modipsl/doc/NEMO_CeCILL.txt)
-   !!----------------------------------------------------------------------
-      !!----------------------------------------------------------------------
-      INTEGER :: kt
-      INTEGER, PARAMETER ::   &
-         jpmois =  12,        &  ! number of months
-         jpriv  = 120,        &  ! maximum number of river
-         jpcoef = 300            ! maximum of river mouth grid-points
-      INTEGER ::   &
-         ji, jj, jl, jm,      &  ! dummy loop indices
-         jn, jr, jc              !    "         "
-      INTEGER ::   &
-         i15, imois, iman,    &  ! temporary integers
-         inbp, idbd, idmeom,  &  !    "          "
-         iipt0, iipt1, ijpt0, &  !    "          "
-         ijpt1, ii, ij           !    "          "
-      INTEGER, DIMENSION(jpriv) ::   &
-         inb                     ! ???
-      INTEGER, DIMENSION(jpcoef,jpriv) ::   &
-         iirnf, ijrnf            ! ???
-
-      REAL(wp) ::   &
-         zconv, zxy,          &  ! temporary scalars
-         zrup1, zrup2, zrup3     !    "         "
-      REAL(wp), DIMENSION(jpmois,jpriv) ::   &
-         zrnfm                   ! ???
-      REAL(wp), DIMENSION(jpi, jpj, 5) ::   &
-         zcoefr                  ! ???
-
-      ! ====
-      ! DATA
-      ! ====
-
-      ! I NORTH AMERICA
-      ! ---------------
-
-      ! 1-Saint Lawrence (Canada) 49N30 67W00               R2
-      ! Old=447 km3/year=14165m3/s; CLIPPER=12430m3/s; ratio=1.14
-      ! --> 13975.249m3/s
-      DATA inb(1)/2/
-      DATA (iirnf(jc,1),jc=1,jpcoef)/ 110, 110, 298*0 /
-      DATA (ijrnf(jc,1),jc=1,jpcoef)/ 110, 111, 298*0 /
-      DATA (zrnfm(jm,1),jm=1,jpmois)/   &
-          11718.,  11971.,  12466.,  14980.,  19257.,  16913.,   &
-          14211.,  13119.,  13367.,  13644.,  13463.,  12594./
-      ! 2-Hudson River (USA) 41N00 77W00                    R1
-      ! Old=12 km3/year=380m3/s; UNESCO(64-83,p255)=408m3/s
-      ! --> 407.416687m3/s
-      DATA inb(2)/1/
-      DATA (iirnf(jc,2),jc=1,jpcoef)/ 106, 299*0 /
-      DATA (ijrnf(jc,2),jc=1,jpcoef)/ 104, 299*0 /
-      DATA (zrnfm(jm,2),jm=1,jpmois)/   &
-          364.,    404.,    648.,    873.,    593.,    301.,   &
-          205.,    171.,    205.,    290.,    379.,    456./
-      ! 3-Missisipi and Atchafalaya (USA) 29N 90W           R1 
-      ! Old=580 km3/year=18380m3/s; CLIPPER=17599m3/s
-      ! --> 17599.334m3/s
-      DATA inb(3)/2/
-      DATA (iirnf(jc,3),jc=1,jpcoef)/  97, 98, 298*0 /
-      DATA (ijrnf(jc,3),jc=1,jpcoef)/  97, 97, 298*0 /
-      DATA (zrnfm(jm,3),jm=1,jpmois)/   &
-          19260.,  20813.,  23026.,  28160.,  26510.,  19650.,   &
-          14181.,  10676.,   9578.,  10167.,  11699.,  17472./
-      ! 4-Brazos (USA) 29N58 95W53                          R1
-      ! Old=7 km3/year=222m3/s; UNESCO(65-84, p298)=200m3/s
-      ! --> 199.666656m3/s
-      DATA inb(4)/1/
-      DATA (iirnf(jc,4),jc=1,jpcoef)/  95, 299*0 /
-      DATA (ijrnf(jc,4),jc=1,jpcoef)/  97, 299*0 /
-      DATA (zrnfm(jm,4),jm=1,jpmois)/   &
-          165.,    224.,    200.,    266.,    457.,    327.,   &
-          143.,     69.,    102.,    132.,    160.,    151./
-      ! 5-Colorado (Mexico) 31N54 114W57                    R5
-      ! Old=20 km3/year=634m3/s
-      ! --> 634.m3/s
-      DATA inb(5)/1/
-      DATA (iirnf(jc,5),jc=1,jpcoef)/  88, 299*0 /
-      DATA (ijrnf(jc,5),jc=1,jpcoef)/  94, 299*0 /
-      DATA (zrnfm(jm,5),jm=1,jpmois)/12*634./ 
-      ! 6-Columbia (USA) 45N60 121W17                       R1
-      ! Old=251 km3/year=7954m3/s; UNESCO(65-84, p312)=5368m3/s
-      ! --> 5368.41602m3/s
-      DATA inb(6)/1/
-      DATA (iirnf(jc,6),jc=1,jpcoef)/  79, 299*0 /
-      DATA (ijrnf(jc,6),jc=1,jpcoef)/ 110, 299*0 /
-      DATA (zrnfm(jm,6),jm=1,jpmois)/   &
-          4752.,    5077.,    5633.,    6075.,    8422.,    9822.,   &
-          6350.,    4112.,    3273.,    3251.,    3548.,    4106./ 
-      ! 7-Fraser (USA)                                      R5
-      ! Old=112 km3/year=3550m3/s;not found in UNESCO
-      ! --> 3549.99976m3/s
-      DATA inb(7)/1/
-      DATA (iirnf(jc,7),jc=1,jpcoef)/  77, 299*0 /
-      DATA (ijrnf(jc,7),jc=1,jpcoef)/ 112, 299*0 /
-      DATA (zrnfm(jm,7),jm=1,jpmois)/12*3550./ 
-      ! 8-Yukon (USA) 63N00 165W00                          R3
-      ! Old=195 km3/year=6179m3/s; UNESCO(65-77+80-84, p321)=4245; ratio=1.46
-      ! --> 6017.24951m3/s
-      DATA inb(8)/1/
-      DATA (iirnf(jc,8),jc=1,jpcoef)/  58, 299*0 /
-      DATA (ijrnf(jc,8),jc=1,jpcoef)/ 122, 299*0 /
-      DATA (zrnfm(jm,8),jm=1,jpmois)/   &
-           1459.,    1213.,    1065.,    1100.,    9067.,   17675.,   &
-          10680.,   11271.,    8658.,    5453.,    2649.,    1917./ 
-     ! 9-Copper (USA) 61N47 144W45                          R2
-     ! Old=39 km3/year=1236m3/s; UNESCO(65-84, p317)=1052m3/s; ratio=1.18
-     ! --> 1236.25m3/s
-      DATA inb(9)/1/
-      DATA (iirnf(jc,9),jc=1,jpcoef)/  67, 299*0 /
-      DATA (ijrnf(jc,9),jc=1,jpcoef)/ 119, 299*0 /
-      DATA (zrnfm(jm,9),jm=1,jpmois)/   &
-           219.,    195.,    181.,    214.,    904.,   2532.,   &
-          4026.,   3627.,   1632.,    697.,    354.,    254./ 
-      ! 10-Susitna (USA)                                    R5
-      ! Old=40 km3/year=1268m3/s
-      ! --> 1268.m3/s
-      DATA inb(10)/1/
-      DATA (iirnf(jc,10),jc=1,jpcoef)/  65, 299*0 /
-      DATA (ijrnf(jc,10),jc=1,jpcoef)/ 118, 299*0 /
-      DATA (zrnfm(jm,10),jm=1,jpmois)/12*1268./ 
-      ! 11-McKenzie (Canada) 70N00 135W00                   R3
-      ! Old=306 km3/year=9697m3/s; UNESCO(67+72+7479+82-84, p168)=8343m3/s;
-      ! ratio=1.16
-      ! --> 9606.66797m3/s
-      DATA inb(11)/1/
-      DATA (iirnf(jc,11),jc=1,jpcoef)/  69, 299*0 /
-      DATA (ijrnf(jc,11),jc=1,jpcoef)/ 131, 299*0 /
-      DATA (zrnfm(jm,11),jm=1,jpmois)/   &
-           3644.,    3400.,    3258.,    3546.,   14188.,   21544.,   &
-          18954.,   15487.,   11284.,   10169.,    5901.,    3905./
-      ! 12-Kuksokwim (USA) 60N00 ~162W30                    R4
-      ! UNESCO(65-84, p318)=1112m3/s (61N87 158E10, Crooked Creek, inland)
-      ! --> 1111.99988m3/s
-      DATA inb(12)/1/
-      DATA (iirnf(jc,12),jc=1,jpcoef)/  60, 299*0 /
-      DATA (ijrnf(jc,12),jc=1,jpcoef)/ 118, 299*0 /
-      DATA (zrnfm(jm,12),jm=1,jpmois)/   &
-          335.,     308.,     295.,     407.,    2415.,    2027.,   &
-          1917.,   2049.,    1585.,    1072.,     545.,     389./     
-
-      DATA inb(13)/0/
-      DATA inb(14)/0/
-      DATA inb(15)/0/
-
-      ! Total run-offs(I)=60974.2539m3/s
-
-      ! II SOUTH AMERICA
-      ! -----------------
-      ! 16-Chira (Peru) 4s54 81w08                          R2
-      ! Old=5 km3/year=158m3/s; UNESCO (65-68, p357)=80.9; ratio=1.95
-      ! --> 158.000015m3/s
-      DATA inb(16)/1/
-      DATA (iirnf(jc,16),jc=1,jpcoef)/ 101, 299*0 /
-      DATA (ijrnf(jc,16),jc=1,jpcoef)/  66, 299*0 /
-      DATA (zrnfm(jm,16),jm=1,jpmois)/   &
-          70.,     105.,     564.,     582.,     230.,      94.,   &
-          87.,      47.,      35.,      32.,      29.,      21./
-      ! 17-Magdalena (Colombia) 11N06 74W51                 R2
-      ! Old=237 km3/year=7510m3/s; UNESCO(71-73+76-79, p348)=6973m3/s;
-      ! ratio=1.08
-      ! --> 7451.58301m3/s
-      DATA inb(17)/1/
-      DATA (iirnf(jc,17),jc=1,jpcoef)/ 104, 299*0 /
-      DATA (ijrnf(jc,17),jc=1,jpcoef)/  86, 299*0 /
-      DATA (zrnfm(jm,17),jm=1,jpmois)/   &
-          6598.,    4441.,    4087.,    5445.,    7380.,    8067.,   &
-          8509.,    7428.,    7631.,    9296.,   10588.,    9949./
-      ! 18-Orinoco (Venezuela) 8N37 62W15!                  R1
-      ! Old=1100 km3/year=34857m3/s; UNESCO(73-75, p364)=28736m3/s;
-      ! HLA=26371m3/s
-      ! Data from HLA: 8.1N - 63.3W   1925/89
-      ! --> 26371.416
-      DATA inb(18)/1/
-      DATA (iirnf(jc,18),jc=1,jpcoef)/ 112, 299*0 /
-      DATA (ijrnf(jc,18),jc=1,jpcoef)/  86, 299*0 /
-      DATA (zrnfm(jm,18),jm=1,jpmois)/   &
-          10670.,    7513.,    6344.,    6887.,   13086.,   26699.,   &
-          43938.,   59957.,   57435.,   39951.,   27410.,   16567./ 
-      ! 19-Amazon (Brazil) 0S05 50W00                       R1
-      ! Old=6300 km3/year=199634m3/s; CLIPPER=166192m3/s~UNESCO;
-      ! HLA=184806m3/s
-      ! UNESCO data are at Obidos which is far inland and therefore does not
-      ! fully represent the actual discharge
-      ! Data from HLA: 1.5N - 55.3W   1970/96
-      ! --> 184845.484m3/s
-      DATA inb(19)/2/
-      DATA (iirnf(jc,19),jc=1,jpcoef)/ 117, 117, 298*0 /
-      DATA (ijrnf(jc,19),jc=1,jpcoef)/  75,  76, 298*0 /
-      DATA (zrnfm(jm,19),jm=1,jpmois)/   &
-          152233., 176415., 199941., 225565., 244016., 244861.,   &
-          230711., 203128., 164531., 128486., 118214., 130045. /
-      ! 20-Sao Francisco (Brazil) 10S30 36W24               R1
-      ! Old=97 km3/year=3074m3/s; UNESCO(77-79, p338)=2645m3/s
-      ! --> 2645.3335m3/s
-      DATA inb(20)/1/
-      DATA (iirnf(jc,20),jc=1,jpcoef)/ 123, 299*0 /
-      DATA (ijrnf(jc,20),jc=1,jpcoef)/  59, 299*0 /
-      DATA (zrnfm(jm,20),jm=1,jpmois)/   &
-          3244.,   4116.,   5139.,   4330.,   2323.,   1800.,   &
-          1587.,   1511.,   1593.,   1740.,   2190.,   2171. /
-      ! 21-Parana+Uruguay(Arg./Uru.) 34S12 58W18            R1 
-      ! Old=470 km3/year=14893m3/s; UNESCO(p323/324)=22263m3/s
-      ! --> 22203.168m3/s 
-      DATA inb(21)/2/
-      DATA (iirnf(jc,21),jc=1,jpcoef)/ 114, 114, 298*0 /
-      DATA (ijrnf(jc,21),jc=1,jpcoef)/  46,  47, 298*0 /
-      DATA (zrnfm(jm,21),jm=1,jpmois)/   &
-          24803.,   26836.,   24585.,   21358.,   20041.,   21689.,   &
-          22647.,   20403.,   19627.,   21400.,   21126.,   21923./
-      ! 22-Negro (Argentina)  41S02  62W47                  R2
-      ! Old=30 km3/year=951m3/s; UNESCO(65-79,p330)=765m3/s; ratio=1.24
-      ! --> 907.666687m3/s
-      DATA inb(22)/1/
-      DATA (iirnf(jc,22),jc=1,jpcoef)/ 110, 299*0 /
-      DATA (ijrnf(jc,22),jc=1,jpcoef)/  43, 299*0 /
-      DATA (zrnfm(jm,22),jm=1,jpmois)/   &
-          882.,     587.,     519.,     482.,     623.,     984.,   &
-          1120.,   1294.,    1092.,    1166.,    1105.,    1038./
-     ! 23- Tocantins :  ~48W00 ~0.5S                        R4
-     ! UNESCO(76-77+80-81, p335)=12167m3/s (station far inland!)
-     ! --> 12173.082m3/s
-      DATA inb(23)/1/
-      DATA (iirnf(jc,23),jc=1,jpcoef)/118, 299*0 /
-      DATA (ijrnf(jc,23),jc=1,jpcoef)/ 73, 299*0 /
-      DATA (zrnfm(jm,23),jm=1,jpmois)/   &
-          15507.,  24855.,  27854.,  24257.,  15860.,   7861.,   &
-          4440.,   3082.,   2728.,   3417.,   5511.,  10705. /
-
-      DATA inb(24)/0/
-      DATA inb(25)/0/
-
-      ! Total run-offs(II)=256755.75
-      ! Total run-offs(I+II)=317730
-
-      ! III EUROPE
-      ! ----------
-      ! 26-Rhone (France) 43N20 4E50                        R1 
-      ! Old=49 km3/year=1553m3/s; UNESCO(65-83, p512)=1229
-      ! --> 1217.16663m3/s
-      DATA inb(26)/1/
-      DATA (iirnf(jc,26),jc=1,jpcoef)/ 150, 299*0 /
-      DATA (ijrnf(jc,26),jc=1,jpcoef)/ 108, 299*0 /
-      DATA (zrnfm(jm,26),jm=1,jpmois)/   &
-          1385.,    1848.,    1568.,    1311.,    1420.,    1191.,   &
-           937.,     796.,     794.,     962.,    1161.,    1233./
-      ! 27-Po (Italy) 11N10 1W09                            R1 
-      ! Old=46 km3/year=1458m3/s; UNESCO(65-84, p527)=1550m3/s
-      ! --> 1550.08337m3/s
-      DATA inb(27)/1/
-      DATA (iirnf(jc,27),jc=1,jpcoef)/ 159, 299*0 /
-      DATA (ijrnf(jc,27),jc=1,jpcoef)/ 110, 299*0 /
-      DATA (zrnfm(jm,27),jm=1,jpmois)/   &
-          1304.,    1605.,    1856.,    1549.,    1958.,    1760.,   &
-          1040.,    1007.,    1615.,    1907.,    1772.,    1228./
-      ! 28-Danube (Roumania) 45N20 29E40                    R1 
-      ! Old=206 km3/year=6527m3/s; UNESCO(65-84, p544)=7011m3/s
-      ! --> 7011.8335m3/s
-      DATA inb(28)/1/
-      DATA (iirnf(jc,28),jc=1,jpcoef)/ 171, 299*0 /
-      DATA (ijrnf(jc,28),jc=1,jpcoef)/ 106, 299*0 /
-      DATA (zrnfm(jm,28),jm=1,jpmois)/   &
-          6625.,    7352.,    7986.,    9279.,    9482.,    8618.,   &
-          7586.,    6001.,    5210.,    4989.,    4921.,    6093./
-
-      DATA inb(29)/0/
-      DATA inb(30)/0/
-
-      ! Total run-offs(III)=9779.08301m3/s
-      ! Total run-offs(I+II+III)=327509.094m3/s
-
-      ! IV EURASIAN ARCTIC
-      ! ------------------
-      ! 31-Iana (Russia)                                    R1
-      ! Old=29 km3/year=919m3/s; UNESCO(65-72+74+76-84, p463)=936m3/s
-      ! --> 0m3/s
-      DATA inb(31)/0/
-      DATA (iirnf(jc,31),jc=1,jpcoef)/ 300*0 /
-      DATA (ijrnf(jc,31),jc=1,jpcoef)/ 300*0 /
-      DATA (zrnfm(jm,31),jm=1,jpmois)/   &
-             2.,       1.,       0.,       0.,     917.,    3514.,   &
-          2968.,    2329.,    1322.,     178.,      36.,       8./ 
-      ! 32-Ob (Russia) 66N45 69E30                          R1 
-      ! Old=385 km3/year=12200m3/s; UNESCO(65-84, p478)=12862m3/s
-      ! --> 12861.5m3/s
-      DATA inb(32)/2/
-      DATA (iirnf(jc,32),jc=1,jpcoef)/ 153, 153, 298*0 /
-      DATA (ijrnf(jc,32),jc=1,jpcoef)/ 146, 147, 298*0 /
-      DATA (zrnfm(jm,32),jm=1,jpmois)/   &
-           4977.,    4150.,    3611.,    3590.,   14723.,   33366.,   &
-          30773.,   22785.,   14701.,    9705.,    6078.,    5879./
-      ! 33-Yenesei (Russia) 71N50 82E40                     R1 
-      ! Old=560 km3/year=17745m3/s; UNESCO(65-84, p472)=17462m3/s
-      ! --> 17461.666m3/s
-      DATA inb(33)/1/
-      DATA (iirnf(jc,33),jc=1,jpcoef)/  32, 299*0 /
-      DATA (ijrnf(jc,33),jc=1,jpcoef)/ 147, 299*0 /
-      DATA (zrnfm(jm,33),jm=1,jpmois)/   &
-           6659.,    7007.,    6895.,    7089.,   24183.,   74965.,   &
-          25325.,   15940.,   16050.,   12678.,    6741.,    6008./
-     ! 34-Severnaya Dvina(Russia)64N32 40E30                R5 
-     ! Old=106 km3/year=3359m3/s; not in UNESCO data base
-     ! --> 3359.00024m3/s
-      DATA inb(34)/1/
-      DATA (iirnf(jc,34),jc=1,jpcoef)/ 153, 299*0 /
-      DATA (ijrnf(jc,34),jc=1,jpcoef)/ 133, 299*0 /
-      DATA (zrnfm(jm,34),jm=1,jpmois)/12*3359./
-      ! 35-Lena (Russia) 72N25 126E40                       R1 
-      ! --> 16192.333m3/s
-      DATA inb(35)/1/
-      DATA (iirnf(jc,35),jc=1,jpcoef)/  36, 299*0 /
-      DATA (ijrnf(jc,35),jc=1,jpcoef)/ 137, 299*0 /
-      DATA (zrnfm(jm,35),jm=1,jpmois)/   &
-           2663.,    2095.,    1629.,    1324.,    7203.,   64230.,   &
-          40810.,   28310.,   25995.,   13299.,    3673.,    3077./
-      ! 36-Kolyma (Russia) 69N30 161E00                     R1
-      ! Old=71 km3/year=2250m3/s; UNESCO(65-84)=2138m3/s
-      ! --> 2138.66675m3/s
-      DATA inb(36)/1/
-      DATA (iirnf(jc,36),jc=1,jpcoef)/  46, 299*0 /
-      DATA (ijrnf(jc,36),jc=1,jpcoef)/ 128, 299*0 /
-      DATA (zrnfm(jm,36),jm=1,jpmois)/   &
-           131.,      92.,      76.,      65.,    2576.,    9143.,   &
-          4670.,    3727.,    3226.,    1119.,     630.,     209./
-     ! 37-Indigirka (russia) 70N48 148E54                   R1 
-     ! Old=55 km3/year=1743m3/s; UNESCO(65-68+70+72-84)=1622.
-     ! --> 1622.00012m3/s
-      DATA inb(37)/1/
-      DATA (iirnf(jc,37),jc=1,jpcoef)/  42, 299*0 /
-      DATA (ijrnf(jc,37),jc=1,jpcoef)/ 131, 299*0 /
-      DATA (zrnfm(jm,37),jm=1,jpmois)/   &
-            32.,      17.,      11.,       7.,     354.,    5629.,   &
-          5210.,    4617.,    2823.,     559.,     142.,      63./
-
-      DATA inb(38)/0/
-      DATA inb(39)/0/
-      DATA inb(40)/0/
-
-      ! Total run-offs(IV)=40773.6641m3/s
-      ! Total run-offs(I+II+III+IV)=368282.75m3/s
-
-      ! V ASIA
-      ! ------
-      ! 41-Amur (Russia) 52N56 141E10                       R2
-      ! Old=325 km3/year=10299m3/s; UNESCO(65-84, p451)=8792m3/s; ratio=1.17
-      ! --> 10299.332m3/s
-      DATA inb(41)/1/
-      DATA (iirnf(jc,41),jc=1,jpcoef)/  33, 299*0 /
-      DATA (ijrnf(jc,41),jc=1,jpcoef)/ 116, 299*0 /
-      DATA (zrnfm(jm,41),jm=1,jpmois)/   &
-           2108.,    1543.,    1252.,    4088.,   15464.,   17336.,   &
-          16100.,   20801.,   21548.,   15658.,    5188.,    2506./
-      ! 42-Liaohe (China) 40N40 122E09                      R5
-      ! Old=6 km3/year=190m3/s; not found in UNESCO data base
-      ! --> 189.999985m3/s
-      DATA inb(42)/1/
-      DATA (iirnf(jc,42),jc=1,jpcoef)/  22, 299*0 /
-      DATA (ijrnf(jc,42),jc=1,jpcoef)/ 104, 299*0 /
-      DATA (zrnfm(jm,42),jm=1,jpmois)/12*190./
-      ! 43-Daling (China) 40N56 121E43                      R5
-      ! Old=1 km3/year=32m3/s; not found in UNESCO data base
-      ! --> 31.9999962m3/s
-      DATA inb(43)/1/
-      DATA (iirnf(jc,43),jc=1,jpcoef)/  24, 299*0 /
-      DATA (ijrnf(jc,43),jc=1,jpcoef)/ 104, 299*0 /
-      DATA (zrnfm(jm,43),jm=1,jpmois)/12*32./
-      ! 44-Haiho (China) 6N17   1E23
-      ! Old=2 km3/year=64m3/s; not found in UNESCO data base
-      ! --> 63.9999924m3/s
-      DATA inb(44)/1/
-      DATA (iirnf(jc,44),jc=1,jpcoef)/  28, 299*0 /
-      DATA (ijrnf(jc,44),jc=1,jpcoef)/ 106, 299*0 /
-      DATA (zrnfm(jm,44),jm=1,jpmois)/12*64./
-      ! 45-Yellow/Huangho (China) 37N32 118E19              R2 
-      ! Old=49 km3/year=1553m3/s; UNESCO(p371)=1210m3/s; ratio=1.28
-      ! --> 1552.5m3/s
-      DATA inb(45)/1/
-      DATA (iirnf(jc,45),jc=1,jpcoef)/  22, 299*0 /
-      DATA (ijrnf(jc,45),jc=1,jpcoef)/ 103, 299*0 /
-      DATA (zrnfm(jm,45),jm=1,jpmois)/   &
-           697.,     674.,    1282.,    1206.,    1109.,     778.,   &
-          1964.,    3270.,    3668.,    2133.,     978.,     871./
-      ! 46-Yangtze (China) 31N48 121E10=Changjiang          R2
-      ! Old=900 km3/year=28520m3/s; UNESCO(p372)=25032m3/s; ratio=1.14
-      ! --> 28519.584m3/s
-      DATA inb(46)/1/
-      DATA (iirnf(jc,46),jc=1,jpcoef)/  23, 299*0 /
-      DATA (ijrnf(jc,46),jc=1,jpcoef)/  98, 299*0 /
-      DATA (zrnfm(jm,46),jm=1,jpmois)/   &
-          10994.,   10548.,   13254.,   24097.,   36972.,   47197.,   &
-          53891.,   42070.,   37627.,   30762.,   22046.,   12777./
-      ! 47-Pearl/Zhujiang (China) 23N36 106E34=Xijiang      R2 
-      ! Old=302 km3/year=9570m3/s; UNESCO(76,82, p373)=6987m3/s; ratio=1.37
-      ! --> 9570.75m3/s
-      DATA inb(47)/1/
-      DATA (iirnf(jc,47),jc=1,jpcoef)/  22, 299*0 /
-      DATA (ijrnf(jc,47),jc=1,jpcoef)/  95, 299*0 /
-      DATA (zrnfm(jm,47),jm=1,jpmois)/   &
-           2357.,    2415.,    2634.,    6446.,   14246.,   18395.,   &
-          19213.,   17558.,   12312.,    7255.,    5864.,    6154./
-      ! 48-Red/Hungho (Vietnam) 20N17 106E34                R5 
-      ! Old=123 km3/year=3897m3/s; not in UNESCO data base
-      ! --> 3897.m3/s
-      DATA inb(48)/1/
-      DATA (iirnf(jc,48),jc=1,jpcoef)/  16, 299*0 /
-      DATA (ijrnf(jc,48),jc=1,jpcoef)/  91, 299*0 /
-      DATA (zrnfm(jm,48),jm=1,jpmois)/12*3897./
-      ! 49-Mekong (Vietnam) 10N33 105E24                    R5
-      ! Old=470 km3/year=14893m3/s; not in UNESCO data base
-      ! --> 14893.0029m3/s
-      DATA inb(49)/2/
-      DATA (iirnf(jc,49),jc=1,jpcoef)/  16,  16, 298*0 /
-      DATA (ijrnf(jc,49),jc=1,jpcoef)/  86,  87, 298*0 /
-      DATA (zrnfm(jm,49),jm=1,jpmois)/12*14893./
-      ! 50-Irrawady (Burma) 15N50  95E06                    R5
-      ! Old=428 km3/year=13563m3/s; not in UNESCO data base
-      ! --> 13563.m3/s
-      DATA inb(50)/1/
-      DATA (iirnf(jc,50),jc=1,jpcoef)/  10, 299*0 /
-      DATA (ijrnf(jc,50),jc=1,jpcoef)/  90, 299*0 /
-      DATA (zrnfm(jm,50),jm=1,jpmois)/12*13563./
-      ! 51-Ganges+Brahmaputra (Beng.) 22N00  91E00          R1
-      ! Old=971 km3/year=30769m3/s; UNESCO(69-70+73-75, p367)=31760m3/s
-      ! --> 32147.498m3/s
-      DATA inb(51)/2/
-      DATA (iirnf(jc,51),jc=1,jpcoef)/   7,  8, 298*0 /
-      DATA (ijrnf(jc,51),jc=1,jpcoef)/  93, 93, 298*0 /
-      DATA (zrnfm(jm,51),jm=1,jpmois)/   &
-           6623.,    6315.,    6432.,    9410.,   17263.,   38302.,   &
-          64688.,   80338.,   84802.,   43387.,   17888.,   10322./
-      ! 52-Mehandi (India)                                  R5
-      ! Old=67 km3/year=2123m3/s; not in UNESCO data base
-      ! --> 0m3/s
-      DATA inb(52)/0/
-      DATA (iirnf(jc,52),jc=1,jpcoef)/ 300*0 /
-      DATA (ijrnf(jc,52),jc=1,jpcoef)/ 300*0 /
-      DATA (zrnfm(jm,52),jm=1,jpmois)/12*2123./
-      ! 53-Damodar (India)                                  R2
-      ! Old=10 km3/year=320m3/s; UNESCO(p386)=173m3/s; ratio=1.85
-      ! --> 432.916656m3/s
-      DATA inb(53)/1/
-      DATA (iirnf(jc,53),jc=1,jpcoef)/   6, 299*0 /
-      DATA (ijrnf(jc,53),jc=1,jpcoef)/  93, 299*0 /
-      DATA (zrnfm(jm,53),jm=1,jpmois)/   &
-           51.,      37.,      42.,      49.,      85.,     296.,   &
-          896.,    1390.,    1591.,     542.,     172.,      44./   
-      ! 54-Godavari (India) 17N00  81E45                    R2
-      ! Old=84 km3/year=2662m3/s; UNESCO(p379)=1916m3/s; ratio=1.39
-      ! --> 3670.5m3/s
-      DATA inb(54)/1/
-      DATA (iirnf(jc,54),jc=1,jpcoef)/   3, 299*0 /
-      DATA (ijrnf(jc,54),jc=1,jpcoef)/  90, 299*0 /
-      DATA (zrnfm(jm,54),jm=1,jpmois)/   &
-           988.,     737.,     280.,     219.,     196.,    1691.,   &
-          8063.,   14571.,   11428.,    3953.,    1221.,     699./   
-     ! 55-Indus (Pakistan) 24N20  67E47                     R2
-     ! Old=238 km3/year=7542m3/s; UNESCO (76-79, p428)=2396m3/s; ratio=3.15
-     ! --> 7562.m3/s
-      DATA inb(55)/2/
-      DATA (iirnf(jc,55),jc=1,jpcoef)/ 176, 176, 298*0 /
-      DATA (ijrnf(jc,55),jc=1,jpcoef)/  93,  94, 298*0 /
-      DATA (zrnfm(jm,55),jm=1,jpmois)/   &
-            988.,     737.,    1904.,    1968.,    3625.,    6143.,   &
-          15969.,    36656.,  17173.,    3215.,    1251.,    1115./   
-      ! 56-Tigris and Euphrates (Iraq) 31N00  47E25         R1
-      ! Old=46 km3/year=1458m3/s; UNESCO(76-79, p428)=2396m3/s
-      ! --> 2248.8335m3/s
-      DATA inb(56)/1/
-      DATA (iirnf(jc,56),jc=1,jpcoef)/ 172, 299*0 /
-      DATA (ijrnf(jc,56),jc=1,jpcoef)/  94, 299*0 /
-      DATA (zrnfm(jm,56),jm=1,jpmois)/   &
-          1872.,    2127.,    2962.,    4944.,    5036.,    3078.,   &
-          1362.,     900.,     786.,     993.,    1184.,    1742./   
-
-      DATA inb(57)/0/
-      DATA inb(58)/0/
-      DATA inb(59)/0/
-      DATA inb(60)/0/
-
-      ! Total run-offs (V)=128642.922m3/s
-      ! Total run-offs (I+II+III+IV+V)=496925.688m3/s
-      
-      ! VI AFRICA
-      ! ----------
-      ! 61-Nile (Egypt) 30N10  31E06                        R1
-      ! Old=30 km3/year=951m3/s; UNESCO(73-84)=1251m3/s
-      ! --> 1251.33337m3/s
-      DATA inb(61)/1/
-      DATA (iirnf(jc,61),jc=1,jpcoef)/ 176, 299*0 /
-      DATA (ijrnf(jc,61),jc=1,jpcoef)/  97, 299*0 /
-      DATA (zrnfm(jm,61),jm=1,jpmois)/   &
-          1239.,    1035.,    1126.,    1110.,    1170.,    1542.,   &
-          1742.,    1559.,    1191.,    1103.,    1075.,    1124./   
-      ! 62-Niger (Nigeria) 5N33   6E33                      R1
-      ! Old=192 km3/year=6084m3/s; HLA=5545m3/s;
-      ! Data from HLA: 4.5N -  6  E   1950/87
-      ! --> 5545.08301m3/s
-      DATA inb(62)/1/
-      DATA (iirnf(jc,62),jc=1,jpcoef)/ 146, 299*0 /
-      DATA (ijrnf(jc,62),jc=1,jpcoef)/  81, 299*0 /
-      DATA (zrnfm(jm,62),jm=1,jpmois)/   &
-          2029.,    1804.,    1573.,    1418.,    1494.,    2561.,   &
-          5272.,    9629.,   14951.,   16227.,    6842.,    2741./   
-      ! 63-Zaire/Congo (Congo) 6S04  12E24                  R1      
-      ! Old=1250 km3/year=39610m3/s; UNESCO(65-83, p149)=41941m3/s; HLA
-      ! =39720m3/s
-      ! Data from HLA: 6S - 12  E   1902/96
-      ! --> 39719.5m3/s
-      DATA inb(63)/2/
-      DATA (iirnf(jc,63),jc=1,jpcoef)/ 148, 148, 298*0 /
-      DATA (ijrnf(jc,63),jc=1,jpcoef)/  63,  64, 298*0 /
-      DATA (zrnfm(jm,63),jm=1,jpmois)/   &
-          47200.,   37692.,   34374.,   36603.,   37862.,   35588.,   &
-          30966.,   30774.,   35981.,   42825.,   51455.,   55314./  
-      ! 64-Orange (Sth Africa) 28S41 16E28                  R5
-      ! Old=11 km3/year=349m3/s
-      ! --> 349.000031m3/s
-      DATA inb(64)/1/
-      DATA (iirnf(jc,64),jc=1,jpcoef)/ 150, 299*0 /
-      DATA (ijrnf(jc,64),jc=1,jpcoef)/  51, 299*0 /
-      DATA (zrnfm(jm,64),jm=1,jpmois)/12*349./
-      ! 65-Zambesi (Mozambique) 18S55 36E04                 R5
-      ! Old=223 km3/year=7066m3/s
-      ! --> 7066.00049m3/s
-      DATA inb(65)/1/
-      DATA (iirnf(jc,65),jc=1,jpcoef)/ 160, 299*0 /
-      DATA (ijrnf(jc,65),jc=1,jpcoef)/  56, 299*0 /
-      DATA (zrnfm(jm,65),jm=1,jpmois)/12*7066./
-      ! 66-Limpopo (mozambique) 25S15 33E30                 R5
-      ! Old=5 km3/year=158m3/s
-      ! --> 158.m3/s
-      DATA inb(66)/1/
-      DATA (iirnf(jc,66),jc=1,jpcoef)/ 159, 299*0 /
-      DATA (ijrnf(jc,66),jc=1,jpcoef)/  52, 299*0 /
-      DATA (zrnfm(jm,66),jm=1,jpmois)/12*158./
-      ! 67-Rufiji (tanzania) 8S00 39E20                     R5
-      ! Old=9 km3/year=285m3/s
-      ! --> 285.m3/s
-      DATA inb(67)/1/
-      DATA (iirnf(jc,67),jc=1,jpcoef)/ 162, 299*0 /
-      DATA (ijrnf(jc,67),jc=1,jpcoef)/  62, 299*0 /
-      DATA (zrnfm(jm,67),jm=1,jpmois)/12*285./
-      ! Other rivers from the Gulf of Guinea
-      ! 68-Kouilou (Congo)  4.5S-11.5E (UNESCO=p51)         R1
-      ! Data from HLA:  4S - 11  E   1953/86 = 919m3/s
-      ! --> 924.25m3/s
-      DATA inb(68)/1/
-      DATA (iirnf(jc,68),jc=1,jpcoef)/ 147, 299*0 /
-      DATA (ijrnf(jc,68),jc=1,jpcoef)/  67, 299*0 /
-      DATA (zrnfm(jm,68),jm=1,jpmois)/   &
-          1115.,    1166.,    1285.,    1462.,    1454.,     725.,   &
-           480.,     378.,     327.,     400.,    1018.,    1281./
-      ! 69-Nyanga (Gabon)  3S-10E (UNESCO=p66)              R1
-      ! Data from HLA: 2.8N - 10  E   1952/82 = 264m3/s
-      ! --> 264.166687m3/s
-      DATA inb(69)/1/
-      DATA (iirnf(jc,69),jc=1,jpcoef)/ 147, 299*0 /
-      DATA (ijrnf(jc,69),jc=1,jpcoef)/  68, 299*0 /
-      DATA (zrnfm(jm,69),jm=1,jpmois)/   &
-          126.,      99.,     143.,     253.,     316.,     300.,   &
-          156.,      93.,     231.,     557.,     592.,     304./
-      ! 70-Ogoue (Gabon) 1S-9E (UNESCO, p69)                R1
-      ! Data from HLA: .75S -  9  E   1930/90 = 4512m3/s
-      ! --> 4511.83301m3/s
-      DATA inb(70)/1/
-      DATA (iirnf(jc,70),jc=1,jpcoef)/ 146, 299*0 /
-      DATA (ijrnf(jc,70),jc=1,jpcoef)/  72, 299*0 /
-      DATA (zrnfm(jm,70),jm=1,jpmois)/   &
-          4514.,    3988.,    4622.,    5403.,    5790.,    4145.,   &
-          2616.,    1961.,    2119.,    4304.,    7870.,    6810./
-      ! 71-Sanaga (Cameroun) 3.5N-9.5E (UNESCO, p31)        R1
-      ! Data from HLA: 3.9N -  9.6E   1944/80 = 1983m3/s
-      ! --> 1982.58337m3/s
-      DATA inb(71)/1/
-      DATA (iirnf(jc,71),jc=1,jpcoef)/ 146, 299*0 /
-      DATA (ijrnf(jc,71),jc=1,jpcoef)/  81, 299*0 /
-      DATA (zrnfm(jm,71),jm=1,jpmois)/   &
-           813.,     578.,     517.,     633.,     909.,    1278.,   &
-          1928.,    2653.,    4361.,    5470.,    3307.,    1344./
-      ! 72-Nyong (Cameroun) 3.1N-9.5E (UNESCO, p36)         R1
-      ! Data from HLA: 3.5N -  9.8E   1951/84 = 432m3/s
-      ! --> 431.166656m3/s      
-      DATA inb(72)/1/
-      DATA (iirnf(jc,72),jc=1,jpcoef)/ 146, 299*0 /
-      DATA (ijrnf(jc,72),jc=1,jpcoef)/  80, 299*0 /
-      DATA (zrnfm(jm,72),jm=1,jpmois)/   &
-          232.,     128.,     150.,     266.,     402.,     459.,   &
-          353.,     306.,     604.,     930.,     878.,     466./
-      ! 73-Ntem (Cameroun) 2.5N-9.5E (UNESCO, p36)          R1
-      ! Data from HLA: 3N -  9.8E   1951/87 = 330m3/s      
-      ! --> 329.75m3/s
-      DATA inb(73)/1/
-      DATA (iirnf(jc,73),jc=1,jpcoef)/ 146, 299*0 /
-      DATA (ijrnf(jc,73),jc=1,jpcoef)/  79, 299*0 /
-      DATA (zrnfm(jm,73),jm=1,jpmois)/   &
-          452.,     455.,     492.,     466.,     482.,     251.,   &
-          151.,     108.,      87.,     133.,     395.,     485./
-      ! 74-Wouri (Cameroun) 4N-9.5E (UNESCO, p34)           R1
-      ! Data from HLA: 4N -  9.5E   1951/84  = 307m3/s
-      ! --> 306.583374m3/s
-      DATA inb(74)/1/
-      DATA (iirnf(jc,74),jc=1,jpcoef)/ 145, 299*0 /
-      DATA (ijrnf(jc,74),jc=1,jpcoef)/  81, 299*0 /
-      DATA (zrnfm(jm,74),jm=1,jpmois)/   &
-           88.,      63.,      68.,      95.,     139.,     222.,   &
-          395.,     633.,     825.,     675.,     328.,     148./
-      ! 75-Oueme (Benin) 6.5N-2.5E (UNESCO, p27)            R1
-      ! Data from HLA:  6.2N -  2.5E   1951/84 =113m3/s
-      ! --> 112.666664m3/s
-      DATA inb(75)/1/
-      DATA (iirnf(jc,75),jc=1,jpcoef)/ 144, 299*0 /
-      DATA (ijrnf(jc,75),jc=1,jpcoef)/  83, 299*0 /
-      DATA (zrnfm(jm,75),jm=1,jpmois)/   &
-           0.,       0.,       0.,       0.,       1.,      18.,   &
-          96.,     324.,     531.,     319.,      58.,       5./
-      ! 76-Mono(Tepetou,Benin/Togo) 6.25N-1.8E (U, p27)     R1
-      ! Data from HLA: 6.2N -  1  E   1951/85  = 96m3/s       
-      ! --> 96.4999924m3/s
-      DATA inb(76)/1/
-      DATA (iirnf(jc,76),jc=1,jpcoef)/ 143, 299*0 /
-      DATA (ijrnf(jc,76),jc=1,jpcoef)/  83, 299*0 /
-      DATA (zrnfm(jm,76),jm=1,jpmois)/   &
-            2.,       1.,       1.,       4.,       8.,      44.,   &
-          168.,     274.,     386.,     215.,      47.,       8./
-      ! 77-Comoe (Ivory Coast) 5N-4W (UNESCO, p56)          R1
-      ! Data from HLA: 5.1N -  4  W   1955/86 = 194m3/s       
-      ! --> 194.25m3/s
-      DATA inb(77)/1/
-      DATA (iirnf(jc,77),jc=1,jpcoef)/ 140, 299*0 /
-      DATA (ijrnf(jc,77),jc=1,jpcoef)/  82, 299*0 /
-      DATA (zrnfm(jm,77),jm=1,jpmois)/   &
-           14.,       7.,       9.,      17.,      29.,      90.,   &
-          145.,     349.,     740.,     688.,     198.,      45./
-      ! 78-Sassandra (Ivory Coast) 5N-6W (UNESCO, p58)      R1
-      ! Data from HLA: 5N -  6  W   1953/86 = 449m3/s
-      ! --> 449.333344m3/s
-      DATA inb(78)/1/
-      DATA (iirnf(jc,78),jc=1,jpcoef)/ 139, 299*0 /
-      DATA (ijrnf(jc,78),jc=1,jpcoef)/  82, 299*0 /
-      DATA (zrnfm(jm,78),jm=1,jpmois)/   &
-          144.,     116.,     119.,     143.,     156.,     334.,   &
-          419.,     692.,    1277.,    1190.,     553.,     249./
-
-      DATA inb(79)/0/
-      DATA inb(80)/0/
-      DATA inb(82)/0/
-      DATA inb(83)/0/
-      DATA inb(84)/0/
-      DATA inb(85)/0/
-      DATA inb(86)/0/
-      DATA inb(87)/0/
-      DATA inb(88)/0/
-      DATA inb(89)/0/
-      DATA inb(90)/0/
-
-      ! Total run-offs (VI)=63977.m3/s
-      ! Total run-offs (I+II+III+IV+V+VI)=560902.688m3/s
-      
-      ! VII OCEANIA
-      ! ------------
-      ! 91-Murray (Austrlia) 35S34 138E54                   R5
-      ! Old=22 km3/year=697m3/s
-      ! --> 696.999939m3/s
-      DATA inb(91)/1/
-      DATA (iirnf(jc,91),jc=1,jpcoef)/  30, 299*0 /
-      DATA (ijrnf(jc,91),jc=1,jpcoef)/  47, 299*0 /
-      DATA (zrnfm(jm,91),jm=1,jpmois)/12*697./
-      ! 92-Haast (N. Zeland) 43S50 169E02                   R5
-      ! Old=6 km3/year=190m3/s
-      ! --> 189.999985m3/s
-      DATA inb(92)/1/
-      DATA (iirnf(jc,92),jc=1,jpcoef)/  48, 299*0 /
-      DATA (ijrnf(jc,92),jc=1,jpcoef)/  41, 299*0 /
-      DATA (zrnfm(jm,92),jm=1,jpmois)/12*190./
-      ! 93-Fly (N. Guinea) 8S30 143E05                      R5
-      ! Old=77 km3/year=2440m3/s
-      ! --> 2440.m3/s
-      DATA inb(93)/1/
-      DATA (iirnf(jc,93),jc=1,jpcoef)/  34, 299*0 /
-      DATA (ijrnf(jc,93),jc=1,jpcoef)/  62, 299*0 /
-      DATA (zrnfm(jm,93),jm=1,jpmois)/12*2440./
-      ! 94-Purari (N. Guinea) 7S25 145E05                   R4
-      ! Old=77 km3/year=2440m3/s; UNESCO(77+83,p595)=2727m3/s
-      ! --> 2727.33325m3/s
-      DATA inb(94)/1/
-      DATA (iirnf(jc,94),jc=1,jpcoef)/  36, 299*0 /
-      DATA (ijrnf(jc,94),jc=1,jpcoef)/  64, 299*0 /
-      DATA (zrnfm(jm,94),jm=1,jpmois)/   &
-          1972.,    3067.,    2332.,    2981.,    2693.,    2876.,   &
-          2867.,    3477.,    2195.,    2967.,    2674.,    2627./
-      ! 95-Choshui ch''i (taiwan) 24N03 120E23              R2
-      ! Old=6 km3/year=190m3/s; UNESCO(65-68, p374)=120m3/s; ratio=1.58 
-      ! --> 0.m3/s
-      DATA inb(95)/0/
-      DATA (iirnf(jc,95),jc=1,jpcoef)/  22, 299*0 /
-      DATA (ijrnf(jc,95),jc=1,jpcoef)/  94, 299*0 /
-      DATA (zrnfm(jm,95),jm=1,jpmois)/   &
-           65.,      63.,      83.,     100.,     177.,     535.,   &
-          332.,     313.,     236.,     193.,     101.,      84./
-      ! 96-Kaoping (Taiwan)                                 R5  
-      ! Old=9 km3/year=285m3/s
-      ! --> 0.m3/s
-      DATA inb(96)/0/
-      DATA (iirnf(jc,96),jc=1,jpcoef)/ 300*0 /
-      DATA (ijrnf(jc,96),jc=1,jpcoef)/ 300*0 /
-      DATA (zrnfm(jm,96),jm=1,jpmois)/12*285./
-      ! 97-Tsengwen (Taiwan)                                R5
-      ! Old=2 km3/year=63m3/s
-      ! --> 0.m3/s
-      DATA inb(97)/0/
-      DATA (iirnf(jc,97),jc=1,jpcoef)/ 300*0 /
-      DATA (ijrnf(jc,97),jc=1,jpcoef)/ 300*0 /
-      DATA (zrnfm(jm,97),jm=1,jpmois)/12*63./
-      ! 98-Hualien (Taiwan)                                 R5
-      ! Old=4 km3/year=127m3/s
-      ! --> 0.m3/s
-      DATA inb(98)/0/
-      DATA (iirnf(jc,98),jc=1,jpcoef)/ 300*0 /
-      DATA (ijrnf(jc,98),jc=1,jpcoef)/ 300*0 /
-      DATA (zrnfm(jm,98),jm=1,jpmois)/12*127./
-      ! 99-Peinan (Taiwan) 22N46 121E10                     R5
-      ! Old=4 km3/year=127m3/s
-      ! --> 127.000008m3/s
-      DATA inb(99)/1/
-      DATA (iirnf(jc,99),jc=1,jpcoef)/  23, 299*0 /
-      DATA (ijrnf(jc,99),jc=1,jpcoef)/  94, 299*0 /
-      DATA (zrnfm(jm,99),jm=1,jpmois)/12*127./
-      ! 100-Hsiukuluan (Taiwan)                             R5
-      ! Old=4 km3/year=127m3/s
-      ! --> 0.m3/s
-      DATA inb(100)/0/
-      DATA (iirnf(jc,100),jc=1,jpcoef)/ 300*0 /
-      DATA (ijrnf(jc,100),jc=1,jpcoef)/ 300*0 /
-      DATA (zrnfm(jm,100),jm=1,jpmois)/12*127./
-
-      DATA inb(101)/0/
-      DATA inb(102)/0/
-
-      !  Total run-offs(VII)=6181.33301m3/s
-      !  Total run-offs(I+II+III+IV+V+VI+VII)= 567084. m3/s
-      !                                      = 1.588668056E-6 kg/m2/s
-      !                                      = 17883.5586 km3/Year
-      !                                      = 0.567083955 Sverdrup
-
-      ! VIII ANTARCTICA
-      ! ---------------
-      ! 103 -Antarcti! ice discharge (2 613 Gt/year) (jacobs et al 1992)
-      ! Data from Jacobs et al 92: =2613 km3/year=82857.686454 m3/s
-      ! --> 82858 m3/s
-      DATA inb(103)/295/
-      DATA (iirnf(jc,103),jc=1,jpcoef)/   &
-            1,    2,    3,    3,    4,    4,    5,    6,    7,    8,   &
-            9,    9,   10,   11,   12,   13,   14,   14,   15,   16,   &
-           17,   18,   19,   20,   21,   22,   23,   24,   25,   26,   &
-           27,   28,   29,   30,   31,   32,   32,   33,   34,   35,   &
-           35,   36,   36,   37,   38,   39,   40,   40,   41,   42,   &
-           42,   43,   43,   44,   44,   44,   44,   44,   44,   44,   &
-           45,   45,   45,   45,   46,   46,   46,   46,   46,   47,   &
-           47,   47,   47,   47,   48,   48,   48,   48,   48,   49,   &
-           50,   51,   51,   52,   53,   54,   55,   56,   57,   58,   &
-           59,   60,   61,   62,   62,   63,   63,   64,   65,   66,   &
-           66,   67,   67,   68,   69,   69,   70,   71,   71,   72,   &
-           73,   73,   74,   75,   76,   76,   77,   78,   78,   79,   &
-           80,   81,   82,   83,   83,   84,   85,   86,   87,   88,   &
-           88,   89,   89,   90,   91,   92,   93,   94,   95,   96,   &
-           97,   98,   99,  100,  101,  102,  102,  103,  104,  104,   &
-          104,  104,  104,  105,  105,  106,  106,  107,  107,  107,   &
-          108,  108,  109,  109,  110,  110,  111,  111,  112,  112,   &
-          112,  112,  112,  112,  112,  112,  112,  112,  112,  112,   &
-          112,  112,  112,  112,  112,  113,  113,  113,  113,  113,   &
-          114,  114,  114,  114,  115,  115,  115,  115,  115,  116,   &
-          116,  117,  117,  118,  118,  119,  119,  120,  121,  122,   &
-          123,  124,  124,  125,  125,  126,  126,  127,  128,  128,   &
-          129,  130,  130,  131,  132,  132,  133,  133,  133,  133,   &
-          133,  134,  135,  135,  136,  136,  136,  137,  137,  138,   &
-          139,  140,  141,  142,  142,  143,  144,  145,  146,  147,   &
-          148,  149,  150,  151,  152,  153,  154,  155,  156,  156,   &
-          157,  158,  159,  160,  161,  161,  162,  163,  163,  164,   &
-          165,  165,  166,  166,  167,  168,  169,  170,  171,  171,   &
-          172,  173,  173,  174,  175,  176,  177,  177,  178,  178,   &
-          179,  180,  180,  181,  181,    0,    0,    0,    0,    0    &
-                                      /
-      DATA (ijrnf(jc,103),jc=1,jpcoef)/   &
-           19,   20,   20,   21,   21,   22,   22,   22,   22,   22,   &
-           22,   23,   23,   23,   23,   23,   22,   23,   22,   22,   &
-           22,   22,   22,   22,   22,   22,   22,   22,   22,   22,   &
-           22,   22,   22,   22,   22,   21,   22,   21,   21,   20,   &
-           21,   19,   20,   19,   19,   19,   18,   19,   18,   17,   &
-           18,   16,   17,    3,    4,    5,    6,    7,    8,   16,   &
-            3,    8,    9,   16,    3,    9,   10,   15,   16,    3,   &
-           10,   11,   14,   15,    3,   11,   12,   13,   14,    3,   &
-            3,    2,    3,    2,    2,    2,    2,    2,    2,    2,   &
-            2,    2,    2,    2,    3,    3,    4,    4,    4,    4,   &
-            5,    5,    6,    6,    6,    7,    7,    7,    8,    8,   &
-            8,    9,    9,    9,    9,   10,   10,   10,   11,   11,   &
-           11,   11,   11,   10,   11,   10,   10,   10,   10,   10,   &
-           11,   11,   12,   12,   12,   12,   12,   12,   12,   12,   &
-           12,   12,   12,   12,   12,   12,   13,   13,   13,   14,   &
-           15,   16,   17,   17,   18,   18,   19,   19,   20,   21,   &
-           21,   22,   22,   23,   23,   24,   24,   25,    8,    9,   &
-           10,   11,   12,   13,   14,   15,   16,   17,   18,   19,   &
-           20,   21,   22,   25,   26,    7,    8,   22,   23,   26,   &
-            7,   23,   24,   26,    6,    7,   24,   25,   26,    5,   &
-            6,    4,    5,    3,    4,    2,    3,    2,    2,    2,   &
-            2,    2,    3,    3,    4,    4,    5,    5,    5,    6,   &
-            6,    6,    7,    7,    7,    8,    8,    9,   10,   11,   &
-           12,   12,   12,   13,   13,   14,   15,   15,   16,   16,   &
-           16,   16,   16,   16,   17,   17,   17,   17,   17,   17,   &
-           17,   17,   17,   17,   17,   17,   17,   17,   17,   18,   &
-           18,   18,   18,   18,   18,   19,   19,   19,   20,   20,   &
-           20,   21,   21,   22,   22,   22,   22,   22,   21,   22,   &
-           21,   20,   21,   20,   20,   20,   19,   20,   18,   19,   &
-           18,   18,   19,   19,   20,    0,    0,    0,    0,    0    &
-                                      /
-      DATA (zrnfm(jm,103),jm=1,jpmois)/12*82858./
-
-      !  Total run-offs(VIII)=82858 m3/s
-      !  Total run-offs(I+II+III+IV+V+VI+VII+VIII)= 649941. m3/s
-      !                                           = 0.649940955 Sverdrup
-
-
-      DATA (inb(jr),jr=104,jpriv)/17*0/
-
-      !!----------------------------------------------------------------------
-      !!  OPA 8.5, LODYC-IPSL (2002)
-      !!----------------------------------------------------------------------
-
-      IF( kt == nit000 ) THEN
-         IF(lwp) WRITE(numout,*)
-         IF(lwp) WRITE(numout,*) 'flxrnf : runoff + orca2'
-         IF(lwp) WRITE(numout,*) '~~~~~~'
-         IF(lwp) WRITE(numout,*) '   nrunoff = ', nrunoff
-         IF(lwp) WRITE(numout,*)
-      ENDIF
-
-      ! 0. initialization and  indexes  modification for zoom version
-      ! --------------------------------------------------------------
-      
-      
-      ! conversion factor m3/s --> kg/m2/s *e1t*e2t
-
-      zconv = rauw
-
-      IF( kt == nit000 ) THEN 
-         ! Set runoffs and upstream coeff to zero
-         runoff (:,:) = 0.e0                                                                                  
-         upsrnfh(:,:) = 0.e0
-         upsrnfz(:)   = 0.e0
-         upsadv (:,:) = 0.e0
-         zcoefr(:,:,:) = 0.e0
-
-         IF( nrunoff == 0 ) THEN
-            IF(lwp) WRITE(numout,*)
-            IF(lwp) WRITE(numout,*) ' NO RUNOFF in this simulation '
-            IF(lwp) WRITE(numout,*)
-         ENDIF
-      ENDIF
-
-      ! 1. initialization of the runoff array
-      ! -------------------------------------
-
-      IF( nrunoff == 1 .OR. nrunoff == 2) THEN
-
-         iman  = jpmois
-         
-         i15 = nday / 16
-         
-         imois = nmonth + i15 - 1
-         IF( imois == 0) imois = iman
-         IF( nn_leap_year == 1 .AND. MOD( nyear, 4 ) == 0 ) THEN
-            ! Number of Day Between Data
-            idbd = nbiss(imois)
-         ELSEIF( nn_leap_year > 1 ) THEN
-            idbd = nn_leap_year
-         ELSE
-            idbd = nobis(imois)
-         ENDIF
-         ! Number of days between imois, 15 and the end of month
-         idmeom = idbd - 15
-         
-         ! 1.1 Calendar computation
-         
-         ! nrnf1 number of the first array record used in the simulation
-         ! nrnf2 number of the last  array record
-         
-         nrnf1 = imois
-         nrnf2 = nrnf1 + 1
-         nrnf1 = MOD( nrnf1, iman )
-         IF( nrnf1 == 0 ) nrnf1 = iman
-         nrnf2 = MOD( nrnf2, iman )
-         IF( nrnf2 == 0 ) nrnf2 = iman
-         
-         IF( kt == nit000 ) THEN
-            IF(lwp) WRITE(numout,*)
-            IF(lwp) WRITE(numout,*) 'flxrnf routine'
-            IF(lwp) WRITE(numout,*) 'first array record used nrnf1 ',nrnf1
-            IF(lwp) WRITE(numout,*) 'last  array record used nrnf2 ',nrnf2
-         ENDIF
-
-         zxy = FLOAT( nday + idmeom - idbd * i15 ) / idbd
-
-         DO jr = 1, jpriv-1
-            inbp = inb(jr)
-            IF( inbp /= 0 ) THEN
-               DO jl = 1, inbp
-                  ii = iirnf(jl,jr) - nimpp + 1 - jpizoom + 1
-                  ij = ijrnf(jl,jr) - njmpp + 1 - jpjzoom + 1
-                  iipt0 = MAX( 1, MIN( ii, jpi   ) )
-                  iipt1 = MAX( 0, MIN( ii, jpi-1 ) )
-                  ijpt0 = MAX( 1, MIN( ij, jpj   ) )
-                  ijpt1 = MAX( 0, MIN( ij, jpj-1 ) )
-                  DO jj = ijpt0, ijpt1
-                     DO ji = iipt0, iipt1
-                        runoff(ji,jj) = runoff(ji,jj)   & 
-                                      + ( (1.-zxy) * zrnfm(nrnf1,jr)   &
-                                      +       zxy  * zrnfm(nrnf2,jr)  ) / FLOAT(inbp)
-                     END DO
-                  END DO
-               END DO
-            ENDIF
-         END DO
-
-
-         ! 2. conversion in mm/day
-         ! -----------------------
-         !     C a u t i o n : upward water flux, runoff is negative
-
-         runoff(:,:) = -zconv * runoff(:,:) / ( e1t(:,:) * e2t(:,:) )
-         
-
-         ! 3. advectve coefficient
-         ! -----------------------
-
-         !  Creates the array coef that contains the coefficient to affect to
-         !  the upstream scheme. advection scheme will be:
-         !  coefr * upstream + (1- coefr) centered
-         !  coefr must be between 0 and 1.
-         !  3 values are used: maximum = zrup1, between zrup2, small zrup3.
-         ! 
-         IF( nrunoff == 2 .AND. kt == nit000 ) THEN
-
-            zrup1 = 0.5
-            zrup2 = 0.25
-            zrup3 = 0.125
-            DO jr = 1, jpriv-1
-               inbp = inb(jr)
-               IF( inbp /= 0) THEN
-                  DO jl = 1, inbp
-                     ! --> zrup1 (originally, = 1 = all upstream).
-                     ! ji,jj
-                     ii = iirnf(jl,jr) - nimpp + 1 - jpizoom - 1
-                     ij = ijrnf(jl,jr) - njmpp + 1 - jpjzoom - 1
-                     iipt0 = MAX( 1, MIN( ii,jpi  ) )
-                     iipt1 = MAX( 0, MIN( ii,jpi-1) )
-                     ijpt0 = MAX( 1, MIN( ij,jpj  ) )
-                     ijpt1 = MAX( 0, MIN( ij,jpj-1) )
-                     DO jj = ijpt0, ijpt1
-                        DO ji = iipt0, iipt1
-                           DO jn = 1, 3
-                              zcoefr(ji,jj,jn) = zrup1
-                           END DO
-                           zcoefr(ji,jj,4) = zrup2
-                           zcoefr(ji,jj,5) = zrup3
-                        END DO
-                     END DO
-                     ! --> zrup2 (originally = 0.5)
-                     ! ji,jj+1
-                     ijpt0 = MAX( 1, MIN( ij+1, jpj  ) )
-                     ijpt1 = MAX( 0, MIN( ij+1, jpj-1) )
-                     DO jj = ijpt0, ijpt1
-                        DO ji = iipt0, iipt1
-                           DO jn = 1, 3
-                              zcoefr(ji,jj,jn) = MAX( zrup2, zcoefr(ji,jj,jn) )
-                           END DO
-                           zcoefr(ji,jj,4) = MAX( zrup3, zcoefr(ji,jj,4) )
-                        END DO
-                     END DO
-                     ! ji,jj-1
-                     ijpt0 = MAX( 1, MIN( ij-1, jpj  ) )
-                     ijpt1 = MAX( 0, MIN( ij-1, jpj-1) )
-                     DO jj = ijpt0, ijpt1
-                        DO ji = iipt0, iipt1
-                           DO jn = 1, 3
-                              zcoefr(ji,jj,jn) = MAX( zrup2, zcoefr(ji,jj,jn) )
-                           END DO
-                           zcoefr(ji,jj,4) = MAX( zrup3, zcoefr(ji,jj,4) )
-                        END DO
-                     END DO
-                     ! ji+1,jj-1
-                     iipt0 = MAX( 1, MIN( ii+1, jpi  ) )
-                     iipt1 = MAX( 0, MIN( ii+1, jpi-1) )
-                     DO jj = ijpt0, ijpt1
-                        DO ji = iipt0, iipt1
-                           DO jn = 1, 3
-                              zcoefr(ji,jj,jn) = MAX( zrup2, zcoefr(ji,jj,jn) )
-                           END DO
-                           zcoefr(ji,jj,4) = MAX( zrup3, zcoefr(ji,jj,4) )
-                        END DO
-                     END DO
-                     ! ji+1,jj
-                     ijpt0 = MAX( 1, MIN( ij, jpj  ) )
-                     ijpt1 = MAX( 0, MIN( ij, jpj-1) )
-                     DO jj = ijpt0, ijpt1
-                        DO ji = iipt0, iipt1
-                           DO jn = 1, 3
-                              zcoefr(ji,jj,jn) = MAX( zrup2, zcoefr(ji,jj,jn) )
-                           END DO
-                           zcoefr(ji,jj,4) = MAX( zrup3, zcoefr(ji,jj,4) )
-                        END DO
-                     END DO
-                     ! ji+1,jj+1
-                     ijpt0 = MAX( 1, MIN( ij+1,jpj  ) )
-                     ijpt1 = MAX( 0, MIN( ij+1,jpj-1) )
-                     DO jj = ijpt0, ijpt1
-                        DO ji = iipt0, iipt1
-                           DO jn = 1, 3
-                              zcoefr(ji,jj,jn) = MAX( zrup2, zcoefr(ji,jj,jn) )
-                           END DO
-                           zcoefr(ji,jj,4) = MAX( zrup3, zcoefr(ji,jj,4) )
-                        END DO
-                     END DO
-                     ! ji-1,jj+1
-                     iipt0 = MAX( 1, MIN( ii-1,jpi  ) )
-                     iipt1 = MAX( 0, MIN( ii-1,jpi-1) )
-                     DO jj = ijpt0, ijpt1
-                        DO ji = iipt0, iipt1
-                           DO jn = 1, 3
-                              zcoefr(ji,jj,jn) = MAX( zrup2, zcoefr(ji,jj,jn) )
-                           END DO
-                           zcoefr(ji,jj,4) = MAX( zrup3, zcoefr(ji,jj,4) )
-                        END DO
-                     END DO
-                     ! ji-1,jj
-                     ijpt0 = MAX( 1, MIN( ij,jpj  ) )
-                     ijpt1 = MAX( 0, MIN( ij,jpj-1) )
-                     DO jj = ijpt0, ijpt1
-                        DO ji = iipt0, iipt1
-                           DO jn = 1, 3
-                              zcoefr(ji,jj,jn) = MAX( zrup2, zcoefr(ji,jj,jn) )
-                           END DO
-                           zcoefr(ji,jj,4) = MAX( zrup3, zcoefr(ji,jj,4) )
-                        END DO
-                     END DO
-                     ! ji-1,jj-1
-                     ijpt0 = MAX( 1, MIN( ij-1,jpj  ) )
-                     ijpt1 = MAX( 0, MIN( ij-1,jpj-1) )
-                     DO jj = ijpt0, ijpt1
-                        DO ji = iipt0, iipt1
-                           DO jn = 1, 3
-                              zcoefr(ji,jj,jn) = MAX( zrup2, zcoefr(ji,jj,jn) )
-                           END DO
-                           zcoefr(ji,jj,4) = MAX( zrup3, zcoefr(ji,jj,4) )
-                        END DO
-                     END DO
-                     ! --> zrup3 (originally = 0.25)
-                     ! ji-1,jj+2
-                     ijpt0 = MAX( 1, MIN( ij+2,jpj  ) )
-                     ijpt1 = MAX( 0, MIN( ij+2,jpj-1) )
-                     DO jj = ijpt0, ijpt1
-                        DO ji = iipt0, iipt1
-                           DO jn = 1, 3
-                              zcoefr(ji,jj,jn) = MAX( zrup3, zcoefr(ji,jj,jn) )
-                           END DO
-                        END DO
-                     END DO
-                     ! ji-1,jj-2
-                     ijpt0 = MAX( 1, MIN( ij-2,jpj  ) )
-                     ijpt1 = MAX( 0, MIN( ij-2,jpj-1) )
-                     DO jj = ijpt0, ijpt1
-                        DO ji = iipt0, iipt1
-                           DO jn = 1, 3 
-                              zcoefr(ji,jj,jn) = MAX( zrup3, zcoefr(ji,jj,jn) )
-                           END DO
-                        END DO
-                     END DO
-                     ! ji,jj-2
-                     iipt0 = MAX( 1, MIN( ii,jpi  ) )
-                     iipt1 = MAX( 0, MIN( ii,jpi-1) )
-                     DO jj = ijpt0, ijpt1
-                        DO ji = iipt0, iipt1
-                           DO jn = 1, 3 
-                              zcoefr(ji,jj,jn) = MAX( zrup3, zcoefr(ji,jj,jn) )
-                           END DO
-                        END DO
-                     END DO
-                     ! ji,jj+2
-                     ijpt0 = MAX( 1, MIN( ij+2,jpj  ) )
-                     ijpt1 = MAX( 0, MIN( ij+2,jpj-1) )
-                     DO jj = ijpt0, ijpt1
-                        DO ji = iipt0, iipt1
-                           DO jn = 1, 3
-                              zcoefr(ji,jj,jn) = MAX( zrup3, zcoefr(ji,jj,jn) )
-                           END DO
-                        END DO
-                     END DO
-                     ! ji+1,jj-2
-                     iipt0 = MAX( 1, MIN( ii+1,jpi  ) )
-                     iipt1 = MAX( 0, MIN( ii+1,jpi-1) )
-                     DO jj = ijpt0, ijpt1
-                        DO ji = iipt0, iipt1
-                           DO jn = 1, 3 
-                              zcoefr(ji,jj,jn) = MAX( zrup3, zcoefr(ji,jj,jn) )
-                           END DO
-                        END DO
-                     END DO
-                     ! ji+1,jj+2
-                     ijpt0 = MAX( 1, MIN( ij+2,jpj  ) )
-                     ijpt1 = MAX( 0, MIN( ij+2,jpj-1) )
-                     DO jj = ijpt0, ijpt1
-                        DO ji = iipt0, iipt1
-                           DO jn = 1, 3 
-                              zcoefr(ji,jj,jn) = MAX( zrup3, zcoefr(ji,jj,jn) )
-                           END DO
-                        END DO
-                     END DO
-                     ! ji+2,jj-2
-                     iipt0 = MAX( 1, MIN( ii+2,jpi  ) )
-                     iipt1 = MAX( 0, MIN( ii+2,jpi-1) )
-                     DO jj = ijpt0, ijpt1
-                        DO ji = iipt0, iipt1
-                           DO jn = 1, 3 
-                              zcoefr(ji,jj,jn) = MAX( zrup3, zcoefr(ji,jj,jn) )
-                           END DO
-                        END DO
-                     END DO
-                     ! ji+2,jj+2
-                     ijpt0 = MAX( 1, MIN( ij+2,jpj  ) )
-                     ijpt1 = MAX( 0, MIN( ij+2,jpj-1) )
-                     DO jj = ijpt0, ijpt1
-                        DO ji = iipt0, iipt1
-                           DO jn = 1, 3  
-                              zcoefr(ji,jj,jn) = MAX( zrup3, zcoefr(ji,jj,jn) )
-                           END DO
-                        END DO
-                     END DO
-                     ! ji+2,jj-1
-                     ijpt0 = MAX( 1, MIN( ij-1,jpj  ) )
-                     ijpt1 = MAX( 0, MIN( ij-1,jpj-1) )
-                     DO jj = ijpt0, ijpt1
-                        DO ji = iipt0, iipt1
-                           DO jn = 1, 3
-                              zcoefr(ji,jj,jn) = MAX( zrup3, zcoefr(ji,jj,jn) )
-                           END DO
-                        END DO
-                     END DO
-                     ! ji+2,jj+1
-                     ijpt0 = MAX( 1, MIN( ij+1,jpj  ) )
-                     ijpt1 = MAX( 0, MIN( ij+1,jpj-1) )
-                     DO jj = ijpt0, ijpt1
-                        DO ji = iipt0, iipt1
-                           DO jn = 1, 3 
-                              zcoefr(ji,jj,jn) = MAX( zrup3, zcoefr(ji,jj,jn) )
-                           END DO
-                        END DO
-                     END DO
-                     ! ji+2,jj
-                     ijpt0 = MAX( 1, MIN( ij,jpj  ) )
-                     ijpt1 = MAX( 0, MIN( ij,jpj-1) )
-                     DO jj = ijpt0, ijpt1
-                        DO ji = iipt0, iipt1
-                           DO jn = 1, 3
-                              zcoefr(ji,jj,jn) = MAX( zrup3, zcoefr(ji,jj,jn) )
-                           END DO
-                        END DO
-                     END DO
-                     ! ji-2,jj
-                     iipt0 = MAX( 1, MIN( ii-2,jpi  ) )
-                     iipt1 = MAX( 0, MIN( ii-2,jpi-1) )
-                     DO jj = ijpt0, ijpt1
-                        DO ji = iipt0, iipt1
-                           DO jn = 1, 3        
-                              zcoefr(ji,jj,jn) = MAX( zrup3, zcoefr(ji,jj,jn) )
-                           END DO
-                        END DO
-                     END DO
-                     ! ji-2,jj-2
-                     ijpt0 = MAX( 1, MIN( ij-2,jpj  ) )
-                     ijpt1 = MAX( 0, MIN( ij-2,jpj-1) )
-                     DO jj = ijpt0, ijpt1
-                        DO ji = iipt0, iipt1
-                           DO jn = 1, 3 
-                              zcoefr(ji,jj,jn) = MAX( zrup3, zcoefr(ji,jj,jn) )
-                           END DO
-                        END DO
-                     END DO
-                     ! ji-2,jj+2
-                     ijpt0 = MAX( 1, MIN( ij+2,jpj  ) )
-                     ijpt1 = MAX( 0, MIN( ij+2,jpj-1) )
-                     DO jj = ijpt0, ijpt1
-                        DO ji = iipt0, iipt1
-                           DO jn = 1, 3
-                              zcoefr(ji,jj,jn) = MAX( zrup3, zcoefr(ji,jj,jn) )
-                           END DO
-                        END DO
-                     END DO
-                     ! ji-2,jj-1
-                     ijpt0 = MAX( 1, MIN( ij-1,jpj  ) )
-                     ijpt1 = MAX( 0, MIN( ij-1,jpj-1) )
-                     DO jj = ijpt0, ijpt1
-                        DO ji = iipt0, iipt1
-                           DO jn = 1, 3 
-                              zcoefr(ji,jj,jn) = MAX( zrup3, zcoefr(ji,jj,jn) )
-                           END DO
-                        END DO
-                     END DO
-                     ! ji-2,jj+1
-                     ijpt0 = MAX( 1, MIN( ij+1,jpj  ) )
-                     ijpt1 = MAX( 0, MIN( ij+1,jpj-1) )
-                     DO jj = ijpt0, ijpt1
-                        DO ji = iipt0, iipt1
-                           DO jn = 1, 3 
-                              zcoefr(ji,jj,jn) = MAX( zrup3, zcoefr(ji,jj,jn) )
-                           END DO
-                        END DO
-                     END DO
-                  END DO
-               ENDIF
-            END DO
-            upsrnfh(:,:) = zcoefr(:,:,1)
-            upsrnfz(:)   = 0.e0
-            upsrnfz(1)   = 1.0
-            upsrnfz(2)   = 1.0
-            upsrnfz(3)   = zrup1
-            upsrnfz(4)   = zrup2
-            upsrnfz(5)   = zrup3
-
-            ! Reduces the sea surface restoring term to Levitus98
-            ! in the Sound Strait
-            DO jj = mj0(116), mj1(116)
-               DO ji = mi0(144), mi1(144)
-                  upsrnfh(ji,jj) = 0.25
-               END DO
-               DO ji = mi0(145), mi1(147)
-                  upsrnfh(ji,jj) = 0.5
-               END DO
-               DO ji = mi0(148), mi1(148)
-                  upsrnfh(ji,jj) = 0.25
-               END DO
-            END DO
-         ENDIF
-      ENDIF
-      
-      ! Upstream and centered scheme in the vicinity of some straits
-
-      IF( kt == nit000 ) THEN
-
-         ! ORCA_R2 configuration :
-         ! Gibraltar Strait
-         DO jj = mj0(102), mj1(102)
-            DO ji = mi0(138), mi1(138)
-               upsadv(ji,jj) = 0.2
-            END DO
-         END DO
-         DO jj = mj0(102), mj1(102)
-            DO ji = mi0(139), mi1(139)
-               upsadv(ji,jj) = 0.4
-            END DO
-         END DO
-         DO jj = mj0(102), mj1(102)
-            DO ji = mi0(140), mi1(140)
-               upsadv(ji,jj) = 0.5
-            END DO
-         END DO
-         DO jj = mj0(101), mj1(102)
-            DO ji = mi0(141), mi1(141)
-               upsadv(ji,jj) = 0.5
-            END DO
-         END DO
-         ! Bab el Mandeb Strait
-         DO jj = mj0(87), mj1(88)
-            DO ji = mi0(164), mi1(164)
-               upsadv(ji,jj) = 0.1
-            END DO
-         END DO
-         DO jj = mj0(88), mj1(88)
-            DO ji = mi0(163), mi1(163)
-               upsadv(ji,jj) = 0.25
-            END DO
-         END DO
-         DO jj = mj0(88), mj1(88)
-            DO ji = mi0(162), mi1(162)
-               upsadv(ji,jj) = 0.4
-            END DO
-         END DO
-         DO jj = mj0(88), mj1(88)
-            DO ji = mi0(160), mi1(161)
-               upsadv(ji,jj) = 0.5
-            END DO
-         END DO
-         DO jj = mj0(89), mj1(89)
-            DO ji = mi0(158), mi1(160)
-               upsadv(ji,jj) = 0.25
-            END DO
-         END DO
-         DO jj = mj0(90), mj1(90)
-            DO ji = mi0(160), mi1(160)
-               upsadv(ji,jj) = 0.25
-            END DO
-         END DO
-         ! Sound Strait
-         DO jj = mj0(116), mj1(116)
-            DO ji = mi0(145), mi1(147)
-               upsadv(ji,jj) = 0.5
-            END DO
-         END DO
-
-      ENDIF
Index: trunk/NEMO/OPA_SRC/SBC/sbcrnf_ORCA_R4.h90
===================================================================
--- trunk/NEMO/OPA_SRC/SBC/sbcrnf_ORCA_R4.h90	(revision 1122)
+++ 	(revision )
@@ -1,1339 +1,0 @@
-   !!----------------------------------------------------------------------
-   !!                    ***  flxrnf_ORCA_R4.h90  ***
-   !!----------------------------------------------------------------------
-
-   !!----------------------------------------------------------------------
-   !!                  OLD ORCA R4 runoff specification
-   !!                      kept here for memory
-   !!                         ( 20/04/2000 )
-   !!----------------------------------------------------------------------
-   !!----------------------------------------------------------------------
-   !!  OPA 9.0 , LOCEAN-IPSL (2005) 
-   !! $ Id: $
-   !! Software governed by the CeCILL licence (modipsl/doc/NEMO_CeCILL.txt)
-   !!----------------------------------------------------------------------
-      !!----------------------------------------------------------------------
-      !! * Local declarations
-      INTEGER  :: kt
-      INTEGER  :: jpmois
-      INTEGER  ::   &
-         jpriv,    & ! jpriv : maximum river number
-         jpcoef      ! jpcoef: maximum point number
-      INTEGER, PARAMETER ::   jpmois=12, jpriv=120, jpcoef=150
-
-      INTEGER  ::   ji, jj, jl, jm, jn, jr, jc
-      INTEGER  ::   i15, imois, iman
-      INTEGER  ::   inbp, idbd, idmeom
-      INTEGER  ::   iipt0, iipt1, ijpt0, ijpt1, ii, ij
-      INTEGER  ::   iirnf(jpcoef,jpriv),ijrnf(jpcoef,jpriv)
-      INTEGER  ::   inb(jpriv)
-
-      REAL(wp) ::   zconv, zxy
-      REAL(wp) ::   zrup1,zrup2,zrup3
-      REAL(wp) ::   zrnfm(jpmois, jpriv)
-      REAL(wp) ::   zcoefr(jpi, jpj, 5)
-      CHARACTER (len=20) ::   criver(jpriv)
-
-      REAL(wp) ::   zbilan , zero
-
-      COMMON/zzzrnf/ zrnfm, iirnf, ijrnf, inb
-
-      ! ====
-      ! DATA
-      ! ====
-
-      zero = 0.0e0
-      IF( kt == nit000 ) THEN 
-          criver = "                    "
-          inb = 0
-          iirnf = 0
-          ijrnf = 0
-          zrnfm  = zero
-          runoff = zero
-         
-          ! I NORTH AMERICA
-          ! ---------------
-          
-          ! 1-Saint Lawrence (Canada) 49N30 67W00           R2
-          ! Old=447 km3/year=14165m3/s; CLIPPER=12430m3/s; ratio=1.14
-          ! --> 13975.249m3/s
-          criver(1) = "Saint Lawrence      "
-          inb(1) = 1
-          iirnf(1:1,1) = 57
-          ijrnf(1:1,1) = 55
-          zrnfm(1:jpmois,1) = (/   &
-              11718.,  11971.,  12466.,  14980.,  19257.,  16913.,   &
-              14211.,  13119.,  13367.,  13644.,  13463.,  12594./)
-          ! 2-Hudson River (USA) 41N00 77W00                R1
-          ! Old=12 km3/year=380m3/s; UNESCO(64-83,p255)=408m3/s
-          ! --> 407.416687m3/s
-          criver(2) = "Hudson River        "
-          inb(2) = 1
-          iirnf(1:1,2) = 54
-          ijrnf(1:1,2) = 53
-          zrnfm(1:jpmois,2) = (/   &
-              364.,    404.,    648.,    873.,    593.,    301.,   &
-              205.,    171.,    205.,    290.,    379.,    456./)
-          ! 3-Missisipi and Atchafalaya (USA) 29N 90W       R1 
-          ! Old=580 km3/year=18380m3/s; CLIPPER=17599m3/s
-          ! --> 17599.334m3/s
-          criver(3) = "Missisipi           "
-          inb(3) = 1
-          iirnf(1:1,3) = 50
-          ijrnf(1:1,3) = 49
-          zrnfm(1:jpmois,3) = (/   &
-              19260.,  20813.,  23026.,  28160.,  26510.,  19650.,   &
-              14181.,  10676.,   9578.,  10167.,  11699.,  17472./)
-          ! 4-Brazos (USA) 29N58 95W53                      R1
-          ! Old=7 km3/year=222m3/s; UNESCO(65-84, p298)=200m3/s
-          ! --> 199.666656m3/s
-          criver(4) = "Brazos              "
-          inb(4) = 1
-          iirnf(1:1,4) = 48
-          ijrnf(1:1,4) = 49
-          zrnfm(1:jpmois,4) = (/   &
-              165.,    224.,    200.,    266.,    457.,    327.,   &
-              143.,     69.,    102.,    132.,    160.,    151. /)
-          ! 5-Colorado (Mexico) 31N54 114W57                R5
-          ! Old=20 km3/year=634m3/s
-          ! --> 634.m3/s
-          criver(5) = "Colorado            "
-          inb(5) = 1
-          iirnf(1:1,5) = 43
-          ijrnf(1:1,5) = 50
-          zrnfm(1:jpmois,5) = 634. 
-          ! 6-Columbia (USA) 45N60 121W17                   R1
-          ! Old=251 km3/year=7954m3/s; UNESCO(65-84, p312)=5368m3/s
-          ! --> 5368.41602m3/s
-          criver(6) = "Columbia            "
-          inb(6) = 1
-          iirnf(1:1,6) = 40
-          ijrnf(1:1,6) = 56
-          zrnfm(1:jpmois,6) = (/   &
-              4752.,    5077.,    5633.,    6075.,    8422.,    9822.,   &
-              6350.,    4112.,    3273.,    3251.,    3548.,    4106. /)
-          ! 7-Fraser (USA) 47N 124W                         R5
-          ! Old=112 km3/year=3550m3/s;not found in UNESCO
-          ! --> 3549.99976m3/s
-          criver(7) = "Fraser              "
-          inb(7) = 1
-          iirnf(1:1,7) = 40
-          ijrnf(1:1,7) = 57
-          zrnfm(1:jpmois,7) = 3550. 
-          ! 8-Yukon (USA) 63N00 165W00                      R3
-          ! Old=195 km3/year=6179m3/s; UNESCO(65-77+80-84, p321)=4245; ratio=1.46
-          ! --> 6017.24951m3/s
-          criver(8) = "Yukon               "
-          inb(8) = 1
-          iirnf(1:1,8) = 30
-          ijrnf(1:1,8) = 62
-          zrnfm(1:jpmois,8) = (/   &
-              1459.,    1213.,    1065.,    1100.,    9067.,   17675.,   &
-              10680.,   11271.,    8658.,    5453.,    2649.,    1917./)
-          ! 9-Copper (USA) 61N47 144W45                     R2
-          ! Old=39 km3/year=1236m3/s; UNESCO(65-84, p317)=1052m3/s; ratio=1.18
-          ! --> 1236.25m3/s
-          criver(9) = "Copper              "
-          inb(9) = 1
-          iirnf(1:1,9) = 35
-          ijrnf(1:1,9) = 60
-          zrnfm(1:jpmois,9) = (/   &
-              219.,    195.,    181.,    214.,    904.,   2532.,   &
-              4026.,   3627.,   1632.,    697.,    354.,    254. /) 
-          ! 10-Susitna (USA)   62N 150E                     R5
-          ! Old=40 km3/year=1268m3/s
-          ! --> 1268.m3/s
-          criver(10) = "Susitna             "
-          inb(10) = 1
-          iirnf(1:1,10) = 34
-          ijrnf(1:1,10) = 60
-          zrnfm(1:jpmois,10) = 1268. 
-          ! 11-McKenzie (Canada) 70N00 135W00               R3
-          ! Old=306 km3/year=9697m3/s; UNESCO(67+72+7479+82-84, p168)=8343m3/s;
-          ! ratio=1.16
-          ! --> 9606.66797m3/s
-          criver(11) = "McKenzie            "
-          inb(11) = 1
-          iirnf(1:1,11) = 35
-          ijrnf(1:1,11) = 66
-          zrnfm(1:jpmois,11) = (/   &
-              3644.,    3400.,    3258.,    3546.,   14188.,   21544.,   &
-              18954.,   15487.,   11284.,   10169.,    5901.,    3905./)
-          ! 12-Kuksokwim (USA) 60N00 ~162W30                R4
-          ! UNESCO(65-84, p318)=1112m3/s (61N87 158E10, Crooked Creek, inland)
-          ! --> 1111.99988m3/s
-          criver(12) = "Kuksokwim           "
-          inb(12) = 1
-          iirnf(1:1,12) = 30
-          ijrnf(1:1,12) = 60
-          zrnfm(1:jpmois,12) = (/   &
-              335.,     308.,     295.,     407.,    2415.,    2027.,   &
-              1917.,   2049.,    1585.,    1072.,     545.,     389./)
-          ! 13-Nelson
-          criver(13) = "Nelson              "
-          inb(13) = 1
-          iirnf(1:1,13) = 51
-          ijrnf(1:1,13) = 65
-          zrnfm(1:jpmois,13) = 0.
-
-          ! Total run-offs(I)=60974.2539m3/s
-
-          ! II SOUTH AMERICA
-          ! -----------------
-          ! 16-Chira (Peru) 4s54 81w08                      R2
-          ! Old=5 km3/year=158m3/s; UNESCO (65-68, p357)=80.9; ratio=1.95
-          ! --> 158.000015m3/s
-          criver(16) = "Chira               "
-          inb(16) = 1
-          iirnf(1:1,16) = 51
-          ijrnf(1:1,16) = 34
-          zrnfm(1:jpmois,16) = (/   &
-              70.,     105.,     564.,     582.,     230.,      94.,   &
-              87.,      47.,      35.,      32.,      29.,      21./)
-          ! 17-Magdalena (Colombia) 11N06 74W51             R2
-          ! Old=237 km3/year=7510m3/s; UNESCO(71-73+76-79, p348)=6973m3/s;
-          ! ratio=1.08
-          ! --> 7451.58301m3/s
-          criver(17) = "Magdelena           "
-          inb(17) = 1
-          iirnf(1:1,17) = 53
-          ijrnf(1:1,17) = 45
-          zrnfm(1:jpmois,17) = (/   &
-              6598.,    4441.,    4087.,    5445.,    7380.,    8067.,   &
-              8509.,    7428.,    7631.,    9296.,   10588.,    9949./)
-          ! 18-Orinoco (Venezuela) 8N37 62W15!              R1
-          ! Old=1100 km3/year=34857m3/s; UNESCO(73-75, p364)=28736m3/s;
-          ! HLA=26371m3/s
-          ! Data from HLA: 8.1N - 63.3W   1925/89
-          ! --> 26371.416
-          criver(18) = "Orinico             "
-          inb(18) = 1
-          iirnf(1:1,18) = 57
-          ijrnf(1:1,18) = 44
-          zrnfm(1:jpmois,18) = (/   &
-              10670.,    7513.,    6344.,    6887.,   13086.,   26699.,   &
-              43938.,   59957.,   57435.,   39951.,   27410.,   16567./)
-          ! 19-Amazon (Brazil) 0S05 50W00                   R1
-          ! Old=6300 km3/year=199634m3/s; CLIPPER=166192m3/s~UNESCO;
-          ! HLA=184806m3/s
-          ! UNESCO data are at Obidos which is far inland and therefore does not
-          ! fully represent the actual discharge
-          ! Data from HLA: 1.5N - 55.3W   1970/96
-          ! --> 184845.484m3/s
-          criver(19) = "Amazon              "
-          inb(19) = 3
-          iirnf(1:3,19) = (/ 60, 60, 60 /)
-          ijrnf(1:3,19) = (/ 39, 40, 41 /)
-          zrnfm(1:jpmois,19) = (/   &
-              152233., 176415., 199941., 225565., 244016., 244861.,   &
-              230711., 203128., 164531., 128486., 118214., 130045. /)
-          ! 20-Sao Francisco (Brazil) 10S30 36W24           R1
-          ! Old=97 km3/year=3074m3/s; UNESCO(77-79, p338)=2645m3/s
-          ! --> 2645.3335m3/s
-          criver(20) = "San Fransisco       "
-          inb(20) = 1
-          iirnf(1:1,20) = 63
-          ijrnf(1:1,20) = 31
-          zrnfm(1:jpmois,20) = (/   &
-              3244.,   4116.,   5139.,   4330.,   2323.,   1800.,     &
-              1587.,   1511.,   1593.,   1740.,   2190.,   2171. /)
-          ! 21-Parana+Uruguay(Arg./Uru.) 34S12 58W18        R1 
-          ! Old=470 km3/year=14893m3/s; UNESCO(p323/324)=22263m3/s
-          ! --> 22203.168m3/s 
-          criver(21) = "Parana              "
-          inb(21) = 1
-          iirnf(1:1,21) = 58
-          ijrnf(1:1,21) = 24
-          zrnfm(1:jpmois,21) = (/   &
-              24803.,   26836.,   24585.,   21358.,   20041.,   21689.,   &
-              22647.,   20403.,   19627.,   21400.,   21126.,   21923./)
-          ! 22-Negro (Argentina)  41S02  62W47              R2
-          ! Old=30 km3/year=951m3/s; UNESCO(65-79,p330)=765m3/s; ratio=1.24
-          ! --> 907.666687m3/s
-          criver(22) = "Rio Negro           "
-          inb(22) = 1
-          iirnf(1:1,22) = 56
-          ijrnf(1:1,22) = 22
-          zrnfm(1:jpmois,22) = (/   &
-              882.,     587.,     519.,     482.,     623.,     984.,   &
-              1120.,   1294.,    1092.,    1166.,    1105.,    1038./)
-          ! 23- Tocantins :  ~48W00 ~0.5S                   R4
-          ! UNESCO(76-77+80-81, p335)=12167m3/s (station far inland!)
-          ! --> 12173.082m3/s
-          criver(23) = "Tocantins           "
-          inb(23) = 1
-          iirnf(1:1,23) = 60
-          ijrnf(1:1,23) = 38
-          zrnfm(1:jpmois,23) = (/   &
-              15507.,  24855.,  27854.,  24257.,  15860.,   7861.,   &
-              4440.,   3082.,   2728.,   3417.,   5511.,  10705. /)
-
-          ! Total run-offs(II)=256755.75
-          ! Total run-offs(I+II)=317730
-
-          ! III EUROPE
-          ! ----------
-          ! 26-Rhone (France) 43N20 4E50                    R1 
-          ! Old=49 km3/year=1553m3/s; UNESCO(65-83, p512)=1229
-          ! --> 1217.16663m3/s
-          criver(26) = "Rhone               "
-          inb(26) = 1
-          iirnf(1:1,26) = 76
-          ijrnf(1:1,26) = 55
-          zrnfm(1:jpmois,26) = (/   &
-              1385.,    1848.,    1568.,    1311.,    1420.,    1191.,   &
-              937.,     796.,     794.,     962.,    1161.,    1233./)
-          ! 27-Po (Italy) 41N10 11E09                       R1 
-          ! Old=46 km3/year=1458m3/s; UNESCO(65-84, p527)=1550m3/s
-          ! --> 1550.08337m3/s
-          criver(27) = "Po                  "
-          inb(27) = 1
-          iirnf(1:1,27) = 79
-          ijrnf(1:1,27) = 54
-          zrnfm(1:jpmois,27) = (/   &
-              1304.,    1605.,    1856.,    1549.,    1958.,    1760.,   &
-              1040.,    1007.,    1615.,    1907.,    1772.,    1228./)
-          ! 28-Danube (Roumania) 45N20 29E40                R1 
-          ! Old=206 km3/year=6527m3/s; UNESCO(65-84, p544)=7011m3/s
-          ! --> 7011.8335m3/s
-          criver(28) = "Danube              "
-          inb(28) = 1
-          iirnf(1:1,28) = 87
-          ijrnf(1:1,28) = 53
-          zrnfm(1:jpmois,28) = (/   &
-              6625.,    7352.,    7986.,    9279.,    9482.,    8618.,   &
-              7586.,    6001.,    5210.,    4989.,    4921.,    6093./)
-
-          ! 29-Volga
-          criver(29) = "Volga               "
-          inb(29) = 1
-          iirnf(1:1,29) = 4
-          ijrnf(1:1,29) = 56
-          zrnfm(1:jpmois,29) = 0.
-          ! 30-Dniepr
-          criver(30) = "Dniepr+Don          "
-          inb(30)= 2
-          iirnf(1:2,30) = (/ 89, 90 /)
-          ijrnf(1:2,30) = (/ 56, 56 /)
-
-          ! Total run-offs(III)=9779.08301m3/s
-          ! Total run-offs(I+II+III)=327509.094m3/s
-
-          ! IV EURASIAN ARCTIC
-          ! ------------------
-          ! 31-Iana (Russia)      126E40 72N25              R1
-          ! Old=29 km3/year=919m3/s; UNESCO(65-72+74+76-84, p463)=936m3/s
-          ! --> 0m3/s
-          criver(31) = "Iana                "
-          inb(31) = 1
-          iirnf(1:1,31) = 20
-          ijrnf(1:1,31) = 68
-          zrnfm(1:jpmois,31) = (/   &
-              2.,       1.,       0.,       0.,     917.,    3514.,   &
-              2968.,    2329.,    1322.,   178.,      36.,       8./)
-          ! 32-Ob (Russia) 66N45 69E30                      R1 
-          ! Old=385 km3/year=12200m3/s; UNESCO(65-84, p478)=12862m3/s
-          ! --> 12861.5m3/s
-          criver(32) = "Ob                  "
-          inb(32) = 1
-          iirnf(1:1,32) = 78
-          ijrnf(1:1,32) = 72
-          zrnfm(1:jpmois,32) = (/   &
-              4977.,    4150.,    3611.,    3590.,   14723.,   33366.,   &
-              30773.,   22785.,   14701.,    9705.,    6078.,    5879./)
-          ! 33-Yenesei (Russia) 71N50 82E40                 R1 
-          ! Old=560 km3/year=17745m3/s; UNESCO(65-84, p472)=17462m3/s
-          ! --> 17461.666m3/s
-          criver(33) = "Yenesei             "
-          inb(33) = 1
-          iirnf(1:1,33) = 77
-          ijrnf(1:1,33) = 75
-          zrnfm(1:jpmois,33) = (/   &
-              6659.,    7007.,    6895.,    7089.,   24183.,   74965.,   &
-              25325.,   15940.,   16050.,   12678.,    6741.,    6008./)
-          ! 34-Severnaya Dvina(Russia)64N32 40E30           R5 
-          ! Old=106 km3/year=3359m3/s; not in UNESCO data base
-          ! --> 3359.00024m3/s
-          criver(34) = "Severnaya Dvina     "
-          inb(34) = 1
-          iirnf(1:1,34) = 78
-          ijrnf(1:1,34) = 68
-          zrnfm(1:jpmois,34) = 3359.
-          ! 35-Lena (Russia) 72N25 126E40                   R1 
-          ! --> 16192.333m3/s
-          criver(35) = "Lena                "
-          inb(35) = 1
-          iirnf(1:1,35) = 20
-          ijrnf(1:1,35) = 68
-          zrnfm(1:jpmois,35) = (/   &
-              2663.,    2095.,    1629.,    1324.,    7203.,   64230.,   &
-              40810.,   28310.,   25995.,   13299.,    3673.,    3077./)
-          ! 36-Kolyma (Russia) 69N30 161E00                 R1
-          ! Old=71 km3/year=2250m3/s; UNESCO(65-84)=2138m3/s
-          ! --> 2138.66675m3/s
-          criver(36) = "Kolyma              "
-          inb(36) = 1
-          iirnf(1:1,36) = 24
-          ijrnf(1:1,36) = 65
-          zrnfm(1:jpmois,36) = (/   &
-              131.,      92.,      76.,      65.,    2576.,    9143.,   &
-              4670.,    3727.,    3226.,    1119.,     630.,     209./)
-          ! 37-Indigirka (russia) 70N48 148E54              R1 
-          ! Old=55 km3/year=1743m3/s; UNESCO(65-68+70+72-84)=1622.
-          ! --> 1622.00012m3/s
-          criver(37) = "Indigirka           "
-          inb(37) = 1
-          iirnf(1:1,37) = 23
-          ijrnf(1:1,37) = 66
-          zrnfm(1:jpmois,37) = (/   &
-              32.,      17.,      11.,       7.,     354.,    5629.,   &
-              5210.,    4617.,    2823.,     559.,     142.,      63./)
-          ! 38-Khatanga
-          criver(38) = "Khatanga             "
-          inb(38) = 1
-          iirnf(1:1,38) = 20
-          ijrnf(1:1,38) = 72
-          zrnfm(1:jpmois,38) = 0.
-
-          ! Total run-offs(IV)=40773.6641m3/s
-          ! Total run-offs(I+II+III+IV)=368282.75m3/s
-
-          ! V ASIA
-          ! ------
-          ! 41-Amur (Russia) 52N56 141E10                   R2
-          ! Old=325 km3/year=10299m3/s; UNESCO(65-84, p451)=8792m3/s; ratio=1.17
-          ! --> 10299.332m3/s
-          criver(41) = "Amur                "
-          inb(41) = 1
-          iirnf(1:1,41) = 18
-          ijrnf(1:1,41) = 59
-          zrnfm(1:jpmois,41) = (/   &
-              2108.,    1543.,    1252.,    4088.,   15464.,   17336.,   &
-              16100.,   20801.,   21548.,   15658.,    5188.,    2506./)
-          ! 42-Liaohe (China) 40N40 122E09                  R5
-          ! Old=6 km3/year=190m3/s; not found in UNESCO data base
-          ! --> 189.999985m3/s
-          criver(42) = "Liaohe              "
-          inb(42) = 1
-          iirnf(1:1,42) = 15
-          ijrnf(1:1,42) = 54
-          zrnfm(1:jpmois,42) = 190.
-          ! 43-Daling (China) 40N56 121E43                  R5
-          ! Old=1 km3/year=32m3/s; not found in UNESCO data base
-          ! --> 31.9999962m3/s
-          criver(43) = "Daling              "
-          inb(43) = 1
-          iirnf(1:1,43) = 15
-          ijrnf(1:1,43) = 54
-          zrnfm(1:jpmois,43) = 32.
-          ! 44-Haiho (China) 6N17   1E23
-          ! Old=2 km3/year=64m3/s; not found in UNESCO data base
-          ! --> 63.9999924m3/s
-          criver(44) = "Haiho               "
-          inb(44) = 0
-          iirnf(1:1,44) = 0
-          ijrnf(1:1,44) = 0
-          zrnfm(1:jpmois,44) = 64.
-          ! 45-Yellow/Huangho (China) 37N32 118E19          R2 
-          ! Old=49 km3/year=1553m3/s; UNESCO(p371)=1210m3/s; ratio=1.28
-          ! --> 1552.5m3/s
-          criver(45) = "Yellow/Huangho      "
-          inb(45) = 1
-          iirnf(1:1,45) = 13
-          ijrnf(1:1,45) = 51
-          zrnfm(1:jpmois,45) = (/   &
-              697.,     674.,    1282.,    1206.,    1109.,     778.,   &
-              1964.,    3270.,    3668.,    2133.,     978.,     871./)
-          ! 46-Yangtze (China) 31N48 121E10=Changjiang      R2
-          ! Old=900 km3/year=28520m3/s; UNESCO(p372)=25032m3/s; ratio=1.14
-          ! --> 28519.584m3/s
-          criver(46) = "Yangtze             "
-          inb(46) = 1
-          iirnf(1:1,46) = 13
-          ijrnf(1:1,46) = 50
-          zrnfm(1:jpmois,46) = (/   &
-              10994.,   10548.,   13254.,   24097.,   36972.,   47197.,   &
-              53891.,   42070.,   37627.,   30762.,   22046.,   12777./)
-          ! 47-Pearl/Zhujiang (China) 23N36 106E34=Xijiang  R2 
-          ! Old=302 km3/year=9570m3/s; UNESCO(76,82, p373)=6987m3/s; ratio=1.37
-          ! --> 9570.75m3/s
-          criver(47) = "Pearl Zhujiang       "
-          inb(47) = 1
-          iirnf(1:1,47) = 12
-          ijrnf(1:1,47) = 48
-          zrnfm(1:jpmois,47) = (/   &
-              2357.,    2415.,    2634.,    6446.,   14246.,   18395.,   &
-              19213.,   17558.,   12312.,    7255.,    5864.,    6154./)
-          ! 48-Red/Hungho (Vietnam) 20N17 106E34            R5 
-          ! Old=123 km3/year=3897m3/s; not in UNESCO data base
-          ! --> 3897.m3/s
-          criver(48) = "Red Hungho          "
-          inb(48) = 1
-          iirnf(1:1,48) = 10
-          ijrnf(1:1,48) = 47
-          zrnfm(1:jpmois,48) = 3897.
-          ! 49-Mekong (Vietnam) 10N33 105E24                R5
-          ! Old=470 km3/year=14893m3/s; not in UNESCO data base
-          ! --> 14893.0029m3/s
-          criver(49) = "Mekong              "
-          inb(49) = 1
-          iirnf(1:1,49) =  8
-          ijrnf(1:1,49) = 44
-          zrnfm(1:jpmois,49) = 14893.
-          ! 50-Irrawady (Burma) 15N50  95E06                R5
-          ! Old=428 km3/year=13563m3/s; not in UNESCO data base
-          ! --> 13563.m3/s
-          criver(50) = "Irrawady            "
-          inb(50) = 1
-          iirnf(1:1,50) =  6
-          ijrnf(1:1,50) = 46
-          zrnfm(1:jpmois,50) = 13563.
-          ! 51-Ganges+Brahmaputra (Beng.) 22N00  91E00      R1
-          ! Old=971 km3/year=30769m3/s; UNESCO(69-70+73-75, p367)=31760m3/s
-          ! --> 32147.498m3/s
-          criver(51) = "Ganges              "
-          inb(51) = 2
-          iirnf(1:2,51) = (/  4,  5 /)
-          ijrnf(1:2,51) = (/ 47, 47 /)
-          zrnfm(1:jpmois,51) = (/   &
-              6623.,    6315.,    6432.,    9410.,   17263.,   38302.,   &
-              64688.,   80338.,   84802.,   43387.,   17888.,   10322./)
-          ! 52-Mehandi (India)                              R5
-          ! Old=67 km3/year=2123m3/s; not in UNESCO data base
-          ! --> 0m3/s
-          criver(52) = "Mehandi             "
-          inb(52) = 0
-          iirnf(1:1,52) = 0
-          ijrnf(1:1,52) = 0
-          zrnfm(1:jpmois,52) = 2123.
-          ! 53-Damodar (India)                              R2
-          ! Old=10 km3/year=320m3/s; UNESCO(p386)=173m3/s; ratio=1.85
-          ! --> 432.916656m3/s
-          criver(53) = "Damodar             "
-          inb(53) = 0
-          iirnf(1:1,53) = 0
-          ijrnf(1:1,53) = 0
-          zrnfm(1:jpmois,53) = (/   &
-              51.,      37.,      42.,      49.,      85.,     296.,   &
-              896.,    1390.,    1591.,     542.,     172.,      44./)  
-          ! 54-Godavari (India) 17N00  81E45                R2
-          ! Old=84 km3/year=2662m3/s; UNESCO(p379)=1916m3/s; ratio=1.39
-          ! --> 3670.5m3/s
-          criver(54) = "Godaravi            "
-          inb(54) = 1
-          iirnf(1:1,54) =  3
-          ijrnf(1:1,54) = 46
-          zrnfm(1:jpmois,54) = (/   &
-              988.,     737.,     280.,     219.,     196.,    1691.,   &
-              8063.,   14571.,   11428.,    3953.,    1221.,     699./)
-          ! 55-Indus (Pakistan) 24N20  67E47                R2
-          ! Old=238 km3/year=7542m3/s; UNESCO (76-79, p428)=2396m3/s; ratio=3.15
-          ! --> 7562.m3/s
-          criver(55) = "Indus               "
-          inb(55) = 1
-          iirnf(1:1,55) = 89
-          ijrnf(1:1,55) = 47
-          zrnfm(1:jpmois,55) = (/   &
-              988.,     737.,    1904.,    1968.,    3625.,    6143.,   &
-              15969.,    36656.,  17173.,    3215.,    1251.,    1115./)
-          ! 56-Tigris and Euphrates (Iraq) 31N00  47E25     R1
-          ! Old=46 km3/year=1458m3/s; UNESCO(76-79, p428)=2396m3/s
-          ! --> 2248.8335m3/s
-          criver(56) = "Tigris              "
-          inb(56) = 1
-          iirnf(1:1,56) = 88
-          ijrnf(1:1,56) = 47
-          zrnfm(1:jpmois,56) = (/   &
-              1872.,    2127.,    2962.,    4944.,    5036.,    3078.,   &
-              1362.,     900.,     786.,     993.,    1184.,    1742./)
-
-          ! Total run-offs (V)=128642.922m3/s
-          ! Total run-offs (I+II+III+IV+V)=496925.688m3/s
-
-          ! VI AFRICA
-          ! ----------
-          ! 61-Nile (Egypt) 30N10  31E06                    R1
-          ! Old=30 km3/year=951m3/s; UNESCO(73-84)=1251m3/s
-          ! --> 1251.33337m3/s
-          criver(61) = "Nile                "
-          inb(61) = 1
-          iirnf(1:1,61) = 90
-          ijrnf(1:1,61) = 50
-          zrnfm(1:jpmois,61) = (/   &
-              1239.,    1035.,    1126.,    1110.,    1170.,    1542.,   &
-              1742.,    1559.,    1191.,    1103.,    1075.,    1124./)
-          ! 62-Niger (Nigeria) 5N33   6E33                  R1
-          ! Old=192 km3/year=6084m3/s; HLA=5545m3/s;
-          ! Data from HLA: 4.5N -  6  E   1950/87
-          ! --> 5545.08301m3/s
-          criver(62) = "Niger               "
-          inb(62) = 1
-          iirnf(1:1,62) = 74
-          ijrnf(1:1,62) = 41
-          zrnfm(1:jpmois,62) = (/   &
-              2029.,    1804.,    1573.,    1418.,    1494.,    2561.,   &
-              5272.,    9629.,   14951.,   16227.,    6842.,    2741./)
-          ! 63-Zaire/Congo (Congo) 6S04  12E24              R1      
-          ! Old=1250 km3/year=39610m3/s; UNESCO(65-83, p149)=41941m3/s; HLA
-          ! =39720m3/s
-          ! Data from HLA: 6S - 12  E   1902/96
-          ! --> 39719.5m3/s
-          criver(63) = "Zaire Congo         "
-          inb(63) = 1
-          iirnf(1:1,63) = 75
-          ijrnf(1:1,63) = 33
-          zrnfm(1:jpmois,63) = (/   &
-              47200.,   37692.,   34374.,   36603.,   37862.,   35588.,   &
-              30966.,   30774.,   35981.,   42825.,   51455.,   55314./)
-          ! 64-Orange (Sth Africa) 28S41 16E28              R5
-          ! Old=11 km3/year=349m3/s
-          ! --> 349.000031m3/s
-          criver(64) = "Orange              "
-          inb(64) = 1
-          iirnf(1:1,64) = 76
-          ijrnf(1:1,64) = 26
-          zrnfm(1:jpmois,64) = 349.
-          ! 65-Zambesi (Mozambique) 18S55 36E04             R5
-          ! Old=223 km3/year=7066m3/s
-          ! --> 7066.00049m3/s
-          criver(65) = "Zambesi             "
-          inb(65) = 1
-          iirnf(1:1,65) = 81
-          ijrnf(1:1,65) = 29
-          zrnfm(1:jpmois,65) = 7066.
-          ! 66-Limpopo (mozambique) 25S15 33E30             R5
-          ! Old=5 km3/year=158m3/s
-          ! --> 158.m3/s
-          criver(66) = "Limpopo             "
-          inb(66) = 1
-          iirnf(1:1,66) = 81
-          ijrnf(1:1,66) = 28
-          zrnfm(1:jpmois,66) = 158.
-          ! 67-Rufiji (tanzania) 8S00 39E20                 R5
-          ! Old=9 km3/year=285m3/s
-          ! --> 285.m3/s
-          criver(67) = "Rufiji             "
-          inb(67) = 1
-          iirnf(1:1,67) = 82
-          ijrnf(1:1,67) = 33
-          zrnfm(1:jpmois,67) = 285.
-          ! Other rivers from the Gulf of Guinea
-          ! 68-Kouilou (Congo)  4.5S-11.5E (UNESCO=p51)     R1
-          ! Data from HLA:  4S - 11  E   1953/86 = 919m3/s
-          ! --> 924.25m3/s
-          criver(68) = "Kouilou            "
-          inb(68) = 1
-          iirnf(1:1,68) = 75
-          ijrnf(1:1,68) = 33
-          zrnfm(1:jpmois,68) = (/   &
-              1115.,    1166.,    1285.,    1462.,    1454.,     725.,   &
-              480.,     378.,     327.,     400.,    1018.,    1281./)
-          ! 69-Nyanga (Gabon)  3S-10E (UNESCO=p66)          R1
-          ! Data from HLA: 2.8N - 10  E   1952/82 = 264m3/s
-          ! --> 264.166687m3/s
-          criver(69) = "Nyanga              "
-          inb(69) = 1
-          iirnf(1:1,69) = 74
-          ijrnf(1:1,69) = 36
-          zrnfm(1:jpmois,69) = (/   &
-              126.,      99.,     143.,     253.,     316.,     300.,   &
-              156.,      93.,     231.,     557.,     592.,     304./)
-          ! 70-Ogoue (Gabon) 1S-9E (UNESCO, p69)            R1
-          ! Data from HLA: .75S -  9  E   1930/90 = 4512m3/s
-          ! --> 4511.83301m3/s
-          criver(70) = "Ogoue               "
-          inb(70) = 1
-          iirnf(1:1,70) = 74
-          ijrnf(1:1,70) = 37
-          zrnfm(1:jpmois,70) = (/   &
-              4514.,    3988.,    4622.,    5403.,    5790.,    4145.,   &
-              2616.,    1961.,    2119.,    4304.,    7870.,    6810./)
-          ! 71-Sanaga (Cameroun) 3.5N-9.5E (UNESCO, p31)    R1
-          ! Data from HLA: 3.9N -  9.6E   1944/80 = 1983m3/s
-          ! --> 1982.58337m3/s
-          criver(71) = "Sanaga              "
-          inb(71) = 1
-          iirnf(1:1,71) = 74
-          ijrnf(1:1,71) = 41
-          zrnfm(1:jpmois,71) = (/   &
-              813.,     578.,     517.,     633.,     909.,    1278.,   &
-              1928.,    2653.,    4361.,    5470.,    3307.,    1344./)
-          ! 72-Nyong (Cameroun) 3.1N-9.5E (UNESCO, p36)     R1
-          ! Data from HLA: 3.5N -  9.8E   1951/84 = 432m3/s
-          ! --> 431.166656m3/s      
-          criver(72) = "Nyong               "
-          inb(72) = 1
-          iirnf(1:1,72) = 74
-          ijrnf(1:1,72) = 41
-          zrnfm(1:jpmois,72) = (/   &
-              232.,     128.,     150.,     266.,     402.,     459.,   &
-              353.,     306.,     604.,     930.,     878.,     466./)
-          ! 73-Ntem (Cameroun) 2.5N-9.5E (UNESCO, p36)      R1
-          ! Data from HLA: 3N -  9.8E   1951/87 = 330m3/s      
-          ! --> 329.75m3/s
-          criver(73) = "Ntem                "
-          inb(73) = 1
-          iirnf(1:1,73) = 74
-          ijrnf(1:1,73) = 41
-          zrnfm(1:jpmois,73) = (/   &
-              452.,     455.,     492.,     466.,     482.,     251.,   &
-              151.,     108.,      87.,     133.,     395.,     485./)
-          ! 74-Wouri (Cameroun) 4N-9.5E (UNESCO, p34)       R1
-          ! Data from HLA: 4N -  9.5E   1951/84  = 307m3/s
-          ! --> 306.583374m3/s
-          criver(74) = "Wouri               "
-          inb(74) = 1
-          iirnf(1:1,74) = 74
-          ijrnf(1:1,74) = 41
-          zrnfm(1:jpmois,74) = (/   &
-              88.,      63.,      68.,      95.,     139.,     222.,   &
-              395.,     633.,     825.,     675.,     328.,     148./)
-          ! 75-Oueme (Benin) 6.5N-2.5E (UNESCO, p27)        R1
-          ! Data from HLA:  6.2N -  2.5E   1951/84 =113m3/s
-          ! --> 112.666664m3/s
-          criver(75) = "Oueme               "
-          inb(75) = 1
-          iirnf(1:1,75) = 73
-          ijrnf(1:1,75) = 42
-          zrnfm(1:jpmois,75) = (/   &
-              0.,       0.,       0.,       0.,       1.,      18.,   &
-              96.,     324.,     531.,     319.,      58.,       5./)
-          ! 76-Mono(Tepetou,Benin/Togo) 6.25N-1.8E (U, p27) R1
-          ! Data from HLA: 6.2N -  1  E   1951/85  = 96m3/s       
-          ! --> 96.4999924m3/s
-          criver(76) = "Mono                "
-          inb(76) = 1
-          iirnf(1:1,76) = 72
-          ijrnf(1:1,76) = 42
-          zrnfm(1:jpmois,76) = (/   &
-              2.,       1.,       1.,       4.,       8.,      44.,   &
-              168.,     274.,     386.,     215.,      47.,       8./)
-          ! 77-Comoe (Ivory Coast) 5N-4W (UNESCO, p56)      R1
-          ! Data from HLA: 5.1N -  4  W   1955/86 = 194m3/s       
-          ! --> 194.25m3/s
-          criver(77) = "Comoe               "
-          inb(77) = 1
-          iirnf(1:1,77) = 71
-          ijrnf(1:1,77) = 41
-          zrnfm(1:jpmois,77) = (/   &
-              14.,       7.,       9.,      17.,      29.,      90.,   &
-              145.,     349.,     740.,     688.,     198.,      45./)
-          ! 78-Sassandra (Ivory Coast) 5N-6W (UNESCO, p58)  R1
-          ! Data from HLA: 5N -  6  W   1953/86 = 449m3/s
-          ! --> 449.333344m3/s
-          criver(78) = "Sassandra           "
-          inb(78) = 1
-          iirnf(1:1,78) = 71
-          ijrnf(1:1,78) = 41
-          zrnfm(1:jpmois,78) = (/   &
-              144.,     116.,     119.,     143.,     156.,     334.,   &
-              419.,     692.,    1277.,    1190.,     553.,     249./)
-          ! 79-Jubba (Somali Coast)
-          criver(79) = "Jubba               "
-          inb(79) = 1
-          iirnf(1:1,79) = 84
-          ijrnf(1:1,79) = 39
-          zrnfm(1:jpmois,79) = 0.
-          ! 80-Senaga
-          criver(80) = "Senegal             "
-          inb(79) = 1
-          iirnf(1:1,80) = 66
-          ijrnf(1:1,80) = 47
-          zrnfm(1:jpmois,79) = 0.
-      
-          ! Total run-offs (VI)=63977.m3/s
-          ! Total run-offs (I+II+III+IV+V+VI)=560902.688m3/s
-          
-          ! VII OCEANIA
-          ! ------------
-          ! 91-Murray (Australia) 35S34 138E54              R5
-          ! Old=22 km3/year=697m3/s
-          ! --> 696.999939m3/s
-          criver(91) = "Murray              "
-          inb(91) = 1
-          iirnf(1:1,91) = 16
-          ijrnf(1:1,91) = 24
-          zrnfm(1:jpmois,91) = 697.
-          ! 92-Haast (N. Zeland) 43S50 169E02               R5
-          ! Old=6 km3/year=190m3/s
-          ! --> 189.999985m3/s
-          criver(92) = "Haast               "
-          inb(92) = 1
-          iirnf(1:1,92) = 24
-          ijrnf(1:1,92) = 22
-          zrnfm(1:jpmois,92) = 190.
-          ! 93-Fly (N. Guinea) 8S30 143E05                  R5
-          ! Old=77 km3/year=2440m3/s
-          ! --> 2440.m3/s
-          criver(93) = "Fly                 "
-          inb(93) = 1
-          iirnf(1:1,93) = 19
-          ijrnf(1:1,93) = 31
-          zrnfm(1:jpmois,93) = 2440.
-          ! 94-Purari (N. Guinea) 7S25 145E05               R4
-          ! Old=77 km3/year=2440m3/s; UNESCO(77+83,p595)=2727m3/s
-          ! --> 2727.33325m3/s
-          criver(94) = "Purari              "
-          inb(94) = 1
-          iirnf(1:1,94) = 19
-          ijrnf(1:1,94) = 33
-          zrnfm(1:jpmois,94) = (/   &
-              1972.,    3067.,    2332.,    2981.,    2693.,    2876.,   &
-              2867.,    3477.,    2195.,    2967.,    2674.,    2627./)
-          ! 95-Choshui ch''i (taiwan) 24N03 120E23          R2
-          ! Old=6 km3/year=190m3/s; UNESCO(65-68, p374)=120m3/s; ratio=1.58 
-          ! --> 0.m3/s
-          criver(95) = "Choshui             "
-          inb(95) = 1
-          iirnf(1:1,95) = 12
-          ijrnf(1:1,95) = 48
-          zrnfm(1:jpmois,95) = (/   &
-              65.,      63.,      83.,     100.,     177.,     535.,   &
-              332.,     313.,     236.,     193.,     101.,      84./)
-          ! 96-Kaoping (Taiwan)                             R5  
-          ! Old=9 km3/year=285m3/s
-          ! --> 0.m3/s
-          criver(96) = "Kaoping             "
-          inb(96) = 0
-          iirnf(1:1,96) = 0
-          ijrnf(1:1,96) = 0
-          zrnfm(1:jpmois,96) = 285.
-          ! 97-Tsengwen (Taiwan)                            R5
-          ! Old=2 km3/year=63m3/s
-          ! --> 0.m3/s
-          criver(97) = "Tsengwen            "
-          inb(97) = 0
-          iirnf(1:1,97) = 0
-          ijrnf(1:1,97) = 0
-          zrnfm(1:jpmois,97) = 63.
-          ! 98-Hualien (Taiwan)                             R5
-          ! Old=4 km3/year=127m3/s
-          ! --> 0.m3/s
-          criver(98) = "Hualien             "
-          inb(98) = 0
-          iirnf(1:1,98) = 0
-          ijrnf(1:1,98) = 0
-          zrnfm(1:jpmois,98) = 127.
-          ! 99-Peinan (Taiwan) 22N46 121E10                 R5
-          ! Old=4 km3/year=127m3/s
-          ! --> 127.000008m3/s
-          criver(99) = "Peinan              "
-          inb(99) = 1
-          iirnf(1:1,99) = 12
-          ijrnf(1:1,99) = 48
-          zrnfm(1:jpmois,99) = 127.
-          ! 100-Hsiukuluan (Taiwan)                         R5
-          ! Old=4 km3/year=127m3/s
-          ! --> 0.m3/s
-          criver(100) = "Hsiukuluan          "
-          inb(100) = 0
-          iirnf(1:1,100) = 0
-          ijrnf(1:1,100) = 0
-          zrnfm(1:jpmois,100) = 127.
-
-          !  Total run-offs(VII)=6181.33301m3/s
-          
-
-          ! VIII ANTARCTICA
-          ! ---------------
-          ! 103 -Antarctic ice discharge (2 613 Gt/year) (jacobs et al 1992)
-          ! Data from Jacobs et al 92: =2613 km3/year=82857.686454 m3/s
-          ! --> 82858 m3/s
-          criver(103) = "Antarctica         "
-          inb(103) = 108
-          iirnf(1:108,103) = (/   &
-               2, 3, 4, 5, 6, 7, 8, 9,10,11,12,13,14,15,16,17,18,19,20,   &
-              21,22,23,24,25,25,24,24,24,25,26,27,28,29,30,31,32,33,34,   &
-              35,36,37,38,39,40,41,42,43,44,45,46,47,48,49,50,51,51,53,   &
-              53,53,53,54,54,54,56,57,57,57,57,57,57,57,57,57,57,58,59,   &
-              60,61,62,63,64,65,66,67,68,69,70,71,72,73,74,75,76,77,78,   &
-              79,80,81,82,83,84,85,86,87,88,89,90,91 /)
-          ijrnf(1:108,103) = (/   &
-              12,12,12,12,12,12,12,12,12,12,12,12,12,12,12,12,12,11,11,   &
-              10,10, 9, 9, 8, 7, 6, 5, 3, 2, 2, 2, 2, 2, 2, 2, 2, 3, 3,   &
-               3, 4, 5, 6, 6, 6, 7, 7, 6, 6, 6, 6, 7, 7, 7, 7, 7, 7, 7,   &
-               8, 9,10,10,11,12,13,12,11,10, 9, 8, 7, 6, 5, 4, 3, 2, 2,   &
-               2, 2, 2, 2, 3, 4, 4, 5, 7, 8, 9, 9, 9, 9, 9, 9,10,10,10,   &
-              10,10,10,10,11,12,12,12,12,11,11,10,10 /)
-          zrnfm(1:jpmois,100) = 82858.0
-                   
-
-          !  Total run-offs(I+II+III+IV+V+VI+VII)= 567084. m3/s
-          !                                      = 1.588668056E-6 kg/m2/s
-          !                                      = 17883.5586 km3/Year
-          !                                      = 0.567083955 Sverdrup
-      ENDIF
-
-      !!----------------------------------------------------------------------
-      !!  OPA 8.5, LODYC-IPSL (2002)
-      !!----------------------------------------------------------------------
-
-      IF( kt == nit000 ) THEN
-         IF(lwp) WRITE(numout,*)
-         IF(lwp) WRITE(numout,*) 'flxrnf : runoff + orca4'
-         IF(lwp) WRITE(numout,*) '~~~~~~'
-         IF(lwp) WRITE(numout,*) '   nrunoff = ', nrunoff
-         IF(lwp) WRITE(numout,*)
-      ENDIF
-
-      ! 0. initialization and  indexes  modification for zoom version
-      ! --------------------------------------------------------------
-      
-      
-      ! conversion factor m3/s --> kg/m2/s *e1t*e2t
-
-      zconv = rauw
-
-      IF( kt == nit000 ) THEN 
-         ! Set runoffs and upstream coeff to zero
-         runoff (:,:) = 0.e0                                                                                  
-         upsrnfh(:,:) = 0.e0
-         upsrnfz(:)   = 0.e0
-         upsadv (:,:) = 0.e0
-         zcoefr(:,:,:) = 0.e0
-
-         IF( nrunoff == 0 ) THEN
-            IF(lwp) WRITE(numout,*)
-            IF(lwp) WRITE(numout,*) ' NO RUNOFF in this simulation '
-            IF(lwp) WRITE(numout,*)
-         ENDIF
-      ENDIF
-
-      IF( kt == nit000 ) THEN
-         IF( nrunoff >= 1 ) THEN 
-            zbilan = 0.0 
-!OCL SCALAR
-            DO jr = 1, jpriv-1
-               inbp = inb(jr)
-               IF( inbp /= 0 ) THEN
-                  ii = mj0(iirnf(1,jr))
-                  ij = mi0(ijrnf(1,jr))
-                  ii = MAX( 2, MIN( ii, jpi-1 ) )
-                  ij = MAX( 2, MIN( ij, jpj-1 ) )
-!OCL SCALAR
-                  IF( tmask(ii,ij,1) <= 0.5 ) THEN 
-                     IF(lwp) WRITE(numout,*) 'Erreur riviere masquee   ', ii, ij, jr, criver(jr)
-                  ENDIF
-                  IF( tmask(ii-1,ij,1)*tmask(ii+1,ij,1)   &
-                     *tmask(ii,ij-1,1)*tmask(ii,ij+1,1) /= 0. )THEN 
-                     IF(lwp) WRITE(numout,*) 'Erreur riviere non cotiere ',ii, ij, jr, criver(jr)
-                  ENDIF
-                  DO nrnf1 = 1, jpmois
-                     zbilan = zbilan + zrnfm(nrnf1,jr) * tmask(ii,ij,1)
-                  END DO
-               ENDIF
-            END DO
-            IF(lwp) WRITE(numout,*) 'zrnfm (runoff) ', zero, zconv, rauw, SUM(zrnfm)
-            IF(lwp) THEN 
-               WRITE(numout,'("Total runoff : ",1E30.15)' ) zbilan
-               WRITE(numout,*) 'Run-off : POSITION and value '
-               WRITE(numout,*) 'ji jj glamt gphit annual mask 12-month '
-!OCL SCALAR
-               DO jr = 1, jpriv-1
-                  inbp = inb(jr)
-                  IF(inbp /= 0) THEN
-!OCL SCALAR
-                     DO jl = 1, inbp
-                        ii = mj0(iirnf(jl,jr))
-                        ij = mi0(ijrnf(jl,jr))
-                        WRITE(numout,'(2I3,1X,1A20,2X,2I4,2F7.1,1I3,1PE10.1,1I2,12E9.1)')   &
-                           jr, jl, criver(jr),   &
-                           ii, ij, glamt(ii,ij), gphit(ii,ij), inbp,   &
-                           SUM(zrnfm(1:jpmois,jr))/float(inbp),   &
-                           NINT(tmask(ii,ij,1)),   &
-                           zrnfm(1:jpmois,jr)/float(inbp)
-                     END DO
-                  ENDIF
-               END DO
-            ENDIF
-            IF(lwp) WRITE(numout,*) ' zrnfm (runoff) ', zero, zconv, rauw, SUM(zrnfm)
-         ENDIF
-      ENDIF
-
-      ! 1. initialization of the runoff array
-      ! -------------------------------------
-
-      IF( nrunoff == 1 .OR. nrunoff == 2) THEN
-
-         iman  = jpmois
-         
-         i15 = nday/16
-         
-         imois = nmonth + i15 - 1
-         IF( imois == 0) imois = iman
-         IF( nn_leap_year == 1 .AND. MOD( nyear, 4 ) == 0 ) THEN
-            ! Number of Day Between Data
-            idbd = nbiss(imois)
-         ELSEIF( nn_leap_year > 1 ) THEN
-            idbd = nn_leap_year
-         ELSE
-            idbd = nobis(imois)
-         ENDIF
-         ! Number of days between imois, 15 and the end of month
-         idmeom = idbd - 15
-         
-         ! 1.1 Calendar computation
-         
-         ! nrnf1 number of the first array record used in the simulation
-         ! nrnf2 number of the last  array record
-         
-         nrnf1 = imois
-         nrnf2 = nrnf1 + 1
-         nrnf1 = MOD( nrnf1, iman )
-         IF( nrnf1 == 0 ) nrnf1 = iman
-         nrnf2 = MOD( nrnf2, iman )
-         IF( nrnf2 == 0 ) nrnf2 = iman
-         
-         IF( kt == nit000 ) THEN
-            IF(lwp) WRITE(numout,*)
-            IF(lwp) WRITE(numout,*) 'flxrnf routine'
-            IF(lwp) WRITE(numout,*) 'first array record used nrnf1 ',nrnf1
-            IF(lwp) WRITE(numout,*) 'last  array record used nrnf2 ',nrnf2
-            IF(lwp) WRITE(numout,*) 'total zrnfm ', SUM(zrnfm)
-            IF(lwp) WRITE(numout,*) 'total surface ', SUM(e1t*e2t*tmask(:,:,1))
-         ENDIF
-
-         zxy = FLOAT( nday + idmeom - idbd * i15 ) / idbd
-
-!OCL SCALAR
-         DO jr = 1, jpriv-1
-            inbp = inb(jr)
-!OCL SCALAR
-            IF( inbp /= 0 ) THEN
-               DO jl = 1, inbp
-                  ii = iirnf(jl,jr) - nimpp + 1 - jpizoom + 1
-                  ij = ijrnf(jl,jr) - njmpp + 1 - jpjzoom + 1
-                  iipt0 = MAX( 1, MIN( ii, jpi   ) )
-                  iipt1 = MAX( 0, MIN( ii, jpi-1 ) )
-                  ijpt0 = MAX( 1, MIN( ij, jpj   ) )
-                  ijpt1 = MAX( 0, MIN( ij, jpj-1 ) )
-                  DO jj = ijpt0, ijpt1
-                     DO ji = iipt0, iipt1
-                        runoff(ji,jj) = runoff(ji,jj)   & 
-                                      + ( (1.-zxy) * zrnfm(nrnf1,jr)   &
-                                      +       zxy  * zrnfm(nrnf2,jr)  ) / FLOAT(inbp)      
-                     END DO
-                  END DO
-               END DO
-            ENDIF
-         END DO
-
-
-         ! 2. conversion in mm/day
-         ! -----------------------
-         !     C a u t i o n : upward water flux, runoff is negative
-
-         runoff(:,:) = -zconv * runoff(:,:) / ( e1t(:,:) * e2t(:,:) )
-         
-
-         ! 3. advectve coefficient
-         ! -----------------------
-
-         !  Creates the array coef that contains the coefficient to affect to
-         !  the upstream scheme. advection scheme will be:
-         !  coefr * upstream + (1- coefr) centered
-         !  coefr must be between 0 and 1.
-         !  3 values are used: maximum = zrup1, between zrup2, small zrup3.
-         ! 
-         IF( nrunoff == 2 .AND. kt == nit000 ) THEN
-
-            zrup1 = 0.5
-            zrup2 = 0.25
-            zrup3 = 0.125
-            DO jr = 1, jpriv-1
-               inbp = inb(jr)
-               IF( inbp /= 0) THEN
-                  DO jl = 1, inbp
-                     ! --> zrup1 (originally, = 1 = all upstream).
-                     ! ji,jj
-                     ii = iirnf(jl,jr) - nimpp + 1 - jpizoom - 1
-                     ij = ijrnf(jl,jr) - njmpp + 1 - jpjzoom - 1
-                     iipt0 = MAX( 1, MIN( ii,jpi  ) )
-                     iipt1 = MAX( 0, MIN( ii,jpi-1) )
-                     ijpt0 = MAX( 1, MIN( ij,jpj  ) )
-                     ijpt1 = MAX( 0, MIN( ij,jpj-1) )
-                     DO jj = ijpt0, ijpt1
-                        DO ji = iipt0, iipt1
-                           DO jn = 1, 3
-                              zcoefr(ji,jj,jn) = zrup1
-                           END DO
-                           zcoefr(ji,jj,4) = zrup2
-                           zcoefr(ji,jj,5) = zrup3
-                        END DO
-                     END DO
-                     ! --> zrup2 (originally = 0.5)
-                     ! ji,jj+1
-                     ijpt0 = MAX( 1, MIN( ij+1, jpj  ) )
-                     ijpt1 = MAX( 0, MIN( ij+1, jpj-1) )
-                     DO jj = ijpt0, ijpt1
-                        DO ji = iipt0, iipt1
-                           DO jn = 1, 3
-                              zcoefr(ji,jj,jn) = MAX( zrup2, zcoefr(ji,jj,jn) )
-                           END DO
-                           zcoefr(ji,jj,4) = MAX( zrup3, zcoefr(ji,jj,4) )
-                        END DO
-                     END DO
-                     ! ji,jj-1
-                     ijpt0 = MAX( 1, MIN( ij-1, jpj  ) )
-                     ijpt1 = MAX( 0, MIN( ij-1, jpj-1) )
-                     DO jj = ijpt0, ijpt1
-                        DO ji = iipt0, iipt1
-                           DO jn = 1, 3
-                              zcoefr(ji,jj,jn) = MAX( zrup2, zcoefr(ji,jj,jn) )
-                           END DO
-                           zcoefr(ji,jj,4) = MAX( zrup3, zcoefr(ji,jj,4) )
-                        END DO
-                     END DO
-                     ! ji+1,jj-1
-                     iipt0 = MAX( 1, MIN( ii+1, jpi  ) )
-                     iipt1 = MAX( 0, MIN( ii+1, jpi-1) )
-                     DO jj = ijpt0, ijpt1
-                        DO ji = iipt0, iipt1
-                           DO jn = 1, 3
-                              zcoefr(ji,jj,jn) = MAX( zrup2, zcoefr(ji,jj,jn) )
-                           END DO
-                           zcoefr(ji,jj,4) = MAX( zrup3, zcoefr(ji,jj,4) )
-                        END DO
-                     END DO
-                     ! ji+1,jj
-                     ijpt0 = MAX( 1, MIN( ij, jpj  ) )
-                     ijpt1 = MAX( 0, MIN( ij, jpj-1) )
-                     DO jj = ijpt0, ijpt1
-                        DO ji = iipt0, iipt1
-                           DO jn = 1, 3
-                              zcoefr(ji,jj,jn) = MAX( zrup2, zcoefr(ji,jj,jn) )
-                           END DO
-                           zcoefr(ji,jj,4) = MAX( zrup3, zcoefr(ji,jj,4) )
-                        END DO
-                     END DO
-                     ! ji+1,jj+1
-                     ijpt0 = MAX( 1, MIN( ij+1,jpj  ) )
-                     ijpt1 = MAX( 0, MIN( ij+1,jpj-1) )
-                     DO jj = ijpt0, ijpt1
-                        DO ji = iipt0, iipt1
-                           DO jn = 1, 3
-                              zcoefr(ji,jj,jn) = MAX( zrup2, zcoefr(ji,jj,jn) )
-                           END DO
-                           zcoefr(ji,jj,4) = MAX( zrup3, zcoefr(ji,jj,4) )
-                        END DO
-                     END DO
-                     ! ji-1,jj+1
-                     iipt0 = MAX( 1, MIN( ii-1,jpi  ) )
-                     iipt1 = MAX( 0, MIN( ii-1,jpi-1) )
-                     DO jj = ijpt0, ijpt1
-                        DO ji = iipt0, iipt1
-                           DO jn = 1, 3
-                              zcoefr(ji,jj,jn) = MAX( zrup2, zcoefr(ji,jj,jn) )
-                           END DO
-                           zcoefr(ji,jj,4) = MAX( zrup3, zcoefr(ji,jj,4) )
-                        END DO
-                     END DO
-                     ! ji-1,jj
-                     ijpt0 = MAX( 1, MIN( ij,jpj  ) )
-                     ijpt1 = MAX( 0, MIN( ij,jpj-1) )
-                     DO jj = ijpt0, ijpt1
-                        DO ji = iipt0, iipt1
-                           DO jn = 1, 3
-                              zcoefr(ji,jj,jn) = MAX( zrup2, zcoefr(ji,jj,jn) )
-                           END DO
-                           zcoefr(ji,jj,4) = MAX( zrup3, zcoefr(ji,jj,4) )
-                        END DO
-                     END DO
-                     ! ji-1,jj-1
-                     ijpt0 = MAX( 1, MIN( ij-1,jpj  ) )
-                     ijpt1 = MAX( 0, MIN( ij-1,jpj-1) )
-                     DO jj = ijpt0, ijpt1
-                        DO ji = iipt0, iipt1
-                           DO jn = 1, 3
-                              zcoefr(ji,jj,jn) = MAX( zrup2, zcoefr(ji,jj,jn) )
-                           END DO
-                           zcoefr(ji,jj,4) = MAX( zrup3, zcoefr(ji,jj,4) )
-                        END DO
-                     END DO
-                     ! --> zrup3 (originally = 0.25)
-                     ! ji-1,jj+2
-                     ijpt0 = MAX( 1, MIN( ij+2,jpj  ) )
-                     ijpt1 = MAX( 0, MIN( ij+2,jpj-1) )
-                     DO jj = ijpt0, ijpt1
-                        DO ji = iipt0, iipt1
-                           DO jn = 1, 3
-                              zcoefr(ji,jj,jn) = MAX( zrup3, zcoefr(ji,jj,jn) )
-                           END DO
-                        END DO
-                     END DO
-                     ! ji-1,jj-2
-                     ijpt0 = MAX( 1, MIN( ij-2,jpj  ) )
-                     ijpt1 = MAX( 0, MIN( ij-2,jpj-1) )
-                     DO jj = ijpt0, ijpt1
-                        DO ji = iipt0, iipt1
-                           DO jn = 1, 3 
-                              zcoefr(ji,jj,jn) = MAX( zrup3, zcoefr(ji,jj,jn) )
-                           END DO
-                        END DO
-                     END DO
-                     ! ji,jj-2
-                     iipt0 = MAX( 1, MIN( ii,jpi  ) )
-                     iipt1 = MAX( 0, MIN( ii,jpi-1) )
-                     DO jj = ijpt0, ijpt1
-                        DO ji = iipt0, iipt1
-                           DO jn = 1, 3 
-                              zcoefr(ji,jj,jn) = MAX( zrup3, zcoefr(ji,jj,jn) )
-                           END DO
-                        END DO
-                     END DO
-                     ! ji,jj+2
-                     ijpt0 = MAX( 1, MIN( ij+2,jpj  ) )
-                     ijpt1 = MAX( 0, MIN( ij+2,jpj-1) )
-                     DO jj = ijpt0, ijpt1
-                        DO ji = iipt0, iipt1
-                           DO jn = 1, 3
-                              zcoefr(ji,jj,jn) = MAX( zrup3, zcoefr(ji,jj,jn) )
-                           END DO
-                        END DO
-                     END DO
-                     ! ji+1,jj-2
-                     iipt0 = MAX( 1, MIN( ii+1,jpi  ) )
-                     iipt1 = MAX( 0, MIN( ii+1,jpi-1) )
-                     DO jj = ijpt0, ijpt1
-                        DO ji = iipt0, iipt1
-                           DO jn = 1, 3 
-                              zcoefr(ji,jj,jn) = MAX( zrup3, zcoefr(ji,jj,jn) )
-                           END DO
-                        END DO
-                     END DO
-                     ! ji+1,jj+2
-                     ijpt0 = MAX( 1, MIN( ij+2,jpj  ) )
-                     ijpt1 = MAX( 0, MIN( ij+2,jpj-1) )
-                     DO jj = ijpt0, ijpt1
-                        DO ji = iipt0, iipt1
-                           DO jn = 1, 3 
-                              zcoefr(ji,jj,jn) = MAX( zrup3, zcoefr(ji,jj,jn) )
-                           END DO
-                        END DO
-                     END DO
-                     ! ji+2,jj-2
-                     iipt0 = MAX( 1, MIN( ii+2,jpi  ) )
-                     iipt1 = MAX( 0, MIN( ii+2,jpi-1) )
-                     DO jj = ijpt0, ijpt1
-                        DO ji = iipt0, iipt1
-                           DO jn = 1, 3 
-                              zcoefr(ji,jj,jn) = MAX( zrup3, zcoefr(ji,jj,jn) )
-                           END DO
-                        END DO
-                     END DO
-                     ! ji+2,jj+2
-                     ijpt0 = MAX( 1, MIN( ij+2,jpj  ) )
-                     ijpt1 = MAX( 0, MIN( ij+2,jpj-1) )
-                     DO jj = ijpt0, ijpt1
-                        DO ji = iipt0, iipt1
-                           DO jn = 1, 3  
-                              zcoefr(ji,jj,jn) = MAX( zrup3, zcoefr(ji,jj,jn) )
-                           END DO
-                        END DO
-                     END DO
-                     ! ji+2,jj-1
-                     ijpt0 = MAX( 1, MIN( ij-1,jpj  ) )
-                     ijpt1 = MAX( 0, MIN( ij-1,jpj-1) )
-                     DO jj = ijpt0, ijpt1
-                        DO ji = iipt0, iipt1
-                           DO jn = 1, 3
-                              zcoefr(ji,jj,jn) = MAX( zrup3, zcoefr(ji,jj,jn) )
-                           END DO
-                        END DO
-                     END DO
-                     ! ji+2,jj+1
-                     ijpt0 = MAX( 1, MIN( ij+1,jpj  ) )
-                     ijpt1 = MAX( 0, MIN( ij+1,jpj-1) )
-                     DO jj = ijpt0, ijpt1
-                        DO ji = iipt0, iipt1
-                           DO jn = 1, 3 
-                              zcoefr(ji,jj,jn) = MAX( zrup3, zcoefr(ji,jj,jn) )
-                           END DO
-                        END DO
-                     END DO
-                     ! ji+2,jj
-                     ijpt0 = MAX( 1, MIN( ij,jpj  ) )
-                     ijpt1 = MAX( 0, MIN( ij,jpj-1) )
-                     DO jj = ijpt0, ijpt1
-                        DO ji = iipt0, iipt1
-                           DO jn = 1, 3
-                              zcoefr(ji,jj,jn) = MAX( zrup3, zcoefr(ji,jj,jn) )
-                           END DO
-                        END DO
-                     END DO
-                     ! ji-2,jj
-                     iipt0 = MAX( 1, MIN( ii-2,jpi  ) )
-                     iipt1 = MAX( 0, MIN( ii-2,jpi-1) )
-                     DO jj = ijpt0, ijpt1
-                        DO ji = iipt0, iipt1
-                           DO jn = 1, 3        
-                              zcoefr(ji,jj,jn) = MAX( zrup3, zcoefr(ji,jj,jn) )
-                           END DO
-                        END DO
-                     END DO
-                     ! ji-2,jj-2
-                     ijpt0 = MAX( 1, MIN( ij-2,jpj  ) )
-                     ijpt1 = MAX( 0, MIN( ij-2,jpj-1) )
-                     DO jj = ijpt0, ijpt1
-                        DO ji = iipt0, iipt1
-                           DO jn = 1, 3 
-                              zcoefr(ji,jj,jn) = MAX( zrup3, zcoefr(ji,jj,jn) )
-                           END DO
-                        END DO
-                     END DO
-                     ! ji-2,jj+2
-                     ijpt0 = MAX( 1, MIN( ij+2,jpj  ) )
-                     ijpt1 = MAX( 0, MIN( ij+2,jpj-1) )
-                     DO jj = ijpt0, ijpt1
-                        DO ji = iipt0, iipt1
-                           DO jn = 1, 3
-                              zcoefr(ji,jj,jn) = MAX( zrup3, zcoefr(ji,jj,jn) )
-                           END DO
-                        END DO
-                     END DO
-                     ! ji-2,jj-1
-                     ijpt0 = MAX( 1, MIN( ij-1,jpj  ) )
-                     ijpt1 = MAX( 0, MIN( ij-1,jpj-1) )
-                     DO jj = ijpt0, ijpt1
-                        DO ji = iipt0, iipt1
-                           DO jn = 1, 3 
-                              zcoefr(ji,jj,jn) = MAX( zrup3, zcoefr(ji,jj,jn) )
-                           END DO
-                        END DO
-                     END DO
-                     ! ji-2,jj+1
-                     ijpt0 = MAX( 1, MIN( ij+1,jpj  ) )
-                     ijpt1 = MAX( 0, MIN( ij+1,jpj-1) )
-                     DO jj = ijpt0, ijpt1
-                        DO ji = iipt0, iipt1
-                           DO jn = 1, 3 
-                              zcoefr(ji,jj,jn) = MAX( zrup3, zcoefr(ji,jj,jn) )
-                           END DO
-                        END DO
-                     END DO
-                  END DO
-               ENDIF
-            END DO
-            upsrnfh(:,:) = zcoefr(:,:,1)
-            upsrnfz(:)   = 0.e0
-            upsrnfz(1)   = 1.0
-            upsrnfz(2)   = 1.0
-            upsrnfz(3)   = zrup1
-            upsrnfz(4)   = zrup2
-            upsrnfz(5)   = zrup3
-         ENDIF
-      ENDIF
-      
-      ! Upstream and centered scheme in the vicinity of some straits
-      
-      IF( kt == nit000 ) THEN
-
-         ! ORCA_R4 configuration :
-         ! Gibraltar Strait
-         DO jj = mj0(52), mj1(53)
-            DO ji = mi0(70), mi1(71)
-               upsadv(ji,jj) = 0.5
-            END DO
-         END DO
-
-      ENDIF
