Ignore:
Timestamp:
2011-11-07T10:25:44+01:00 (10 years ago)
Author:
cbricaud
Message:

add changes from dev_r2787_MERCATOR2_tidalharm in dev_MERCATOR_2011_MERGE

Location:
branches/2011/dev_MERCATOR_2011_MERGE
Files:
5 edited
5 copied

Legend:

Unmodified
Added
Removed
  • branches/2011/dev_MERCATOR_2011_MERGE/DOC/TexFiles/Chapters/Chap_DIA.tex

    r2541 r3046  
    685685 
    686686% ------------------------------------------------------------------------------------------------------------- 
     687%       Harmonic analysis of tidal constituents 
     688% ------------------------------------------------------------------------------------------------------------- 
     689\section{Harmonic analysis of tidal constituents (\key{diaharm}) } 
     690\label{DIA_diag_harm} 
     691 
     692A module is available to compute the amplitude and phase for tidal waves.  
     693This diagnostic is actived with \key{diaharm}. 
     694 
     695%------------------------------------------namdia_harm---------------------------------------------------- 
     696\namdisplay{namdia_harm} 
     697%---------------------------------------------------------------------------------------------------------- 
     698 
     699Concerning the on-line Harmonic analysis, some parameters are available in namelist: 
     700 
     701- \texttt{nit000\_han} is the first time step used for harmonic analysis 
     702 
     703- \texttt{nitend\_han} is the last time step used for harmonic analysis 
     704 
     705- \texttt{nstep\_han} is the time step frequency for harmonic analysis 
     706 
     707- \texttt{nb\_ana} is the number of harmonics to analyse 
     708 
     709- \texttt{tname} is an array with names of tidal constituents to analyse 
     710 
     711\texttt{nit000\_han} and \texttt{nitend\_han} must be between \texttt{nit000} and \texttt{nitend} of the simulation. 
     712The restart capability is not implemented. 
     713 
     714The Harmonic analysis solve this equation: 
     715\begin{equation} 
     716h_{i} - A_{0} + \sum^{nb\_ana}_{j=1}[A_{j}cos(\nu_{j}t_{j}-\phi_{j})] = e_{i} 
     717\end{equation} 
     718 
     719With $A_{j}$,$\nu_{j}$,$\phi_{j}$, the amplitude, frequency and phase for each wave and $e_{i}$ the error. 
     720$h_{i}$ is the sea level for the time $t_{i}$ and $A_{0}$ is the mean sea level. \\ 
     721We can rewrite this equation: 
     722\begin{equation} 
     723h_{i} - A_{0} + \sum^{nb\_ana}_{j=1}[C_{j}cos(\nu_{j}t_{j})+S_{j}sin(\nu_{j}t_{j})] = e_{i} 
     724\end{equation} 
     725with $A_{j}=\sqrt{C^{2}_{j}+S^{2}_{j}}$ et $\phi_{j}=arctan(S_{j}/C_{j})$. 
     726 
     727We obtain in output $C_{j}$ and $S_{j}$ for each tidal wave. 
     728 
     729% ------------------------------------------------------------------------------------------------------------- 
    687730%       Other Diagnostics 
    688731% ------------------------------------------------------------------------------------------------------------- 
     
    735778 
    736779 
    737  
    738780% ================================================================ 
    739781% Steric effect in sea surface height 
  • branches/2011/dev_MERCATOR_2011_MERGE/NEMOGCM/CONFIG/ORCA2_LIM/EXP00/iodef.xml

    r2729 r3046  
    129129   <field id="ticemel_cea"  description="Rate of Melt at Upper Surface of Sea Ice (cell average)"      unit="kg/m2/s"  /> 
    130130 
     131 
     132        <!-- tidal composante --> 
     133        <field id="M2x"          description="M2 Elevation harmonic real part "                             unit="m"        /> 
     134        <field id="M2y"          description="M2 Elevation harmonic imaginary part"                         unit="m"        /> 
     135        <field id="S2x"          description="M2 Elevation harmonic real part "                             unit="m"        /> 
     136        <field id="S2y"          description="M2 Elevation harmonic imaginary part"                         unit="m"        /> 
     137        <field id="N2x"          description="M2 Elevation harmonic real part "                             unit="m"        /> 
     138        <field id="N2y"          description="M2 Elevation harmonic imaginary part"                         unit="m"        /> 
     139        <field id="K1x"          description="M2 Elevation harmonic real part "                             unit="m"        /> 
     140        <field id="K1y"          description="M2 Elevation harmonic imaginary part"                         unit="m"        /> 
     141        <field id="O1x"          description="M2 Elevation harmonic real part "                             unit="m"        /> 
     142        <field id="O1y"          description="M2 Elevation harmonic imaginary part"                         unit="m"        /> 
     143        <field id="Q1x"          description="M2 Elevation harmonic real part "                             unit="m"        /> 
     144        <field id="Q1y"          description="M2 Elevation harmonic imaginary part"                         unit="m"        /> 
     145        <field id="M4x"          description="M2 Elevation harmonic real part "                             unit="m"        /> 
     146        <field id="M4y"          description="M2 Elevation harmonic imaginary part"                         unit="m"        /> 
     147        <field id="K2x"          description="M2 Elevation harmonic real part "                             unit="m"        /> 
     148        <field id="K2y"          description="M2 Elevation harmonic imaginary part"                         unit="m"        /> 
     149        <field id="P1x"          description="M2 Elevation harmonic real part "                             unit="m"        /> 
     150        <field id="P1y"          description="M2 Elevation harmonic imaginary part"                         unit="m"        /> 
     151        <field id="Mfx"          description="M2 Elevation harmonic real part "                             unit="m"        /> 
     152        <field id="Mfy"          description="M2 Elevation harmonic imaginary part"                         unit="m"        /> 
     153        <field id="Mmx"          description="M2 Elevation harmonic real part "                             unit="m"        /> 
     154        <field id="Mmy"          description="M2 Elevation harmonic imaginary part"                         unit="m"        /> 
     155 
     156        <field id="M2x_u"          description="M2 current barotrope along i-axis harmonic real part "      unit="m/s"        /> 
     157        <field id="M2y_u"          description="M2 current barotrope along i-axis harmonic imaginary part " unit="m/s"        /> 
     158        <field id="S2x_u"          description="S2 current barotrope along i-axis harmonic real part "      unit="m/s"        /> 
     159        <field id="S2y_u"          description="S2 current barotrope along i-axis harmonic imaginary part " unit="m/s"        /> 
     160        <field id="N2x_u"          description="N2 current barotrope along i-axis harmonic real part "      unit="m/s"        /> 
     161        <field id="N2y_u"          description="N2 current barotrope along i-axis harmonic imaginary part " unit="m/s"        /> 
     162        <field id="K1x_u"          description="K1 current barotrope along i-axis harmonic real part "      unit="m/s"        /> 
     163        <field id="K1y_u"          description="K1 current barotrope along i-axis harmonic imaginary part " unit="m/s"        /> 
     164        <field id="O1x_u"          description="O1 current barotrope along i-axis harmonic real part "      unit="m/s"        /> 
     165        <field id="O1y_u"          description="O1 current barotrope along i-axis harmonic imaginary part " unit="m/s"        /> 
     166        <field id="Q1x_u"          description="Q1 current barotrope along i-axis harmonic real part "      unit="m/s"        /> 
     167        <field id="Q1y_u"          description="Q1 current barotrope along i-axis harmonic imaginary part " unit="m/s"        /> 
     168        <field id="M4x_u"          description="M4 current barotrope along i-axis harmonic real part "      unit="m/s"        /> 
     169        <field id="M4y_u"          description="M4 current barotrope along i-axis harmonic imaginary part " unit="m/s"        /> 
     170        <field id="K2x_u"          description="K2 current barotrope along i-axis harmonic real part "      unit="m/s"        /> 
     171        <field id="K2y_u"          description="K2 current barotrope along i-axis harmonic imaginary part " unit="m/s"        /> 
     172        <field id="P1x_u"          description="P1 current barotrope along i-axis harmonic real part "      unit="m/s"        /> 
     173        <field id="P1y_u"          description="P1 current barotrope along i-axis harmonic imaginary part " unit="m/s"        /> 
     174        <field id="Mfx_u"          description="Mf current barotrope along i-axis harmonic real part "      unit="m/s"        /> 
     175        <field id="Mfy_u"          description="Mf current barotrope along i-axis harmonic imaginary part " unit="m/s"        /> 
     176        <field id="Mmx_u"          description="Mm current barotrope along i-axis harmonic real part "      unit="m/s"        /> 
     177        <field id="Mmy_u"          description="Mm current barotrope along i-axis harmonic imaginary part " unit="m/s"        /> 
     178 
     179        <field id="M2x_v"          description="M2 current barotrope along j-axis harmonic real part "      unit="m/s"        /> 
     180        <field id="M2y_v"          description="M2 current barotrope along j-axis harmonic imaginary part " unit="m/s"        /> 
     181        <field id="S2x_v"          description="S2 current barotrope along j-axis harmonic real part "      unit="m/s"        /> 
     182        <field id="S2y_v"          description="S2 current barotrope along j-axis harmonic imaginary part " unit="m/s"        /> 
     183        <field id="N2x_v"          description="N2 current barotrope along j-axis harmonic real part "      unit="m/s"        /> 
     184        <field id="N2y_v"          description="N2 current barotrope along j-axis harmonic imaginary part " unit="m/s"        /> 
     185        <field id="K1x_v"          description="K1 current barotrope along j-axis harmonic real part "      unit="m/s"        /> 
     186        <field id="K1y_v"          description="K1 current barotrope along j-axis harmonic imaginary part " unit="m/s"        /> 
     187        <field id="O1x_v"          description="O1 current barotrope along j-axis harmonic real part "      unit="m/s"        /> 
     188        <field id="O1y_v"          description="O1 current barotrope along j-axis harmonic imaginary part " unit="m/s"        /> 
     189        <field id="Q1x_v"          description="Q1 current barotrope along j-axis harmonic real part "      unit="m/s"        /> 
     190        <field id="Q1y_v"          description="Q1 current barotrope along j-axis harmonic imaginary part " unit="m/s"        /> 
     191        <field id="M4x_v"          description="M4 current barotrope along j-axis harmonic real part "      unit="m/s"        /> 
     192        <field id="M4y_v"          description="M4 current barotrope along j-axis harmonic imaginary part " unit="m/s"        /> 
     193        <field id="K2x_v"          description="K2 current barotrope along j-axis harmonic real part "      unit="m/s"        /> 
     194        <field id="K2y_v"          description="K2 current barotrope along j-axis harmonic imaginary part " unit="m/s"        /> 
     195        <field id="P1x_v"          description="P1 current barotrope along j-axis harmonic real part "      unit="m/s"        /> 
     196        <field id="P1y_v"          description="P1 current barotrope along j-axis harmonic imaginary part " unit="m/s"        /> 
     197        <field id="Mfx_v"          description="Mf current barotrope along j-axis harmonic real part "      unit="m/s"        /> 
     198        <field id="Mfy_v"          description="Mf current barotrope along j-axis harmonic imaginary part " unit="m/s"        /> 
     199        <field id="Mmx_v"          description="Mm current barotrope along j-axis harmonic real part "      unit="m/s"        /> 
     200        <field id="Mmy_v"          description="Mm current barotrope along j-axis harmonic imaginary part " unit="m/s"        /> 
    131201      </group> 
    132202 
     
    261331       <field ref="sst"    name="sst_1d_ave"                       />     <!-- mean --> 
    262332       <field ref="sst"    name="sst_1d_max"  operation="t_max(X)" />     <!-- max --> 
     333            <field ref="M2x"     name="M2_x_elev"  /> 
     334            <field ref="M2y"     name="M2_y_elev"  /> 
    263335     </file> 
    264336   </group> 
     337 
     338 
     339 
    265340 
    266341      </group> 
  • branches/2011/dev_MERCATOR_2011_MERGE/NEMOGCM/CONFIG/ORCA2_LIM/EXP00/namelist

    r2715 r3046  
    863863    salfixmin = -9999      !  Minimum salinity after applying the increments 
    864864/ 
     865!----------------------------------------------------------------------- 
     866&nam_diaharm   !   Harmonic analysis of tidal constituents ('key_diaharm') 
     867!----------------------------------------------------------------------- 
     868    nit000_han=76           ! First time step used for harmonic analysis 
     869    nitend_han=5475         ! Last time step used for harmonic analysis 
     870    nstep_han=15            ! Time step frequency for harmonic analysis 
     871    nb_ana=2               ! Number of harmonics to analyse 
     872    tname(1)='M2'          ! Name of tidal constituents 
     873    tname(2)='K1' 
     874/ 
  • branches/2011/dev_MERCATOR_2011_MERGE/NEMOGCM/NEMO/OPA_SRC/step.F90

    r2715 r3046  
    165165      IF( ln_diaptr  )   CALL dia_ptr( kstp )         ! Poleward TRansports diagnostics 
    166166      IF( lk_diaar5  )   CALL dia_ar5( kstp )         ! ar5 diag 
     167      IF( lk_diaharm )   CALL dia_harm( kstp )        ! Tidal harmonic analysis 
    167168                         CALL dia_wri( kstp )         ! ocean model: outputs 
    168169 
  • branches/2011/dev_MERCATOR_2011_MERGE/NEMOGCM/NEMO/OPA_SRC/step_oce.F90

    r2528 r3046  
    8585   USE diafwb           ! freshwater budget                (dia_fwb routine) 
    8686   USE diahsb           ! heat, salt and volume budgets    (dia_hsb routine) 
     87   USE diaharm 
    8788   USE flo_oce          ! floats variables 
    8889   USE floats           ! floats computation               (flo_stp routine) 
Note: See TracChangeset for help on using the changeset viewer.