Changeset 4463
- Timestamp:
- 2014-02-04T13:17:48+01:00 (10 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2011/DEV_r2739_STFC_dCSE/NEMOGCM/NEMO/OPA_SRC/BDY/bdydta.F90
r3432 r4463 340 340 ii = nbi(ib,igrd) 341 341 ij = nbj(ib,igrd) 342 DO ik = 1, jpkm1342 DO ik = 1, mbkmax(ii,ij)-1 ! jpkm1 343 343 tbdy(ib,ik) = tn(ii,ij,ik) 344 344 sbdy(ib,ik) = sn(ii,ij,ik) … … 352 352 ii = nbi(ib,igrd) 353 353 ij = nbj(ib,igrd) 354 DO ik = 1, jpkm1354 DO ik = 1, mbkmax(ii,ij)-1 ! jpkm1 355 355 ubdy(ib,ik) = un(ii, ij, ik) 356 356 END DO … … 361 361 ii = nbi(ib,igrd) 362 362 ij = nbj(ib,igrd) 363 DO ik = 1, jpkm1363 DO ik = 1, mbkmax(ii,ij)-1 ! jpkm1 364 364 vbdy(ib,ik) = vn(ii, ij, ik) 365 365 END DO … … 409 409 ! 410 410 DO ib = 1, nblen(igrd) 411 DO ik = 1, jpkm1 411 412 ii = nbi(ib,igrd) 413 ij = nbj(ib,igrd) 414 415 DO ik = 1, mbkmax(ii,ij)-1 ! jpkm1 412 416 tbdydta(ib,ik,2) = zdta(nbmap(ib,igrd),1,ik) 413 417 END DO … … 427 431 ! 428 432 DO ib = 1, nblen(igrd) 429 DO ik = 1, jpkm1 433 434 ii = nbi(ib,igrd) 435 ij = nbj(ib,igrd) 436 437 DO ik = 1, mbkmax(ii,ij)-1 ! jpkm1 430 438 sbdydta(ib,ik,2) = zdta(nbmap(ib,igrd),1,ik) 431 439 END DO … … 610 618 zdta(1:ipi,1:ipj,1:ipk), nbdy_a, lzfirst=.FALSE. ) 611 619 DO ib = 1, nblen(igrd) 612 DO ik = 1, jpkm1 620 621 ii = nbi(ib,igrd) 622 ij = nbj(ib,igrd) 623 624 DO ik = 1, mbkmax(ii,ij)-1 ! jpkm1 613 625 tbdydta(ib,ik,2) = zdta(nbmap(ib,igrd),1,ik) 614 626 END DO … … 622 634 zdta(1:ipi,1:ipj,1:ipk), nbdy_a, lzfirst=.FALSE. ) 623 635 DO ib = 1, nblen(igrd) 624 DO ik = 1, jpkm1 636 ii = nbi(ib,igrd) 637 ij = nbj(ib,igrd) 638 DO ik = 1, mbkmax(ii,ij)-1 ! jpkm1 625 639 sbdydta(ib,ik,2) = zdta(nbmap(ib,igrd),1,ik) 626 640 END DO … … 698 712 igrd = 1 ! temperature & salinity 699 713 DO ib = 1, nblen(igrd) 700 DO ik = 1, jpkm1 714 ii = nbi(ib,igrd) 715 ij = nbj(ib,igrd) 716 DO ik = 1, mbkmax(ii,ij)-1 ! jpkm1 701 717 tbdy(ib,ik) = zxy * tbdydta(ib,ik,2) + (1.-zxy) * tbdydta(ib,ik,1) 702 718 sbdy(ib,ik) = zxy * sbdydta(ib,ik,2) + (1.-zxy) * sbdydta(ib,ik,1) … … 708 724 igrd = 2 ! u-velocity 709 725 DO ib = 1, nblen(igrd) 710 DO ik = 1, jpkm1 726 ii = nbi(ib,igrd) 727 ij = nbj(ib,igrd) 728 DO ik = 1, mbkmax(ii,ij)-1 ! jpkm1 711 729 ubdy(ib,ik) = zxy * ubdydta(ib,ik,2) + (1.-zxy) * ubdydta(ib,ik,1) 712 730 END DO … … 715 733 igrd = 3 ! v-velocity 716 734 DO ib = 1, nblen(igrd) 717 DO ik = 1, jpkm1 735 ii = nbi(ib,igrd) 736 ij = nbj(ib,igrd) 737 DO ik = 1, mbkmax(ii,ij)-1 ! jpkm1 718 738 vbdy(ib,ik) = zxy * vbdydta(ib,ik,2) + (1.-zxy) * vbdydta(ib,ik,1) 719 739 END DO
Note: See TracChangeset
for help on using the changeset viewer.