Changeset 3047
- Timestamp:
- 2011-11-07T12:13:35+01:00 (13 years ago)
- Location:
- branches/2011/dev_MERCATOR_2011_MERGE
- Files:
-
- 7 edited
- 3 copied
Legend:
- Unmodified
- Added
- Removed
-
branches/2011/dev_MERCATOR_2011_MERGE/DOC/TexFiles/Chapters/Chap_SBC.tex
r2541 r3047 645 645 646 646 % ================================================================ 647 % Tidal Potential 648 % ================================================================ 649 \section [Tidal Potential (\textit{sbctide})] 650 {Tidal Potential (\mdl{sbctide})} 651 \label{SBC_tide} 652 653 A module is available to use the tidal potential forcing and is activated with with \key{tide}. 654 655 656 %------------------------------------------nam_tide---------------------------------------------------- 657 \namdisplay{nam_tide} 658 %------------------------------------------------------------------------------------------------------------- 659 660 Concerning the tidal potential, some parameters are available in namelist: 661 662 - \texttt{ln\_tide\_pot} activate the tidal potential forcing 663 664 - \texttt{nb\_harmo} is the number of constituent used 665 666 - \texttt{clname} is the name of constituent 667 668 669 The tide is generated by the forces of gravity ot the Earth-Moon and Earth-Sun sytem; 670 they are expressed as the gradient of the astronomical potential ($\vec{\nabla}\Pi_{a}$). \\ 671 672 The potential astronomical expressed, for the three types of tidal frequencies 673 following, by : \\ 674 Tide long period : 675 \begin{equation} 676 \Pi_{a}=gA_{k}(\frac{1}{2}-\frac{3}{2}sin^{2}\phi)cos(\omega_{k}t+V_{0k}) 677 \end{equation} 678 diurnal Tide : 679 \begin{equation} 680 \Pi_{a}=gA_{k}(sin 2\phi)cos(\omega_{k}t+\lambda+V_{0k}) 681 \end{equation} 682 Semi-diurnal tide: 683 \begin{equation} 684 \Pi_{a}=gA_{k}(cos^{2}\phi)cos(\omega_{k}t+2\lambda+V_{0k}) 685 \end{equation} 686 687 688 $A_{k}$ is the amplitude of the wave k, $\omega_{k}$ the pulsation of the wave k, $V_{0k}$ the astronomical phase of the wave 689 $k$ to Greenwich. 690 691 We make corrections to the astronomical potential. 692 We obtain : 693 \begin{equation} 694 \Pi-g\delta = (1+k-h) \Pi_{A}(\lambda,\phi) 695 \end{equation} 696 with $k$ a number of Love estimated to 0.6 which parametrized the astronomical tidal land, 697 and $h$ a number of Love to 0.3 which parametrized the parametrization due to the astronomical tidal land. 698 699 % ================================================================ 647 700 % River runoffs 648 701 % ================================================================ -
branches/2011/dev_MERCATOR_2011_MERGE/NEMOGCM/CONFIG/GYRE/EXP00/namelist
r2715 r3047 365 365 rn_sponge_tra = 2880. ! coefficient for tracer sponge layer [m2/s] 366 366 rn_sponge_dyn = 2880. ! coefficient for dynamics sponge layer [m2/s] 367 / 368 !----------------------------------------------------------------------- 369 ! nam_tide tide parameters (#ifdef key_tide) 370 !----------------------------------------------------------------------- 371 ! ln_tide_pot = use tidal potential forcing 372 ! nb_harmo = number of constituents used 373 ! name(1) = 'M2', 'K1', etc name of constituent 374 375 &nam_tide 376 ln_tide_pot = .true. 377 nb_harmo = 11 378 clname(1) = 'M2' 379 clname(2) = 'S2' 380 clname(3) = 'N2' 381 clname(4) = 'K1' 382 clname(5) = 'O1' 383 clname(6) = 'Q1' 384 clname(7) = 'M4' 385 clname(8) = 'K2' 386 clname(9) = 'P1' 387 clname(10) = 'Mf' 388 clname(11) = 'Mm' 367 389 / 368 390 !----------------------------------------------------------------------- -
branches/2011/dev_MERCATOR_2011_MERGE/NEMOGCM/CONFIG/ORCA2_LIM/EXP00/namelist
r3046 r3047 365 365 rn_sponge_tra = 2880. ! coefficient for tracer sponge layer [m2/s] 366 366 rn_sponge_dyn = 2880. ! coefficient for dynamics sponge layer [m2/s] 367 / 368 !----------------------------------------------------------------------- 369 ! nam_tide tide parameters (#ifdef key_tide) 370 !----------------------------------------------------------------------- 371 ! ln_tide_pot = use tidal potential forcing 372 ! nb_harmo = number of constituents used 373 ! name(1) = 'M2', 'K1', etc name of constituent 374 375 &nam_tide 376 ln_tide_pot = .true. 377 nb_harmo = 11 378 clname(1) = 'M2' 379 clname(2) = 'S2' 380 clname(3) = 'N2' 381 clname(4) = 'K1' 382 clname(5) = 'O1' 383 clname(6) = 'Q1' 384 clname(7) = 'M4' 385 clname(8) = 'K2' 386 clname(9) = 'P1' 387 clname(10) = 'Mf' 388 clname(11) = 'Mm' 367 389 / 368 390 !----------------------------------------------------------------------- -
branches/2011/dev_MERCATOR_2011_MERGE/NEMOGCM/CONFIG/POMME/EXP00/namelist
r2650 r3047 365 365 rn_sponge_tra = 2880. ! coefficient for tracer sponge layer [s] 366 366 rn_sponge_dyn = 2880. ! coefficient for dynamics sponge layer [s] 367 / 368 !----------------------------------------------------------------------- 369 ! nam_tide tide parameters (#ifdef key_tide) 370 !----------------------------------------------------------------------- 371 ! ln_tide_pot = use tidal potential forcing 372 ! nb_harmo = number of constituents used 373 ! name(1) = 'M2', 'K1', etc name of constituent 374 375 &nam_tide 376 ln_tide_pot = .true. 377 nb_harmo = 11 378 clname(1) = 'M2' 379 clname(2) = 'S2' 380 clname(3) = 'N2' 381 clname(4) = 'K1' 382 clname(5) = 'O1' 383 clname(6) = 'Q1' 384 clname(7) = 'M4' 385 clname(8) = 'K2' 386 clname(9) = 'P1' 387 clname(10) = 'Mf' 388 clname(11) = 'Mm' 367 389 / 368 390 !----------------------------------------------------------------------- -
branches/2011/dev_MERCATOR_2011_MERGE/NEMOGCM/NEMO/OPA_SRC/DYN/dynspg_ts.F90
r2724 r3047 35 35 USE bdydyn ! unstructured open boundaries 36 36 USE bdytides ! tidal forcing at unstructured open boundaries. 37 USE sbctide 38 USE updtide 37 39 USE lib_mpp ! distributed memory computing library 38 40 USE lbclnk ! ocean lateral boundary conditions (or mpp link) … … 371 373 IF( lk_obc ) CALL obc_dta_bt ( kt, jn ) 372 374 IF( lk_bdy ) CALL bdy_dta_fla( kt, jn+1, icycle ) 375 IF ( ln_tide_pot ) CALL upd_tide( kt, jn ) 373 376 374 377 ! !* after ssh_e … … 418 421 zu_spg = -grav * ( sshn_e(ji+1,jj) - sshn_e(ji,jj) ) / e1u(ji,jj) 419 422 zv_spg = -grav * ( sshn_e(ji,jj+1) - sshn_e(ji,jj) ) / e2v(ji,jj) 423 ENDIF 424 ! add tidal astronomical forcing 425 IF ( ln_tide_pot ) THEN 426 zu_spg = zu_spg + grav * ( pot_astro(ji+1,jj) - pot_astro(ji,jj) ) / e1u(ji,jj) 427 zv_spg = zv_spg + grav * ( pot_astro(ji,jj+1) - pot_astro(ji,jj) ) / e2v(ji,jj) 420 428 ENDIF 421 429 ! energy conserving formulation for planetary vorticity term … … 447 455 zv_spg = -grav * ( sshn_e(ji,jj+1) - sshn_e(ji,jj) ) / e2v(ji,jj) 448 456 ENDIF 457 ! add tidal astronomical forcing 458 IF ( ln_tide_pot ) THEN 459 zu_spg = zu_spg + grav * ( pot_astro(ji+1,jj) - pot_astro(ji,jj) ) / e1u(ji,jj) 460 zv_spg = zv_spg + grav * ( pot_astro(ji,jj+1) - pot_astro(ji,jj) ) / e2v(ji,jj) 461 ENDIF 449 462 ! enstrophy conserving formulation for planetary vorticity term 450 463 zy1 = z1_8 * ( zwy(ji ,jj-1) + zwy(ji+1,jj-1) + zwy(ji,jj) + zwy(ji+1,jj ) ) / e1u(ji,jj) … … 472 485 zu_spg = -grav * ( sshn_e(ji+1,jj) - sshn_e(ji,jj) ) / e1u(ji,jj) 473 486 zv_spg = -grav * ( sshn_e(ji,jj+1) - sshn_e(ji,jj) ) / e2v(ji,jj) 487 ENDIF 488 ! add tidal astronomical forcing 489 IF ( ln_tide_pot ) THEN 490 zu_spg = zu_spg + grav * ( pot_astro(ji+1,jj) - pot_astro(ji,jj) ) / e1u(ji,jj) 491 zv_spg = zv_spg + grav * ( pot_astro(ji,jj+1) - pot_astro(ji,jj) ) / e2v(ji,jj) 474 492 ENDIF 475 493 ! energy/enstrophy conserving formulation for planetary vorticity term -
branches/2011/dev_MERCATOR_2011_MERGE/NEMOGCM/NEMO/OPA_SRC/step.F90
r3046 r3047 97 97 IF( lk_dtasal ) CALL dta_sal( kstp ) ! update 3D salinity data 98 98 CALL sbc ( kstp ) ! Sea Boundary Condition (including sea-ice) 99 IF( lk_tide ) CALL sbc_tide( kstp ) 99 100 IF( lk_obc ) CALL obc_dta( kstp ) ! update dynamic and tracer data at open boundaries 100 101 IF( lk_obc ) CALL obc_rad( kstp ) ! compute phase velocities at open boundaries -
branches/2011/dev_MERCATOR_2011_MERGE/NEMOGCM/NEMO/OPA_SRC/step_oce.F90
r3046 r3047 23 23 USE sbccpl ! surface boundary condition: coupled formulation (call send at end of step) 24 24 USE cpl_oasis3, ONLY : lk_cpl 25 USE sbctide ! Tide initialisation 25 26 26 27 USE traqsr ! solar radiation penetration (tra_qsr routine)
Note: See TracChangeset
for help on using the changeset viewer.