- Timestamp:
- 2021-03-03T16:04:57+01:00 (3 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/branches/2021/dev_r14273_HPC-02_Daley_Tiling/src/OCE/LDF/ldfslp.F90
r13497 r14574 229 229 !!gm end modif 230 230 END_3D 231 CALL lbc_lnk _multi( 'ldfslp', zwz, 'U', -1.0_wp, zww, 'V', -1.0_wp ) ! lateral boundary conditions231 CALL lbc_lnk( 'ldfslp', zwz, 'U', -1.0_wp, zww, 'V', -1.0_wp ) ! lateral boundary conditions 232 232 ! 233 233 ! !* horizontal Shapiro filter … … 245 245 & + 4.* zww(ji,jj ,jk) ) 246 246 END_2D 247 DO jj = 3, jpj-2 ! other rows248 DO ji = 2, jpim1 ! vector opt.249 uslp(ji,jj,jk) = z1_16 * ( zwz(ji-1,jj-1,jk) + zwz(ji+1,jj-1,jk) &250 & + zwz(ji-1,jj+1,jk) + zwz(ji+1,jj+1,jk) &251 & + 2.*( zwz(ji ,jj-1,jk) + zwz(ji-1,jj ,jk) &252 & + zwz(ji+1,jj ,jk) + zwz(ji ,jj+1,jk) ) &253 & + 4.* zwz(ji ,jj ,jk) )254 vslp(ji,jj,jk) = z1_16 * ( zww(ji-1,jj-1,jk) + zww(ji+1,jj-1,jk) &255 & + zww(ji-1,jj+1,jk) + zww(ji+1,jj+1,jk) &256 & + 2.*( zww(ji ,jj-1,jk) + zww(ji-1,jj ,jk) &257 & + zww(ji+1,jj ,jk) + zww(ji ,jj+1,jk) ) &258 & + 4.* zww(ji,jj ,jk) )259 END DO260 END DO261 247 ! !* decrease along coastal boundaries 262 248 DO_2D( 0, 0, 0, 0 ) … … 303 289 !!gm end modif 304 290 END_3D 305 CALL lbc_lnk _multi( 'ldfslp', zwz, 'T', -1.0_wp, zww, 'T', -1.0_wp ) ! lateral boundary conditions291 CALL lbc_lnk( 'ldfslp', zwz, 'T', -1.0_wp, zww, 'T', -1.0_wp ) ! lateral boundary conditions 306 292 ! 307 293 ! !* horizontal Shapiro filter … … 321 307 & + 4.* zww(ji ,jj ,jk) ) * zcofw 322 308 END_2D 323 DO jj = 3, jpj-2 ! other rows324 DO ji = 2, jpim1 ! vector opt.325 zcofw = wmask(ji,jj,jk) * z1_16326 wslpi(ji,jj,jk) = ( zwz(ji-1,jj-1,jk) + zwz(ji+1,jj-1,jk) &327 & + zwz(ji-1,jj+1,jk) + zwz(ji+1,jj+1,jk) &328 & + 2.*( zwz(ji ,jj-1,jk) + zwz(ji-1,jj ,jk) &329 & + zwz(ji+1,jj ,jk) + zwz(ji ,jj+1,jk) ) &330 & + 4.* zwz(ji ,jj ,jk) ) * zcofw331 332 wslpj(ji,jj,jk) = ( zww(ji-1,jj-1,jk) + zww(ji+1,jj-1,jk) &333 & + zww(ji-1,jj+1,jk) + zww(ji+1,jj+1,jk) &334 & + 2.*( zww(ji ,jj-1,jk) + zww(ji-1,jj ,jk) &335 & + zww(ji+1,jj ,jk) + zww(ji ,jj+1,jk) ) &336 & + 4.* zww(ji ,jj ,jk) ) * zcofw337 END DO338 END DO339 309 ! !* decrease in vicinity of topography 340 310 DO_2D( 0, 0, 0, 0 ) … … 348 318 ! IV. Lateral boundary conditions 349 319 ! =============================== 350 CALL lbc_lnk _multi( 'ldfslp', uslp , 'U', -1.0_wp , vslp , 'V', -1.0_wp , wslpi, 'W', -1.0_wp, wslpj, 'W', -1.0_wp )320 CALL lbc_lnk( 'ldfslp', uslp , 'U', -1.0_wp , vslp , 'V', -1.0_wp , wslpi, 'W', -1.0_wp, wslpj, 'W', -1.0_wp ) 351 321 352 322 IF(sn_cfctl%l_prtctl) THEN … … 689 659 END_2D 690 660 !!gm this lbc_lnk should be useless.... 691 CALL lbc_lnk _multi( 'ldfslp', uslpml , 'U', -1.0_wp , vslpml , 'V', -1.0_wp , wslpiml, 'W', -1.0_wp , wslpjml, 'W', -1.0_wp )661 CALL lbc_lnk( 'ldfslp', uslpml , 'U', -1.0_wp , vslpml , 'V', -1.0_wp , wslpiml, 'W', -1.0_wp , wslpjml, 'W', -1.0_wp ) 692 662 ! 693 663 END SUBROUTINE ldf_slp_mxl … … 757 727 ! END DO 758 728 ! END DO 759 ! CALL lbc_lnk _multi( 'ldfslp', uslp , 'U', -1. ; CALL lbc_lnk( 'ldfslp', vslp , 'V', -1., wslpi, 'W', -1., wslpj, 'W', -1. )729 ! CALL lbc_lnk( 'ldfslp', uslp , 'U', -1. ; CALL lbc_lnk( 'ldfslp', vslp , 'V', -1., wslpi, 'W', -1., wslpj, 'W', -1. ) 760 730 !!gm ENDIF 761 731 ENDIF
Note: See TracChangeset
for help on using the changeset viewer.