Changeset 15278 for NEMO/branches/UKMO/tools_r4.0-HEAD_dev_MEs/DOMAINcfg/src/mes.F90

Timestamp:

2021-09-22T18:44:26+02:00 (18 months ago)

Author:

dbruciaferri

Message:

reorganising code in more general structure and adding partial-steps for MEs

File:

• NEMO/branches/UKMO/tools_r4.0-HEAD_dev_MEs/DOMAINcfg/src/mes.F90 (modified) (4 diffs)

NEMO/branches/UKMO/tools_r4.0-HEAD_dev_MEs/DOMAINcfg/src/mes.F90

@@ -44 +44 @@
    REAL(wp)           :: rn_e_al(max_nn_env)     ! Array specifing stretching parameter rn_alpha of SF12 for
                                                  ! each vertical sub-zone
-   LOGICAL, PUBLIC    :: ln_loc_mes              ! To use localised MEs (.TRUE.) or not (.FALSE.
+   LOGICAL, PUBLIC    :: ln_pst_mes              ! To apply partial steps to MEs (.TRUE.) or not (.FALSE.).
+   LOGICAL, PUBLIC    :: ln_pst_l2g              ! To apply partial steps to the transition zone (.TRUE.) 
+                                                 ! or not (.FALSE.) when ln_loc_zgr = .TRUE. 
+   REAL, PUBLIC       :: rn_e3pst_min            ! partial steps parameters (require ln_pst_mes = .TRUE.)
+   REAL, PUBLIC       :: rn_e3pst_rat            ! partial steps parameters (require ln_pst_mes = .TRUE.)
    !
    REAL(wp), POINTER, DIMENSION(:,:,:) :: envlt  ! array for the envelopes
@@ -669 +673 @@
       REAL(wp)                     :: rn_bot_max ! max depth of ocean bottom (= ocean depth) (>0) (m)
 
-      NAMELIST/namzgr_mes/rn_bot_min , rn_bot_max , ln_envl, &
-                          nn_strt    , nn_slev    , rn_e_hc, &
-                          rn_e_th    , rn_e_bb    , rn_e_ba, &
-                          rn_e_al    , ln_loc_mes
+      NAMELIST/namzgr_mes/rn_bot_min  , rn_bot_max  , ln_envl   , &
+          &               nn_strt     , nn_slev     , rn_e_hc   , &
+          &               rn_e_th     , rn_e_bb     , rn_e_ba   , &
+          &               rn_e_al     , ln_pst_mes  , ln_pst_l2g, &
+          &               rn_e3pst_min, rn_e3pst_rat
       !!----------------------------------------------------------------------
       !
@@ -720 +725 @@
       WRITE(ctlmes,*) 'TOT number of levels (jpk) IS DIFFERENT from sum over nn_slev(:)'
       IF ( SUM(nn_slev(:)) /= jpk ) CALL ctl_stop( ctlmes )
+      !
+      IF( .NOT.(ln_loc_zgr) .AND. ln_pst_l2g) ln_pst_l2g = .FALSE.
       !
       ! 3) Reading Bathymetry and envelopes
@@ -810 +817 @@
             WRITE(numout,*) '          ------------------------------------------------------------------'
          ENDDO
+         IF( ln_loc_zgr) THEN
+           WRITE(numout,*) ''
+           WRITE(numout,*) '          LOCALISED MEs '
+           WRITE(numout,*) '          ============= '
+         ENDIF
+         IF( ln_pst_mes .OR. ln_pst_l2g) THEN
+           WRITE(numout,*) ''
+           WRITE(numout,*) '          APPLYING PARTIAL-STEPS to '
+           WRITE(numout,*) '          ========================= '
+           WRITE(numout,*) ''
+           WRITE(numout,*) '            a) MEs area       : ', ln_pst_mes
+           IF( ln_loc_zgr) WRITE(numout,*) '            b) Transition area: ', ln_pst_l2g
+           WRITE(numout,*) ''
+           WRITE(numout,*) '            Partial step thickness is set larger than the minimum of'
+           WRITE(numout,*) ''
+           WRITE(numout,*) '                  rn_e3pst_min = ', rn_e3pst_min
+           WRITE(numout,*) '            and'
+           WRITE(numout,*) '                  rn_e3pst_rat = ', rn_e3pst_rat
+           WRITE(numout,*) ''
+           WRITE(numout,*) '            with 0 < rn_e3pst_rat < 1'
+         ENDIF
       ENDIF