Changeset 9117
 Timestamp:
 20171218T16:17:39+01:00 (5 years ago)
 File:

 1 edited
Legend:
 Unmodified
 Added
 Removed

branches/2017/dev_merge_2017/NEMOGCM/NEMO/OPA_SRC/SBC/sbcwave.F90
r9115 r9117 170 170 END DO 171 171 END DO 172 zsqrtpi = SQRT(rpi) 173 z_two_thirds = 2.0_wp / 3.0_wp 172 174 DO jk = 1, jpkm1 173 175 DO jj = 2, jpjm1 174 176 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 depth178 zkb_v = 2.0_wp * zk_v(ji,jj) * zbot_v ! 2k * bottom depth177 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 179 181 ! 180 182 zke3_u = MAX(1.e8_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.