Changeset 217 for trunk/NEMO/OPA_SRC/DYN

Timestamp:

2005-03-18T15:57:51+01:00 (19 years ago)

Author:

opalod

Message:

CT : UPDATE152 : optimization: reduce the MPI communications in doing only one CALL lbc_lnk() with a 3D array instead of jpk CALL lbc_lnk() with a 2D array; or suppressing useless CALL lbc_lnk()

File:

• trunk/NEMO/OPA_SRC/DYN/dynspg_fsc.F90 (modified) (4 diffs)

trunk/NEMO/OPA_SRC/DYN/dynspg_fsc.F90

@@ -145 +145 @@
       znugdt =  rnu * grav * z2dt
       znurau =  znugdt * zraur
-      IF( lk_mpp ) THEN
-         ! Mpp : export boundary values of to neighboring processors
-         !!bug :  I don t understand why this only in mpp???? ==> Can be suppressed, no?
-         CALL lbc_lnk( ua, 'U', -1. )
-         CALL lbc_lnk( va, 'V', -1. )
-      ENDIF
 
       ! 1. Surface pressure gradient (now)
@@ -351 +345 @@
       ! Euler (forward) time stepping, no time filter
       IF( neuler == 0 .AND. kt == nit000 ) THEN
-         DO jj = 2, jpjm1            !!bug  if 1, jpj  the lbc_lnk call can be suppress
+         DO jj = 1, jpj
             DO ji = 1, jpi
                ! after free surface elevation
@@ -362 +356 @@
       ELSE
          ! Leap-frog time stepping and time filter
-         DO jj = 2, jpjm1            !!bug  if 1, jpj  the lbc_lnk call can be suppress
+         DO jj = 1, jpj
             DO ji = 1, jpi
                ! after free surface elevation
@@ -373 +367 @@
       ENDIF
 
-      ! Boundary conditions on sshn
-      CALL lbc_lnk( sshn, 'T', 1. )
 
       IF(l_ctl) THEN         ! print sum trends (used for debugging)
          WRITE(numout,*) ' spg - ssh:', SUM( sshn(2:nictl,2:njctl)*tmask(2:nictl,2:njctl,1) )
       ENDIF
-
 
    END SUBROUTINE dyn_spg_fsc