Changeset 8306 for branches/2017/dev_r8183_ICEMODEL/NEMOGCM/NEMO/OPA_SRC/BDY/bdyice_lim.F90

Timestamp:

2017-07-10T12:18:03+02:00 (7 years ago)

Author:

clem

Message:

step1: remove LIM2 from the code

File:

• branches/2017/dev_r8183_ICEMODEL/NEMOGCM/NEMO/OPA_SRC/BDY/bdyice_lim.F90 (modified) (6 diffs)

branches/2017/dev_r8183_ICEMODEL/NEMOGCM/NEMO/OPA_SRC/BDY/bdyice_lim.F90

@@ -2 +2 @@
    !!======================================================================
    !!                       ***  MODULE  bdyice_lim  ***
-   !! Unstructured Open Boundary Cond. :  Open boundary conditions for sea-ice (LIM2 and LIM3)
+   !! Unstructured Open Boundary Cond. :  Open boundary conditions for sea-ice (LIM3)
    !!======================================================================
    !!  History :  3.3  !  2010-09 (D. Storkey)  Original code
@@ -8 +8 @@
    !!              -   !  2012-01 (C. Rousset)  add lim3 and remove useless jk loop 
    !!----------------------------------------------------------------------
-#if defined key_lim2 || defined key_lim3
-   !!----------------------------------------------------------------------
-   !!   'key_lim2'                                                 LIM-2 sea ice model
+#if defined key_lim3
+   !!----------------------------------------------------------------------
    !!   'key_lim3'                                                 LIM-3 sea ice model
    !!----------------------------------------------------------------------
@@ -20 +19 @@
    USE eosbn2          ! equation of state
    USE oce             ! ocean dynamics and tracers variables
-#if defined key_lim2
-   USE par_ice_2
-   USE ice_2           ! LIM_2 ice variables
-   USE dom_ice_2       ! sea-ice domain
-#elif defined key_lim3
+#if defined key_lim3
    USE ice             ! LIM_3 ice variables
    USE limvar
@@ -55 +50 @@
       !!                  ***  SUBROUTINE bdy_ice_lim  ***
       !!
-      !! ** Purpose : - Apply open boundary conditions for ice (LIM2 and LIM3)
+      !! ** Purpose : - Apply open boundary conditions for ice (LIM3)
       !!
       !!----------------------------------------------------------------------
@@ -110 +105 @@
       REAL(wp) ::   zwgt, zwgt1        ! local scalar
       REAL(wp) ::   ztmelts, zdh
-#if  defined key_lim2 && ! defined key_lim2_vp && defined key_agrif
-     USE ice_2, vt_s => hsnm
-     USE ice_2, vt_i => hicm
-#endif
       !!------------------------------------------------------------------------------
       !
@@ -120 +111 @@
       jgrd = 1      ! Everything is at T-points here
       !
-#if defined key_lim2
-      DO jb = 1, idx%nblenrim(jgrd)
-         ji    = idx%nbi(jb,jgrd)
-         jj    = idx%nbj(jb,jgrd)
-         zwgt  = idx%nbw(jb,jgrd)
-         zwgt1 = 1.e0 - idx%nbw(jb,jgrd)
-         frld (ji,jj) = ( frld (ji,jj) * zwgt1 + dta%frld (jb) * zwgt ) * tmask(ji,jj,1)     ! Leads fraction 
-         hicif(ji,jj) = ( hicif(ji,jj) * zwgt1 + dta%hicif(jb) * zwgt ) * tmask(ji,jj,1)     ! Ice depth 
-         hsnif(ji,jj) = ( hsnif(ji,jj) * zwgt1 + dta%hsnif(jb) * zwgt ) * tmask(ji,jj,1)     ! Snow depth
-      END DO 
-
-      CALL lbc_bdy_lnk( frld,  'T', 1., ib_bdy )                                         ! lateral boundary conditions
-      CALL lbc_bdy_lnk( hicif, 'T', 1., ib_bdy )
-      CALL lbc_bdy_lnk( hsnif, 'T', 1., ib_bdy )
-
-      vt_i(:,:) = hicif(:,:) * frld(:,:)
-      vt_s(:,:) = hsnif(:,:) * frld(:,:)
-      !
-#elif defined key_lim3
+#if defined key_lim3
 
       DO jl = 1, jpl