Index: NEMO/trunk/src/ICE/ice.F90
===================================================================
--- NEMO/trunk/src/ICE/ice.F90	(revision 10527)
+++ NEMO/trunk/src/ICE/ice.F90	(revision 10531)
@@ -114,5 +114,5 @@
    INTEGER           , PUBLIC ::   nlay_i           !: number of ice  layers 
    INTEGER           , PUBLIC ::   nlay_s           !: number of snow layers 
-   INTEGER           , PUBLIC ::   nn_virtual_itd   !: virtual ITD mono-category parameterizations (1-3) or not (0)
+   LOGICAL           , PUBLIC ::   ln_virtual_itd   !: virtual ITD mono-category parameterization (T) or not (F)
    LOGICAL           , PUBLIC ::   ln_icedyn        !: flag for ice dynamics (T) or not (F)
    LOGICAL           , PUBLIC ::   ln_icethd        !: flag for ice thermo   (T) or not (F)
Index: NEMO/trunk/src/ICE/icedyn_rdgrft.F90
===================================================================
--- NEMO/trunk/src/ICE/icedyn_rdgrft.F90	(revision 10527)
+++ NEMO/trunk/src/ICE/icedyn_rdgrft.F90	(revision 10531)
@@ -586,12 +586,4 @@
                   &                                      + rhos * vsrft(ji) * ( 1._wp - rn_fsnwrft ) ) * r1_rdtice
 
-               ! Put the melt pond water into the ocean
-               ! clem: I think the following lines must be commented since there
-               !       is no net mass flux between melt ponds and the ocean (see icethd_pnd.F90 for ex.)
-               !IF ( ln_pnd_fwb ) THEN
-               !   wfx_pnd_1d(ji) = wfx_pnd_1d(ji) + ( rhow * vprdg(ji) * ( 1._wp - rn_fpndrdg )   &        ! fresh water source for ocean
-               !      &                              + rhow * vprft(ji) * ( 1._wp - rn_fpndrft ) ) * r1_rdtice
-               !ENDIF
-
                ! virtual salt flux to keep salinity constant
                IF( nn_icesal /= 2 )  THEN
Index: NEMO/trunk/src/ICE/iceforcing.F90
===================================================================
--- NEMO/trunk/src/ICE/iceforcing.F90	(revision 10527)
+++ NEMO/trunk/src/ICE/iceforcing.F90	(revision 10531)
@@ -145,5 +145,5 @@
          SELECT CASE( nice_jules )
          CASE( np_jules_ACTIVE )  !    compute conduction flux and surface temperature (as in Jules surface module)
-                                  CALL blk_ice_qcn    ( nn_virtual_itd, t_su, t_bo, h_s, h_i )
+                                  CALL blk_ice_qcn    ( ln_virtual_itd, t_su, t_bo, h_s, h_i )
          END SELECT
       CASE ( jp_purecpl )         !--- coupled formulation
Index: NEMO/trunk/src/ICE/icestp.F90
===================================================================
--- NEMO/trunk/src/ICE/icestp.F90	(revision 10527)
+++ NEMO/trunk/src/ICE/icestp.F90	(revision 10531)
@@ -297,5 +297,5 @@
       INTEGER  ::   ios                 ! Local integer
       !!
-      NAMELIST/nampar/ jpl, nlay_i, nlay_s, nn_virtual_itd, ln_icedyn, ln_icethd, rn_amax_n, rn_amax_s,  &
+      NAMELIST/nampar/ jpl, nlay_i, nlay_s, ln_virtual_itd, ln_icedyn, ln_icethd, rn_amax_n, rn_amax_s,  &
          &             cn_icerst_in, cn_icerst_indir, cn_icerst_out, cn_icerst_outdir
       !!-------------------------------------------------------------------
@@ -317,5 +317,5 @@
          WRITE(numout,*) '         number of ice  layers                               nlay_i    = ', nlay_i
          WRITE(numout,*) '         number of snow layers                               nlay_s    = ', nlay_s
-         WRITE(numout,*) '         virtual ITD param for jpl=1 (1-3) or not (0)   nn_virtual_itd = ', nn_virtual_itd
+         WRITE(numout,*) '         virtual ITD param for jpl=1 (T) or not (F)     ln_virtual_itd = ', ln_virtual_itd
          WRITE(numout,*) '         Ice dynamics       (T) or not (F)                   ln_icedyn = ', ln_icedyn
          WRITE(numout,*) '         Ice thermodynamics (T) or not (F)                   ln_icethd = ', ln_icethd
@@ -324,8 +324,8 @@
       ENDIF
       !                                        !--- check consistency
-      IF ( jpl > 1 .AND. nn_virtual_itd == 1 ) THEN
-         nn_virtual_itd = 0
+      IF ( jpl > 1 .AND. ln_virtual_itd ) THEN
+         ln_virtual_itd = .FALSE.
          IF(lwp) WRITE(numout,*)
-         IF(lwp) WRITE(numout,*) '   nn_virtual_itd forced to 0 as jpl>1, no need with multiple categories to emulate them'
+         IF(lwp) WRITE(numout,*) '   ln_virtual_itd forced to false as jpl>1, no need with multiple categories to emulate them'
       ENDIF
       !
Index: NEMO/trunk/src/ICE/icethd.F90
===================================================================
--- NEMO/trunk/src/ICE/icethd.F90	(revision 10527)
+++ NEMO/trunk/src/ICE/icethd.F90	(revision 10531)
@@ -233,13 +233,6 @@
                               CALL ice_thd_temp                     ! --- temperature update --- !
             !
-!!gm please create a new logical (l_thd_mono or a better explicit name) set one for all in icestp.F90 module
-!!gm        l_thd_mono = ln_icedH .AND. ( ( nn_virtual_itd == 1 .OR. nn_virtual_itd == 4 ) .AND. jpl == 1 )
-!!gm        by the way, the different options associated with nn_virtual_itd =1 to 4  are quite impossible to identify
-!!gm        more comment to add when ready the namelist, with an explicit print in the ocean.output
-            IF( ln_icedH ) THEN
-               IF ( ( nn_virtual_itd == 1 .OR. nn_virtual_itd == 3 ) .AND. jpl == 1 ) THEN
-                              CALL ice_thd_mono                      ! --- extra lateral melting if virtual_itd --- !
-               END IF
-            END IF
+            IF( ln_icedH .AND. ln_virtual_itd ) &
+               &              CALL ice_thd_mono                     ! --- extra lateral melting if virtual_itd --- !
             !
             IF( ln_icedA )    CALL ice_thd_da                       ! --- lateral melting --- !
Index: NEMO/trunk/src/ICE/icethd_zdf_bl99.F90
===================================================================
--- NEMO/trunk/src/ICE/icethd_zdf_bl99.F90	(revision 10527)
+++ NEMO/trunk/src/ICE/icethd_zdf_bl99.F90	(revision 10531)
@@ -265,7 +265,5 @@
          zghe(1:npti) = 1._wp
          !
-         SELECT CASE ( nn_virtual_itd )
-         !
-         CASE ( 1 , 2 )
+         IF( ln_virtual_itd ) THEN
             !
             zepsilon = 0.1_wp
@@ -277,5 +275,5 @@
             END DO
             !
-         END SELECT
+         ENDIF
          !
          !-----------------
