Changeset 13433
- Timestamp:
- 2020-08-24T15:33:43+02:00 (4 years ago)
- 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 62 62 !! Reference : Madec, G. and M. Imbard, 1996: Climate Dynamics, 12, 381-388. 63 63 !!---------------------------------------------------------------------- 64 USE ldftra_oce, ONLY: aht0 65 64 66 INTEGER :: ioptio ! ??? 65 67 LOGICAL :: ll_print = .FALSE. ! Logical flag for printing viscosity coef. … … 67 69 NAMELIST/namdyn_ldf/ ln_dynldf_lap , ln_dynldf_bilap, & 68 70 & 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 70 73 !!---------------------------------------------------------------------- 71 74 72 75 REWIND( numnam ) ! Read Namelist namdyn_ldf : Lateral physics 73 76 READ ( numnam, namdyn_ldf ) 77 78 IF ( rn_ahm_eq_orca < 0.0_wp ) rn_ahm_eq_orca = aht0 74 79 75 80 IF(lwp) THEN ! Parameter print … … 86 91 WRITE(numout,*) ' background viscosity rn_ahmb_0 = ', rn_ahmb_0 87 92 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 88 94 ENDIF 89 95 -
NEMO/branches/NERC/dev_release-3.4_NEMOTAM_consolidated/NEMOGCM/NEMO/OPA_SRC/LDF/ldfdyn_c2d.h90
r3294 r13433 217 217 ! symmetric in the south hemisphere) 218 218 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 221 225 DO jj = 1, jpj 222 226 DO ji = 1, jpi … … 300 304 !! 301 305 !!---------------------------------------------------------------------- 302 USE ldftra_oce, ONLY: aht0303 306 ! 304 307 LOGICAL, INTENT (in) :: ld_print ! If true, output arrays on numout … … 382 385 ! symmetric in the south hemisphere) 383 386 384 zahmeq = aht0387 zahmeq = rn_ahm_eq_orca 385 388 zam20s = ahm0*COS( rad * 20. ) 386 389 -
NEMO/branches/NERC/dev_release-3.4_NEMOTAM_consolidated/NEMOGCM/NEMO/OPA_SRC/LDF/ldfdyn_c3d.h90
r3294 r13433 300 300 301 301 IF( jp_cfg == 2 ) THEN 302 zahmeq = aht0302 zahmeq = rn_ahm_eq_orca 303 303 zahmm = ahm0 304 304 zahm0(:,:) = ahm0 … … 306 306 307 307 IF( jp_cfg == 1 ) THEN 308 zahmeq = aht0 ! reduced to aht0 on equator; set to ahm0 if no tropical reduction is required308 zahmeq = rn_ahm_eq_orca 309 309 zahmm = ahm0 310 310 zahm0(:,:) = ahm0 -
NEMO/branches/NERC/dev_release-3.4_NEMOTAM_consolidated/NEMOGCM/NEMO/OPA_SRC/LDF/ldfdyn_oce.F90
r2715 r13433 22 22 REAL(wp), PUBLIC :: rn_ahmb_0 = 0._wp !: lateral laplacian background eddy viscosity (m2/s) 23 23 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 24 27 REAL(wp), PUBLIC :: ahm0, ahmb0, ahm0_blp !: OLD namelist names 25 28
Note: See TracChangeset
for help on using the changeset viewer.