Changeset 8378 for branches/2017/dev_r8183_ICEMODEL/NEMOGCM/NEMO/LIM_SRC_3/limthd.F90

Timestamp:

2017-07-26T15:55:59+02:00 (7 years ago)

Author:

clem

Message:

some cleaning

File:

• branches/2017/dev_r8183_ICEMODEL/NEMOGCM/NEMO/LIM_SRC_3/limthd.F90 (modified) (3 diffs)

branches/2017/dev_r8183_ICEMODEL/NEMOGCM/NEMO/LIM_SRC_3/limthd.F90

@@ -42 +42 @@
    USE lbclnk         ! lateral boundary condition - MPP links
    USE lib_mpp        ! MPP library
-   USE wrk_nemo       ! work arrays
    USE lib_fortran    ! Fortran utilities (allows no signed zero when 'key_nosignedzero' defined)  
    USE timing         ! Timing
@@ -234 +233 @@
          IF( nidx > 0 ) THEN  ! If there is no ice, do nothing.
             !                                                                
-            s_i_new   (1:nidx) = 0._wp ; dh_s_tot (1:nidx) = 0._wp  ! --- some init --- !
-            dh_i_surf (1:nidx) = 0._wp ; dh_i_bott(1:nidx) = 0._wp
-            dh_snowice(1:nidx) = 0._wp ; dh_i_sub (1:nidx) = 0._wp
-
                               CALL lim_thd_1d2d( jl, 1 )            ! --- Move to 1D arrays --- !
             !
@@ -246 +241 @@
                WHERE( ht_s_1d(1:nidx)>0._wp ) e_s_1d(1:nidx,jk) = e_s_1d(1:nidx,jk) / (ht_s_1d(1:nidx) * a_i_1d(1:nidx)) * nlay_s
             ENDDO
+            !
+            s_i_new   (1:nidx) = 0._wp ; dh_s_tot (1:nidx) = 0._wp  ! --- some init --- !  (important to have them here) 
+            dh_i_surf (1:nidx) = 0._wp ; dh_i_bott(1:nidx) = 0._wp
+            dh_snowice(1:nidx) = 0._wp ; dh_i_sub (1:nidx) = 0._wp
             !
             IF( ln_limdH )    CALL lim_thd_dif                      ! --- Ice/Snow Temperature profile --- !