Changeset 834 for trunk/NEMO/LIM_SRC_3/limmsh.F90

Timestamp:

2008-03-07T18:11:35+01:00 (16 years ago)

Author:

ctlod

Message:

Clean comments and useless lines, see ticket:#72

File:

• trunk/NEMO/LIM_SRC_3/limmsh.F90 (modified) (6 diffs)

trunk/NEMO/LIM_SRC_3/limmsh.F90

@@ -6 +6 @@
 #if defined key_lim3
    !!----------------------------------------------------------------------
-   !!   'key_lim3'                                     LIM sea-ice model
+   !!   'key_lim3'                                      LIM3 sea-ice model
    !!----------------------------------------------------------------------
    !!   lim_msh   : definition of the ice mesh
@@ -24 +24 @@
 
    !!----------------------------------------------------------------------
-   !!   LIM 2.0,  UCL-LOCEAN-IPSL (2005) 
+   !!   LIM 3.0,  UCL-ASTR-LOCEAN-IPSL (2008) 
    !! $Header: /home/opalod/NEMOCVSROOT/NEMO/LIM_SRC/limmsh.F90,v 1.5 2005/03/27 18:34:42 opalod Exp $ 
    !! This software is governed by the CeCILL licence see modipsl/doc/NEMO_CeCILL.txt 
@@ -75 +75 @@
       fcor(:,:) = 2. * omega * SIN( gphit(:,:) * rad )   !  coriolis factor
  
-!i    DO jj = 1, jpj
-!i       zmsk(jj) = SUM( tmask(:,jj,:) )   ! = 0          if land  everywhere on a j-line
-!!ii     write(numout,*) jj, zind(jj)
-!i    END DO
-
       IF( fcor(1,1) * fcor(1,nlcj) < 0.e0 ) THEN   ! local domain include both hemisphere
          l_jeq = .TRUE.
@@ -118 +113 @@
       tmv(:,:) = 0.0e0 ! CGrid EVP
 !!i
-      
-      
       ! metric coefficients for sea ice dynamic
       !----------------------------------------
@@ -169 +162 @@
                &   + e2t(ji-1,jj-1) * wght(ji,jj,1,1)
 
-! better written but change the last digit and thus solver in less than 100 timestep
-!           zh1p  = e1t(ji-1,jj  ) * wght(ji,jj,1,2) + e1t(ji,jj  ) * wght(ji,jj,2,2)   &
-!              &  + e1t(ji-1,jj-1) * wght(ji,jj,1,1) + e1t(ji,jj-1) * wght(ji,jj,2,1) 
-
-!           zh2p  = e2t(ji-1,jj  ) * wght(ji,jj,1,2) + e2t(ji,jj  ) * wght(ji,jj,2,2)   &
-!              &  + e2t(ji-1,jj-1) * wght(ji,jj,1,1) + e2t(ji,jj-1) * wght(ji,jj,2,1)
-
-!!ibug =0   zusden = 1.0 / ( zh1p * zh2p * 4.e0 )
             zusden = 1.0 / MAX( zh1p * zh2p * 4.e0 , 1.e-20 )
             zusden2 = zusden * 2.0 
@@ -231 +216 @@
       tms(:,:) = tmask(:,:,1)      ! ice T-point  : use surface tmask
 
-!i here we can use umask with a i and j shift of -1,-1
       tmu(:,1) = 0.e0
       tmu(1,:) = 0.e0
       tmv(:,1) = 0.e0
       tmv(1,:) = 0.e0
-!     DO jj = 2, jpj               ! ice U.V-point: computed from ice T-point mask
-!        DO ji = 2, jpim1
-! bug
+
       DO jj = 1, jpj - 1
          DO ji = 2 , jpi - 1