- Timestamp:
- 2019-11-22T15:29:17+01:00 (4 years ago)
- Location:
- NEMO/branches/2019/dev_r11943_MERGE_2019/src
- Files:
-
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/branches/2019/dev_r11943_MERGE_2019/src
- Property svn:mergeinfo deleted
-
NEMO/branches/2019/dev_r11943_MERGE_2019/src/OCE/SBC/sbcwave.F90
r11536 r11949 80 80 CONTAINS 81 81 82 SUBROUTINE sbc_stokes( )82 SUBROUTINE sbc_stokes( Kmm ) 83 83 !!--------------------------------------------------------------------- 84 84 !! *** ROUTINE sbc_stokes *** … … 92 92 !! ** action 93 93 !!--------------------------------------------------------------------- 94 INTEGER, INTENT(in) :: Kmm ! ocean time level index 94 95 INTEGER :: jj, ji, jk ! dummy loop argument 95 96 INTEGER :: ik ! local integer … … 152 153 DO jj = 2, jpjm1 153 154 DO ji = 2, jpim1 154 zdep_u = 0.5_wp * ( gdept _n(ji,jj,jk) + gdept_n(ji+1,jj,jk) )155 zdep_v = 0.5_wp * ( gdept _n(ji,jj,jk) + gdept_n(ji,jj+1,jk) )155 zdep_u = 0.5_wp * ( gdept(ji,jj,jk,Kmm) + gdept(ji+1,jj,jk,Kmm) ) 156 zdep_v = 0.5_wp * ( gdept(ji,jj,jk,Kmm) + gdept(ji,jj+1,jk,Kmm) ) 156 157 ! 157 158 zkh_u = zk_u(ji,jj) * zdep_u ! k * depth … … 179 180 DO jj = 2, jpjm1 180 181 DO ji = 2, jpim1 181 zbot_u = ( gdepw _n(ji,jj,jk+1) + gdepw_n(ji+1,jj,jk+1) ) ! 2 * bottom depth182 zbot_v = ( gdepw _n(ji,jj,jk+1) + gdepw_n(ji,jj+1,jk+1) ) ! 2 * bottom depth182 zbot_u = ( gdepw(ji,jj,jk+1,Kmm) + gdepw(ji+1,jj,jk+1,Kmm) ) ! 2 * bottom depth 183 zbot_v = ( gdepw(ji,jj,jk+1,Kmm) + gdepw(ji,jj+1,jk+1,Kmm) ) ! 2 * bottom depth 183 184 zkb_u = zk_u(ji,jj) * zbot_u ! 2 * k * bottom depth 184 185 zkb_v = zk_v(ji,jj) * zbot_v ! 2 * k * bottom depth 185 186 ! 186 zke3_u = MAX(1.e-8_wp, 2.0_wp * zk_u(ji,jj) * e3u _n(ji,jj,jk)) ! 2k * thickness187 zke3_v = MAX(1.e-8_wp, 2.0_wp * zk_v(ji,jj) * e3v _n(ji,jj,jk)) ! 2k * thickness187 zke3_u = MAX(1.e-8_wp, 2.0_wp * zk_u(ji,jj) * e3u(ji,jj,jk,Kmm)) ! 2k * thickness 188 zke3_v = MAX(1.e-8_wp, 2.0_wp * zk_v(ji,jj) * e3v(ji,jj,jk,Kmm)) ! 2k * thickness 188 189 189 190 ! Depth attenuation .... do u component first.. … … 223 224 DO jj = 2, jpj 224 225 DO ji = fs_2, jpi 225 ze3divh(ji,jj,jk) = ( e2u(ji ,jj) * e3u _n(ji ,jj,jk) * usd(ji ,jj,jk) &226 & - e2u(ji-1,jj) * e3u _n(ji-1,jj,jk) * usd(ji-1,jj,jk) &227 & + e1v(ji,jj ) * e3v _n(ji,jj ,jk) * vsd(ji,jj ,jk) &228 & - e1v(ji,jj-1) * e3v _n(ji,jj-1,jk) * vsd(ji,jj-1,jk) ) * r1_e1e2t(ji,jj)226 ze3divh(ji,jj,jk) = ( e2u(ji ,jj) * e3u(ji ,jj,jk,Kmm) * usd(ji ,jj,jk) & 227 & - e2u(ji-1,jj) * e3u(ji-1,jj,jk,Kmm) * usd(ji-1,jj,jk) & 228 & + e1v(ji,jj ) * e3v(ji,jj ,jk,Kmm) * vsd(ji,jj ,jk) & 229 & - e1v(ji,jj-1) * e3v(ji,jj-1,jk,Kmm) * vsd(ji,jj-1,jk) ) * r1_e1e2t(ji,jj) 229 230 END DO 230 231 END DO … … 307 308 308 309 309 SUBROUTINE sbc_wave( kt )310 SUBROUTINE sbc_wave( kt, Kmm ) 310 311 !!--------------------------------------------------------------------- 311 312 !! *** ROUTINE sbc_wave *** … … 322 323 !!--------------------------------------------------------------------- 323 324 INTEGER, INTENT(in ) :: kt ! ocean time step 325 INTEGER, INTENT(in ) :: Kmm ! ocean time index 324 326 !!--------------------------------------------------------------------- 325 327 ! … … 361 363 ! 362 364 IF( ( ll_st_bv_li .AND. jp_hsw>0 .AND. jp_wmp>0 .AND. jp_usd>0 .AND. jp_vsd>0 ) .OR. & 363 & ( ll_st_peakfr .AND. jp_wfr>0 .AND. jp_usd>0 .AND. jp_vsd>0 ) ) CALL sbc_stokes( )365 & ( ll_st_peakfr .AND. jp_wfr>0 .AND. jp_usd>0 .AND. jp_vsd>0 ) ) CALL sbc_stokes( Kmm ) 364 366 ! 365 367 ENDIF
Note: See TracChangeset
for help on using the changeset viewer.