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 15150 for NEMO/trunk/src/OCE/DYN/divhor.F90 – NEMO

Ignore:
Timestamp:
2021-07-27T12:38:24+02:00 (3 years ago)
Author:
smasson
Message:

trunk: remove some of the changes done in [15145]

File:
1 edited

Legend:

Unmodified
Added
Removed

  • NEMO/trunk/src/OCE/DYN/divhor.F90

    r15145 r15150  
    7979      ENDIF 
    8080      ! 
    81       DO_3D_OVR( nn_hls-1, nn_hls-1, nn_hls-1, nn_hls-1, 1, jpkm1 )                                      !==  Horizontal divergence  ==! 
     81      DO_3D_OVR( nn_hls-1, nn_hls, nn_hls-1, nn_hls, 1, jpkm1 )                                          !==  Horizontal divergence  ==! 
    8282         ! round brackets added to fix the order of floating point operations 
    8383         ! needed to ensure halo 1 - halo 2 compatibility 
    84          hdiv(ji,jj,jk) = (  ( e2u(ji  ,jj) * e3u(ji  ,jj,jk,Kmm) * uu(ji  ,jj,jk,Kmm)     & ! Warning: with qco, e3u uses r3u that 
    85             &                - e2u(ji-1,jj) * e3u(ji-1,jj,jk,Kmm) * uu(ji-1,jj,jk,Kmm)     & !          is not defined in jpi 
     84         hdiv(ji,jj,jk) = (  ( e2u(ji  ,jj) * e3u(ji  ,jj,jk,Kmm) * uu(ji  ,jj,jk,Kmm)     & 
     85            &                - e2u(ji-1,jj) * e3u(ji-1,jj,jk,Kmm) * uu(ji-1,jj,jk,Kmm)     & 
    8686            &                )                                                             & ! bracket for halo 1 - halo 2 compatibility 
    87             &              + ( e1v(ji,jj  ) * e3v(ji,jj  ,jk,Kmm) * vv(ji,jj  ,jk,Kmm)     & ! Warning: with qco, e3v uses r3v that 
    88             &                - e1v(ji,jj-1) * e3v(ji,jj-1,jk,Kmm) * vv(ji,jj-1,jk,Kmm)     & !          is not defined in jpj 
     87            &              + ( e1v(ji,jj  ) * e3v(ji,jj  ,jk,Kmm) * vv(ji,jj  ,jk,Kmm)     & 
     88            &                - e1v(ji,jj-1) * e3v(ji,jj-1,jk,Kmm) * vv(ji,jj-1,jk,Kmm)     & 
    8989            &                )                                                             & ! bracket for halo 1 - halo 2 compatibility 
    9090            &             )  * r1_e1e2t(ji,jj) / e3t(ji,jj,jk,Kmm) 
Note: See TracChangeset for help on using the changeset viewer.