Changeset 4573 for branches

Timestamp:

2014-03-26T09:27:21+01:00 (10 years ago)

Author:

pabouttier

Message:

Fix the case where jpi is odd, in northfold treatment for T-point in div_cur_adj (divcur_tam.F90 module), see Ticket #1270

File:

• branches/2012/dev_v3_4_STABLE_2012/NEMOGCM/NEMO/OPATAM_SRC/DYN/divcur_tam.F90 (modified) (2 diffs)

branches/2012/dev_v3_4_STABLE_2012/NEMOGCM/NEMO/OPATAM_SRC/DYN/divcur_tam.F90

@@ -355 +355 @@
       REAL(KIND=wp), POINTER, DIMENSION(:,:) :: &
          & zwv    ! Workspace
+      REAL(wp) :: &
+         & ztmp
       !!----------------------------------------------------------------------
       IF( nn_timing == 1 )  CALL timing_start('div_cur_adj')
@@ -451 +453 @@
          ELSEIF( nperio == 5 .OR. nperio == 6 ) THEN
             ! north fold ( Grid defined with a F-point pivot) ORCA 0.5 degree
-            DO ji = jpi+2, jpi/2+1, -1
+            DO ji = jpi+2, jpi/2+1, -1 ! if jpi is odd, you can get ji=ijt hence the use of ztmp
                ijt = jpi - ji + 1
-               zwv(ijt,jpjm1) = zwv(ijt,jpjm1) - zwv(ji,jpjm1)
-               zwv(ji ,jpjm1) = 0.0_wp
+               ztmp = zwv(ji, jpjm1)
+               zwv(ji, jpjm1) = 0.0_wp
+               zwv(ijt,jpjm1) = zwv(ijt,jpjm1) - ztmp
             END DO
             DO ji = jpi+2, -1, -1