Index: branches/2017/dev_r8183_ICEMODEL/NEMOGCM/NEMO/LIM_SRC_3/icestp.F90
===================================================================
--- branches/2017/dev_r8183_ICEMODEL/NEMOGCM/NEMO/LIM_SRC_3/icestp.F90	(revision 8516)
+++ branches/2017/dev_r8183_ICEMODEL/NEMOGCM/NEMO/LIM_SRC_3/icestp.F90	(revision 8517)
@@ -258,16 +258,16 @@
       CALL ice_itd_init                ! ice thickness distribution initialization
       !
-      CALL ice_dyn_init                ! set ice dynamics parameters
-      !
-      CALL ice_rdgrft_init             ! set ice ridging/rafting parameters
-      !
-      CALL ice_rhg_init                ! set ice rheology parameters
-      !
-      CALL ice_adv_init                ! set ice advection parameters
-      !
-      CALL ice_thd_init                ! set ice thermodynics parameters
-      !
-      CALL ice_thd_sal_init            ! set ice salinity parameters
-       
+      IF( ln_icedyn )   THEN
+         CALL ice_dyn_init             ! set ice dynamics parameters
+         CALL ice_rdgrft_init          ! set ice ridging/rafting parameters
+         CALL ice_rhg_init             ! set ice rheology parameters
+         CALL ice_adv_init             ! set ice advection parameters
+      ENDIF
+
+      IF( ln_icethd ) THEN
+         CALL ice_thd_init             ! set ice thermodynics parameters
+         CALL ice_thd_sal_init         ! set ice salinity parameters
+      ENDIF
+   
       ! MV MP 2016
       CALL lim_mp_init                 ! set melt ponds parameters
@@ -283,13 +283,14 @@
       CALL ice_var_glo2eqv
       !
-      CALL ice_update_init                 ! ice surface boundary condition
-      !
-      CALL ice_alb_init                    ! ice surface albedo
-      !
-      CALL ice_dia_init                    ! initialization for diags
-      !
-      fr_i  (:,:)   = at_i(:,:)         ! initialisation of sea-ice fraction
-      tn_ice(:,:,:) = t_su(:,:,:)       ! initialisation of surface temp for coupled simu
-      !
+      CALL ice_update_init             ! ice surface boundary condition
+      !
+      CALL ice_alb_init                ! ice surface albedo
+      !
+      CALL ice_dia_init                ! initialization for diags
+      !
+      fr_i  (:,:)   = at_i(:,:)        ! initialisation of sea-ice fraction
+      tn_ice(:,:,:) = t_su(:,:,:)      ! initialisation of surface temp for coupled simu
+      !
+      !                                ! set max concentration in both hemispheres
       WHERE( gphit(:,:) > 0._wp )   ;   rn_amax_2d(:,:) = rn_amax_n  ! NH
       ELSEWHERE                     ;   rn_amax_2d(:,:) = rn_amax_s  ! SH
@@ -306,11 +307,11 @@
       !!
       !! ** Method  :   Read the namice_run namelist and check the parameter
-      !!              values called at the first timestep (nit000)
+      !!                values called at the first timestep (nit000)
       !!
       !! ** input   :   Namelist namice_run
       !!-------------------------------------------------------------------
       INTEGER  ::   ios                 ! Local integer output status for namelist read
-      NAMELIST/namice_run/ jpl, nlay_i, nlay_s, nn_monocat, rn_amax_n, rn_amax_s, cn_icerst_in, cn_icerst_indir,   &
-         &                cn_icerst_out, cn_icerst_outdir
+      NAMELIST/namice_run/ jpl, nlay_i, nlay_s, nn_monocat, ln_icedyn, ln_icethd, rn_amax_n, rn_amax_s,  &
+         &                 cn_icerst_in, cn_icerst_indir, cn_icerst_out, cn_icerst_outdir
       !!-------------------------------------------------------------------
       !
@@ -333,9 +334,11 @@
          WRITE(numout,*) '      number of snow layers                                  nlay_s = ', nlay_s
          WRITE(numout,*) '      virtual ITD mono-category param (1-4) or not (0)   nn_monocat = ', nn_monocat
+         WRITE(numout,*) '      Ice dynamics       (T) or not (F)                  ln_icedyn  = ', ln_icedyn
+         WRITE(numout,*) '      Ice thermodynamics (T) or not (F)                  ln_icethd  = ', ln_icethd
          WRITE(numout,*) '      maximum ice concentration for NH                              = ', rn_amax_n 
          WRITE(numout,*) '      maximum ice concentration for SH                              = ', rn_amax_s
       ENDIF
       !
-      !                                         !--- check consistency
+      !                                        !--- check consistency
       IF ( jpl > 1 .AND. nn_monocat == 1 ) THEN
          nn_monocat = 0
@@ -349,9 +352,9 @@
       IF( ln_bdy .AND. ln_icediachk )   CALL ctl_warn('ice_run_init: online conservation check does not work with BDY')
       !
-      rdt_ice   = REAL(nn_fsbc) * rdt           !--- sea-ice timestep and inverse
+      rdt_ice   = REAL(nn_fsbc) * rdt          !--- sea-ice timestep and inverse
       r1_rdtice = 1._wp / rdt_ice
-      IF( lwp ) WRITE(numout,*) '   ice timestep rdt_ice  = ', rdt_ice
-      !
-      r1_nlay_i = 1._wp / REAL( nlay_i, wp )    !--- inverse of nlay_i and nlay_s
+      IF( lwp ) WRITE(numout,*) '   ice timestep rdt_ice = ', rdt_ice
+      !
+      r1_nlay_i = 1._wp / REAL( nlay_i, wp )   !--- inverse of nlay_i and nlay_s
       r1_nlay_s = 1._wp / REAL( nlay_s, wp )
       !
