Changeset 7027 for branches/UKMO/dev_r5518_nn_etau_revision/NEMOGCM/NEMO/OPA_SRC/ZDF/zdftke.F90

Timestamp:

2016-10-13T18:19:43+02:00 (8 years ago)

Author:

hadcv

Message:

Rewrite case 4 of nn_htau for tinkering.

File:

• branches/UKMO/dev_r5518_nn_etau_revision/NEMOGCM/NEMO/OPA_SRC/ZDF/zdftke.F90 (modified) (2 diffs)

branches/UKMO/dev_r5518_nn_etau_revision/NEMOGCM/NEMO/OPA_SRC/ZDF/zdftke.F90

@@ -757 +757 @@
       INTEGER ::   ji, jj, jk   ! dummy loop indices
       INTEGER ::   ios
+      REAL(wp) :: zmin, zmax_h_n, zmax_h_s, zmax_s, zmax_scale
       !!
       NAMELIST/namzdf_tke/ rn_ediff, rn_ediss , rn_ebb , rn_emin  ,   &
@@ -829 +830 @@
          CASE( 3 )                                 ! F(latitude) : 0.5m to 15m poleward of 20 degrees
             htau(:,:) = MAX(  0.5_wp, MIN( 15._wp, 45._wp* ABS( SIN( rpi/180._wp * gphit(:,:) ) ) )   )
-         CASE( 4 )                                 ! F(latitude) : 0.5m to 10m/30m poleward of 13/40 degrees north/south
-            DO jj = 2, jpjm1
-               DO ji = fs_2, fs_jpim1   ! vector opt.
-                  IF( gphit(ji,jj) <= 0._wp ) THEN
-                     htau(ji,jj) = MAX(  0.5_wp, MIN( 30._wp, 45._wp* ABS( SIN( rpi/180._wp * gphit(ji,jj) ) ) )   )
-                  ELSE
-                     htau(ji,jj) = MAX(  0.5_wp, MIN( 10._wp, 45._wp* ABS( SIN( rpi/180._wp * gphit(ji,jj) ) ) )   )
-                  ENDIF
-               END DO
-            END DO
+         CASE( 4 )                                 ! F(latitude) : single sine curve per hemisphere
+            zmax_scale = 1._wp                     ! Amount to scale the maximum value by
+            zmin = 0.5_wp                          ! The minimum value
+            zmax_h_n = 10._wp * zmax_scale         ! The hard maximum value (cutoff imposed on sine curve)
+            zmax_h_s = 30._wp * zmax_scale
+            zmax_s = 45._wp * zmax_scale           ! The soft maximum value (peak of sine curve)
+
+            htau(:,:) =                   MAX( zmin, MIN( zmax_h_n, zmax_s * ABS( SIN( rpi/180._wp * gphit(:,:) ) ) ) )
+            WHERE( gphit < 0._wp ) htau = MAX( zmin, MIN( zmax_h_s, zmax_s * ABS( SIN( rpi/180._wp * gphit      ) ) ) )
          CASE ( 5 )                                ! F(latitude) : 0.5m to 10m poleward of 13 degrees north/south,
             DO jj = 2, jpjm1                       !               10m to 30m between 30/45 degrees south