Changeset 1103 for trunk/NEMO/LIM_SRC_3/limdyn.F90

Timestamp:

2008-06-10T15:11:32+02:00 (16 years ago)

Author:

rblod

Message:

Fixes in kil3, change results, see ticket #206

File:

• trunk/NEMO/LIM_SRC_3/limdyn.F90 (modified) (3 diffs)

trunk/NEMO/LIM_SRC_3/limdyn.F90

@@ -115 +115 @@
                   i_j1 = i_j1 + 1
                END DO
-               i_j1 = MAX( 1, i_j1-1 )
+               i_j1 = MAX( 1, i_j1-2 )
                IF(ln_ctl) CALL prt_ctl_info( 'lim_dyn  : NH  i_j1 = ', ivar1=i_j1, clinfo2=' ij_jpj = ', ivar2=i_jpj )
                CALL lim_rhg( i_j1, i_jpj )
@@ -125 +125 @@
                   i_jpj = i_jpj - 1
                END DO
-               i_jpj = MIN( jpj, i_jpj+2 )
+               i_jpj = MIN( jpj, i_jpj+1 )
                IF(ln_ctl) CALL prt_ctl_info( 'lim_dyn  : SH  i_j1 = ', ivar1=i_j1, clinfo2=' ij_jpj = ', ivar2=i_jpj )
 
-               CALL lim_rhg( i_j1, i_jpj )
-
-            ELSE                                 ! local domain extends over one hemisphere only
-               !                                 ! Rheology is computed only over the ice cover
-               !                                 ! latitude strip
-               i_j1  = 1
+       CALL lim_rhg( i_j1, i_jpj )
+
+    ELSE                                 ! local domain extends over one hemisphere only
+       !                                 ! Rheology is computed only over the ice cover
+       !                                 ! latitude strip
+       i_j1  = 1
                DO WHILE ( i_j1 <= jpj .AND. zind(i_j1) == FLOAT(jpi) .AND. zmsk(i_j1) /=0 )
                   i_j1 = i_j1 + 1
                END DO
-               i_j1 = MAX( 1, i_j1-1 )
+               i_j1 = MAX( 1, i_j1-2 )
 
                i_jpj  = jpj
@@ -143 +143 @@
                   i_jpj = i_jpj - 1
                END DO
-               i_jpj = MIN( jpj, i_jpj+2)
+               i_jpj = MIN( jpj, i_jpj+1)
 
                IF(ln_ctl) CALL prt_ctl_info( 'lim_dyn  : one hemisphere:  i_j1 = ', ivar1=i_j1, clinfo2=' ij_jpj = ', ivar2=i_jpj )