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 5460 – NEMO

Changeset 5460


Ignore:
Timestamp:
2015-06-22T14:20:11+02:00 (9 years ago)
Author:
jchanut
Message:

Correct output for bottom velocities introduced in ticket #1474

File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/NEMOGCM/NEMO/OPA_SRC/DIA/diawri.F90

    r5407 r5460  
    125125      !! 
    126126      INTEGER                      ::   ji, jj, jk              ! dummy loop indices 
     127      INTEGER                      ::   jkbot                   ! 
    127128      REAL(wp)                     ::   zztmp, zztmpx, zztmpy   !  
    128129      !! 
     
    154155         DO jj = 1, jpj 
    155156            DO ji = 1, jpi 
    156                z2d(ji,jj) = tsn(ji,jj,MAX(mbathy(ji,jj),1),jp_tem) 
     157               jkbot = mbkt(ji,jj) 
     158               z2d(ji,jj) = tsn(ji,jj,jkbot,jp_tem) 
    157159            END DO 
    158160         END DO 
     
    165167         DO jj = 1, jpj 
    166168            DO ji = 1, jpi 
    167                z2d(ji,jj) = tsn(ji,jj,MAX(mbathy(ji,jj),1),jp_sal) 
     169               jkbot = mbkt(ji,jj) 
     170               z2d(ji,jj) = tsn(ji,jj,jkbot,jp_sal) 
    168171            END DO 
    169172         END DO 
     
    176179         DO jj = 1, jpj 
    177180            DO ji = 1, jpi 
    178                z2d(ji,jj) = un(ji,jj,MAX(mbathy(ji,jj),1)) 
     181               jkbot = mbku(ji,jj) 
     182               z2d(ji,jj) = un(ji,jj,jkbot) 
    179183            END DO 
    180184         END DO 
     
    187191         DO jj = 1, jpj 
    188192            DO ji = 1, jpi 
    189                z2d(ji,jj) = vn(ji,jj,MAX(mbathy(ji,jj),1)) 
     193               jkbot = mbkv(ji,jj) 
     194               z2d(ji,jj) = vn(ji,jj,jkbot) 
    190195            END DO 
    191196         END DO 
Note: See TracChangeset for help on using the changeset viewer.