Changeset 13433


Ignore:
Timestamp:
2020-08-24T15:33:43+02:00 (5 months 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:
4 edited

Legend:

Unmodified
Added
Removed
  • NEMO/branches/NERC/dev_release-3.4_NEMOTAM_consolidated/NEMOGCM/NEMO/OPA_SRC/LDF/ldfdyn.F90

    r3294 r13433  
    6262      !! Reference :   Madec, G. and M. Imbard, 1996: Climate Dynamics, 12, 381-388. 
    6363      !!---------------------------------------------------------------------- 
     64      USE ldftra_oce, ONLY: aht0 
     65 
    6466      INTEGER ::   ioptio         ! ??? 
    6567      LOGICAL ::   ll_print = .FALSE.    ! Logical flag for printing viscosity coef. 
     
    6769      NAMELIST/namdyn_ldf/ ln_dynldf_lap  , ln_dynldf_bilap,                  & 
    6870         &                 ln_dynldf_level, ln_dynldf_hor  , ln_dynldf_iso,   & 
    69          &                 rn_ahm_0_lap   , rn_ahmb_0      , rn_ahm_0_blp 
     71         &                 rn_ahm_0_lap   , rn_ahmb_0      , rn_ahm_0_blp,    & 
     72         &                 rn_ahm_eq_orca 
    7073      !!---------------------------------------------------------------------- 
    7174 
    7275      REWIND( numnam )                  ! Read Namelist namdyn_ldf : Lateral physics 
    7376      READ  ( numnam, namdyn_ldf ) 
     77 
     78      IF ( rn_ahm_eq_orca < 0.0_wp ) rn_ahm_eq_orca = aht0 
    7479 
    7580      IF(lwp) THEN                      ! Parameter print 
     
    8691         WRITE(numout,*) '      background viscosity                    rn_ahmb_0       = ', rn_ahmb_0 
    8792         WRITE(numout,*) '      horizontal bilaplacian eddy viscosity   rn_ahm_0_blp    = ', rn_ahm_0_blp 
     93         WRITE(numout,*) '      hor. lap. visc. near eq. in ORCA1/ORCA2 rn_ahm_eq_orca  = ', rn_ahm_eq_orca 
    8894      ENDIF 
    8995 
  • NEMO/branches/NERC/dev_release-3.4_NEMOTAM_consolidated/NEMOGCM/NEMO/OPA_SRC/LDF/ldfdyn_c2d.h90

    r3294 r13433  
    217217      ! symmetric in the south hemisphere) 
    218218 
    219       zahmeq = aht0 
    220        
     219      IF (jp_cfg == 2) THEN 
     220          zahmeq = rn_ahm_eq_orca 
     221      ELSE 
     222         zahmeq = aht0 
     223      END IF 
     224 
    221225      DO jj = 1, jpj 
    222226         DO ji = 1, jpi 
     
    300304      !! 
    301305      !!---------------------------------------------------------------------- 
    302       USE ldftra_oce, ONLY:   aht0 
    303306      ! 
    304307      LOGICAL, INTENT (in) ::   ld_print   ! If true, output arrays on numout 
     
    382385      ! symmetric in the south hemisphere) 
    383386 
    384       zahmeq = aht0 
     387      zahmeq = rn_ahm_eq_orca 
    385388      zam20s = ahm0*COS( rad * 20. ) 
    386389       
  • NEMO/branches/NERC/dev_release-3.4_NEMOTAM_consolidated/NEMOGCM/NEMO/OPA_SRC/LDF/ldfdyn_c3d.h90

    r3294 r13433  
    300300 
    301301      IF( jp_cfg == 2 )   THEN 
    302          zahmeq     = aht0 
     302         zahmeq     = rn_ahm_eq_orca 
    303303         zahmm      = ahm0 
    304304         zahm0(:,:) = ahm0 
     
    306306 
    307307      IF( jp_cfg == 1 )   THEN 
    308          zahmeq     = aht0  ! reduced to aht0 on equator; set to ahm0 if no tropical reduction is required 
     308         zahmeq     = rn_ahm_eq_orca 
    309309         zahmm      = ahm0 
    310310         zahm0(:,:) = ahm0 
  • NEMO/branches/NERC/dev_release-3.4_NEMOTAM_consolidated/NEMOGCM/NEMO/OPA_SRC/LDF/ldfdyn_oce.F90

    r2715 r13433  
    2222   REAL(wp), PUBLIC ::   rn_ahmb_0       =     0._wp   !: lateral laplacian background eddy viscosity (m2/s) 
    2323   REAL(wp), PUBLIC ::   rn_ahm_0_blp    =     0._wp   !: lateral bilaplacian eddy viscosity (m4/s) 
     24   REAL(wp), PUBLIC ::   rn_ahm_eq_orca  =   -1.0_wp   !: lateral laplacian eddy viscosity near the Equator 
     25                                                       !  in ORCA1/ORCA2 configurations (m2/s); negative values 
     26                                                       !  will be replaced by the value of aht0 
    2427   REAL(wp), PUBLIC ::   ahm0, ahmb0, ahm0_blp         !: OLD namelist names 
    2528 
Note: See TracChangeset for help on using the changeset viewer.