Changeset 11664
- Timestamp:
- 2019-10-08T19:32:15+02:00 (5 years ago)
- Location:
- NEMO/branches/2019/dev_r10742_ENHANCE-12_SimonM-Tides
- Files:
-
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/branches/2019/dev_r10742_ENHANCE-12_SimonM-Tides/cfgs/SHARED/namelist_ref
r11663 r11664 565 565 !----------------------------------------------------------------------- 566 566 ln_tide = .false. ! Activate tides 567 nn_tide_var = 0 ! Variant of tidal parameter set and tide-potential computation 567 nn_tide_var = 1 ! Variant of tidal parameter set and tide-potential computation 568 ! ! (1: default; 0: compatibility with previous versions) 568 569 ln_tide_pot = .true. ! use tidal potential forcing 569 570 rn_tide_gamma = 0.7 ! Tidal tilt factor -
NEMO/branches/2019/dev_r10742_ENHANCE-12_SimonM-Tides/src/OCE/TDE/tide_mod.F90
r11663 r11664 111 111 WRITE(numout,*) ' Namelist nam_tide' 112 112 WRITE(numout,*) ' Use tidal components ln_tide = ', ln_tide 113 WRITE(numout,*) ' Variant ( 0: legacy option)nn_tide_var = ', nn_tide_var113 WRITE(numout,*) ' Variant (1: default; 0: legacy option) nn_tide_var = ', nn_tide_var 114 114 WRITE(numout,*) ' Apply astronomical potential ln_tide_pot = ', ln_tide_pot 115 115 WRITE(numout,*) ' Tidal tilt factor rn_tide_gamma = ', rn_tide_gamma … … 280 280 ztmp = tide_harmonics(jk)%v0 + tide_harmonics(jk)%u + tide_components(jk)%nutide * zlon 281 281 ! le potentiel est composé des effets des astres: 282 IF ( tide_components(jk)%nutide == 1 ) THEN ; zcs = zcons * SIN( 2._wp*zlat ) 283 ELSEIF( tide_components(jk)%nutide == 2 ) THEN ; zcs = zcons * COS( zlat )**2 284 ELSE ; zcs = 0._wp 285 ENDIF 282 SELECT CASE( tide_components(jk)%nutide ) 283 CASE( 0 ) ! long-periodic tidal constituents (included unless 284 zcs = zcons * ( 0.5_wp - 1.5_wp * SIN( zlat )**2 ) ! compatibility with original formulation is requested) 285 IF ( nn_tide_var < 1 ) zcs = 0.0_wp 286 CASE( 1 ) ! diurnal tidal constituents 287 zcs = zcons * SIN( 2.0_wp*zlat ) 288 CASE( 2 ) ! semi-diurnal tidal constituents 289 zcs = zcons * COS( zlat )**2 290 CASE DEFAULT ! constituents of higher frequency are not included 291 zcs = 0.0_wp 292 END SELECT 286 293 ztmp1 = ztmp1 + zcs * COS( ztmp ) 287 294 ztmp2 = ztmp2 - zcs * SIN( ztmp )
Note: See TracChangeset
for help on using the changeset viewer.