Changeset 2182 for branches/dev_r1821_mixed_ldfdyn/NEMO/OPA_SRC/LDF
- Timestamp:
- 2010-10-07T17:01:49+02:00 (14 years ago)
- Location:
- branches/dev_r1821_mixed_ldfdyn/NEMO/OPA_SRC/LDF
- Files:
-
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/dev_r1821_mixed_ldfdyn/NEMO/OPA_SRC/LDF/ldfdyn.F90
r1839 r2182 67 67 NAMELIST/namdyn_ldf/ ln_dynldf_lap , ln_dynldf_bilap, & 68 68 & ln_dynldf_level, ln_dynldf_hor , ln_dynldf_iso, & 69 & rn_ahm_0 69 & rn_ahm_0_lap , rn_ahmb_0 , rn_ahm_0_blp 70 70 !!---------------------------------------------------------------------- 71 71 … … 83 83 WRITE(numout,*) ' horizontal (geopotential) ln_dynldf_hor = ', ln_dynldf_hor 84 84 WRITE(numout,*) ' iso-neutral ln_dynldf_iso = ', ln_dynldf_iso 85 WRITE(numout,*) ' horizontal laplacian eddy viscosity rn_ahm_0 = ', rn_ahm_085 WRITE(numout,*) ' horizontal laplacian eddy viscosity rn_ahm_0_lap = ', rn_ahm_0_lap 86 86 WRITE(numout,*) ' background viscosity rn_ahmb_0 = ', rn_ahmb_0 87 WRITE(numout,*) ' horizontal bilaplacian eddy viscosity rn_ahm_0 = ', rn_ahm_088 89 ENDIF 90 91 ahm0 = rn_ahm_0 87 WRITE(numout,*) ' horizontal bilaplacian eddy viscosity rn_ahm_0_blp = ', rn_ahm_0_blp 88 89 ENDIF 90 91 ahm0 = rn_ahm_0_lap ! OLD namelist variables defined from DOCTOR namelist variables 92 92 ahmb0 = rn_ahmb_0 93 93 ahm0_blp = rn_ahm_0_blp … … 120 120 IF( ln_dynldf_bilap ) THEN 121 121 IF(lwp) WRITE(numout,*) ' biharmonic momentum diffusion' 122 IF( ahm0_blp == 0.0 ) ahm0_blp = ahm0 ! Old namelist method: bilap specified with ahm0123 122 IF( .NOT. ln_dynldf_lap ) ahm0 = ahm0_blp ! Allow spatially varying coefs, which use ahm0 as input 124 123 IF( ahm0_blp > 0 .AND. .NOT. lk_esopa ) CALL ctl_stop( 'The horizontal viscosity coef. ahm0 must be negative' ) -
branches/dev_r1821_mixed_ldfdyn/NEMO/OPA_SRC/LDF/ldfdyn_oce.F90
r1839 r2182 17 17 LOGICAL , PUBLIC :: ln_dynldf_hor = .TRUE. !: horizontal (geopotential) direction 18 18 LOGICAL , PUBLIC :: ln_dynldf_iso = .FALSE. !: iso-neutral direction 19 REAL(wp), PUBLIC :: rn_ahm_0 19 REAL(wp), PUBLIC :: rn_ahm_0_lap = 40000._wp !: lateral laplacian eddy viscosity (m2/s) 20 20 REAL(wp), PUBLIC :: rn_ahmb_0 = 0._wp !: lateral laplacian background eddy viscosity (m2/s) 21 21 REAL(wp), PUBLIC :: rn_ahm_0_blp = 0._wp !: lateral bilaplacian eddy viscosity (m4/s) 22 22 23 REAL(wp), PUBLIC :: ahm0, ahmb0, ahm0_blp ! OLD namelist names23 REAL(wp), PUBLIC :: ahm0, ahmb0, ahm0_blp ! OLD namelist names 24 24 25 25 #if defined key_dynldf_c3d
Note: See TracChangeset
for help on using the changeset viewer.