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

Changeset 9117


Ignore:
Timestamp:
2017-12-18T16:17:39+01:00 (6 years ago)
Author:
acc
Message:

Branch 2017/dev_merge_2017. Fix uninitialised variables in sbcwave.F90 + 1 efficiency change

File:
1 edited

Legend:

Unmodified
Added
Removed
  • branches/2017/dev_merge_2017/NEMOGCM/NEMO/OPA_SRC/SBC/sbcwave.F90

    r9115 r9117  
    170170            END DO 
    171171         END DO 
     172         zsqrtpi = SQRT(rpi) 
     173         z_two_thirds = 2.0_wp / 3.0_wp 
    172174         DO jk = 1, jpkm1 
    173175            DO jj = 2, jpjm1 
    174176               DO ji = 2, jpim1 
    175                   zbot_u = 0.5_wp * ( gdepw_n(ji,jj,jk+1) + gdepw_n(ji+1,jj,jk+1) ) 
    176                   zbot_v = 0.5_wp * ( gdepw_n(ji,jj,jk+1) + gdepw_n(ji,jj+1,jk+1) ) 
    177                   zkb_u  = 2.0_wp * zk_u(ji,jj) * zbot_u                             ! 2k * bottom depth 
    178                   zkb_v  = 2.0_wp * zk_v(ji,jj) * zbot_v                             ! 2k * bottom depth 
     177                  zbot_u = ( gdepw_n(ji,jj,jk+1) + gdepw_n(ji+1,jj,jk+1) )  ! 2 * bottom depth 
     178                  zbot_v = ( gdepw_n(ji,jj,jk+1) + gdepw_n(ji,jj+1,jk+1) )  ! 2 * bottom depth 
     179                  zkb_u  = zk_u(ji,jj) * zbot_u                             ! 2 * k * bottom depth 
     180                  zkb_v  = zk_v(ji,jj) * zbot_v                             ! 2 * k * bottom depth 
    179181                  ! 
    180182                  zke3_u = MAX(1.e-8_wp, 2.0_wp * zk_u(ji,jj) * e3u_n(ji,jj,jk))     ! 2k * thickness 
Note: See TracChangeset for help on using the changeset viewer.