Changeset 9097 for branches/2017/dev_merge_2017/NEMOGCM/NEMO/OPA_SRC/TRD/trddyn.F90

Timestamp:

2017-12-15T18:36:32+01:00 (6 years ago)

Author:

cetlod

Message:

Use of lbclnk_multi in subdir USR & TRD

File:

• branches/2017/dev_merge_2017/NEMOGCM/NEMO/OPA_SRC/TRD/trddyn.F90 (modified) (3 diffs)

branches/2017/dev_merge_2017/NEMOGCM/NEMO/OPA_SRC/TRD/trddyn.F90

@@ -121 +121 @@
                               z3dx(:,:,:) = 0._wp                  ! U.dxU & V.dyV (approximation)
                               z3dy(:,:,:) = 0._wp
-                              DO jk = 1, jpkm1                 ! no mask as un,vn are masked
+                              DO jk = 1, jpkm1   ! no mask as un,vn are masked
                                  DO jj = 2, jpjm1
                                     DO ji = 2, jpim1
@@ -129 +129 @@
                                  END DO
                               END DO
-                              CALL lbc_lnk( z3dx, 'U', -1. )
-                              CALL lbc_lnk( z3dy, 'V', -1. )
+                              CALL lbc_lnk_multi( z3dx, 'U', -1., z3dy, 'V', -1. )
                               CALL iom_put( "utrd_udx", z3dx  )
                               CALL iom_put( "vtrd_vdy", z3dy  )
@@ -151 +150 @@
                               IF( ln_drgimp ) THEN
                                  ALLOCATE( z3dx(jpi,jpj,jpk) , z3dy(jpi,jpj,jpk) )
-                                 z3dx(:,:,:) = 0._wp   ;   z3dy(:,:,:) = 0._wp  ! after velocity known (now filed at this stage)
-                                 DO jk = 1, jpkm1
+                             z3dx(:,:,:) = 0._wp   ;   z3dy(:,:,:) = 0._wp  ! after velocity known (now filed at this stage)
+                            DO jk = 1, jpkm1
                                     DO jj = 2, jpjm1
                                        DO ji = 2, jpim1
-                                          ikbu = mbku(ji,jj)          ! deepest ocean u- & v-levels
+                                      ikbu = mbku(ji,jj)          ! deepest ocean u- & v-levels
                                           ikbv = mbkv(ji,jj)
                                           z3dx(ji,jj,jk) = 0.5*( rCdU_bot(ji+1,jj)+rCdU_bot(ji,jj) )*un(ji,jj,ikbu)/e3u_n(ji,jj,ikbu)
                                           z3dy(ji,jj,jk) = 0.5*( rCdU_bot(ji,jj+1)+rCdU_bot(ji,jj) )*vn(ji,jj,ikbv)/e3v_n(ji,jj,ikbv)
-                                       END DO
-                                    END DO
-                                 END DO
-                                 CALL lbc_lnk( z3dx, 'U', -1. )   ;   CALL lbc_lnk( z3dy, 'V', -1. )
-                                 CALL iom_put( "utrd_bfr", z3dx )
-                                 CALL iom_put( "vtrd_bfr", z3dy )
+                                    END DO
+                                 END DO
+                              END DO
+                              CALL lbc_lnk_multi( z3dx, 'U', -1., z3dy, 'V', -1. )
+                              CALL iom_put( "utrd_bfr", z3dx )
+                              CALL iom_put( "vtrd_bfr", z3dy )
                                  DEALLOCATE( z3dx , z3dy )
                               ENDIF