- Timestamp:
- 2011-12-11T16:00:26+01:00 (12 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2011/DEV_r2739_STFC_dCSE/NEMOGCM/NEMO/OPA_SRC/OBC/obcdyn.F90
r2528 r3211 45 45 REAL(wp) :: rtaue , rtauw , rtaun , rtaus , & 46 46 rtauein, rtauwin, rtaunin, rtausin 47 48 !! * Control permutation of array indices 49 # include "oce_ftrans.h90" 50 # include "dom_oce_ftrans.h90" 51 # include "obc_oce_ftrans.h90" 47 52 48 53 !!--------------------------------------------------------------------------------- … … 147 152 ! 1.1 U zonal velocity 148 153 ! -------------------- 154 #if defined key_z_first 155 DO jj = 1, jpj 156 DO ji = nie0, nie1 157 DO jk = 1, jpkm1 158 #else 149 159 DO ji = nie0, nie1 150 160 DO jk = 1, jpkm1 151 161 DO jj = 1, jpj 162 #endif 152 163 ua(ji,jj,jk) = ua(ji,jj,jk) * (1.-uemsk(jj,jk)) + & 153 164 uemsk(jj,jk)*ufoe(jj,jk) … … 158 169 ! 1.2 V meridional velocity 159 170 ! ------------------------- 171 #if defined key_z_first 172 DO jj = 1, jpj 173 DO ji = nie0+1, nie1+1 174 DO jk = 1, jpkm1 175 #else 160 176 DO ji = nie0+1, nie1+1 161 177 DO jk = 1, jpkm1 162 178 DO jj = 1, jpj 179 #endif 163 180 va(ji,jj,jk) = va(ji,jj,jk) * (1.-vemsk(jj,jk)) + & 164 181 vfoe(jj,jk)*vemsk(jj,jk) … … 191 208 ! ... radiative conditions on the total part + relaxation toward climatology 192 209 ! ... (jpjedp1, jpjefm1),jpieob 210 #if defined key_z_first 211 DO jj = 1, jpj 212 DO ji = nie0, nie1 213 DO jk = 1, jpkm1 214 #else 193 215 DO ji = nie0, nie1 194 216 DO jk = 1, jpkm1 195 217 DO jj = 1, jpj 218 #endif 196 219 z05cx = u_cxebnd(jj,jk) 197 220 z05cx = z05cx / e1t(ji,jj) … … 229 252 ! ... radiative condition 230 253 ! ... (jpjedp1, jpjefm1), jpieob+1 254 #if defined key_z_first 255 DO jj = 1, jpj 256 DO ji = nie0+1, nie1+1 257 DO jk = 1, jpkm1 258 #else 231 259 DO ji = nie0+1, nie1+1 232 260 DO jk = 1, jpkm1 233 261 DO jj = 1, jpj 262 #endif 234 263 z05cx = v_cxebnd(jj,jk) 235 264 z05cx = z05cx / e1f(ji-1,jj) … … 289 318 ! 1.1 U zonal velocity 290 319 ! --------------------- 320 #if defined key_z_first 321 DO jj = 1, jpj 322 DO ji = niw0, niw1 323 DO jk = 1, jpkm1 324 #else 291 325 DO ji = niw0, niw1 292 326 DO jk = 1, jpkm1 293 327 DO jj = 1, jpj 328 #endif 294 329 ua(ji,jj,jk) = ua(ji,jj,jk) * (1.-uwmsk(jj,jk)) + & 295 330 uwmsk(jj,jk)*ufow(jj,jk) … … 300 335 ! 1.2 V meridional velocity 301 336 ! ------------------------- 337 #if defined key_z_first 338 DO jj = 1, jpj 339 DO ji = niw0, niw1 340 DO jk = 1, jpkm1 341 #else 302 342 DO ji = niw0, niw1 303 343 DO jk = 1, jpkm1 304 344 DO jj = 1, jpj 345 #endif 305 346 va(ji,jj,jk) = va(ji,jj,jk) * (1.-vwmsk(jj,jk)) + & 306 347 vfow(jj,jk)*vwmsk(jj,jk) … … 333 374 ! ... radiative conditions on the total part + relaxation toward climatology 334 375 ! ... (jpjwdp1, jpjwfm1), jpiwob 376 #if defined key_z_first 377 DO jj = 1, jpj 378 DO ji = niw0, niw1 379 DO jk = 1, jpkm1 380 #else 335 381 DO ji = niw0, niw1 336 382 DO jk = 1, jpkm1 337 383 DO jj = 1, jpj 384 #endif 338 385 z05cx = u_cxwbnd(jj,jk) 339 386 z05cx = z05cx / e1t(ji+1,jj) … … 370 417 ! ... radiative condition plus Raymond-Kuo 371 418 ! ... (jpjwdp1, jpjwfm1),jpiwob 419 #if defined key_z_first 420 DO jj = 1, jpj 421 DO ji = niw0, niw1 422 DO jk = 1, jpkm1 423 #else 372 424 DO ji = niw0, niw1 373 425 DO jk = 1, jpkm1 374 426 DO jj = 1, jpj 427 #endif 375 428 z05cx = v_cxwbnd(jj,jk) 376 429 z05cx = z05cx / e1f(ji,jj) … … 429 482 ! -------------------- 430 483 DO jj = njn0+1, njn1+1 484 #if defined key_z_first 485 DO ji = 1, jpi 486 DO jk = 1, jpkm1 487 #else 431 488 DO jk = 1, jpkm1 432 489 DO ji = 1, jpi 490 #endif 433 491 ua(ji,jj,jk)= ua(ji,jj,jk) * (1.-unmsk(ji,jk)) + & 434 492 ufon(ji,jk)*unmsk(ji,jk) … … 440 498 ! ------------------------- 441 499 DO jj = njn0, njn1 500 #if defined key_z_first 501 DO ji = 1, jpi 502 DO jk = 1, jpkm1 503 #else 442 504 DO jk = 1, jpkm1 443 505 DO ji = 1, jpi 506 #endif 444 507 va(ji,jj,jk)= va(ji,jj,jk) * (1.-vnmsk(ji,jk)) + & 445 508 vfon(ji,jk)*vnmsk(ji,jk) … … 474 537 ! ... jpjnob+1,(jpindp1, jpinfm1) 475 538 DO jj = njn0+1, njn1+1 539 #if defined key_z_first 540 DO ji = 1, jpi 541 DO jk = 1, jpkm1 542 #else 476 543 DO jk = 1, jpkm1 477 544 DO ji = 1, jpi 545 #endif 478 546 z05cx= u_cynbnd(ji,jk) 479 547 z05cx = z05cx / e2f(ji, jj-1) … … 518 586 ! ... jpjnob,(jpindp1, jpinfm1) 519 587 DO jj = njn0, njn1 588 #if defined key_z_first 589 DO ji = 1, jpi 590 DO jk = 1, jpkm1 591 #else 520 592 DO jk = 1, jpkm1 521 593 DO ji = 1, jpi 594 #endif 522 595 ! ... 2* gradj(v) (T-point i=nibm, time mean) 523 596 z05cx = v_cynbnd(ji,jk) … … 580 653 ! -------------------- 581 654 DO jj = njs0, njs1 655 #if defined key_z_first 656 DO ji = 1, jpi 657 DO jk = 1, jpkm1 658 #else 582 659 DO jk = 1, jpkm1 583 660 DO ji = 1, jpi 661 #endif 584 662 ua(ji,jj,jk)= ua(ji,jj,jk) * (1.-usmsk(ji,jk)) + & 585 663 usmsk(ji,jk) * ufos(ji,jk) … … 591 669 ! ------------------------- 592 670 DO jj = njs0, njs1 671 #if defined key_z_first 672 DO ji = 1, jpi 673 DO jk = 1, jpkm1 674 #else 593 675 DO jk = 1, jpkm1 594 676 DO ji = 1, jpi 677 #endif 595 678 va(ji,jj,jk)= va(ji,jj,jk) * (1.-vsmsk(ji,jk)) + & 596 679 vsmsk(ji,jk) * vfos(ji,jk) … … 624 707 ! ... jpjsob,(jpisdp1, jpisfm1) 625 708 DO jj = njs0, njs1 709 #if defined key_z_first 710 DO ji = 1, jpi 711 DO jk = 1, jpkm1 712 #else 626 713 DO jk = 1, jpkm1 627 714 DO ji = 1, jpi 715 #endif 628 716 z05cx= u_cysbnd(ji,jk) 629 717 z05cx = z05cx / e2f(ji, jj) … … 665 753 ! ... jpjsob,(jpisdp1,jpisfm1) 666 754 DO jj = njs0, njs1 755 #if defined key_z_first 756 DO ji = 1, jpi 757 DO jk = 1, jpkm1 758 #else 667 759 DO jk = 1, jpkm1 668 760 DO ji = 1, jpi 761 #endif 669 762 z05cx = v_cysbnd(ji,jk) 670 763 z05cx = z05cx / e2t(ji,jj+1)
Note: See TracChangeset
for help on using the changeset viewer.