Changeset 4963

Timestamp:

2014-12-02T17:00:20+01:00 (9 years ago)

Author:

clem

Message:

LIM3 debug

File:

• branches/2014/dev_MERGE_2014/NEMOGCM/NEMO/LIM_SRC_3/limtrp.F90 (modified) (3 diffs)

branches/2014/dev_MERGE_2014/NEMOGCM/NEMO/LIM_SRC_3/limtrp.F90

@@ -393 +393 @@
                      zsmv = smv_i(ji,jj,jl)
                      zes  = e_s  (ji,jj,1,jl)
-                     zei  = SUM( e_i(ji,jj,:,jl) )
+                     zei  = SUM( e_i(ji,jj,1:nlay_i,jl) )
                      zdv  = v_i(ji,jj,jl) - zviold(ji,jj,jl)   
                      !zda = a_i(ji,jj,jl) - zaiold(ji,jj,jl)   
@@ -414 +414 @@
                      smv_i(ji,jj,jl) = zindh * smv_i(ji,jj,jl)
                      e_s(ji,jj,1,jl) = zindh * e_s(ji,jj,1,jl)
-                     e_i(ji,jj,:,jl) = zindh * e_i(ji,jj,:,jl)
+                     e_i(ji,jj,1:nlay_i,jl) = zindh * e_i(ji,jj,1:nlay_i,jl)
 
                      ! Update mass fluxes
@@ -421 +421 @@
                      sfx_res(ji,jj) = sfx_res(ji,jj) - ( smv_i(ji,jj,jl) - zsmv ) * rhoic * r1_rdtice 
                      hfx_res(ji,jj) = hfx_res(ji,jj) + ( e_s(ji,jj,1,jl) - zes ) * unit_fac / area(ji,jj) * r1_rdtice ! W.m-2 <0
-                     hfx_res(ji,jj) = hfx_res(ji,jj) + ( SUM( e_i(ji,jj,:,jl) ) - zei ) * unit_fac / area(ji,jj) * r1_rdtice ! W.m-2 <0
+                     hfx_res(ji,jj) = hfx_res(ji,jj) + ( SUM( e_i(ji,jj,1:nlay_i,jl) ) - zei ) * unit_fac / area(ji,jj) * r1_rdtice ! W.m-2 <0
 
                   ENDIF