Changeset 1857 for branches/DEV_r1837_mass_heat_salt_fluxes/NEMO/LIM_SRC_2/limmsh_2.F90

Timestamp:

2010-05-03T13:59:46+02:00 (14 years ago)

Author:

gm

Message:

ticket:#665 Reverting previous commit and going back to revision 1850

File:

• branches/DEV_r1837_mass_heat_salt_fluxes/NEMO/LIM_SRC_2/limmsh_2.F90 (modified) (7 diffs)

branches/DEV_r1837_mass_heat_salt_fluxes/NEMO/LIM_SRC_2/limmsh_2.F90

@@ -4 +4 @@
    !! LIM 2.0 ice model :   definition of the ice mesh parameters
    !!======================================================================
-   !! ** History :  LIM  ! 2001-04 (UCL)  original code
-   !!               2.0  ! 2002-08 (C. Ethe, G. Madec) F90, module
-   !!----------------------------------------------------------------------
 #if defined key_lim2
    !!----------------------------------------------------------------------
@@ -13 +10 @@
    !!   lim_msh_2   : definition of the ice mesh
    !!----------------------------------------------------------------------
+   !! * Modules used
    USE phycst
    USE dom_oce
@@ -22 +20 @@
    PRIVATE
 
-   PUBLIC   lim_msh_2   ! routine called by ice_ini_2.F90
-
-   !!----------------------------------------------------------------------
-   !! NEMO/LIM 3.3,  UCL-LOCEAN-IPSL (2010) 
+   !! * Accessibility
+   PUBLIC lim_msh_2      ! routine called by ice_ini_2.F90
+
+   !!----------------------------------------------------------------------
+   !!   LIM 2.0,  UCL-LOCEAN-IPSL (2005) 
    !! $Id$
-   !! Software governed by the CeCILL licence (modipsl/doc/NEMO_CeCILL.txt)
+   !! This software is governed by the CeCILL licence see modipsl/doc/NEMO_CeCILL.txt 
    !!----------------------------------------------------------------------
 
@@ -43 +42 @@
       !!              - Initialization of the ice masks (tmsk, umsk)
       !! 
-      !! Reference  : Deleersnijder et al. Ocean Modelling 100, 7-10 
+      !! ** Refer.  : Deleersnijder et al. Ocean Modelling 100, 7-10 
+      !!
+      !! ** History :
+      !!         original    : 01-04 (LIM)
+      !!         addition    : 02-08 (C. Ethe, G. Madec)
       !!--------------------------------------------------------------------- 
-      INTEGER  ::   ji, jj   ! dummy loop indices
-      REAL(wp) ::   zh1p, zd1d2p, zusden    ! temporary scalars
-      REAL(wp) ::   zh2p, zd2d1p, zusden2   !    -         -
-      REAL(wp), DIMENSION(jpi,jpj) ::   zd2d1 , zd1d2   ! 2D workspace
+      !! * Local variables
+      INTEGER :: ji, jj      ! dummy loop indices
+
+      REAL(wp), DIMENSION(jpi,jpj) ::  &
+         zd2d1 , zd1d2       ! Derivative of zh2 (resp. zh1) in the x direction
+         !                   ! (resp. y direction) (defined at the center)
+      REAL(wp) ::         &
+         zh1p  , zh2p   , &  ! Idem zh1, zh2 for the bottom left corner of the grid
+         zd2d1p, zd1d2p , &  ! Idem zd2d1, zd1d2 for the bottom left corner of the grid
+         zusden, zusden2     ! temporary scalars
       !!---------------------------------------------------------------------
 
@@ -103 +112 @@
       !-------------------
 !!ibug ???
-      akappa(:,:,:,:)     = 0.e0
-      wght  (:,:,:,:)     = 0.e0
+      akappa(:,:,:,:) = 0.e0
+      wght(:,:,:,:) = 0.e0
       alambd(:,:,:,:,:,:) = 0.e0
-      tmu   (:,:)         = 0.e0
+      tmu(:,:) = 0.e0
 !!i
       
@@ -229 +238 @@
          END DO
       END DO
-      CALL lbc_lnk( tmu(:,:), 'I', 1. )      !--lateral boundary conditions    
-      
-      area(:,:) = e1t(:,:) * e2t(:,:)        ! unmasked and masked area of T-grid cell
+      
+      !--lateral boundary conditions    
+      CALL lbc_lnk( tmu(:,:), 'I', 1. )
+      
+      ! unmasked and masked area of T-grid cell
+      area(:,:) = e1t(:,:) * e2t(:,:)
       
    END SUBROUTINE lim_msh_2
@@ -239 +251 @@
    !!   Default option            Dummy Module         NO LIM sea-ice model
    !!----------------------------------------------------------------------
+CONTAINS
+   SUBROUTINE lim_msh_2           ! Dummy routine
+   END SUBROUTINE lim_msh_2
 #endif