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

Timestamp:

2017-05-24T10:02:23+02:00 (7 years ago)

Author:

vancop

Message:

Quick commit on empirical melt ponds

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

@@ -417 +417 @@
       !--------------------------------------------------------------------
       !
-      IF ( ln_limMP ) THEN
-         a_ip(:,:,:) = 0.1_wp
-         v_ip(:,:,:) = 0.1_wp
-         h_ip(:,:,:) = 0._wp
-         a_ip_frac(:,:,:) = 0._wp
-
-         at_ip(:,:) = 0._wp
-         vt_ip(:,:) = 0._wp
-      ENDIF
+      !IF ( ln_pnd   ) THEN
+      SELECT CASE ( nn_pnd_scheme )
+
+         z1_jpl =  1 / REAL(jpl)
+
+         CASE ( 0 )           !--- Prescribed melt ponds
+
+            DO jl = 1, jpl
+
+               a_ip(:,:,jl) = rn_apnd * z1_jpl * zswitch(:,:)
+               h_ip(:,:,jl) = 0.1 * zswitch(:,:)
+
+            END DO
+
+         CASE ( 1, 2 )        !--- Prognostic melt ponds
+
+            DO jl = 1, jpl
+
+               a_ip(:,:,jl) = 0.1 * zswitch(:,:)
+               h_ip(:,:,jl) = 0.1 * zswitch(:,:)
+
+            END DO
+
+      END SELECT
+
+      v_ip(:,:,:)      = a_ip(:,:,:)  * h_i_p(:,:,:)
+      a_ip_frac(:,:,:) = a_ip(:,:,:) / a_i(:,:,:)
 
       ! END MV MP 2016
@@ -473 +491 @@
 
       ! MV MP 2016
-      IF ( ln_limMP ) THEN
+      IF ( nn_pnd_scheme > 1 ) THEN
          sxap  (:,:,:) = 0._wp    ; sxvp  (:,:,:) = 0._wp 
          syap  (:,:,:) = 0._wp    ; syvp  (:,:,:) = 0._wp