Changeset 869 for trunk/NEMO/LIM_SRC_3/limrhg.F90

Timestamp:

2008-03-26T10:21:54+01:00 (16 years ago)

Author:

rblod

Message:

Parallelisation of LIM3. This commit seems to ensure the reproducibility mono/mpp. See ticket #77.

File:

• trunk/NEMO/LIM_SRC_3/limrhg.F90 (modified) (5 diffs)

trunk/NEMO/LIM_SRC_3/limrhg.F90

@@ -314 +314 @@
 
       !-Initialise stress tensor 
-      zs1(:,:)  = stress1_i(:,:)
+      zs1(:,:)  = stress1_i(:,:) 
       zs2(:,:)  = stress2_i(:,:)
       zs12(:,:) = stress12_i(:,:)
@@ -387 +387 @@
               END DO
            END DO
+           CALL lbc_lnk( v_ice1(:,:), 'U', -1. )
+           CALL lbc_lnk( u_ice2(:,:), 'V', -1. )
 
 !CDIR NOVERRCHK
@@ -454 +456 @@
                     &        - dtotel*( (1.0-alphaevp)*ecc2*zs12(ji,jj) - zds(ji,jj) / &
                     &          ( 2.0*deltac(ji,jj) ) * zpreshc(ji,jj))) &
-                    &         / ( 1.0 + alphaevp*ecc2*dtotel )
+                    &         / ( 1.0 + alphaevp*ecc2*dtotel ) 
 
               END DO ! ji
@@ -625 +627 @@
       END DO
 
+      CALL lbc_lnk( u_ice(:,:), 'U', -1. ) 
+      CALL lbc_lnk( v_ice(:,:), 'V', -1. ) 
+
       DO jj = k_j1+1, k_jpj-1 
          DO ji = fs_2, fs_jpim1
@@ -641 +646 @@
       END DO
 
-      CALL lbc_lnk( u_ice2(:,:), 'U', -1. ) 
-      CALL lbc_lnk( v_ice1(:,:), 'V', -1. )
+      CALL lbc_lnk( u_ice2(:,:), 'V', -1. ) 
+      CALL lbc_lnk( v_ice1(:,:), 'U', -1. )
 
       ! Recompute delta, shear and div, inputs for mechanical redistribution