Changeset 8085 for branches/2016/dev_r6859_LIM3_meltponds/NEMOGCM/NEMO/LIM_SRC_3/limvar.F90

Timestamp:

2017-05-29T16:58:42+02:00 (7 years ago)

Author:

vancop

Message:

More on melt ponds in LIM

File:

• branches/2016/dev_r6859_LIM3_meltponds/NEMOGCM/NEMO/LIM_SRC_3/limvar.F90 (modified) (3 diffs)

branches/2016/dev_r6859_LIM3_meltponds/NEMOGCM/NEMO/LIM_SRC_3/limvar.F90

@@ -526 +526 @@
       !!-------------------------------------------------------------------
       INTEGER  ::   ji, jj, jl, jk   ! dummy loop indices
-      REAL(wp) ::   zsal, zvi, zvs, zei, zes
+      REAL(wp) ::   zsal, zvi, zvs, zei, zes, zvp
       !!-------------------------------------------------------------------
       at_i (:,:) = 0._wp
@@ -566 +566 @@
                zvs  = v_s  (ji,jj,  jl)
                zes  = e_s  (ji,jj,1,jl)
+               zvp  = v_ip (ji,jj  ,jl)
                !-----------------------------------------------------------------
                ! Zap snow energy 
@@ -582 +583 @@
                oa_i (ji,jj,jl) = oa_i (ji,jj,jl) * rswitch
                smv_i(ji,jj,jl) = smv_i(ji,jj,jl) * rswitch
+
+               ! MV MP 2016
+               a_ip (ji,jj,jl) = a_ip (ji,jj,jl) * rswitch
+               v_ip (ji,jj,jl) = v_ip (ji,jj,jl) * rswitch
+               wfx_res(ji,jj)  = wfx_res(ji,jj) - ( v_ip(ji,jj,jl)  - zvp  ) * rhofw * r1_rdtice
+               ! END MV MP 2016
 
                ! update exchanges with ocean