New URL for NEMO forge!   http://forge.nemo-ocean.eu

Since March 2022 along with NEMO 4.2 release, the code development moved to a self-hosted GitLab.
This present forge is now archived and remained online for history.
Changeset 3047 – NEMO

Changeset 3047


Ignore:
Timestamp:
2011-11-07T12:13:35+01:00 (12 years ago)
Author:
cbricaud
Message:

add changes from dev_r2787_MERCATOR3_tidalpot

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  
    645645 
    646646% ================================================================ 
     647%        Tidal Potential 
     648% ================================================================ 
     649\section   [Tidal Potential (\textit{sbctide})] 
     650                        {Tidal Potential (\mdl{sbctide})} 
     651\label{SBC_tide} 
     652 
     653A 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 
     660Concerning 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 
     669The tide is generated by the forces of gravity ot the Earth-Moon and Earth-Sun sytem; 
     670they are expressed as the gradient of the astronomical potential ($\vec{\nabla}\Pi_{a}$). \\ 
     671 
     672The potential astronomical expressed, for the three types of tidal frequencies 
     673following, by : \\ 
     674Tide 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} 
     678diurnal Tide : 
     679\begin{equation} 
     680\Pi_{a}=gA_{k}(sin 2\phi)cos(\omega_{k}t+\lambda+V_{0k}) 
     681\end{equation} 
     682Semi-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 
     691We make corrections to the astronomical potential. 
     692We obtain :  
     693\begin{equation} 
     694\Pi-g\delta = (1+k-h) \Pi_{A}(\lambda,\phi) 
     695\end{equation} 
     696with $k$ a number of Love estimated to 0.6 which parametrized the astronomical tidal land, 
     697and $h$ a number of Love to 0.3 which parametrized the parametrization due to the astronomical tidal land. 
     698 
     699% ================================================================ 
    647700%        River runoffs 
    648701% ================================================================ 

  • branches/2011/dev_MERCATOR_2011_MERGE/NEMOGCM/CONFIG/GYRE/EXP00/namelist

    r2715 r3047  
    365365   rn_sponge_tra = 2880.   !  coefficient for tracer   sponge layer [m2/s] 
    366366   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' 
    367389/ 
    368390!----------------------------------------------------------------------- 

  • branches/2011/dev_MERCATOR_2011_MERGE/NEMOGCM/CONFIG/ORCA2_LIM/EXP00/namelist

    r3046 r3047  
    365365   rn_sponge_tra = 2880.   !  coefficient for tracer   sponge layer [m2/s] 
    366366   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' 
    367389/ 
    368390!----------------------------------------------------------------------- 

  • branches/2011/dev_MERCATOR_2011_MERGE/NEMOGCM/CONFIG/POMME/EXP00/namelist

    r2650 r3047  
    365365   rn_sponge_tra = 2880.   !  coefficient for tracer   sponge layer [s] 
    366366   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' 
    367389/ 
    368390!----------------------------------------------------------------------- 

  • branches/2011/dev_MERCATOR_2011_MERGE/NEMOGCM/NEMO/OPA_SRC/DYN/dynspg_ts.F90

    r2724 r3047  
    3535   USE bdydyn          ! unstructured open boundaries 
    3636   USE bdytides        ! tidal forcing at unstructured open boundaries. 
     37   USE sbctide 
     38   USE updtide 
    3739   USE lib_mpp         ! distributed memory computing library 
    3840   USE lbclnk          ! ocean lateral boundary conditions (or mpp link) 
    … …  
    371373         IF( lk_obc )   CALL obc_dta_bt ( kt, jn   ) 
    372374         IF( lk_bdy )   CALL bdy_dta_fla( kt, jn+1, icycle ) 
     375         IF ( ln_tide_pot ) CALL upd_tide( kt, jn ) 
    373376 
    374377         !                                                !* after ssh_e 
    … …  
    418421                     zu_spg = -grav * ( sshn_e(ji+1,jj) - sshn_e(ji,jj) ) / e1u(ji,jj) 
    419422                     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) 
    420428                  ENDIF 
    421429                  ! energy conserving formulation for planetary vorticity term 
    … …  
    447455                     zv_spg = -grav * ( sshn_e(ji,jj+1) - sshn_e(ji,jj) ) / e2v(ji,jj) 
    448456                  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 
    449462                  ! enstrophy conserving formulation for planetary vorticity term 
    450463                  zy1 =   z1_8 * ( zwy(ji  ,jj-1) + zwy(ji+1,jj-1) + zwy(ji,jj) + zwy(ji+1,jj  ) ) / e1u(ji,jj) 
    … …  
    472485                     zu_spg = -grav * ( sshn_e(ji+1,jj) - sshn_e(ji,jj) ) / e1u(ji,jj) 
    473486                     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) 
    474492                  ENDIF 
    475493                  ! energy/enstrophy conserving formulation for planetary vorticity term 

  • branches/2011/dev_MERCATOR_2011_MERGE/NEMOGCM/NEMO/OPA_SRC/step.F90

    r3046 r3047  
    9797      IF( lk_dtasal  )   CALL dta_sal( kstp )         ! update 3D salinity data 
    9898                         CALL sbc    ( kstp )         ! Sea Boundary Condition (including sea-ice) 
     99      IF( lk_tide    )   CALL sbc_tide( kstp ) 
    99100      IF( lk_obc     )   CALL obc_dta( kstp )         ! update dynamic and tracer data at open boundaries 
    100101      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  
    2323   USE sbccpl           ! surface boundary condition: coupled formulation (call send at end of step) 
    2424   USE cpl_oasis3, ONLY : lk_cpl 
     25   USE sbctide          ! Tide initialisation 
    2526 
    2627   USE traqsr           ! solar radiation penetration      (tra_qsr routine) 
Note: See TracChangeset for help on using the changeset viewer.