Changeset 15579

Timestamp:

2021-12-07T12:48:08+01:00 (2 years ago)

Author:

emmafiedler

Message:

Namelist updates

File:

• NEMO/branches/UKMO/NEMO_4.0.4_ice_strength/src/ICE/icedyn_rdgrft.F90 (modified) (13 diffs)

NEMO/branches/UKMO/NEMO_4.0.4_ice_strength/src/ICE/icedyn_rdgrft.F90

@@ -62 +62 @@
    ! ** namelist (namdyn_rdgrft) **
    LOGICAL  ::   ln_str_H79       ! ice strength parameterization (Hibler79)
-   REAL(wp) ::   rn_pstar         ! determines ice strength, Hibler JPO79
+   REAL(wp) ::   rn_pstar         !    determines ice strength, Hibler JPO79
    LOGICAL  ::   ln_str_R75       ! ice strength parameterization (Rothrock75)
    REAL(wp) ::   rn_Cf            !    ratio of ridging work to PE change in ridging, R75
-   INTEGER  ::   rdg_redistr      ! redistribution of ridged ice (0 = uniform (Hibler 1980), 1 = exponential)
-   REAL(wp) ::   rn_murdg         !    gives e-folding scale of ridged ice (m^.5) for rdg_redistr = 1
+   INTEGER  ::   ismooth          ! smoothing the resistance to deformation (strength)
+   LOGICAL  ::   ln_redist_ufm    ! uniform redistribution of ridged ice (Hibler, 1980) 
+   LOGICAL  ::   ln_redist_exp    ! exponential redistribution of ridged ice
+   REAL(wp) ::   rn_murdg         !    gives e-folding scale of ridged ice (m^.5) for ln_redist_exp
    REAL(wp) ::   rn_csrdg         ! fraction of shearing energy contributing to ridging            
    LOGICAL  ::   ln_partf_lin     ! participation function linear (Thorndike et al. (1975))
@@ -454 +456 @@
                hraft  (ji,jl) = zhi(ji,jl) * zfac
                !
-               IF (rdg_redistr == 0) THEN
+               IF ( ln_redist_ufm ) THEN
                  hrmax  (ji,jl) = 2._wp * zhmean - hrmin(ji,jl)
                  hi_hrdg(ji,jl) = zhi(ji,jl) / zhmean
-               ELSE IF (rdg_redistr == 1) THEN
+               ELSE IF ( ln_redist_exp ) THEN
                  hrexp  (ji,jl) = rn_murdg * SQRT( zhi(ji,jl) )               
                  hi_hrdg(ji,jl) = zhi(ji,jl) / ( hrmin(ji,jl) + hrexp(ji,jl) )     
@@ -703 +705 @@
                  ! Compute the fraction of ridged ice area and volume going to thickness category jl2
                  !
-                 IF (rdg_redistr == 0) THEN ! Hibler 1980 uniform formulation
+                 IF ( ln_redist_ufm ) THEN ! Hibler 1980 uniform formulation
                  !
                   IF( hrmin(ji,jl1) <= hi_max(jl2) .AND. hrmax(ji,jl1) > hi_max(jl2-1) ) THEN
@@ -717 +719 @@
                   ENDIF
                   !
-                 ELSE IF (rdg_redistr == 1) THEN ! exponential formulation
+                 ELSE IF ( ln_redist_exp ) THEN ! exponential formulation
                   !
                   IF (jl2 < jpl) THEN
@@ -743 +745 @@
                   END IF            ! jl2 < jpl
                   ! 
-                 END IF             ! rdg_redistr
+                 END IF             ! ridge redistribution
 
                   ! Compute the fraction of rafted ice area and volume going to thickness category jl2
@@ -826 +828 @@
       !!----------------------------------------------------------------------
       INTEGER             ::   ji, jj, jl  ! dummy loop indices
-      INTEGER             ::   ismooth     ! smoothing the resistance to deformation
       INTEGER             ::   itframe     ! number of time steps for the P smoothing
       REAL(wp)            ::   zp, z1_3    ! local scalars
@@ -842 +843 @@
       !                              !--------------------------------------------------!
         strength(:,:) = rn_pstar * SUM( v_i(:,:,:), dim=3 ) * EXP( -rn_crhg * ( 1._wp - SUM( a_i(:,:,:), dim=3 ) ) )
-        ismooth = 1
+      ! Recommend ismooth = 1
       !                              !--------------------------------------------------!
       ELSE IF( ln_str_R75 ) THEN     ! Ice strength => Rothrock (1975) method           !
@@ -849 +850 @@
         ! Initialise the strength field at each timestep
         strength(:,:) = 0._wp
-        ismooth = 0
+      ! Recommend ismooth = 0
 
         !--------------------------------
@@ -891 +892 @@
               IF ( apartf(ji,jl) > 0._wp ) THEN
                 !
-                IF (rdg_redistr == 0) THEN   ! Uniform distribution of ridged ice               
+                IF ( ln_redist_ufm ) THEN   ! Uniform distribution of ridged ice               
                   h2rdg = (1._wp/3._wp) * (hrmax(ji,jl)**3._wp - hrmin(ji,jl)**3._wp)  &
                            / (hrmax(ji,jl) - hrmin(ji,jl))
                 !
-                ELSE IF (rdg_redistr == 1) THEN   ! Exponential distribution of ridged ice
+                ELSE IF ( ln_redist_exp ) THEN   ! Exponential distribution of ridged ice
         h2rdg = hrmin(ji,jl) * hrmin(ji,jl) & 
                 + 2._wp * hrmin(ji,jl) * hrexp(ji,jl) & 
@@ -939 +940 @@
          !                           !--------------------------------------------------!
          strength(:,:) = 0._wp
-         ismooth = 0
+         !
       ENDIF  ! Ice strength formulation
       !      
@@ -1086 +1087 @@
       INTEGER :: ios                 ! Local integer output status for namelist read
       !!
-      NAMELIST/namdyn_rdgrft/ ln_str_H79, rn_pstar,          &
-         &                    ln_str_R75, rn_Cf,             &
-         &                    rdg_redistr, rn_murdg,         &
+      NAMELIST/namdyn_rdgrft/ ismooth, ln_str_H79, rn_pstar,          &
+         &                    ln_str_R75, rn_Cf,                      &
+         &                    ln_redist_ufm, ln_redist_exp, rn_murdg, &
          &                    rn_crhg, rn_csrdg,             &
          &                    ln_partf_lin, rn_gstar,        &
@@ -1114 +1115 @@
          WRITE(numout,*) '      ice strength parameterization Rothrock (1975)            ln_str_R75   = ', ln_str_R75
          WRITE(numout,*) '            ratio of ridging work to PE change in ridging      rn_Cf        = ', rn_Cf
-         WRITE(numout,*) '      redistribution of ridged ice                             rdg_redistr  = ', rdg_redistr
-         WRITE(numout,*) '            (0 = uniform (Hibler 1980), 1 = exponential)                      '
-         WRITE(numout,*) '            e-folding scale of ridged ice for rdg_redistr=1    rn_murdg     = ', rn_murdg
+         WRITE(numout,*) '      smoothing the resistance to deformation (strength)       ismooth      = ', ismooth
+         WRITE(numout,*) '      uniform redistribution of ridged ice (Hibler, 1980)     ln_redist_ufm = ', ln_redist_ufm
+         WRITE(numout,*) '      exponential redistribution of ridged ice                ln_redist_exp = ', ln_redist_exp
+         WRITE(numout,*) '            e-folding scale of ridged ice for ln_redist_exp    rn_murdg     = ', rn_murdg
          WRITE(numout,*) '      Fraction of shear energy contributing to ridging         rn_csrdg     = ', rn_csrdg 
          WRITE(numout,*) '      linear ridging participation function                    ln_partf_lin = ', ln_partf_lin
@@ -1142 +1144 @@
       ENDIF
       !
-!      IF ( ( rdg_redistr .NE. 0 ) .OR. ( rdg_redistr .NE. 1 ) ) THEN  ! why doesn't this work?
-!         CALL ctl_stop( 'ice_dyn_rdgrft_init: ice redistribution function rdg_redistr must be 0 (uniform) or 1 (exponential)' )
-!      ENDIF
-      !
+      IF ( ( ln_redist_ufm .AND. ln_redist_exp ) .OR. ( .NOT.ln_redist_ufm .AND. .NOT.ln_redist_exp ) ) THEN
+         CALL ctl_stop( 'ice_dyn_rdgrft_init: choose one and only one ice redistribution function (ln_redist_ufm or ln_redist_exp)' )
+      ENDIF
+      !
+      !! Add check for ismooth too, =0, 1 or 2
+
       IF( .NOT. ln_icethd ) THEN
          rn_porordg = 0._wp