Changeset 13433 for NEMO/branches/NERC

Timestamp:

2020-08-24T15:33:43+02:00 (4 years ago)

Author:

smueller

Message:

Addition of configurability of the equatorial viscosity reduction in ORCA1 and ORCA2 configurations

Location:

NEMO/branches/NERC/dev_release-3.4_NEMOTAM_consolidated/NEMOGCM/NEMO/OPA_SRC/LDF

Files:

• ldfdyn.F90 (modified) (3 diffs)
• ldfdyn_c2d.h90 (modified) (3 diffs)
• ldfdyn_c3d.h90 (modified) (2 diffs)
• ldfdyn_oce.F90 (modified) (1 diff)

NEMO/branches/NERC/dev_release-3.4_NEMOTAM_consolidated/NEMOGCM/NEMO/OPA_SRC/LDF/ldfdyn.F90

@@ -62 +62 @@
       !! Reference :   Madec, G. and M. Imbard, 1996: Climate Dynamics, 12, 381-388.
       !!----------------------------------------------------------------------
+      USE ldftra_oce, ONLY: aht0
+
       INTEGER ::   ioptio         ! ???
       LOGICAL ::   ll_print = .FALSE.    ! Logical flag for printing viscosity coef.
@@ -67 +69 @@
       NAMELIST/namdyn_ldf/ ln_dynldf_lap  , ln_dynldf_bilap,                  &
          &                 ln_dynldf_level, ln_dynldf_hor  , ln_dynldf_iso,   &
-         &                 rn_ahm_0_lap   , rn_ahmb_0      , rn_ahm_0_blp
+         &                 rn_ahm_0_lap   , rn_ahmb_0      , rn_ahm_0_blp,    &
+         &                 rn_ahm_eq_orca
       !!----------------------------------------------------------------------
 
       REWIND( numnam )                  ! Read Namelist namdyn_ldf : Lateral physics
       READ  ( numnam, namdyn_ldf )
+
+      IF ( rn_ahm_eq_orca < 0.0_wp ) rn_ahm_eq_orca = aht0
 
       IF(lwp) THEN                      ! Parameter print
@@ -86 +91 @@
          WRITE(numout,*) '      background viscosity                    rn_ahmb_0       = ', rn_ahmb_0
          WRITE(numout,*) '      horizontal bilaplacian eddy viscosity   rn_ahm_0_blp    = ', rn_ahm_0_blp
+         WRITE(numout,*) '      hor. lap. visc. near eq. in ORCA1/ORCA2 rn_ahm_eq_orca  = ', rn_ahm_eq_orca
       ENDIF
 

NEMO/branches/NERC/dev_release-3.4_NEMOTAM_consolidated/NEMOGCM/NEMO/OPA_SRC/LDF/ldfdyn_c2d.h90

@@ -217 +217 @@
       ! symmetric in the south hemisphere)
 
-      zahmeq = aht0
-      
+      IF (jp_cfg == 2) THEN
+          zahmeq = rn_ahm_eq_orca
+      ELSE
+         zahmeq = aht0
+      END IF
+
       DO jj = 1, jpj
          DO ji = 1, jpi
@@ -300 +304 @@
       !!
       !!----------------------------------------------------------------------
-      USE ldftra_oce, ONLY:   aht0
       !
       LOGICAL, INTENT (in) ::   ld_print   ! If true, output arrays on numout
@@ -382 +385 @@
       ! symmetric in the south hemisphere)
 
-      zahmeq = aht0
+      zahmeq = rn_ahm_eq_orca
       zam20s = ahm0*COS( rad * 20. )
       

NEMO/branches/NERC/dev_release-3.4_NEMOTAM_consolidated/NEMOGCM/NEMO/OPA_SRC/LDF/ldfdyn_c3d.h90

@@ -300 +300 @@
 
       IF( jp_cfg == 2 )   THEN
-         zahmeq     = aht0
+         zahmeq     = rn_ahm_eq_orca
          zahmm      = ahm0
          zahm0(:,:) = ahm0
@@ -306 +306 @@
 
       IF( jp_cfg == 1 )   THEN
-         zahmeq     = aht0  ! reduced to aht0 on equator; set to ahm0 if no tropical reduction is required
+         zahmeq     = rn_ahm_eq_orca
          zahmm      = ahm0
          zahm0(:,:) = ahm0

NEMO/branches/NERC/dev_release-3.4_NEMOTAM_consolidated/NEMOGCM/NEMO/OPA_SRC/LDF/ldfdyn_oce.F90

@@ -22 +22 @@
    REAL(wp), PUBLIC ::   rn_ahmb_0       =     0._wp   !: lateral laplacian background eddy viscosity (m2/s)
    REAL(wp), PUBLIC ::   rn_ahm_0_blp    =     0._wp   !: lateral bilaplacian eddy viscosity (m4/s)
+   REAL(wp), PUBLIC ::   rn_ahm_eq_orca  =   -1.0_wp   !: lateral laplacian eddy viscosity near the Equator
+                                                       !  in ORCA1/ORCA2 configurations (m2/s); negative values
+                                                       !  will be replaced by the value of aht0
    REAL(wp), PUBLIC ::   ahm0, ahmb0, ahm0_blp         !: OLD namelist names