Changeset 3816

Timestamp:

2013-02-20T16:30:12+01:00 (11 years ago)

Author:

clevy

Message:

Bugfix on articifial ice source, see ticket:#1054

File:

• branches/2012/dev_v3_4_STABLE_2012/NEMOGCM/NEMO/LIM_SRC_3/limupdate.F90 (modified) (2 diffs)

branches/2012/dev_v3_4_STABLE_2012/NEMOGCM/NEMO/LIM_SRC_3/limupdate.F90

@@ -81 +81 @@
       INTEGER ::   i_ice_switch
       INTEGER ::   ind_im, layer      ! indices for internal melt
-      REAL(wp) ::   zweight, zesum, zhimax, z_da_i, z_dv_i
+      REAL(wp) ::   zweight, zesum, zhimax, z_da_i
       REAL(wp) ::   zindb, zindsn, zindic, zacrith
       REAL(wp) ::   zrtt, zindg, zh, zdvres, zviold
@@ -823 +823 @@
                zindb   =  MAX( rzero, SIGN( rone, v_i(ji,jj,jl) ) ) 
                z_da_i = a_i(ji,jj,jl) * z_da_ex / MAX( at_i(ji,jj), epsi06 ) * zindb
-               z_dv_i = v_i(ji,jj,jl) * z_da_i  / MAX( at_i(ji,jj), epsi06 )
                a_i(ji,jj,jl) = a_i(ji,jj,jl) - z_da_i
-               v_i(ji,jj,jl) = v_i(ji,jj,jl) + z_dv_i
             END DO