Changeset 6717 for branches/2016/dev_r6409_SIMPLIF_2_usrdef/NEMOGCM/NEMO/LIM_SRC_2/limdmp_2.F90

Timestamp:

2016-06-17T12:00:46+02:00 (8 years ago)

Author:

gm

Message:

#1692 - branch SIMPLIF_2_usrdef: numerous improvement in the user defined interface

File:

• branches/2016/dev_r6409_SIMPLIF_2_usrdef/NEMOGCM/NEMO/LIM_SRC_2/limdmp_2.F90 (modified) (2 diffs)

branches/2016/dev_r6409_SIMPLIF_2_usrdef/NEMOGCM/NEMO/LIM_SRC_2/limdmp_2.F90

@@ -144 +144 @@
       
          resto_ice(:,:,:) = 0._wp
-         !      Re-calculate the North and South boundary restoring term
-         !      because those boundaries may change with the prescribed zoom area.
          !
          irelax  = 16                     ! width of buffer zone with respect to close boundary
@@ -156 +154 @@
          ! REM: if there is no ice in the model and in the data, 
          !      no restoring even with non zero resto_ice
-         DO jj = mj0(jpjzoom - 1 + 1), mj1(jpjzoom -1 + irelax)
-            zreltim = zdmpmin + zfactor * ( mjg(jj) - jpjzoom + 1 )
+         DO jj = mj0(1), mj1( irelax)
+            zreltim = zdmpmin + zfactor * mjg(jj)
             resto_ice(:,jj,:) = 1._wp / ( zreltim * 86400._wp )
          END DO
 
          ! North boundary restoring term
-         DO jj =  mj0(jpjzoom -1 + jpjglo - irelax), mj1(jpjzoom - 1 + jpjglo)
-            zreltim = zdmpmin + zfactor * (jpjglo - ( mjg(jj) - jpjzoom + 1 ))
+         DO jj =  mj0(jpjglo - irelax), mj1(jpjglo)
+            zreltim = zdmpmin + zfactor * (jpjglo - mjg(jj))
             resto_ice(:,jj,:) = 1.e0 / ( zreltim * 86400 )
          END DO