New URL for NEMO forge!   http://forge.nemo-ocean.eu

Since March 2022 along with NEMO 4.2 release, the code development moved to a self-hosted GitLab.
This present forge is now archived and remained online for history.
Changeset 4463 – NEMO

Changeset 4463


Ignore:
Timestamp:
2014-02-04T13:17:48+01:00 (10 years ago)
Author:
trackstand2
Message:

Use of mbkmax in bdy_dta_fra

File:
1 edited

Legend:

Unmodified
Added
Removed
  • branches/2011/DEV_r2739_STFC_dCSE/NEMOGCM/NEMO/OPA_SRC/BDY/bdydta.F90

    r3432 r4463  
    340340                  ii = nbi(ib,igrd) 
    341341                  ij = nbj(ib,igrd) 
    342                   DO ik = 1, jpkm1 
     342                  DO ik = 1, mbkmax(ii,ij)-1 ! jpkm1 
    343343                     tbdy(ib,ik) = tn(ii,ij,ik) 
    344344                     sbdy(ib,ik) = sn(ii,ij,ik) 
     
    352352                  ii = nbi(ib,igrd) 
    353353                  ij = nbj(ib,igrd) 
    354                   DO ik = 1, jpkm1 
     354                  DO ik = 1, mbkmax(ii,ij)-1 ! jpkm1 
    355355                     ubdy(ib,ik) = un(ii, ij, ik) 
    356356                  END DO 
     
    361361                  ii = nbi(ib,igrd) 
    362362                  ij = nbj(ib,igrd) 
    363                   DO ik = 1, jpkm1 
     363                  DO ik = 1, mbkmax(ii,ij)-1 ! jpkm1 
    364364                     vbdy(ib,ik) = vn(ii, ij, ik) 
    365365                  END DO 
     
    409409               ! 
    410410               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 
    412416                     tbdydta(ib,ik,2) =  zdta(nbmap(ib,igrd),1,ik) 
    413417                  END DO 
     
    427431               ! 
    428432               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 
    430438                     sbdydta(ib,ik,2) =  zdta(nbmap(ib,igrd),1,ik) 
    431439                  END DO 
     
    610618                             zdta(1:ipi,1:ipj,1:ipk), nbdy_a, lzfirst=.FALSE. ) 
    611619              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 
    613625                    tbdydta(ib,ik,2) = zdta(nbmap(ib,igrd),1,ik) 
    614626                 END DO 
     
    622634                             zdta(1:ipi,1:ipj,1:ipk), nbdy_a, lzfirst=.FALSE. ) 
    623635              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 
    625639                    sbdydta(ib,ik,2) = zdta(nbmap(ib,igrd),1,ik) 
    626640                 END DO 
     
    698712             igrd = 1                                   ! temperature & salinity 
    699713             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 
    701717                 tbdy(ib,ik) = zxy * tbdydta(ib,ik,2) + (1.-zxy) * tbdydta(ib,ik,1) 
    702718                 sbdy(ib,ik) = zxy * sbdydta(ib,ik,2) + (1.-zxy) * sbdydta(ib,ik,1) 
     
    708724             igrd = 2                                   ! u-velocity 
    709725             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 
    711729                 ubdy(ib,ik) = zxy * ubdydta(ib,ik,2) + (1.-zxy) * ubdydta(ib,ik,1)    
    712730               END DO 
     
    715733             igrd = 3                                   ! v-velocity 
    716734             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 
    718738                 vbdy(ib,ik) = zxy * vbdydta(ib,ik,2) + (1.-zxy) * vbdydta(ib,ik,1)    
    719739               END DO 
Note: See TracChangeset for help on using the changeset viewer.