Changeset 3046
- Timestamp:
- 2011-11-07T10:25:44+01:00 (13 years ago)
- 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 685 685 686 686 % ------------------------------------------------------------------------------------------------------------- 687 % Harmonic analysis of tidal constituents 688 % ------------------------------------------------------------------------------------------------------------- 689 \section{Harmonic analysis of tidal constituents (\key{diaharm}) } 690 \label{DIA_diag_harm} 691 692 A module is available to compute the amplitude and phase for tidal waves. 693 This diagnostic is actived with \key{diaharm}. 694 695 %------------------------------------------namdia_harm---------------------------------------------------- 696 \namdisplay{namdia_harm} 697 %---------------------------------------------------------------------------------------------------------- 698 699 Concerning 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. 712 The restart capability is not implemented. 713 714 The Harmonic analysis solve this equation: 715 \begin{equation} 716 h_{i} - A_{0} + \sum^{nb\_ana}_{j=1}[A_{j}cos(\nu_{j}t_{j}-\phi_{j})] = e_{i} 717 \end{equation} 718 719 With $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. \\ 721 We can rewrite this equation: 722 \begin{equation} 723 h_{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} 725 with $A_{j}=\sqrt{C^{2}_{j}+S^{2}_{j}}$ et $\phi_{j}=arctan(S_{j}/C_{j})$. 726 727 We obtain in output $C_{j}$ and $S_{j}$ for each tidal wave. 728 729 % ------------------------------------------------------------------------------------------------------------- 687 730 % Other Diagnostics 688 731 % ------------------------------------------------------------------------------------------------------------- … … 735 778 736 779 737 738 780 % ================================================================ 739 781 % Steric effect in sea surface height -
branches/2011/dev_MERCATOR_2011_MERGE/NEMOGCM/CONFIG/ORCA2_LIM/EXP00/iodef.xml
r2729 r3046 129 129 <field id="ticemel_cea" description="Rate of Melt at Upper Surface of Sea Ice (cell average)" unit="kg/m2/s" /> 130 130 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" /> 131 201 </group> 132 202 … … 261 331 <field ref="sst" name="sst_1d_ave" /> <!-- mean --> 262 332 <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" /> 263 335 </file> 264 336 </group> 337 338 339 265 340 266 341 </group> -
branches/2011/dev_MERCATOR_2011_MERGE/NEMOGCM/CONFIG/ORCA2_LIM/EXP00/namelist
r2715 r3046 863 863 salfixmin = -9999 ! Minimum salinity after applying the increments 864 864 / 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 165 165 IF( ln_diaptr ) CALL dia_ptr( kstp ) ! Poleward TRansports diagnostics 166 166 IF( lk_diaar5 ) CALL dia_ar5( kstp ) ! ar5 diag 167 IF( lk_diaharm ) CALL dia_harm( kstp ) ! Tidal harmonic analysis 167 168 CALL dia_wri( kstp ) ! ocean model: outputs 168 169 -
branches/2011/dev_MERCATOR_2011_MERGE/NEMOGCM/NEMO/OPA_SRC/step_oce.F90
r2528 r3046 85 85 USE diafwb ! freshwater budget (dia_fwb routine) 86 86 USE diahsb ! heat, salt and volume budgets (dia_hsb routine) 87 USE diaharm 87 88 USE flo_oce ! floats variables 88 89 USE floats ! floats computation (flo_stp routine)
Note: See TracChangeset
for help on using the changeset viewer.