Changeset 8106 for branches/2016/dev_r6859_LIM3_meltponds/NEMOGCM/NEMO/LIM_SRC_3/limistate.F90

Timestamp:

2017-06-01T09:54:36+02:00 (7 years ago)

Author:

vancop

Message:

Follow up on melt ponds and albedo

File:

• branches/2016/dev_r6859_LIM3_meltponds/NEMOGCM/NEMO/LIM_SRC_3/limistate.F90 (modified) (2 diffs)

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

@@ -383 +383 @@
 
             DO jl = 1, jpl
-               a_ip(:,:,jl) = 0.1 * zswitch(:,:)
-               v_ip(:,:,jl) = 0.1 * zswitch(:,:) * a_ip(:,:,jl)
+
+               a_ip_frac(:,:,jl) = 0.2 *  zswitch(:,:)
+               h_ip     (:,:,jl) = 0.05 * zswitch(:,:)
+               a_ip(:,:,jl)      = a_ip_frac(:,:,jl) * a_i (:,:,jl) 
+               v_ip(:,:,jl)      = h_ip     (:,:,jl) * a_ip(:,:,jl)
+
             END DO
 
          ELSE
 
-            a_ip = 0._wp
-            v_ip = 0._wp
+            a_ip(:,:,:)      = 0._wp
+            v_ip(:,:,:)      = 0._wp
+            a_ip_frac(:,:,:) = 0._wp
+            h_ip     (:,:,:) = 0._wp
 
          ENDIF
@@ -418 +424 @@
          END DO
 
-         ! MV MP 2016
          a_ip(:,:,:)      = 0._wp
          v_ip(:,:,:)      = 0._wp
-         ! END MV MP 2016
+         a_ip_frac(:,:,:) = 0._wp
+         h_ip     (:,:,:) = 0._wp
 
       ENDIF ! ln_limini