Changeset 4072 for branches/2013/dev_r4028_CNRS_LIM3/NEMOGCM/NEMO/LIM_SRC_3/limupdate2.F90

Timestamp:

2013-10-17T15:10:59+02:00 (11 years ago)

Author:

clem

Message:

few rewritings

File:

• branches/2013/dev_r4028_CNRS_LIM3/NEMOGCM/NEMO/LIM_SRC_3/limupdate2.F90 (modified) (5 diffs)

branches/2013/dev_r4028_CNRS_LIM3/NEMOGCM/NEMO/LIM_SRC_3/limupdate2.F90

@@ -39 +39 @@
    USE wrk_nemo         ! work arrays
    USE lib_fortran     ! glob_sum
+   USE timing          ! Timing
 
    IMPLICIT NONE
@@ -95 +96 @@
       REAL(wp), POINTER, DIMENSION(:,:,:) ::   zviold, zvsold, zsmvold   ! old ice volume...
       !!-------------------------------------------------------------------
+      IF( nn_timing == 1 )  CALL timing_start('limupdate2')
 
       CALL wrk_alloc( jpi,jpj,jpl, internal_melt )   ! integer
@@ -558 +560 @@
       !------------
       DO jj = 2, jpjm1
-         DO ji = fs_2, fs_jpim1
+         DO ji = 2, jpim1
             IF ( at_i(ji,jj) .EQ. 0.0 ) THEN ! what to do if there is no ice
                IF ( at_i(ji+1,jj) .EQ. 0.0 ) u_ice(ji,jj)   = 0.0 ! right side
@@ -567 +569 @@
          END DO
       END DO
+      !lateral boundary conditions
+      CALL lbc_lnk( u_ice(:,:), 'U', -1. )
+      CALL lbc_lnk( v_ice(:,:), 'V', -1. )
       !mask velocities
       u_ice(:,:) = u_ice(:,:) * tmu(:,:)
       v_ice(:,:) = v_ice(:,:) * tmv(:,:)
-      !lateral boundary conditions
-      CALL lbc_lnk( u_ice(:,:), 'U', -1. )
-      CALL lbc_lnk( v_ice(:,:), 'V', -1. )
  
       !--------------------------------
@@ -687 +689 @@
       CALL wrk_dealloc( jpi,jpj,jpl,zviold, zvsold, zsmvold )   ! clem
 
+      IF( nn_timing == 1 )  CALL timing_stop('limupdate2')
    END SUBROUTINE lim_update2
 #else