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

Timestamp:

2008-05-13T10:28:52+02:00 (16 years ago)

Author:

rblod

Message:

Correct indentation and print for debug in LIM3, see ticket #134, step I

File:

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

trunk/NEMO/LIM_SRC_3/limmsh.F90

@@ -53 +53 @@
       REAL(wp), DIMENSION(jpi,jpj) ::  &
          zd2d1 , zd1d2       ! Derivative of zh2 (resp. zh1) in the x direction
-         !                   ! (resp. y direction) (defined at the center)
+      !                   ! (resp. y direction) (defined at the center)
       REAL(wp) ::         &
          zh1p  , zh2p   , &  ! Idem zh1, zh2 for the bottom left corner of the grid
@@ -65 +65 @@
          WRITE(numout,*) '~~~~~~~'
       ENDIF
-      
+
       !----------------------------------------------------------                          
       !    Initialization of local and some global (common) variables 
       !------------------------------------------------------------------ 
-      
+
       njeq   = INT( jpj / 2 )   !i bug mpp potentiel
       njeqm1 = njeq - 1 
 
       fcor(:,:) = 2. * omega * SIN( gphit(:,:) * rad )   !  coriolis factor
- 
+
       IF( fcor(1,1) * fcor(1,nlcj) < 0.e0 ) THEN   ! local domain include both hemisphere
          l_jeq = .TRUE.
@@ -97 +97 @@
       !   For each grid, definition of geometric tables 
       !------------------------------------------------------------------
-      
+
       !-------------------
       ! Conventions :    !
@@ -106 +106 @@
       !  3 = corner SW x(i-1/2),y(j-1/2)
       !-------------------
-!!ibug ???
+      !!ibug ???
       akappa(:,:,:,:) = 0.e0
       wght(:,:,:,:) = 0.e0
@@ -112 +112 @@
       tmu(:,:) = 0.e0
       tmv(:,:) = 0.0e0 ! CGrid EVP
-!!i
+      !!i
       ! metric coefficients for sea ice dynamic
       !----------------------------------------
@@ -130 +130 @@
       akappa(:,:,2,1) = zd2d1(:,:) / ( 4.0 * e1t(:,:) * e2t(:,:) )
       akappa(:,:,2,2) =        1.0 / ( 2.0 * e2t(:,:) )
-      
+
       !                                                      ! weights (wght)
       DO jj = 2, jpj
@@ -146 +146 @@
       CALL lbc_lnk( wght(:,:,2,1), 'I', 1. )      ! but it is never used
       CALL lbc_lnk( wght(:,:,2,2), 'I', 1. )
-    
+
       ! Coefficients for divergence of the stress tensor
       !-------------------------------------------------
@@ -209 +209 @@
       CALL lbc_lnk( alambd(:,:,2,1,1,1), 'I', 1. )      !
       CALL lbc_lnk( alambd(:,:,2,1,1,2), 'I', 1. )      !
-            
+
 
       ! Initialization of ice masks
       !----------------------------
-      
+
       tms(:,:) = tmask(:,:,1)      ! ice T-point  : use surface tmask
 
-!      tmu(:,1) = 0.e0
-!      tmu(1,:) = 0.e0
-!      tmv(:,1) = 0.e0
-!      tmv(1,:) = 0.e0
+      !      tmu(:,1) = 0.e0
+      !      tmu(1,:) = 0.e0
+      !      tmv(:,1) = 0.e0
+      !      tmv(1,:) = 0.e0
 
       DO jj = 1, jpj - 1
@@ -226 +226 @@
             tmv(ji,jj) =  tms(ji,jj) * tms(ji,jj+1)
             tmf(ji,jj) =  tms(ji,jj) * tms(ji+1,jj) * tms(ji,jj+1) * &
-                          tms(ji+1,jj+1)
-         END DO
-      END DO
-      
+               tms(ji+1,jj+1)
+         END DO
+      END DO
+
       !--lateral boundary conditions    
       CALL lbc_lnk( tmu(:,:), 'U', 1. )
       CALL lbc_lnk( tmv(:,:), 'V', 1. )
       CALL lbc_lnk( tmf(:,:), 'F', 1. )
-      
+
       ! unmasked and masked area of T-grid cell
       area(:,:) = e1t(:,:) * e2t(:,:)
-      
+
    END SUBROUTINE lim_msh