Changeset 8125 for branches/2016/dev_r6859_LIM3_meltponds/NEMOGCM/NEMO/LIM_SRC_3/limmp.F90

Timestamp:

2017-06-02T08:53:53+02:00 (7 years ago)

Author:

vancop

Message:

Melt pond podifications

File:

• branches/2016/dev_r6859_LIM3_meltponds/NEMOGCM/NEMO/LIM_SRC_3/limmp.F90 (modified) (6 diffs)

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

@@ -83 +83 @@
       !!-------------------------------------------------------------------
       INTEGER  ::   ios                 ! Local integer output status for namelist read
-      NAMELIST/namicemp/  ln_pnd, nn_pnd_scheme, nn_pnd_cpl, rn_apnd, rn_hpnd
+      NAMELIST/namicemp/  ln_pnd, ln_pnd_rad, ln_pnd_fw, nn_pnd_scheme, rn_apnd, rn_hpnd
       !!-------------------------------------------------------------------
 
@@ -95 +95 @@
       IF(lwm) WRITE ( numoni, namicemp )
       
-      ! Shortcut for radiatively active melt ponds
-      ln_pnd_rad = ln_pnd .AND. ( ( nn_pnd_cpl == 1 ) .OR. ( nn_pnd_cpl == 2 ) )
-
       IF(lwp) THEN                        ! control print
          WRITE(numout,*)
          WRITE(numout,*) 'lim_mp_init : ice parameters for melt ponds'
          WRITE(numout,*) '~~~~~~~~~~~~'
-         WRITE(numout,*)'    Activate melt ponds                                         ln_pnd        = ', ln_pnd
-         WRITE(numout,*)'    Type of melt pond scheme =0 presc, =1 empirical = 2 topo    nn_pnd_scheme = ', nn_pnd_scheme
-         WRITE(numout,*)'    Type of melt pond coupling =0 pass., =1 full, =2 rad, 3=fw  nn_pnd_cpl    = ', nn_pnd_cpl
-         WRITE(numout,*)'    Prescribed pond fraction                                    rn_apnd       = ', rn_apnd
-         WRITE(numout,*)'    Prescribed pond depth                                       rn_hpnd       = ', rn_hpnd
-         WRITE(numout,*)'    Melt ponds radiatively active                               ln_pnd_rad    = ', ln_pnd_rad
+         WRITE(numout,*) '   Active melt ponds                                           ln_pnd        = ', ln_pnd
+         WRITE(numout,*) '   Active melt ponds radiative coupling                        ln_pnd_rad    = ', ln_pnd_rad
+         WRITE(numout,*) '   Active melt ponds freshwater coupling                       ln_pnd_fw     = ', ln_pnd_fw
+         WRITE(numout,*) '   Type of melt pond scheme =0 presc, =1 empirical = 2 topo    nn_pnd_scheme = ', nn_pnd_scheme
+         WRITE(numout,*) '   Prescribed pond fraction                                    rn_apnd       = ', rn_apnd
+         WRITE(numout,*) '   Prescribed pond depth                                       rn_hpnd       = ', rn_hpnd
       ENDIF
 
@@ -113 +110 @@
          WRITE(numout,*)
          WRITE(numout,*) ' Melt ponds are not activated '
-         WRITE(numout,*) ' nn_pnd_scheme, nn_pnd_cpl, rn_apnd and rn_hpnd are set to zero '
+         WRITE(numout,*) ' ln_pnd_rad and ln_pnd_fw set to .FALSE. '
+         WRITE(numout,*) ' nn_pnd_scheme, rn_apnd, rn_hpnd set to zero '
+         ln_pnd_rad    = .FALSE.
+         ln_pnd_fw     = .FALSE.
          nn_pnd_scheme = 0
-         nn_pnd_cpl    = 0
          rn_apnd       = 0._wp
          rn_hpnd       = 0._wp
+
+         IF(lwp) THEN                     ! control print
+            WRITE(numout,*) '    Active melt ponds radiative coupling                        ln_pnd_rad    = ', ln_pnd_rad
+            WRITE(numout,*) '    Active melt ponds freshwater coupling                       ln_pnd_fw     = ', ln_pnd_fw
+            WRITE(numout,*) '    Type of melt pond scheme =0 presc, =1 empirical = 2 topo    nn_pnd_scheme = ', nn_pnd_scheme
+            WRITE(numout,*) '    Prescribed pond fraction                                    rn_apnd       = ', rn_apnd
+            WRITE(numout,*) '    Prescribed pond depth                                       rn_hpnd       = ', rn_hpnd
+         ENDIF
+      ENDIF
+
+      IF ( ln_pnd .AND. ( nn_pnd_scheme == 2 ) .AND. ( jpl == 1 ) ) THEN
+         WRITE(numout,*) ' Topographic melt ponds are incompatible with jpl = 1 '
+         WRITE(numout,*) ' Run aborted '
+         CALL ctl_stop( 'STOP', 'lim_mp_init: uncompatible options, reset namelist_ice_ref ' )
       ENDIF
 
@@ -314 +327 @@
                h_ip(ji,jj,jl)      = 0._wp
 
-               IF ( ( nn_pnd_cpl == 1 ) .OR. ( nn_pnd_cpl == 3 ) ) & !--- Give freshwater to the ocean
+               IF ( ln_pnd_fw ) & !--- Give freshwater to the ocean
                   wfx_pnd(ji,jj)   = wfx_pnd(ji,jj) + v_ip(ji,jj,jl) 
 
@@ -332 +345 @@
                !--- Dump meltwater due to refreezing ( of course this is wrong
                !--- but this parameterization is too simple )
-               IF ( ( nn_pnd_cpl == 1 ) .OR. ( nn_pnd_cpl == 3 ) ) THEN
-
+               IF ( ln_pnd_fw ) &
                   wfx_pnd(ji,jj)   = wfx_pnd(ji,jj) + rhofw * ( v_ip(ji,jj,jl) - zvpold ) * r1_rdtice
-
-               ENDIF
 
                a_ip_frac(ji,jj,jl) = MIN( 1._wp , SQRT( v_ip(ji,jj,jl) * z1_zpnd_aspect / a_i(ji,jj,jl) ) )
@@ -359 +369 @@
        !--- Remove retained meltwater from surface fluxes 
 
-       IF ( ( nn_pnd_cpl .EQ. 1 ) .OR. ( nn_pnd_cpl .EQ. 3 ) ) THEN
+       IF ( ln_pnd_fw ) THEN
 
            wfx_snw(:,:) = wfx_snw(:,:) *  ( 1. - zrmin - ( zrmax - zrmin ) * at_i(:,:) )