Changeset 1222 for trunk/NEMO/C1D_SRC/dynnxt_c1d.F90

Timestamp:

2008-10-31T14:55:47+01:00 (16 years ago)

Author:

ctlod

Message:

add a missing jk indice in a DO jk loop, see ticket: #275

File:

• trunk/NEMO/C1D_SRC/dynnxt_c1d.F90 (modified) (1 diff)

trunk/NEMO/C1D_SRC/dynnxt_c1d.F90

@@ -79 +79 @@
       DO jk = 1, jpkm1      ! Time filter and swap of dynamics arrays
          IF( neuler == 0 .AND. kt == nit000 ) THEN      ! Euler (forward) time stepping
-             ub(:,:,:) = un(:,:,:)
-             vb(:,:,:) = vn(:,:,:)
-             un(:,:,:) = ua(:,:,:)
-             vn(:,:,:) = va(:,:,:)
+             ub(:,:,jk) = un(:,:,jk)
+             vb(:,:,jk) = vn(:,:,jk)
+             un(:,:,jk) = ua(:,:,jk)
+             vn(:,:,jk) = va(:,:,jk)
          ELSE                                           ! Leap-frog time stepping
-             ub(:,:,:) = atfp * ( ub(:,:,:) + ua(:,:,:) ) + atfp1 * un(:,:,:)
-             vb(:,:,:) = atfp * ( vb(:,:,:) + va(:,:,:) ) + atfp1 * vn(:,:,:)
-             un(:,:,:) = ua(:,:,:)
-             vn(:,:,:) = va(:,:,:)
+             ub(:,:,jk) = atfp * ( ub(:,:,jk) + ua(:,:,jk) ) + atfp1 * un(:,:,jk)
+             vb(:,:,jk) = atfp * ( vb(:,:,jk) + va(:,:,jk) ) + atfp1 * vn(:,:,jk)
+             un(:,:,jk) = ua(:,:,jk)
+             vn(:,:,jk) = va(:,:,jk)
          ENDIF
       END DO