Changeset 9118 for branches/2017/dev_merge_2017/NEMOGCM/NEMO/LIM_SRC_3/icevar.F90

Timestamp:

2017-12-18T16:17:44+01:00 (6 years ago)

Author:

clem

Message:

debug ice test cases with thermo or dynamics only

File:

• branches/2017/dev_merge_2017/NEMOGCM/NEMO/LIM_SRC_3/icevar.F90 (modified) (3 diffs)

branches/2017/dev_merge_2017/NEMOGCM/NEMO/LIM_SRC_3/icevar.F90

@@ -124 +124 @@
          tm_si(:,:) = SUM( t_si(:,:,:) * a_i(:,:,:) , dim=3 ) * z1_at_i(:,:)
          om_i (:,:) = SUM( oa_i(:,:,:)              , dim=3 ) * z1_at_i(:,:)
+         sm_i (:,:) = SUM( sv_i(:,:,:)              , dim=3 ) * z1_vt_i(:,:)
          !
          tm_i(:,:) = 0._wp
-         sm_i(:,:) = 0._wp
          DO jl = 1, jpl
             DO jk = 1, nlay_i
                tm_i(:,:) = tm_i(:,:) + r1_nlay_i * t_i (:,:,jk,jl) * v_i(:,:,jl) * z1_vt_i(:,:)
-               sm_i(:,:) = sm_i(:,:) + r1_nlay_i * sz_i(:,:,jk,jl) * v_i(:,:,jl) * z1_vt_i(:,:)
             END DO
          END DO
@@ -298 +297 @@
          !            !---------------------------------------!
          sz_i(:,:,:,:) = rn_icesal
-         s_i(:,:,:)   = rn_icesal
+         s_i (:,:,:)   = rn_icesal
          !
          !            !---------------------------------------------!
@@ -313 +312 @@
          !                                      ! Slope of the linear profile 
          WHERE( h_i(:,:,:) > epsi20 )   ;   z_slope_s(:,:,:) = 2._wp * s_i(:,:,:) / h_i(:,:,:)
-         ELSEWHERE                       ;   z_slope_s(:,:,:) = 0._wp
+         ELSEWHERE                      ;   z_slope_s(:,:,:) = 0._wp
          END WHERE
          !