- Timestamp:
- 2014-03-26T09:27:21+01:00 (10 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2012/dev_v3_4_STABLE_2012/NEMOGCM/NEMO/OPATAM_SRC/DYN/divcur_tam.F90
r3627 r4573 355 355 REAL(KIND=wp), POINTER, DIMENSION(:,:) :: & 356 356 & zwv ! Workspace 357 REAL(wp) :: & 358 & ztmp 357 359 !!---------------------------------------------------------------------- 358 360 IF( nn_timing == 1 ) CALL timing_start('div_cur_adj') … … 451 453 ELSEIF( nperio == 5 .OR. nperio == 6 ) THEN 452 454 ! north fold ( Grid defined with a F-point pivot) ORCA 0.5 degree 453 DO ji = jpi+2, jpi/2+1, -1 455 DO ji = jpi+2, jpi/2+1, -1 ! if jpi is odd, you can get ji=ijt hence the use of ztmp 454 456 ijt = jpi - ji + 1 455 zwv(ijt,jpjm1) = zwv(ijt,jpjm1) - zwv(ji,jpjm1) 456 zwv(ji ,jpjm1) = 0.0_wp 457 ztmp = zwv(ji, jpjm1) 458 zwv(ji, jpjm1) = 0.0_wp 459 zwv(ijt,jpjm1) = zwv(ijt,jpjm1) - ztmp 457 460 END DO 458 461 DO ji = jpi+2, -1, -1
Note: See TracChangeset
for help on using the changeset viewer.