Changeset 5058 for branches/2015/dev_r5044_CNRS_LIM3CLEAN/NEMOGCM/NEMO/OPA_SRC/SBC/sbcice_lim.F90

Timestamp:

2015-02-04T12:25:57+01:00 (9 years ago)

Author:

clem

Message:

LIM3: add a namelist parameter to change the shape of ice thickness categories

File:

• branches/2015/dev_r5044_CNRS_LIM3CLEAN/NEMOGCM/NEMO/OPA_SRC/SBC/sbcice_lim.F90 (modified) (3 diffs)

branches/2015/dev_r5044_CNRS_LIM3CLEAN/NEMOGCM/NEMO/OPA_SRC/SBC/sbcice_lim.F90

@@ -423 +423 @@
       !! ** Purpose :   Initializes the ice thickness distribution
       !! ** Method  :   ...
-      !!    Note    : hi_max(jpl) is here set up to a value close to 7 m for
-      !!              limistate (only) and is changed to 99 m in sbc_lim_init
-      !!------------------------------------------------------------------
+      !! ** input   :   Namelist namiceitd
+      !!-------------------------------------------------------------------
+      INTEGER  ::   ios                 ! Local integer output status for namelist read
+      NAMELIST/namiceitd/ nn_itdshp, rn_itmean
+      !
       INTEGER  ::   jl                   ! dummy loop index
       REAL(wp) ::   zc1, zc2, zc3, zx1   ! local scalars
       REAL(wp) ::   zhmax, znum, zden, zalpha !
       !!------------------------------------------------------------------
-
-      IF(lwp) WRITE(numout,*)
-      IF(lwp) WRITE(numout,*) 'lim_itd_init : Initialization of ice thickness distribution '
-      IF(lwp) WRITE(numout,*) '~~~~~~~~~~~~'
-
-      !------------------------------------------------------------------------------!
-      ! 1) Ice thickness distribution parameters initialization    
-      !------------------------------------------------------------------------------!
-      IF(lwp) WRITE(numout,*) ' Number of ice categories jpl = ', jpl
-
+      !
+      REWIND( numnam_ice_ref )              ! Namelist namiceitd in reference namelist : Parameters for ice
+      READ  ( numnam_ice_ref, namiceitd, IOSTAT = ios, ERR = 903)
+903   IF( ios /= 0 ) CALL ctl_nam ( ios , 'namiceitd in reference namelist', lwp )
+
+      REWIND( numnam_ice_cfg )              ! Namelist namiceitd in configuration namelist : Parameters for ice
+      READ  ( numnam_ice_cfg, namiceitd, IOSTAT = ios, ERR = 904 )
+904   IF( ios /= 0 ) CALL ctl_nam ( ios , 'namiceitd in configuration namelist', lwp )
+      IF(lwm) WRITE ( numoni, namiceitd )
+      !
+      !
+      IF(lwp) THEN                        ! control print
+         WRITE(numout,*)
+         WRITE(numout,*) 'ice_itd : ice cat distribution'
+         WRITE(numout,*) ' ~~~~~~'
+         WRITE(numout,*) '   shape of ice categories distribution                          nn_itdshp = ', nn_itdshp
+         WRITE(numout,*) '   mean ice thickness in the domain (only active if nn_itdshp=2) rn_itmean = ', rn_itmean
+      ENDIF
+
+      !----------------------------------
       !- Thickness categories boundaries 
       !----------------------------------
-      !  Clem - je sais pas encore dans quelle namelist les mettre, ca depend des chgts liés à bdy
-      nn_hibnd  = 2                !  thickness category boundaries: tanh function (1) h^(-alpha) (2)
-      rn_hibnd  = 2.5              !  mean thickness of the domain (used to compute category boundaries, nn_hibnd = 2 only)
+      IF(lwp) WRITE(numout,*)
+      IF(lwp) WRITE(numout,*) 'lim_itd_init : Initialization of ice cat distribution '
+      IF(lwp) WRITE(numout,*) '~~~~~~~~~~~~'
 
       hi_max(:) = 0._wp
 
-      SELECT CASE ( nn_hibnd  )       
+      SELECT CASE ( nn_itdshp  )       
                                    !----------------------
          CASE (1)                  ! tanh function (CICE)
@@ -466 +478 @@
          zalpha = 0.05             ! exponent of the transform function
 
-         zhmax  = 3.*rn_hibnd
+         zhmax  = 3.*rn_itmean
 
          DO jl = 1, jpl 
@@ -479 +491 @@
          hi_mean(jl) = ( hi_max(jl) + hi_max(jl-1) ) * 0.5_wp
       END DO
+
       ! Set hi_max(jpl) to a big value to ensure that all ice is thinner than hi_max(jpl)
       hi_max(jpl) = 99._wp