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