Changeset 12205
- Timestamp:
- 2019-12-12T11:52:50+01:00 (5 years ago)
- Location:
- NEMO/branches/2019/dev_r11943_MERGE_2019
- Files:
-
- 7 deleted
- 44 edited
- 3 copied
Legend:
- Unmodified
- Added
- Removed
-
NEMO/branches/2019/dev_r11943_MERGE_2019/cfgs/AGRIF_DEMO/EXPREF/1_namelist_cfg
r12182 r12205 381 381 !! namdiu Cool skin and warm layer models (default: OFF) 382 382 !! namflo float parameters (default: OFF) 383 !! nam_diaharm Harmonic analysis of tidal constituents (default: OFF)384 383 !! nam_diadct transports through some sections (default: OFF) 385 384 !! nam_diatmb Top Middle Bottom Output (default: OFF) -
NEMO/branches/2019/dev_r11943_MERGE_2019/cfgs/AGRIF_DEMO/EXPREF/2_namelist_cfg
r12182 r12205 330 330 !! namdiu Cool skin and warm layer models (default: OFF) 331 331 !! namflo float parameters (default: OFF) 332 !! nam_diaharm Harmonic analysis of tidal constituents (default: OFF)333 332 !! nam_diadct transports through some sections (default: OFF) 334 333 !! nam_diatmb Top Middle Bottom Output (default: OFF) -
NEMO/branches/2019/dev_r11943_MERGE_2019/cfgs/AGRIF_DEMO/EXPREF/3_namelist_cfg
r12182 r12205 312 312 !! namdiu Cool skin and warm layer models (default: OFF) 313 313 !! namflo float parameters (default: OFF) 314 !! nam_diaharm Harmonic analysis of tidal constituents (default: OFF)315 314 !! nam_diadct transports through some sections (default: OFF) 316 315 !! nam_diatmb Top Middle Bottom Output (default: OFF) -
NEMO/branches/2019/dev_r11943_MERGE_2019/cfgs/AGRIF_DEMO/EXPREF/namelist_cfg
r12182 r12205 381 381 !! namdiu Cool skin and warm layer models (default: OFF) 382 382 !! namflo float parameters (default: OFF) 383 !! nam_diaharm Harmonic analysis of tidal constituents (default: OFF)384 383 !! nam_diadct transports through some sections (default: OFF) 385 384 !! nam_diatmb Top Middle Bottom Output (default: OFF) -
NEMO/branches/2019/dev_r11943_MERGE_2019/cfgs/AMM12/EXPREF/context_nemo.xml
r12193 r12205 22 22 <field_definition src="./field_def_nemo-oce.xml"/> <!-- NEMO ocean dynamics --> 23 23 24 <!-- Override field definitions for multiple-linear-regression analysis (diamlr) --> 25 <field_definition level="1" prec="4" operation="average" enabled=".TRUE." default_value="1.e20" > 26 <field_group id="diamlr_fields"> 27 <!-- Time --> 28 <field id="diamlr_time" grid_ref="diamlr_grid_T_2D" prec="8" /> 29 <!-- Regressors for tidal harmonic analysis --> 30 <field id="diamlr_r001" field_ref="diamlr_time" expr="sin( __TDE_M2_omega__ * diamlr_time )" enabled=".TRUE." comment="harmonic:sin:M2" /> 31 <field id="diamlr_r002" field_ref="diamlr_time" expr="cos( __TDE_M2_omega__ * diamlr_time )" enabled=".TRUE." comment="harmonic:cos:M2" /> 32 <field id="diamlr_r003" field_ref="diamlr_time" expr="sin( __TDE_K1_omega__ * diamlr_time )" enabled=".TRUE." comment="harmonic:sin:K1" /> 33 <field id="diamlr_r004" field_ref="diamlr_time" expr="cos( __TDE_K1_omega__ * diamlr_time )" enabled=".TRUE." comment="harmonic:cos:K1" /> 34 <field id="diamlr_r005" enabled=".FALSE." /> 35 <field id="diamlr_r006" enabled=".FALSE." /> 36 <field id="diamlr_r007" enabled=".FALSE." /> 37 <field id="diamlr_r008" enabled=".FALSE." /> 38 <field id="diamlr_r009" enabled=".FALSE." /> 39 <field id="diamlr_r010" enabled=".FALSE." /> 40 <field id="diamlr_r011" enabled=".FALSE." /> 41 <field id="diamlr_r012" enabled=".FALSE." /> 42 <field id="diamlr_r013" enabled=".FALSE." /> 43 <field id="diamlr_r014" enabled=".FALSE." /> 44 <field id="diamlr_r015" enabled=".FALSE." /> 45 <field id="diamlr_r016" enabled=".FALSE." /> 46 <field id="diamlr_r017" enabled=".FALSE." /> 47 <field id="diamlr_r018" enabled=".FALSE." /> 48 <field id="diamlr_r019" enabled=".FALSE." /> 49 <field id="diamlr_r020" enabled=".FALSE." /> 50 <field id="diamlr_r021" enabled=".FALSE." /> 51 <field id="diamlr_r022" enabled=".FALSE." /> 52 <field id="diamlr_r023" enabled=".FALSE." /> 53 <field id="diamlr_r024" enabled=".FALSE." /> 54 <field id="diamlr_r025" enabled=".FALSE." /> 55 <field id="diamlr_r026" enabled=".FALSE." /> 56 <field id="diamlr_r027" enabled=".FALSE." /> 57 <field id="diamlr_r028" enabled=".FALSE." /> 58 <field id="diamlr_r029" enabled=".FALSE." /> 59 <field id="diamlr_r030" enabled=".FALSE." /> 60 <field id="diamlr_r031" enabled=".FALSE." /> 61 <field id="diamlr_r032" enabled=".FALSE." /> 62 <field id="diamlr_r033" enabled=".FALSE." /> 63 <field id="diamlr_r034" enabled=".FALSE." /> 64 <field id="diamlr_r035" enabled=".FALSE." /> 65 <field id="diamlr_r036" enabled=".FALSE." /> 66 <field id="diamlr_r037" enabled=".FALSE." /> 67 <field id="diamlr_r038" enabled=".FALSE." /> 68 <field id="diamlr_r101" field_ref="diamlr_time" expr="diamlr_time^0.0" enabled=".TRUE." comment="mean" /> 69 <!-- Fields selected for regression analysis --> 70 <field id="diamlr_f001" field_ref="ssh" enabled=".TRUE." /> 71 <field id="diamlr_f002" field_ref="uoce" enabled=".TRUE." /> 72 <field id="diamlr_f003" field_ref="voce" enabled=".TRUE." /> 73 <field id="diamlr_f004" field_ref="toce" enabled=".TRUE." /> 74 </field_group> 75 </field_definition> 76 24 77 <!-- Files definition --> 25 78 <file_definition src="./file_def_nemo-oce.xml"/> <!-- NEMO ocean dynamics --> 26 79 80 <file_definition type="multiple_file" name="@expname@_@freq@_@startdate@_@enddate@" sync_freq="10d" min_digits="4"> 81 82 <!-- Activation of intermediate output for multiple-linear-regression analysis (diamlr) --> 83 <file_group id="diamlr_files" output_freq="1d" output_level="10" enabled=".TRUE."/> 84 85 <!-- Activation and selection of daily detided model diagnostics (diadetide) --> 86 <file_group id="diadetide_files" output_freq="1d" output_level="10" enabled=".TRUE."> 87 <file id="file22" name_suffix="_M2detided_grid_T" description="M2-detided ocean T-grid variables"> 88 <field id="diadetide_ssh" field_ref="diadetide_weight_grid_T_2D" operation="accumulate"> this * ssh </field> 89 </file> 90 <file id="file23" name_suffix="_M2detided_grid_U" description="M2-detided ocean U-grid variables"> 91 <field id="diadetide_uoce" field_ref="diadetide_weight_grid_U_3D" operation="accumulate"> this * uoce </field> 92 </file> 93 <file id="file24" name_suffix="_M2detided_grid_V" description="M2-detided ocean V-grid variables"> 94 <field id="diadetide_voce" field_ref="diadetide_weight_grid_V_3D" operation="accumulate"> this * voce </field> 95 </file> 96 </field_group> 97 98 </file_definition> 99 27 100 <!-- Axis definition --> 28 101 <axis_definition src="./axis_def_nemo.xml"/> 29 102 30 103 <!-- Domain definition --> 31 104 <domain_definition src="./domain_def_nemo.xml"/> -
NEMO/branches/2019/dev_r11943_MERGE_2019/cfgs/AMM12/EXPREF/file_def_nemo-oce.xml
r10641 r12205 49 49 50 50 <file_group id="25h_mean" output_freq="1d" output_level="10" enabled=".TRUE."> 51 52 51 53 52 <file id="file12" name_suffix="_25hourm_grid_T" description="ocean T grid variables, 25h meaned" enabled=".TRUE." > -
NEMO/branches/2019/dev_r11943_MERGE_2019/cfgs/AMM12/EXPREF/namelist_cfg
r11536 r12205 154 154 &nam_tide ! tide parameters (default: OFF) 155 155 !----------------------------------------------------------------------- 156 ln_tide = .true. ! Activate tides 157 ln_tide_pot = .true. ! use tidal potential forcing 158 159 clname(1) = 'Q1' ! name of constituent 160 clname(2) = 'O1' 161 clname(3) = 'P1' 162 clname(4) = 'S1' 163 clname(5) = 'K1' 164 clname(6) = '2N2' 165 clname(7) = 'MU2' 166 clname(8) = 'N2' 167 clname(9) = 'NU2' 168 clname(10) = 'M2' 169 clname(11) = 'L2' 170 clname(12) = 'T2' 171 clname(13) = 'S2' 172 clname(14) = 'K2' 173 clname(15) = 'M4' 156 ln_tide = .true. ! Activate tides 157 ln_tide_pot = .true. ! use tidal potential forcing 158 sn_tide_cnames(1) = 'Q1' ! name of constituent 159 sn_tide_cnames(2) = 'O1' 160 sn_tide_cnames(3) = 'P1' 161 sn_tide_cnames(4) = 'S1' 162 sn_tide_cnames(5) = 'K1' 163 sn_tide_cnames(6) = '2N2' 164 sn_tide_cnames(7) = 'MU2' 165 sn_tide_cnames(8) = 'N2' 166 sn_tide_cnames(9) = 'NU2' 167 sn_tide_cnames(10) = 'M2' 168 sn_tide_cnames(11) = 'L2' 169 sn_tide_cnames(12) = 'T2' 170 sn_tide_cnames(13) = 'S2' 171 sn_tide_cnames(14) = 'K2' 172 sn_tide_cnames(15) = 'M4' 174 173 / 175 174 !----------------------------------------------------------------------- … … 351 350 !! namdiu Cool skin and warm layer models (default: OFF) 352 351 !! namflo float parameters (default: OFF) 353 !! nam_diaharm Harmonic analysis of tidal constituents (default: OFF)354 352 !! nam_diadct transports through some sections (default: OFF) 355 353 !! nam_diatmb Top Middle Bottom Output (default: OFF) -
NEMO/branches/2019/dev_r11943_MERGE_2019/cfgs/C1D_PAPA/EXPREF/namelist_cfg
r12182 r12205 420 420 !! namdiu Cool skin and warm layer models (default: OFF) 421 421 !! namflo float parameters (default: OFF) 422 !! nam_diaharm Harmonic analysis of tidal constituents (default: OFF)423 422 !! nam_diadct transports through some sections (default: OFF) 424 423 !! nam_diatmb Top Middle Bottom Output (default: OFF) … … 448 447 / 449 448 !----------------------------------------------------------------------- 450 &nam_diaharm ! Harmonic analysis of tidal constituents (default: OFF)451 !-----------------------------------------------------------------------452 /453 !-----------------------------------------------------------------------454 449 &nam_diadct ! transports through some sections (default: OFF) 455 450 !----------------------------------------------------------------------- -
NEMO/branches/2019/dev_r11943_MERGE_2019/cfgs/GYRE_BFM/EXPREF/namelist_cfg
r11536 r12205 228 228 !! namdiu Cool skin and warm layer models (default: OFF) 229 229 !! namflo float parameters (default: OFF) 230 !! nam_diaharm Harmonic analysis of tidal constituents (default: OFF)231 230 !! nam_diadct transports through some sections (default: OFF) 232 231 !! nam_diatmb Top Middle Bottom Output (default: OFF) -
NEMO/branches/2019/dev_r11943_MERGE_2019/cfgs/GYRE_PISCES/EXPREF/namelist_cfg
r11536 r12205 222 222 !! namdiu Cool skin and warm layer models (default: OFF) 223 223 !! namflo float parameters (default: OFF) 224 !! nam_diaharm Harmonic analysis of tidal constituents (default: OFF)225 224 !! nam_diadct transports through some sections (default: OFF) 226 225 !! nam_diatmb Top Middle Bottom Output (default: OFF) -
NEMO/branches/2019/dev_r11943_MERGE_2019/cfgs/ORCA2_ICE_PISCES/EXPREF/namelist_cfg
r12182 r12205 391 391 !! namdiu Cool skin and warm layer models (default: OFF) 392 392 !! namflo float parameters (default: OFF) 393 !! nam_diaharm Harmonic analysis of tidal constituents (default: OFF)394 393 !! nam_diadct transports through some sections (default: OFF) 395 394 !! nam_diatmb Top Middle Bottom Output (default: OFF) -
NEMO/branches/2019/dev_r11943_MERGE_2019/cfgs/ORCA2_OFF_PISCES/EXPREF/namelist_cfg
r11536 r12205 377 377 !! namdiu Cool skin and warm layer models (default: OFF) 378 378 !! namflo float parameters (default: OFF) 379 !! nam_diaharm Harmonic analysis of tidal constituents (default: OFF)380 379 !! nam_diadct transports through some sections (default: OFF) 381 380 !! nam_diatmb Top Middle Bottom Output (default: OFF) … … 405 404 / 406 405 !----------------------------------------------------------------------- 407 &nam_diaharm ! Harmonic analysis of tidal constituents (default: OFF)408 !-----------------------------------------------------------------------409 /410 !-----------------------------------------------------------------------411 406 &nam_diadct ! transports through some sections (default: OFF) 412 407 !----------------------------------------------------------------------- -
NEMO/branches/2019/dev_r11943_MERGE_2019/cfgs/ORCA2_OFF_TRC/EXPREF/namelist_cfg
r11536 r12205 375 375 !! namdiu Cool skin and warm layer models (default: OFF) 376 376 !! namflo float parameters (default: OFF) 377 !! nam_diaharm Harmonic analysis of tidal constituents ("key_diaharm")378 377 !! namdct transports through some sections ("key_diadct") 379 378 !! nam_diatmb Top Middle Bottom Output (default: OFF) … … 403 402 / 404 403 !----------------------------------------------------------------------- 405 &nam_diaharm ! Harmonic analysis of tidal constituents ("key_diaharm")406 !-----------------------------------------------------------------------407 /408 !-----------------------------------------------------------------------409 404 &namdct ! transports through some sections ("key_diadct") 410 405 !----------------------------------------------------------------------- -
NEMO/branches/2019/dev_r11943_MERGE_2019/cfgs/ORCA2_SAS_ICE/EXPREF/namelist_cfg
r12182 r12205 196 196 !! namdiu Cool skin and warm layer models (default: OFF) 197 197 !! namflo float parameters (default: OFF) 198 !! nam_diaharm Harmonic analysis of tidal constituents (default: OFF)199 198 !! nam_diadct transports through some sections (default: OFF) 200 199 !! nam_diatmb Top Middle Bottom Output (default: OFF) -
NEMO/branches/2019/dev_r11943_MERGE_2019/cfgs/SHARED/field_def_nemo-oce.xml
r12193 r12205 1 1 <?xml version="1.0"?> 2 2 <!-- $id$ --> 3 4 <field_definition level="1" prec="4" operation="average" enabled=".TRUE." default_value="1.e20" > <!-- time step automaticaly defined --> 5 6 <!-- 7 ===================================================================================================== 8 = Configurable diagnostics = 9 ===================================================================================================== 10 --> 11 12 <field_group id="diamlr_fields"> 13 14 <!-- 15 ===================================================================================================== 16 Configuration of multiple-linear-regression analysis (diamlr) 17 ===================================================================================================== 18 19 This field group configures diamlr for tidal harmonic analysis of field 20 ssh: in addition to a regressor for fitting the mean value (diamlr_r101), 21 it includes the regressors for the analysis of the tidal constituents 22 that are available in the tidal-forcing implementation (see 23 ./src/OCE/SBC/tide.h90). 24 25 --> 26 27 <!-- Time --> 28 <field id="diamlr_time" grid_ref="diamlr_grid_T_2D" prec="8" /> 29 30 <!-- Regressors for tidal harmonic analysis --> 31 <field id="diamlr_r001" field_ref="diamlr_time" expr="sin( __TDE_M2_omega__ * diamlr_time )" enabled=".TRUE." comment="harmonic:sin:M2" /> 32 <field id="diamlr_r002" field_ref="diamlr_time" expr="cos( __TDE_M2_omega__ * diamlr_time )" enabled=".TRUE." comment="harmonic:cos:M2" /> 33 <field id="diamlr_r003" field_ref="diamlr_time" expr="sin( __TDE_N2_omega__ * diamlr_time )" enabled=".TRUE." comment="harmonic:sin:N2" /> 34 <field id="diamlr_r004" field_ref="diamlr_time" expr="cos( __TDE_N2_omega__ * diamlr_time )" enabled=".TRUE." comment="harmonic:cos:N2" /> 35 <field id="diamlr_r005" field_ref="diamlr_time" expr="sin( __TDE_2N2_omega__ * diamlr_time )" enabled=".TRUE." comment="harmonic:sin:2N2" /> 36 <field id="diamlr_r006" field_ref="diamlr_time" expr="cos( __TDE_2N2_omega__ * diamlr_time )" enabled=".TRUE." comment="harmonic:cos:2N2" /> 37 <field id="diamlr_r007" field_ref="diamlr_time" expr="sin( __TDE_S2_omega__ * diamlr_time )" enabled=".TRUE." comment="harmonic:sin:S2" /> 38 <field id="diamlr_r008" field_ref="diamlr_time" expr="cos( __TDE_S2_omega__ * diamlr_time )" enabled=".TRUE." comment="harmonic:cos:S2" /> 39 <field id="diamlr_r009" field_ref="diamlr_time" expr="sin( __TDE_K2_omega__ * diamlr_time )" enabled=".TRUE." comment="harmonic:sin:K2" /> 40 <field id="diamlr_r010" field_ref="diamlr_time" expr="cos( __TDE_K2_omega__ * diamlr_time )" enabled=".TRUE." comment="harmonic:cos:K2" /> 41 <field id="diamlr_r011" field_ref="diamlr_time" expr="sin( __TDE_K1_omega__ * diamlr_time )" enabled=".TRUE." comment="harmonic:sin:K1" /> 42 <field id="diamlr_r012" field_ref="diamlr_time" expr="cos( __TDE_K1_omega__ * diamlr_time )" enabled=".TRUE." comment="harmonic:cos:K1" /> 43 <field id="diamlr_r013" field_ref="diamlr_time" expr="sin( __TDE_O1_omega__ * diamlr_time )" enabled=".TRUE." comment="harmonic:sin:O1" /> 44 <field id="diamlr_r014" field_ref="diamlr_time" expr="cos( __TDE_O1_omega__ * diamlr_time )" enabled=".TRUE." comment="harmonic:cos:O1" /> 45 <field id="diamlr_r015" field_ref="diamlr_time" expr="sin( __TDE_Q1_omega__ * diamlr_time )" enabled=".TRUE." comment="harmonic:sin:Q1" /> 46 <field id="diamlr_r016" field_ref="diamlr_time" expr="cos( __TDE_Q1_omega__ * diamlr_time )" enabled=".TRUE." comment="harmonic:cos:Q1" /> 47 <field id="diamlr_r017" field_ref="diamlr_time" expr="sin( __TDE_P1_omega__ * diamlr_time )" enabled=".TRUE." comment="harmonic:sin:P1" /> 48 <field id="diamlr_r018" field_ref="diamlr_time" expr="cos( __TDE_P1_omega__ * diamlr_time )" enabled=".TRUE." comment="harmonic:cos:P1" /> 49 <field id="diamlr_r019" field_ref="diamlr_time" expr="sin( __TDE_M4_omega__ * diamlr_time )" enabled=".TRUE." comment="harmonic:sin:M4" /> 50 <field id="diamlr_r020" field_ref="diamlr_time" expr="cos( __TDE_M4_omega__ * diamlr_time )" enabled=".TRUE." comment="harmonic:cos:M4" /> 51 <field id="diamlr_r021" field_ref="diamlr_time" expr="sin( __TDE_Mf_omega__ * diamlr_time )" enabled=".TRUE." comment="harmonic:sin:Mf" /> 52 <field id="diamlr_r022" field_ref="diamlr_time" expr="cos( __TDE_Mf_omega__ * diamlr_time )" enabled=".TRUE." comment="harmonic:cos:Mf" /> 53 <field id="diamlr_r023" field_ref="diamlr_time" expr="sin( __TDE_Mm_omega__ * diamlr_time )" enabled=".TRUE." comment="harmonic:sin:Mm" /> 54 <field id="diamlr_r024" field_ref="diamlr_time" expr="cos( __TDE_Mm_omega__ * diamlr_time )" enabled=".TRUE." comment="harmonic:cos:Mm" /> 55 <field id="diamlr_r025" field_ref="diamlr_time" expr="sin( __TDE_Msqm_omega__ * diamlr_time )" enabled=".TRUE." comment="harmonic:sin:Msqm" /> 56 <field id="diamlr_r026" field_ref="diamlr_time" expr="cos( __TDE_Msqm_omega__ * diamlr_time )" enabled=".TRUE." comment="harmonic:cos:Msqm" /> 57 <field id="diamlr_r027" field_ref="diamlr_time" expr="sin( __TDE_Mtm_omega__ * diamlr_time )" enabled=".TRUE." comment="harmonic:sin:Mtm" /> 58 <field id="diamlr_r028" field_ref="diamlr_time" expr="cos( __TDE_Mtm_omega__ * diamlr_time )" enabled=".TRUE." comment="harmonic:cos:Mtm" /> 59 <field id="diamlr_r029" field_ref="diamlr_time" expr="sin( __TDE_S1_omega__ * diamlr_time )" enabled=".TRUE." comment="harmonic:sin:S1" /> 60 <field id="diamlr_r030" field_ref="diamlr_time" expr="cos( __TDE_S1_omega__ * diamlr_time )" enabled=".TRUE." comment="harmonic:cos:S1" /> 61 <field id="diamlr_r031" field_ref="diamlr_time" expr="sin( __TDE_MU2_omega__ * diamlr_time )" enabled=".TRUE." comment="harmonic:sin:MU2" /> 62 <field id="diamlr_r032" field_ref="diamlr_time" expr="cos( __TDE_MU2_omega__ * diamlr_time )" enabled=".TRUE." comment="harmonic:cos:MU2" /> 63 <field id="diamlr_r033" field_ref="diamlr_time" expr="sin( __TDE_NU2_omega__ * diamlr_time )" enabled=".TRUE." comment="harmonic:sin:NU2" /> 64 <field id="diamlr_r034" field_ref="diamlr_time" expr="cos( __TDE_NU2_omega__ * diamlr_time )" enabled=".TRUE." comment="harmonic:cos:NU2" /> 65 <field id="diamlr_r035" field_ref="diamlr_time" expr="sin( __TDE_L2_omega__ * diamlr_time )" enabled=".TRUE." comment="harmonic:sin:L2" /> 66 <field id="diamlr_r036" field_ref="diamlr_time" expr="cos( __TDE_L2_omega__ * diamlr_time )" enabled=".TRUE." comment="harmonic:cos:L2" /> 67 <field id="diamlr_r037" field_ref="diamlr_time" expr="sin( __TDE_T2_omega__ * diamlr_time )" enabled=".TRUE." comment="harmonic:sin:T2" /> 68 <field id="diamlr_r038" field_ref="diamlr_time" expr="cos( __TDE_T2_omega__ * diamlr_time )" enabled=".TRUE." comment="harmonic:cos:T2" /> 69 <field id="diamlr_r101" field_ref="diamlr_time" expr="diamlr_time^0.0" enabled=".TRUE." comment="mean" /> 70 71 <!-- Fields selected for regression analysis --> 72 <field id="diamlr_f001" field_ref="ssh" enabled=".TRUE." /> 73 74 </field_group> 3 75 4 76 <!-- … … 8 80 ============================================================================================================ 9 81 --> 10 <field_definition level="1" prec="4" operation="average" enabled=".TRUE." default_value="1.e20" > <!-- time step automaticaly defined --> 82 83 <field_group id="diadetide_fields"> 84 85 <!-- 86 ===================================================================================================== 87 Weight fields for the computation of daily detided model diagnostics (diadetide) 88 ===================================================================================================== 89 90 --> 91 92 <field id="diadetide_weight" grid_ref="diadetide_grid_T_2D" enabled=".TRUE." /> 93 <field id="diadetide_weight_grid_T_2D" field_ref="diadetide_weight" grid_ref="diadetide_grid_T_2D" enabled=".TRUE." > this </field> 94 <field id="diadetide_weight_grid_U_2D" field_ref="diadetide_weight" grid_ref="diadetide_grid_U_2D" enabled=".TRUE." > this </field> 95 <field id="diadetide_weight_grid_V_2D" field_ref="diadetide_weight" grid_ref="diadetide_grid_V_2D" enabled=".TRUE." > this </field> 96 <field id="diadetide_weight_grid_T_3D" field_ref="diadetide_weight" grid_ref="diadetide_grid_2D_to_grid_T_3D" enabled=".TRUE." > this </field> 97 <field id="diadetide_weight_grid_U_3D" field_ref="diadetide_weight" grid_ref="diadetide_grid_2D_to_grid_U_3D" enabled=".TRUE." > this </field> 98 <field id="diadetide_weight_grid_V_3D" field_ref="diadetide_weight" grid_ref="diadetide_grid_2D_to_grid_V_3D" enabled=".TRUE." > this </field> 99 <field id="diadetide_weight_grid_W_3D" field_ref="diadetide_weight" grid_ref="diadetide_grid_2D_to_grid_W_3D" enabled=".TRUE." > this </field> 100 101 </field_group> 11 102 12 103 <!-- … … 138 229 139 230 <field_group id="Tides_T" grid_ref="grid_T_2D" operation="once" > 140 <!-- tidal composante --> 141 <field id="M2x" long_name="M2 Elevation harmonic real part " unit="m" /> 142 <field id="M2y" long_name="M2 Elevation harmonic imaginary part" unit="m" /> 143 <field id="S2x" long_name="S2 Elevation harmonic real part " unit="m" /> 144 <field id="S2y" long_name="S2 Elevation harmonic imaginary part" unit="m" /> 145 <field id="N2x" long_name="N2 Elevation harmonic real part " unit="m" /> 146 <field id="N2y" long_name="N2 Elevation harmonic imaginary part" unit="m" /> 147 <field id="K1x" long_name="K1 Elevation harmonic real part " unit="m" /> 148 <field id="K1y" long_name="K1 Elevation harmonic imaginary part" unit="m" /> 149 <field id="O1x" long_name="O1 Elevation harmonic real part " unit="m" /> 150 <field id="O1y" long_name="O1 Elevation harmonic imaginary part" unit="m" /> 151 <field id="Q1x" long_name="Q1 Elevation harmonic real part " unit="m" /> 152 <field id="Q1y" long_name="Q1 Elevation harmonic imaginary part" unit="m" /> 153 <field id="M4x" long_name="M4 Elevation harmonic real part " unit="m" /> 154 <field id="M4y" long_name="M4 Elevation harmonic imaginary part" unit="m" /> 155 <field id="K2x" long_name="K2 Elevation harmonic real part " unit="m" /> 156 <field id="K2y" long_name="K2 Elevation harmonic imaginary part" unit="m" /> 157 <field id="P1x" long_name="P1 Elevation harmonic real part " unit="m" /> 158 <field id="P1y" long_name="P1 Elevation harmonic imaginary part" unit="m" /> 159 <field id="Mfx" long_name="Mf Elevation harmonic real part " unit="m" /> 160 <field id="Mfy" long_name="Mf Elevation harmonic imaginary part" unit="m" /> 161 <field id="Mmx" long_name="Mm Elevation harmonic real part " unit="m" /> 162 <field id="Mmy" long_name="Mm Elevation harmonic imaginary part" unit="m" /> 163 </field_group> 164 165 <field_group id="Tides_U" grid_ref="grid_U_2D" operation="once" > 166 <field id="M2x_u" long_name="M2 current barotrope along i-axis harmonic real part " unit="m/s" /> 167 <field id="M2y_u" long_name="M2 current barotrope along i-axis harmonic imaginary part " unit="m/s" /> 168 <field id="S2x_u" long_name="S2 current barotrope along i-axis harmonic real part " unit="m/s" /> 169 <field id="S2y_u" long_name="S2 current barotrope along i-axis harmonic imaginary part " unit="m/s" /> 170 <field id="N2x_u" long_name="N2 current barotrope along i-axis harmonic real part " unit="m/s" /> 171 <field id="N2y_u" long_name="N2 current barotrope along i-axis harmonic imaginary part " unit="m/s" /> 172 <field id="K1x_u" long_name="K1 current barotrope along i-axis harmonic real part " unit="m/s" /> 173 <field id="K1y_u" long_name="K1 current barotrope along i-axis harmonic imaginary part " unit="m/s" /> 174 <field id="O1x_u" long_name="O1 current barotrope along i-axis harmonic real part " unit="m/s" /> 175 <field id="O1y_u" long_name="O1 current barotrope along i-axis harmonic imaginary part " unit="m/s" /> 176 <field id="Q1x_u" long_name="Q1 current barotrope along i-axis harmonic real part " unit="m/s" /> 177 <field id="Q1y_u" long_name="Q1 current barotrope along i-axis harmonic imaginary part " unit="m/s" /> 178 <field id="M4x_u" long_name="M4 current barotrope along i-axis harmonic real part " unit="m/s" /> 179 <field id="M4y_u" long_name="M4 current barotrope along i-axis harmonic imaginary part " unit="m/s" /> 180 <field id="K2x_u" long_name="K2 current barotrope along i-axis harmonic real part " unit="m/s" /> 181 <field id="K2y_u" long_name="K2 current barotrope along i-axis harmonic imaginary part " unit="m/s" /> 182 <field id="P1x_u" long_name="P1 current barotrope along i-axis harmonic real part " unit="m/s" /> 183 <field id="P1y_u" long_name="P1 current barotrope along i-axis harmonic imaginary part " unit="m/s" /> 184 <field id="Mfx_u" long_name="Mf current barotrope along i-axis harmonic real part " unit="m/s" /> 185 <field id="Mfy_u" long_name="Mf current barotrope along i-axis harmonic imaginary part " unit="m/s" /> 186 <field id="Mmx_u" long_name="Mm current barotrope along i-axis harmonic real part " unit="m/s" /> 187 <field id="Mmy_u" long_name="Mm current barotrope along i-axis harmonic imaginary part " unit="m/s" /> 188 </field_group> 189 190 <field_group id="Tides_V" grid_ref="grid_V_2D" operation="once" > 191 <field id="M2x_v" long_name="M2 current barotrope along j-axis harmonic real part " unit="m/s" /> 192 <field id="M2y_v" long_name="M2 current barotrope along j-axis harmonic imaginary part " unit="m/s" /> 193 <field id="S2x_v" long_name="S2 current barotrope along j-axis harmonic real part " unit="m/s" /> 194 <field id="S2y_v" long_name="S2 current barotrope along j-axis harmonic imaginary part " unit="m/s" /> 195 <field id="N2x_v" long_name="N2 current barotrope along j-axis harmonic real part " unit="m/s" /> 196 <field id="N2y_v" long_name="N2 current barotrope along j-axis harmonic imaginary part " unit="m/s" /> 197 <field id="K1x_v" long_name="K1 current barotrope along j-axis harmonic real part " unit="m/s" /> 198 <field id="K1y_v" long_name="K1 current barotrope along j-axis harmonic imaginary part " unit="m/s" /> 199 <field id="O1x_v" long_name="O1 current barotrope along j-axis harmonic real part " unit="m/s" /> 200 <field id="O1y_v" long_name="O1 current barotrope along j-axis harmonic imaginary part " unit="m/s" /> 201 <field id="Q1x_v" long_name="Q1 current barotrope along j-axis harmonic real part " unit="m/s" /> 202 <field id="Q1y_v" long_name="Q1 current barotrope along j-axis harmonic imaginary part " unit="m/s" /> 203 <field id="M4x_v" long_name="M4 current barotrope along j-axis harmonic real part " unit="m/s" /> 204 <field id="M4y_v" long_name="M4 current barotrope along j-axis harmonic imaginary part " unit="m/s" /> 205 <field id="K2x_v" long_name="K2 current barotrope along j-axis harmonic real part " unit="m/s" /> 206 <field id="K2y_v" long_name="K2 current barotrope along j-axis harmonic imaginary part " unit="m/s" /> 207 <field id="P1x_v" long_name="P1 current barotrope along j-axis harmonic real part " unit="m/s" /> 208 <field id="P1y_v" long_name="P1 current barotrope along j-axis harmonic imaginary part " unit="m/s" /> 209 <field id="Mfx_v" long_name="Mf current barotrope along j-axis harmonic real part " unit="m/s" /> 210 <field id="Mfy_v" long_name="Mf current barotrope along j-axis harmonic imaginary part " unit="m/s" /> 211 <field id="Mmx_v" long_name="Mm current barotrope along j-axis harmonic real part " unit="m/s" /> 212 <field id="Mmy_v" long_name="Mm current barotrope along j-axis harmonic imaginary part " unit="m/s" /> 231 <!-- Tidal potential --> 232 <field id="tide_pot" long_name="Total tidal potential" unit="m" /> 233 <field id="tide_pot_M2" long_name="M2 tidal potential" unit="m" /> 234 <field id="tide_pot_N2" long_name="N2 tidal potential" unit="m" /> 235 <field id="tide_pot_2N2" long_name="2N2 tidal potential" unit="m" /> 236 <field id="tide_pot_S2" long_name="S2 tidal potential" unit="m" /> 237 <field id="tide_pot_K2" long_name="K2 tidal potential" unit="m" /> 238 <field id="tide_pot_K1" long_name="K1 tidal potential" unit="m" /> 239 <field id="tide_pot_O1" long_name="O1 tidal potential" unit="m" /> 240 <field id="tide_pot_Q1" long_name="Q1 tidal potential" unit="m" /> 241 <field id="tide_pot_P1" long_name="P1 tidal potential" unit="m" /> 242 <field id="tide_pot_M4" long_name="M4 tidal potential" unit="m" /> 243 <field id="tide_pot_Mf" long_name="Mf tidal potential" unit="m" /> 244 <field id="tide_pot_Mm" long_name="Mm tidal potential" unit="m" /> 245 <field id="tide_pot_Msqm" long_name="Msqm tidal potential" unit="m" /> 246 <field id="tide_pot_Mtm" long_name="Mtm tidal potential" unit="m" /> 247 <field id="tide_pot_S1" long_name="S1 tidal potential" unit="m" /> 248 <field id="tide_pot_MU2" long_name="MU2 tidal potential" unit="m" /> 249 <field id="tide_pot_NU2" long_name="NU2 tidal potential" unit="m" /> 250 <field id="tide_pot_L2" long_name="L2 tidal potential" unit="m" /> 251 <field id="tide_pot_T2" long_name="T2 tidal potential" unit="m" /> 213 252 </field_group> 214 253 -
NEMO/branches/2019/dev_r11943_MERGE_2019/cfgs/SHARED/grid_def_nemo.xml
r12193 r12205 1 <?xml version="1.0"?> 1 2 <!-- 2 3 ============================================================================================================ … … 4 5 ============================================================================================================ 5 6 --> 6 7 7 8 <grid_definition> 8 9 … … 59 60 <!-- --> 60 61 <grid id="grid_T_nfloat" > 61 <domain id="grid_T" /> 62 <axis id="nfloat" /> 63 </grid> 64 65 <!-- ABL grid definition --> 66 <grid id="grid_TA_2D"> 67 <domain domain_ref="grid_T" /> 68 </grid> 69 <grid id="grid_TA_3D"> 70 <domain domain_ref="grid_T" /> 71 <axis id="ght_abl" /> 72 </grid> 73 <grid id="grid_WA_3D"> 74 <domain domain_ref="grid_T" /> 75 <axis id="ghw_abl" /> 76 </grid> 77 <!-- --> 62 <domain domain_ref="grid_T" /> 63 <axis axis_ref="nfloat" /> 64 </grid> 65 <!-- --> 78 66 <grid id="grid_EqT" > 79 67 <domain domain_ref="EqT" /> … … 237 225 </grid> 238 226 239 <grid id="grid_V_4strait_ice_hsum"> 240 <scalar > 241 <reduce_domain operation="sum" local="true"/> 242 <reduce_scalar operation="sum" /> 243 </scalar> 227 <grid id="grid_V_4strait_ice_hsum"> 228 <scalar > 229 <reduce_domain operation="sum" local="true"/> 230 <reduce_scalar operation="sum" /> 231 </scalar> 232 <axis axis_ref="section_ice"/> 233 </grid> 234 235 <grid id="grid_4strait_ice"> 244 236 <axis axis_ref="section_ice"/> 245 </grid> 246 247 <grid id="grid_4strait_ice"> 248 <axis axis_ref="section_ice"/> 249 </grid> 237 </grid> 250 238 251 239 <!-- scalars --> … … 253 241 <scalar/> 254 242 </grid> 255 256 </grid_definition> 257 243 244 <!-- ABL grid definition --> 245 <grid id="grid_TA_2D"> 246 <domain domain_ref="grid_T" /> 247 </grid> 248 <grid id="grid_TA_3D"> 249 <domain domain_ref="grid_T" /> 250 <axis id="ght_abl" /> 251 </grid> 252 <grid id="grid_WA_3D"> 253 <domain domain_ref="grid_T" /> 254 <axis id="ghw_abl" /> 255 </grid> 256 <!-- --> 257 258 <!-- grid definitions for multiple-linear-regression analysis (diamlr) --> 259 <grid id="diamlr_grid_scalar" > 260 <scalar /> 261 <scalar /> 262 </grid> 263 <grid id="diamlr_grid_T_2D" > 264 <domain domain_ref="grid_T" /> 265 <scalar /> 266 </grid> 267 <grid id="diamlr_grid_U_2D" > 268 <domain domain_ref="grid_U" /> 269 <scalar /> 270 </grid> 271 <grid id="diamlr_grid_V_2D" > 272 <domain domain_ref="grid_V" /> 273 <scalar /> 274 </grid> 275 <grid id="diamlr_grid_W_2D" > 276 <domain domain_ref="grid_W" /> 277 <scalar /> 278 </grid> 279 <grid id="diamlr_grid_2D_to_grid_T_3D" > 280 <domain domain_ref="grid_T" /> 281 <axis axis_ref="deptht"> 282 <duplicate_scalar /> 283 </axis> 284 </grid> 285 <grid id="diamlr_grid_2D_to_grid_U_3D" > 286 <domain domain_ref="grid_U" /> 287 <axis axis_ref="depthu"> 288 <duplicate_scalar /> 289 </axis> 290 </grid> 291 <grid id="diamlr_grid_2D_to_grid_V_3D" > 292 <domain domain_ref="grid_V" /> 293 <axis axis_ref="depthv"> 294 <duplicate_scalar /> 295 </axis> 296 </grid> 297 <grid id="diamlr_grid_2D_to_grid_W_3D" > 298 <domain domain_ref="grid_W" /> 299 <axis axis_ref="depthw"> 300 <duplicate_scalar /> 301 </axis> 302 </grid> 303 <grid id="diamlr_grid_2D_to_scalar" > 304 <scalar> 305 <reduce_domain operation="average" /> 306 </scalar> 307 <scalar /> 308 </grid> 309 <!-- grid definitions for the computation of daily detided model diagnostics (diadetide) --> 310 <grid id="diadetide_grid_T_2D" > 311 <domain domain_ref="grid_T" /> 312 <scalar /> 313 </grid> 314 <grid id="diadetide_grid_U_2D" > 315 <domain domain_ref="grid_U" /> 316 <scalar /> 317 </grid> 318 <grid id="diadetide_grid_V_2D" > 319 <domain domain_ref="grid_V" /> 320 <scalar /> 321 </grid> 322 <grid id="diadetide_grid_2D_to_grid_T_3D" > 323 <domain domain_ref="grid_T" /> 324 <axis axis_ref="deptht"> 325 <duplicate_scalar /> 326 </axis> 327 </grid> 328 <grid id="diadetide_grid_2D_to_grid_U_3D" > 329 <domain domain_ref="grid_U" /> 330 <axis axis_ref="depthu"> 331 <duplicate_scalar /> 332 </axis> 333 </grid> 334 <grid id="diadetide_grid_2D_to_grid_V_3D" > 335 <domain domain_ref="grid_V" /> 336 <axis axis_ref="depthv"> 337 <duplicate_scalar /> 338 </axis> 339 </grid> 340 <grid id="diadetide_grid_2D_to_grid_W_3D" > 341 <domain domain_ref="grid_W" /> 342 <axis axis_ref="depthw"> 343 <duplicate_scalar /> 344 </axis> 345 </grid> 346 347 </grid_definition> -
NEMO/branches/2019/dev_r11943_MERGE_2019/cfgs/SHARED/namelist_ref
r12193 r12205 642 642 !----------------------------------------------------------------------- 643 643 ln_tide = .false. ! Activate tides 644 nn_tide_var = 1 ! Variant of tidal parameter set and tide-potential computation 645 ! ! (1: default; 0: compatibility with previous versions) 646 ln_tide_dia = .false. ! Enable tidal diagnostic output 644 647 ln_tide_pot = .false. ! use tidal potential forcing 648 rn_tide_gamma = 0.7 ! Tidal tilt factor 645 649 ln_scal_load = .false. ! Use scalar approximation for 646 650 rn_scal_load = 0.094 ! load potential … … 649 653 ! 650 654 ln_tide_ramp = .false. ! Use linear ramp for tides at startup 651 r dttideramp = 0.! ramp duration in days652 clname(1)= 'DUMMY' ! name of constituent - all tidal components must be set in namelist_cfg655 rn_tide_ramp_dt = 0. ! ramp duration in days 656 sn_tide_cnames(1) = 'DUMMY' ! name of constituent - all tidal components must be set in namelist_cfg 653 657 / 654 658 !----------------------------------------------------------------------- … … 726 730 filtide = 'bdydta/amm12_bdytide_' ! file name root of tidal forcing files 727 731 ln_bdytide_2ddta = .false. ! 728 ln_bdytide_conj = .false. !729 732 / 730 733 … … 1192 1195 !! namdiu Cool skin and warm layer models (default: OFF) 1193 1196 !! namflo float parameters (default: OFF) 1194 !! nam_diaharm Harmonic analysis of tidal constituents (default: OFF)1195 1197 !! nam_diadct transports through some sections (default: OFF) 1196 1198 !! nam_diatmb Top Middle Bottom Output (default: OFF) … … 1245 1247 ln_ariane = .true. ! Input with Ariane tool convention(T) 1246 1248 ln_flo_ascii= .true. ! Output with Ariane tool netcdf convention(F) or ascii file (T) 1247 /1248 !-----------------------------------------------------------------------1249 &nam_diaharm ! Harmonic analysis of tidal constituents (default: OFF)1250 !-----------------------------------------------------------------------1251 ln_diaharm = .false. ! Choose tidal harmonic output or not1252 nit000_han = 1 ! First time step used for harmonic analysis1253 nitend_han = 75 ! Last time step used for harmonic analysis1254 nstep_han = 15 ! Time step frequency for harmonic analysis1255 tname(1) = 'M2' ! Name of tidal constituents1256 tname(2) = 'K1' ! ---1257 1249 / 1258 1250 !----------------------------------------------------------------------- -
NEMO/branches/2019/dev_r11943_MERGE_2019/cfgs/SPITZ12/EXPREF/namelist_cfg
r12182 r12205 138 138 &nam_tide ! tide parameters (default: OFF) 139 139 !----------------------------------------------------------------------- 140 ln_tide = .true. ! Activate tides141 ln_tide_pot = .false.! use tidal potential forcing142 clname(1) = 'M2'! name of constituent143 clname(2)= 'S2'144 clname(3)= 'N2'145 clname(4)= 'K2'146 clname(5)= 'K1'147 clname(6)= 'O1'148 clname(7)= 'P1'149 clname(8)= 'Q1'140 ln_tide = .true. ! Activate tides 141 ln_tide_pot = .false. ! use tidal potential forcing 142 sn_tide_cnames(1) = 'M2' ! name of constituent 143 sn_tide_cnames(2) = 'S2' 144 sn_tide_cnames(3) = 'N2' 145 sn_tide_cnames(4) = 'K2' 146 sn_tide_cnames(5) = 'K1' 147 sn_tide_cnames(6) = 'O1' 148 sn_tide_cnames(7) = 'P1' 149 sn_tide_cnames(8) = 'Q1' 150 150 / 151 151 !----------------------------------------------------------------------- … … 348 348 !! namdiu Cool skin and warm layer models (default: OFF) 349 349 !! namflo float parameters (default: OFF) 350 !! nam_diaharm Harmonic analysis of tidal constituents (default: OFF)351 350 !! nam_diadct transports through some sections (default: OFF) 352 351 !! nam_diatmb Top Middle Bottom Output (default: OFF) -
NEMO/branches/2019/dev_r11943_MERGE_2019/doc/latex/NEMO/subfiles/chap_DIA.tex
r11693 r12205 1580 1580 1581 1581 %% ================================================================================================= 1582 \section[Harmonic analysis of tidal constituents (\texttt{\textbf{key\_diaharm}})]{Harmonic analysis of tidal constituents (\protect\key{diaharm})}1583 \label{sec:DIA_diag_harm}1584 1585 \begin{listing}1586 \nlst{nam_diaharm}1587 \caption{\forcode{&nam_diaharm}}1588 \label{lst:nam_diaharm}1589 \end{listing}1590 1591 A module is available to compute the amplitude and phase of tidal waves.1592 This on-line Harmonic analysis is actived with \key{diaharm}.1593 1594 Some parameters are available in namelist \nam{_diaharm}{\_diaharm}:1595 1596 - \np{nit000_han}{nit000\_han} is the first time step used for harmonic analysis1597 1598 - \np{nitend_han}{nitend\_han} is the last time step used for harmonic analysis1599 1600 - \np{nstep_han}{nstep\_han} is the time step frequency for harmonic analysis1601 1602 % - \np{nb_ana}{nb\_ana} is the number of harmonics to analyse1603 1604 - \np{tname}{tname} is an array with names of tidal constituents to analyse1605 1606 \np{nit000_han}{nit000\_han} and \np{nitend_han}{nitend\_han} must be between \np{nit000}{nit000} and \np{nitend}{nitend} of the simulation.1607 The restart capability is not implemented.1608 1609 The Harmonic analysis solve the following equation:1610 1611 \[1612 h_{i} - A_{0} + \sum^{nb\_ana}_{j=1}[A_{j}cos(\nu_{j}t_{j}-\phi_{j})] = e_{i}1613 \]1614 1615 With $A_{j}$, $\nu_{j}$, $\phi_{j}$, the amplitude, frequency and phase for each wave and $e_{i}$ the error.1616 $h_{i}$ is the sea level for the time $t_{i}$ and $A_{0}$ is the mean sea level. \\1617 We can rewrite this equation:1618 1619 \[1620 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}1621 \]1622 1623 with $A_{j}=\sqrt{C^{2}_{j}+S^{2}_{j}}$ and $\phi_{j}=arctan(S_{j}/C_{j})$.1624 1625 We obtain in output $C_{j}$ and $S_{j}$ for each tidal wave.1626 1627 %% =================================================================================================1628 1582 \section[Transports across sections (\texttt{\textbf{key\_diadct}})]{Transports across sections (\protect\key{diadct})} 1629 1583 \label{sec:DIA_diag_dct} -
NEMO/branches/2019/dev_r11943_MERGE_2019/src/OCE/BDY/bdydta.F90
r12145 r12205 23 23 USE phycst ! physical constants 24 24 USE sbcapr ! atmospheric pressure forcing 25 USE sbctide ! Tidal forcing or not25 USE tide_mod, ONLY: ln_tide ! tidal forcing 26 26 USE bdy_oce ! ocean open boundary conditions 27 27 USE bdytides ! tidal forcing at boundaries … … 75 75 CONTAINS 76 76 77 SUBROUTINE bdy_dta( kt, Kmm, k it, kt_offset )77 SUBROUTINE bdy_dta( kt, Kmm, kt_offset ) 78 78 !!---------------------------------------------------------------------- 79 79 !! *** SUBROUTINE bdy_dta *** … … 86 86 INTEGER, INTENT(in) :: kt ! ocean time-step index 87 87 INTEGER, INTENT(in) :: Kmm ! ocean time level index 88 INTEGER, INTENT(in), OPTIONAL :: kit ! subcycle time-step index (for timesplitting option) 89 INTEGER, INTENT(in), OPTIONAL :: kt_offset ! time offset in units of timesteps. NB. if kit 88 INTEGER, INTENT(in), OPTIONAL :: kt_offset ! time offset in units of timesteps 90 89 ! ! is present then units = subcycle timesteps. 91 90 ! ! kt_offset = 0 => get data at "now" time level … … 106 105 ! Initialise data arrays once for all from initial conditions where required 107 106 !--------------------------------------------------------------------------- 108 IF( kt == nit000 .AND. .NOT.PRESENT(kit)) THEN107 IF( kt == nit000 ) THEN 109 108 110 109 ! Calculate depth-mean currents … … 217 216 ! read/update all bdy data 218 217 ! ------------------------ 219 CALL fld_read( kt, 1, bf_alias, k it = kit, kt_offset = kt_offset )218 CALL fld_read( kt, 1, bf_alias, kt_offset = kt_offset ) 220 219 ! apply some corrections in some specific cases... 221 220 ! -------------------------------------------------- … … 275 274 END DO 276 275 ENDIF ! ltotvel 277 278 ! update tidal harmonic forcing279 IF( PRESENT(kit) .AND. nn_dyn2d_dta(jbdy) .GE. 2 ) THEN280 CALL bdytide_update( kt = kt, idx = idx_bdy(jbdy), dta = dta_alias, td = tides(jbdy), &281 & kit = kit, kt_offset = kt_offset )282 ENDIF283 276 284 277 ! atm surface pressure : add inverted barometer effect to ssh if it was read -
NEMO/branches/2019/dev_r11943_MERGE_2019/src/OCE/BDY/bdyini.F90
r11960 r12205 22 22 USE bdydta ! open boundary cond. setting (bdy_dta_init routine) 23 23 USE bdytides ! open boundary cond. setting (bdytide_init routine) 24 USE sbctide ! Tidal forcing or not24 USE tide_mod, ONLY: ln_tide ! tidal forcing 25 25 USE phycst , ONLY: rday 26 26 ! -
NEMO/branches/2019/dev_r11943_MERGE_2019/src/OCE/BDY/bdytides.F90
r11960 r12205 18 18 USE phycst ! physical constants 19 19 USE bdy_oce ! ocean open boundary conditions 20 USE tide ini!20 USE tide_mod ! 21 21 USE daymod ! calendar 22 22 ! … … 30 30 31 31 PUBLIC bdytide_init ! routine called in bdy_init 32 PUBLIC bdytide_update ! routine called in bdy_dta33 32 PUBLIC bdy_dta_tides ! routine called in dyn_spg_ts 34 33 … … 45 44 TYPE(OBC_DATA) , PUBLIC, DIMENSION(jp_bdy) :: dta_bdy_s !: bdy external data (slow component) 46 45 46 INTEGER :: kt_tide 47 47 48 !!---------------------------------------------------------------------- 48 49 !! NEMO/OCE 4.0 , NEMO Consortium (2018) … … 64 65 CHARACTER(len=80) :: filtide !: Filename root for tidal input files 65 66 LOGICAL :: ln_bdytide_2ddta !: If true, read 2d harmonic data 66 LOGICAL :: ln_bdytide_conj !: If true, assume complex conjugate tidal data67 67 !! 68 68 INTEGER :: ib_bdy, itide, ib !: dummy loop indices … … 79 79 TYPE(TIDES_DATA), POINTER :: td !: local short cut 80 80 !! 81 NAMELIST/nambdy_tide/filtide, ln_bdytide_2ddta , ln_bdytide_conj81 NAMELIST/nambdy_tide/filtide, ln_bdytide_2ddta 82 82 !!---------------------------------------------------------------------- 83 83 ! … … 119 119 IF(lwp) WRITE(numout,*) ' Namelist nambdy_tide : tidal harmonic forcing at open boundaries' 120 120 IF(lwp) WRITE(numout,*) ' read tidal data in 2d files: ', ln_bdytide_2ddta 121 IF(lwp) WRITE(numout,*) ' assume complex conjugate : ', ln_bdytide_conj122 121 IF(lwp) WRITE(numout,*) ' Number of tidal components to read: ', nb_harmo 123 122 IF(lwp) THEN 124 123 WRITE(numout,*) ' Tidal components: ' 125 124 DO itide = 1, nb_harmo 126 WRITE(numout,*) ' ', Wave(ntide(itide))%cname_tide125 WRITE(numout,*) ' ', tide_harmonics(itide)%cname_tide 127 126 END DO 128 127 ENDIF … … 165 164 igrd = 1 ! Everything is at T-points here 166 165 DO itide = 1, nb_harmo 167 CALL iom_get( inum, jpdom_autoglo, TRIM( Wave(ntide(itide))%cname_tide)//'_z1', ztr(:,:) )168 CALL iom_get( inum, jpdom_autoglo, TRIM( Wave(ntide(itide))%cname_tide)//'_z2', zti(:,:) )166 CALL iom_get( inum, jpdom_autoglo, TRIM(tide_harmonics(itide)%cname_tide)//'_z1', ztr(:,:) ) 167 CALL iom_get( inum, jpdom_autoglo, TRIM(tide_harmonics(itide)%cname_tide)//'_z2', zti(:,:) ) 169 168 DO ib = 1, ilen0(igrd) 170 169 ii = idx_bdy(ib_bdy)%nbi(ib,igrd) … … 182 181 igrd = 2 ! Everything is at U-points here 183 182 DO itide = 1, nb_harmo 184 CALL iom_get ( inum, jpdom_autoglo, TRIM( Wave(ntide(itide))%cname_tide)//'_u1', ztr(:,:) )185 CALL iom_get ( inum, jpdom_autoglo, TRIM( Wave(ntide(itide))%cname_tide)//'_u2', zti(:,:) )183 CALL iom_get ( inum, jpdom_autoglo, TRIM(tide_harmonics(itide)%cname_tide)//'_u1', ztr(:,:) ) 184 CALL iom_get ( inum, jpdom_autoglo, TRIM(tide_harmonics(itide)%cname_tide)//'_u2', zti(:,:) ) 186 185 DO ib = 1, ilen0(igrd) 187 186 ii = idx_bdy(ib_bdy)%nbi(ib,igrd) … … 199 198 igrd = 3 ! Everything is at V-points here 200 199 DO itide = 1, nb_harmo 201 CALL iom_get ( inum, jpdom_autoglo, TRIM( Wave(ntide(itide))%cname_tide)//'_v1', ztr(:,:) )202 CALL iom_get ( inum, jpdom_autoglo, TRIM( Wave(ntide(itide))%cname_tide)//'_v2', zti(:,:) )200 CALL iom_get ( inum, jpdom_autoglo, TRIM(tide_harmonics(itide)%cname_tide)//'_v1', ztr(:,:) ) 201 CALL iom_get ( inum, jpdom_autoglo, TRIM(tide_harmonics(itide)%cname_tide)//'_v2', zti(:,:) ) 203 202 DO ib = 1, ilen0(igrd) 204 203 ii = idx_bdy(ib_bdy)%nbi(ib,igrd) … … 224 223 DO itide = 1, nb_harmo 225 224 ! ! SSH fields 226 clfile = TRIM(filtide)//TRIM( Wave(ntide(itide))%cname_tide)//'_grid_T.nc'225 clfile = TRIM(filtide)//TRIM(tide_harmonics(itide)%cname_tide)//'_grid_T.nc' 227 226 CALL iom_open( clfile, inum ) 228 227 CALL fld_map( inum, 'z1' , dta_read(1:ilen0(1),1:1,1:1) , 1, idx_bdy(ib_bdy)%nbmap(:,1) ) … … 232 231 CALL iom_close( inum ) 233 232 ! ! U fields 234 clfile = TRIM(filtide)//TRIM( Wave(ntide(itide))%cname_tide)//'_grid_U.nc'233 clfile = TRIM(filtide)//TRIM(tide_harmonics(itide)%cname_tide)//'_grid_U.nc' 235 234 CALL iom_open( clfile, inum ) 236 235 CALL fld_map( inum, 'u1' , dta_read(1:ilen0(2),1:1,1:1) , 1, idx_bdy(ib_bdy)%nbmap(:,2) ) … … 240 239 CALL iom_close( inum ) 241 240 ! ! V fields 242 clfile = TRIM(filtide)//TRIM( Wave(ntide(itide))%cname_tide)//'_grid_V.nc'241 clfile = TRIM(filtide)//TRIM(tide_harmonics(itide)%cname_tide)//'_grid_V.nc' 243 242 CALL iom_open( clfile, inum ) 244 243 CALL fld_map( inum, 'v1' , dta_read(1:ilen0(3),1:1,1:1) , 1, idx_bdy(ib_bdy)%nbmap(:,3) ) … … 254 253 ENDIF ! ln_bdytide_2ddta=.true. 255 254 ! 256 IF( ln_bdytide_conj ) THEN ! assume complex conjugate in data files257 td%ssh0(:,:,2) = - td%ssh0(:,:,2)258 td%u0 (:,:,2) = - td%u0 (:,:,2)259 td%v0 (:,:,2) = - td%v0 (:,:,2)260 ENDIF261 !262 255 ! Allocate slow varying data in the case of time splitting: 263 256 ! Do it anyway because at this stage knowledge of free surface scheme is unknown … … 274 267 ! 275 268 END SUBROUTINE bdytide_init 276 277 278 SUBROUTINE bdytide_update( kt, idx, dta, td, kit, kt_offset )279 !!----------------------------------------------------------------------280 !! *** SUBROUTINE bdytide_update ***281 !!282 !! ** Purpose : - Add tidal forcing to ssh, u2d and v2d OBC data arrays.283 !!284 !!----------------------------------------------------------------------285 INTEGER , INTENT(in ) :: kt ! Main timestep counter286 TYPE(OBC_INDEX) , INTENT(in ) :: idx ! OBC indices287 TYPE(OBC_DATA) , INTENT(inout) :: dta ! OBC external data288 TYPE(TIDES_DATA) , INTENT(inout) :: td ! tidal harmonics data289 INTEGER, OPTIONAL, INTENT(in ) :: kit ! Barotropic timestep counter (for timesplitting option)290 INTEGER, OPTIONAL, INTENT(in ) :: kt_offset ! time offset in units of timesteps. NB. if kit291 ! ! is present then units = subcycle timesteps.292 ! ! kt_offset = 0 => get data at "now" time level293 ! ! kt_offset = -1 => get data at "before" time level294 ! ! kt_offset = +1 => get data at "after" time level295 ! ! etc.296 !297 INTEGER :: itide, igrd, ib ! dummy loop indices298 INTEGER :: time_add ! time offset in units of timesteps299 INTEGER, DIMENSION(3) :: ilen0 ! length of boundary data (from OBC arrays)300 REAL(wp) :: z_arg, z_sarg, zflag, zramp ! local scalars301 REAL(wp), DIMENSION(jpmax_harmo) :: z_sist, z_cost302 !!----------------------------------------------------------------------303 !304 ilen0(1) = SIZE(td%ssh(:,1,1))305 ilen0(2) = SIZE(td%u(:,1,1))306 ilen0(3) = SIZE(td%v(:,1,1))307 308 zflag=1309 IF ( PRESENT(kit) ) THEN310 IF ( kit /= 1 ) zflag=0311 ENDIF312 313 IF ( (nsec_day == NINT(0.5_wp * rdt) .OR. kt==nit000) .AND. zflag==1 ) THEN314 !315 kt_tide = kt - (nsec_day - 0.5_wp * rdt)/rdt316 !317 IF(lwp) THEN318 WRITE(numout,*)319 WRITE(numout,*) 'bdytide_update : (re)Initialization of the tidal bdy forcing at kt=',kt320 WRITE(numout,*) '~~~~~~~~~~~~~~ '321 ENDIF322 !323 CALL tide_init_elevation ( idx, td )324 CALL tide_init_velocities( idx, td )325 !326 ENDIF327 328 time_add = 0329 IF( PRESENT(kt_offset) ) THEN330 time_add = kt_offset331 ENDIF332 333 IF( PRESENT(kit) ) THEN334 z_arg = ((kt-kt_tide) * rdt + (kit+0.5_wp*(time_add-1)) * rdt / REAL(nn_baro,wp) )335 ELSE336 z_arg = ((kt-kt_tide)+time_add) * rdt337 ENDIF338 339 ! Linear ramp on tidal component at open boundaries340 zramp = 1._wp341 IF (ln_tide_ramp) zramp = MIN(MAX( (z_arg + (kt_tide-nit000)*rdt)/(rdttideramp*rday),0._wp),1._wp)342 343 DO itide = 1, nb_harmo344 z_sarg = z_arg * omega_tide(itide)345 z_cost(itide) = COS( z_sarg )346 z_sist(itide) = SIN( z_sarg )347 END DO348 349 DO itide = 1, nb_harmo350 igrd=1 ! SSH on tracer grid351 DO ib = 1, ilen0(igrd)352 dta%ssh(ib) = dta%ssh(ib) + zramp*(td%ssh(ib,itide,1)*z_cost(itide) + td%ssh(ib,itide,2)*z_sist(itide))353 END DO354 igrd=2 ! U grid355 DO ib = 1, ilen0(igrd)356 dta%u2d(ib) = dta%u2d(ib) + zramp*(td%u (ib,itide,1)*z_cost(itide) + td%u (ib,itide,2)*z_sist(itide))357 END DO358 igrd=3 ! V grid359 DO ib = 1, ilen0(igrd)360 dta%v2d(ib) = dta%v2d(ib) + zramp*(td%v (ib,itide,1)*z_cost(itide) + td%v (ib,itide,2)*z_sist(itide))361 END DO362 END DO363 !364 END SUBROUTINE bdytide_update365 269 366 270 … … 406 310 ! Linear ramp on tidal component at open boundaries 407 311 zramp = 1. 408 IF (ln_tide_ramp) zramp = MIN(MAX( (z_arg - nit000*rdt)/(r dttideramp*rday),0.),1.)312 IF (ln_tide_ramp) zramp = MIN(MAX( (z_arg - nit000*rdt)/(rn_tide_ramp_dt*rday),0.),1.) 409 313 410 314 DO ib_bdy = 1,nb_bdy … … 447 351 DO itide = 1, nb_harmo 448 352 ! 449 z_sarg = (z_arg + zoff) * omega_tide(itide)353 z_sarg = (z_arg + zoff) * tide_harmonics(itide)%omega 450 354 z_cost = zramp * COS( z_sarg ) 451 355 z_sist = zramp * SIN( z_sarg ) … … 505 409 END DO 506 410 DO ib = 1 , ilen0(igrd) 507 mod_tide(ib)=mod_tide(ib)* ftide(itide)508 phi_tide(ib)=phi_tide(ib)+ v0tide(itide)+utide(itide)411 mod_tide(ib)=mod_tide(ib)*tide_harmonics(itide)%f 412 phi_tide(ib)=phi_tide(ib)+tide_harmonics(itide)%v0+tide_harmonics(itide)%u 509 413 ENDDO 510 414 DO ib = 1 , ilen0(igrd) … … 544 448 END DO 545 449 DO ib = 1, ilen0(igrd) 546 mod_tide(ib)=mod_tide(ib)* ftide(itide)547 phi_tide(ib)=phi_tide(ib)+ v0tide(itide)+utide(itide)450 mod_tide(ib)=mod_tide(ib)*tide_harmonics(itide)%f 451 phi_tide(ib)=phi_tide(ib)+tide_harmonics(itide)%v0 + tide_harmonics(itide)%u 548 452 ENDDO 549 453 DO ib = 1, ilen0(igrd) … … 565 469 END DO 566 470 DO ib = 1, ilen0(igrd) 567 mod_tide(ib)=mod_tide(ib)* ftide(itide)568 phi_tide(ib)=phi_tide(ib)+ v0tide(itide)+utide(itide)471 mod_tide(ib)=mod_tide(ib)*tide_harmonics(itide)%f 472 phi_tide(ib)=phi_tide(ib)+tide_harmonics(itide)%v0 + tide_harmonics(itide)%u 569 473 ENDDO 570 474 DO ib = 1, ilen0(igrd) -
NEMO/branches/2019/dev_r11943_MERGE_2019/src/OCE/DYN/dynspg.F90
r11960 r12205 21 21 USE dynspg_exp ! surface pressure gradient (dyn_spg_exp routine) 22 22 USE dynspg_ts ! surface pressure gradient (dyn_spg_ts routine) 23 USE sbctide ! 24 USE updtide ! 23 USE tide_mod ! 25 24 USE trd_oce ! trends: ocean variables 26 25 USE trddyn ! trend manager: dynamics … … 43 42 INTEGER, PARAMETER :: np_EXP = 0 ! explicit time stepping 44 43 INTEGER, PARAMETER :: np_NO =-1 ! no surface pressure gradient, no scheme 44 ! 45 REAL(wp) :: zt0step ! Time of day at the beginning of the time step 45 46 46 47 !! * Substitutions … … 116 117 IF( .NOT.ln_dynspg_ts .AND. ( ln_tide_pot .AND. ln_tide ) ) THEN ! N.B. added directly at sub-time-step in ts-case 117 118 ! 118 CALL upd_tide( kt ) ! update tide potential 119 ! Update tide potential at the beginning of current time step 120 zt0step = REAL(nsec_day, wp)-0.5_wp*rdt 121 CALL upd_tide(zt0step, Kmm) 119 122 ! 120 123 DO jj = 2, jpjm1 ! add tide potential forcing -
NEMO/branches/2019/dev_r11943_MERGE_2019/src/OCE/DYN/dynspg_ts.F90
r12150 r12205 44 44 USE bdytides ! open boundary condition data 45 45 USE bdydyn2d ! open boundary conditions on barotropic variables 46 USE sbctide ! tides 47 USE updtide ! tide potential 46 USE tide_mod ! 48 47 USE sbcwave ! surface wave 49 48 USE diatmb ! Top,middle,bottom output … … 173 172 REAL(wp), ALLOCATABLE, DIMENSION(:,:) :: ztwdmask, zuwdmask, zvwdmask ! ROMS wetting and drying masks at t,u,v points 174 173 REAL(wp), ALLOCATABLE, DIMENSION(:,:) :: zuwdav2, zvwdav2 ! averages over the sub-steps of zuwdmask and zvwdmask 174 REAL(wp) :: zt0substep ! Time of day at the beginning of the time substep 175 175 !!---------------------------------------------------------------------- 176 176 ! … … 444 444 ! 445 445 IF( ln_bdy .AND. ln_tide ) CALL bdy_dta_tides( kt, kit=jn, kt_offset= noffset+1 ) 446 IF( ln_tide_pot .AND. ln_tide ) CALL upd_tide ( kt, kit=jn, kt_offset= noffset ) 446 ! Update tide potential at the beginning of current time substep 447 IF( ln_tide_pot .AND. ln_tide ) THEN 448 zt0substep = REAL(nsec_day, wp) - 0.5_wp*rdt + (jn + noffset - 1) * rdt / REAL(nn_baro, wp) 449 CALL upd_tide(zt0substep, Kmm) 450 END IF 447 451 ! 448 452 ! !== extrapolation at mid-step ==! (jn+1/2) -
NEMO/branches/2019/dev_r11943_MERGE_2019/src/OCE/IOM/iom.F90
r12193 r12205 55 55 LOGICAL, PUBLIC, PARAMETER :: lk_iomput = .FALSE. !: iom_put flag 56 56 #endif 57 PUBLIC iom_init, iom_ swap, iom_open, iom_close, iom_setkt, iom_varid, iom_get, iom_get_var57 PUBLIC iom_init, iom_init_closedef, iom_swap, iom_open, iom_close, iom_setkt, iom_varid, iom_get, iom_get_var 58 58 PUBLIC iom_chkatt, iom_getatt, iom_putatt, iom_getszuld, iom_rstput, iom_delay_rst, iom_put 59 PUBLIC iom_use, iom_context_finalize, iom_ miss_val59 PUBLIC iom_use, iom_context_finalize, iom_update_file_name, iom_miss_val 60 60 61 61 PRIVATE iom_rp0d, iom_rp1d, iom_rp2d, iom_rp3d … … 64 64 #if defined key_iomput 65 65 PRIVATE iom_set_domain_attr, iom_set_axis_attr, iom_set_field_attr, iom_set_file_attr, iom_get_file_attr, iom_set_grid_attr 66 PRIVATE set_grid, set_grid_bounds, set_scalar, set_xmlatt, set_mooring, iom_ update_file_name, iom_sdate66 PRIVATE set_grid, set_grid_bounds, set_scalar, set_xmlatt, set_mooring, iom_sdate 67 67 PRIVATE iom_set_rst_context, iom_set_rstw_active, iom_set_rstr_active 68 68 # endif … … 92 92 CONTAINS 93 93 94 SUBROUTINE iom_init( cdname, fname, ld_tmppatch )94 SUBROUTINE iom_init( cdname, fname, ld_tmppatch, ld_closedef ) 95 95 !!---------------------------------------------------------------------- 96 96 !! *** ROUTINE *** … … 102 102 CHARACTER(len=*), OPTIONAL, INTENT(in) :: fname 103 103 LOGICAL , OPTIONAL, INTENT(in) :: ld_tmppatch 104 LOGICAL , OPTIONAL, INTENT(in) :: ld_closedef 104 105 #if defined key_iomput 105 106 ! … … 116 117 INTEGER :: nldi_save, nlei_save !: and close boundaries in output files 117 118 INTEGER :: nldj_save, nlej_save !: 119 LOGICAL :: ll_closedef = .TRUE. 118 120 !!---------------------------------------------------------------------- 119 121 ! … … 130 132 IF( njmpp + jpj - 1 == jpjglo ) nlej = jpj 131 133 ENDIF 134 IF ( PRESENT(ld_closedef) ) ll_closedef = ld_closedef 132 135 ! 133 136 ALLOCATE( zt_bnds(2,jpk), zw_bnds(2,jpk) ) … … 268 271 ENDIF 269 272 ! 270 ! end file definition273 ! set time step length 271 274 dtime%second = rdt 272 275 CALL xios_set_timestep( dtime ) 273 CALL xios_close_context_definition() 274 CALL xios_update_calendar( 0 ) 276 ! 277 ! conditional closure of context definition 278 IF ( ll_closedef ) CALL iom_init_closedef 275 279 ! 276 280 DEALLOCATE( zt_bnds, zw_bnds ) … … 283 287 ! 284 288 END SUBROUTINE iom_init 289 290 SUBROUTINE iom_init_closedef 291 !!---------------------------------------------------------------------- 292 !! *** SUBROUTINE iom_init_closedef *** 293 !!---------------------------------------------------------------------- 294 !! 295 !! ** Purpose : Closure of context definition 296 !! 297 !!---------------------------------------------------------------------- 298 299 CALL xios_close_context_definition() 300 CALL xios_update_calendar( 0 ) 301 302 END SUBROUTINE iom_init_closedef 285 303 286 304 SUBROUTINE iom_set_rstw_var_active(field) -
NEMO/branches/2019/dev_r11943_MERGE_2019/src/OCE/TDE/tide_mod.F90
r12065 r12205 17 17 !!---------------------------------------------------------------------- 18 18 19 USE oce, ONLY : ssh n! sea-surface height19 USE oce, ONLY : ssh ! sea-surface height 20 20 USE par_oce ! ocean parameters 21 21 USE phycst, ONLY : rpi ! pi … … 132 132 sn_tide_cnames(:) = '' 133 133 ! Read Namelist nam_tide 134 REWIND( numnam_ref ) ! Namelist nam_tide in reference namelist : Tides135 134 READ ( numnam_ref, nam_tide, IOSTAT = ios, ERR = 901) 136 135 901 IF( ios /= 0 ) CALL ctl_nam ( ios , 'nam_tide in reference namelist' ) 137 136 ! 138 REWIND( numnam_cfg ) ! Namelist nam_tide in configuration namelist : Tides139 137 READ ( numnam_cfg, nam_tide, IOSTAT = ios, ERR = 902 ) 140 138 902 IF( ios > 0 ) CALL ctl_nam ( ios , 'nam_tide in configuration namelist' ) … … 715 713 716 714 717 SUBROUTINE upd_tide(pdelta )715 SUBROUTINE upd_tide(pdelta, Kmm) 718 716 !!---------------------------------------------------------------------- 719 717 !! *** ROUTINE upd_tide *** … … 726 724 !!---------------------------------------------------------------------- 727 725 REAL, INTENT(in) :: pdelta ! Temporal offset in seconds 726 INTEGER, INTENT(IN) :: Kmm ! Time level index 728 727 INTEGER :: jk ! Dummy loop index 729 728 REAL(wp) :: zt, zramp ! Local scalars … … 755 754 ! 756 755 IF( ln_tide_dia ) THEN ! Output total tidal potential (incl. load potential) 757 IF ( iom_use( "tide_pot" ) ) CALL iom_put( "tide_pot", pot_astro(:,:) + rn_scal_load * ssh n(:,:) )756 IF ( iom_use( "tide_pot" ) ) CALL iom_put( "tide_pot", pot_astro(:,:) + rn_scal_load * ssh(:,:,Kmm) ) 758 757 END IF 759 758 ! -
NEMO/branches/2019/dev_r11943_MERGE_2019/src/OCE/nemogcm.F90
r12193 r12205 46 46 USE closea ! treatment of closed seas (for ln_closea) 47 47 USE usrdef_nam ! user defined configuration 48 USE tide ini ! tidal components initialization (tide_iniroutine)48 USE tide_mod, ONLY : tide_init ! tidal components initialization (tide_init routine) 49 49 USE bdy_oce, ONLY : ln_bdy 50 50 USE bdyini ! open boundary cond. setting (bdy_init routine) … … 59 59 USE diaobs ! Observation diagnostics (dia_obs_init routine) 60 60 USE diacfl ! CFL diagnostics (dia_cfl_init routine) 61 USE dia harm ! tidal harmonics diagnostics (dia_harm_init routine)61 USE diamlr ! IOM context management for multiple-linear-regression analysis 62 62 USE step ! NEMO time-stepping (stp routine) 63 63 USE isfstp ! ice shelf (isf_stp_init routine) … … 75 75 USE diatmb ! Top,middle,bottom output 76 76 USE dia25h ! 25h mean output 77 USE diadetide ! Weights computation for daily detiding of model diagnostics 77 78 USE sbc_oce , ONLY : lk_oasis 78 79 USE wet_dry ! Wetting and drying setting (wad_init routine) … … 492 493 CALL dia_tmb_init ! TMB outputs 493 494 CALL dia_25h_init( Nbb ) ! 25h mean outputs 494 CALL dia_ harm_init ! tidal harmonics outputs495 CALL dia_detide_init ! Weights computation for daily detiding of model diagnostics 495 496 IF( ln_diaobs ) CALL dia_obs( nit000-1, Nnn ) ! Observation operator for restart 497 CALL dia_mlr_init ! Initialisation of IOM context management for multiple-linear-regression analysis 496 498 497 499 ! ! Assimilation increments -
NEMO/branches/2019/dev_r11943_MERGE_2019/src/OCE/step.F90
r12193 r12205 107 107 108 108 IF( kstp == nit000 ) THEN ! initialize IOM context (must be done after nemo_init for AGRIF+XIOS+OASIS) 109 CALL iom_init( cxios_context ) ! for model grid (including passible AGRIF zoom) 109 CALL iom_init( cxios_context, ld_closedef=.FALSE. ) ! for model grid (including passible AGRIF zoom) 110 IF( lk_diamlr ) CALL dia_mlr_iom_init ! with additional setup for multiple-linear-regression analysis 111 CALL iom_init_closedef 110 112 IF( ln_crs ) CALL iom_init( TRIM(cxios_context)//"_crs" ) ! for coarse grid 111 113 ENDIF … … 117 119 ! Update external forcing (tides, open boundaries, ice shelf interaction and surface boundary condition (including sea-ice) 118 120 !<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<< 119 IF( ln_tide ) CALL sbc_tide( kstp )! update tide potential121 IF( ln_tide ) CALL tide_update( kstp ) ! update tide potential 120 122 IF( ln_apr_dyn ) CALL sbc_apr ( kstp ) ! atmospheric pressure (NB: call before bdy_dta which needs ssh_ib) 121 123 IF( ln_bdy ) CALL bdy_dta ( kstp, Nnn, kt_offset = +1 ) ! update dynamic & tracer data at open boundaries … … 211 213 ! diagnostics and outputs 212 214 !<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<< 213 IF( ln_floats ) CALL flo_stp ( kstp, Nbb, Nnn )! drifting Floats214 IF( ln_diacfl ) CALL dia_cfl ( kstp, Nnn )! Courant number diagnostics215 CALL dia_hth ( kstp, Nnn )! Thermocline depth (20 degres isotherm depth)216 IF( ln_diadct ) CALL dia_dct ( kstp, Nnn )! Transports217 CALL dia_ar5 ( kstp, Nnn )! ar5 diag218 CALL dia_ptr ( kstp, Nnn )! Poleward adv/ldf TRansports diagnostics219 IF( ln_diaharm ) CALL dia_harm( kstp, Nnn ) ! Tidal harmonic analysis220 CALL dia_wri ( kstp, Nnn ) ! ocean model: outputs221 !222 IF( l n_crs ) CALL crs_fld ( kstp, Nnn ) ! ocean model: online field coarsening & output215 IF( ln_floats ) CALL flo_stp ( kstp, Nbb, Nnn ) ! drifting Floats 216 IF( ln_diacfl ) CALL dia_cfl ( kstp, Nnn ) ! Courant number diagnostics 217 CALL dia_hth ( kstp, Nnn ) ! Thermocline depth (20 degres isotherm depth) 218 IF( ln_diadct ) CALL dia_dct ( kstp, Nnn ) ! Transports 219 CALL dia_ar5 ( kstp, Nnn ) ! ar5 diag 220 CALL dia_ptr ( kstp, Nnn ) ! Poleward adv/ldf TRansports diagnostics 221 CALL dia_wri ( kstp, Nnn ) ! ocean model: outputs 222 IF( ln_crs ) CALL crs_fld ( kstp, Nnn ) ! ocean model: online field coarsening & output 223 IF( lk_diadetide ) CALL dia_detide( kstp ) ! Weights computation for daily detiding of model diagnostics 224 IF( lk_diamlr ) CALL dia_mlr ! Update time used in multiple-linear-regression analysis 223 225 224 226 #if defined key_top -
NEMO/branches/2019/dev_r11943_MERGE_2019/src/OCE/step_oce.F90
r12150 r12205 19 19 USE sbccpl ! surface boundary condition: coupled formulation (call send at end of step) 20 20 USE sbcapr ! surface boundary condition: atmospheric pressure 21 USE sbctide ! Tide initialisation21 USE tide_mod, ONLY : ln_tide, tide_update 22 22 USE sbcwave ! Wave intialisation 23 23 … … 81 81 USE diahth ! thermocline depth (dia_hth routine) 82 82 USE diahsb ! heat, salt and volume budgets (dia_hsb routine) 83 USE diaharm84 83 USE diacfl 85 84 USE diaobs ! Observation operator 85 USE diadetide ! Weights computation for daily detiding of model diagnostics 86 USE diamlr ! IOM context management for multiple-linear-regression analysis 86 87 USE flo_oce ! floats variables 87 88 USE floats ! floats computation (flo_stp routine) -
NEMO/branches/2019/dev_r11943_MERGE_2019/tests/CANAL/EXPREF/field_def_nemo-oce.xml
r9572 r12205 110 110 <field id="tpt_dep" long_name="T-point depth" standard_name="depth_below_geoid" unit="m" grid_ref="grid_T_3D" /> 111 111 <field id="e3tdef" long_name="T-cell thickness deformation" unit="%" grid_ref="grid_T_3D" /> 112 </field_group>113 114 <!-- Tides -->115 116 <field_group id="Tides_T" grid_ref="grid_T_2D" operation="once" >117 <!-- tidal composante -->118 <field id="M2x" long_name="M2 Elevation harmonic real part " unit="m" />119 <field id="M2y" long_name="M2 Elevation harmonic imaginary part" unit="m" />120 <field id="S2x" long_name="S2 Elevation harmonic real part " unit="m" />121 <field id="S2y" long_name="S2 Elevation harmonic imaginary part" unit="m" />122 <field id="N2x" long_name="N2 Elevation harmonic real part " unit="m" />123 <field id="N2y" long_name="N2 Elevation harmonic imaginary part" unit="m" />124 <field id="K1x" long_name="K1 Elevation harmonic real part " unit="m" />125 <field id="K1y" long_name="K1 Elevation harmonic imaginary part" unit="m" />126 <field id="O1x" long_name="O1 Elevation harmonic real part " unit="m" />127 <field id="O1y" long_name="O1 Elevation harmonic imaginary part" unit="m" />128 <field id="Q1x" long_name="Q1 Elevation harmonic real part " unit="m" />129 <field id="Q1y" long_name="Q1 Elevation harmonic imaginary part" unit="m" />130 <field id="M4x" long_name="M4 Elevation harmonic real part " unit="m" />131 <field id="M4y" long_name="M4 Elevation harmonic imaginary part" unit="m" />132 <field id="K2x" long_name="K2 Elevation harmonic real part " unit="m" />133 <field id="K2y" long_name="K2 Elevation harmonic imaginary part" unit="m" />134 <field id="P1x" long_name="P1 Elevation harmonic real part " unit="m" />135 <field id="P1y" long_name="P1 Elevation harmonic imaginary part" unit="m" />136 <field id="Mfx" long_name="Mf Elevation harmonic real part " unit="m" />137 <field id="Mfy" long_name="Mf Elevation harmonic imaginary part" unit="m" />138 <field id="Mmx" long_name="Mm Elevation harmonic real part " unit="m" />139 <field id="Mmy" long_name="Mm Elevation harmonic imaginary part" unit="m" />140 </field_group>141 142 <field_group id="Tides_U" grid_ref="grid_U_2D" operation="once" >143 <field id="M2x_u" long_name="M2 current barotrope along i-axis harmonic real part " unit="m/s" />144 <field id="M2y_u" long_name="M2 current barotrope along i-axis harmonic imaginary part " unit="m/s" />145 <field id="S2x_u" long_name="S2 current barotrope along i-axis harmonic real part " unit="m/s" />146 <field id="S2y_u" long_name="S2 current barotrope along i-axis harmonic imaginary part " unit="m/s" />147 <field id="N2x_u" long_name="N2 current barotrope along i-axis harmonic real part " unit="m/s" />148 <field id="N2y_u" long_name="N2 current barotrope along i-axis harmonic imaginary part " unit="m/s" />149 <field id="K1x_u" long_name="K1 current barotrope along i-axis harmonic real part " unit="m/s" />150 <field id="K1y_u" long_name="K1 current barotrope along i-axis harmonic imaginary part " unit="m/s" />151 <field id="O1x_u" long_name="O1 current barotrope along i-axis harmonic real part " unit="m/s" />152 <field id="O1y_u" long_name="O1 current barotrope along i-axis harmonic imaginary part " unit="m/s" />153 <field id="Q1x_u" long_name="Q1 current barotrope along i-axis harmonic real part " unit="m/s" />154 <field id="Q1y_u" long_name="Q1 current barotrope along i-axis harmonic imaginary part " unit="m/s" />155 <field id="M4x_u" long_name="M4 current barotrope along i-axis harmonic real part " unit="m/s" />156 <field id="M4y_u" long_name="M4 current barotrope along i-axis harmonic imaginary part " unit="m/s" />157 <field id="K2x_u" long_name="K2 current barotrope along i-axis harmonic real part " unit="m/s" />158 <field id="K2y_u" long_name="K2 current barotrope along i-axis harmonic imaginary part " unit="m/s" />159 <field id="P1x_u" long_name="P1 current barotrope along i-axis harmonic real part " unit="m/s" />160 <field id="P1y_u" long_name="P1 current barotrope along i-axis harmonic imaginary part " unit="m/s" />161 <field id="Mfx_u" long_name="Mf current barotrope along i-axis harmonic real part " unit="m/s" />162 <field id="Mfy_u" long_name="Mf current barotrope along i-axis harmonic imaginary part " unit="m/s" />163 <field id="Mmx_u" long_name="Mm current barotrope along i-axis harmonic real part " unit="m/s" />164 <field id="Mmy_u" long_name="Mm current barotrope along i-axis harmonic imaginary part " unit="m/s" />165 </field_group>166 167 <field_group id="Tides_V" grid_ref="grid_V_2D" operation="once" >168 <field id="M2x_v" long_name="M2 current barotrope along j-axis harmonic real part " unit="m/s" />169 <field id="M2y_v" long_name="M2 current barotrope along j-axis harmonic imaginary part " unit="m/s" />170 <field id="S2x_v" long_name="S2 current barotrope along j-axis harmonic real part " unit="m/s" />171 <field id="S2y_v" long_name="S2 current barotrope along j-axis harmonic imaginary part " unit="m/s" />172 <field id="N2x_v" long_name="N2 current barotrope along j-axis harmonic real part " unit="m/s" />173 <field id="N2y_v" long_name="N2 current barotrope along j-axis harmonic imaginary part " unit="m/s" />174 <field id="K1x_v" long_name="K1 current barotrope along j-axis harmonic real part " unit="m/s" />175 <field id="K1y_v" long_name="K1 current barotrope along j-axis harmonic imaginary part " unit="m/s" />176 <field id="O1x_v" long_name="O1 current barotrope along j-axis harmonic real part " unit="m/s" />177 <field id="O1y_v" long_name="O1 current barotrope along j-axis harmonic imaginary part " unit="m/s" />178 <field id="Q1x_v" long_name="Q1 current barotrope along j-axis harmonic real part " unit="m/s" />179 <field id="Q1y_v" long_name="Q1 current barotrope along j-axis harmonic imaginary part " unit="m/s" />180 <field id="M4x_v" long_name="M4 current barotrope along j-axis harmonic real part " unit="m/s" />181 <field id="M4y_v" long_name="M4 current barotrope along j-axis harmonic imaginary part " unit="m/s" />182 <field id="K2x_v" long_name="K2 current barotrope along j-axis harmonic real part " unit="m/s" />183 <field id="K2y_v" long_name="K2 current barotrope along j-axis harmonic imaginary part " unit="m/s" />184 <field id="P1x_v" long_name="P1 current barotrope along j-axis harmonic real part " unit="m/s" />185 <field id="P1y_v" long_name="P1 current barotrope along j-axis harmonic imaginary part " unit="m/s" />186 <field id="Mfx_v" long_name="Mf current barotrope along j-axis harmonic real part " unit="m/s" />187 <field id="Mfy_v" long_name="Mf current barotrope along j-axis harmonic imaginary part " unit="m/s" />188 <field id="Mmx_v" long_name="Mm current barotrope along j-axis harmonic real part " unit="m/s" />189 <field id="Mmy_v" long_name="Mm current barotrope along j-axis harmonic imaginary part " unit="m/s" />190 112 </field_group> 191 113 -
NEMO/branches/2019/dev_r11943_MERGE_2019/tests/CANAL/EXPREF/namelist_cfg
r11536 r12205 276 276 !! namdiu Cool skin and warm layer models (default: OFF) 277 277 !! namflo float parameters (default: OFF) 278 !! nam_diaharm Harmonic analysis of tidal constituents (default: OFF)279 278 !! nam_diadct transports through some sections (default: OFF) 280 279 !! nam_diatmb Top Middle Bottom Output (default: OFF) -
NEMO/branches/2019/dev_r11943_MERGE_2019/tests/ICE_ADV1D/EXPREF/namelist_cfg
r11536 r12205 202 202 !! namdiu Cool skin and warm layer models (default: OFF) 203 203 !! namflo float parameters (default: OFF) 204 !! nam_diaharm Harmonic analysis of tidal constituents (default: OFF)205 204 !! nam_diadct transports through some sections (default: OFF) 206 205 !! nam_diatmb Top Middle Bottom Output (default: OFF) -
NEMO/branches/2019/dev_r11943_MERGE_2019/tests/ICE_ADV1D/EXPREF/namelist_cfg_120pts
r11536 r12205 202 202 !! namdiu Cool skin and warm layer models (default: OFF) 203 203 !! namflo float parameters (default: OFF) 204 !! nam_diaharm Harmonic analysis of tidal constituents (default: OFF)205 204 !! nam_diadct transports through some sections (default: OFF) 206 205 !! nam_diatmb Top Middle Bottom Output (default: OFF) -
NEMO/branches/2019/dev_r11943_MERGE_2019/tests/ICE_ADV1D/EXPREF/namelist_cfg_240pts
r11536 r12205 202 202 !! namdiu Cool skin and warm layer models (default: OFF) 203 203 !! namflo float parameters (default: OFF) 204 !! nam_diaharm Harmonic analysis of tidal constituents (default: OFF)205 204 !! nam_diadct transports through some sections (default: OFF) 206 205 !! nam_diatmb Top Middle Bottom Output (default: OFF) -
NEMO/branches/2019/dev_r11943_MERGE_2019/tests/ICE_ADV1D/EXPREF/namelist_cfg_60pts
r11536 r12205 202 202 !! namdiu Cool skin and warm layer models (default: OFF) 203 203 !! namflo float parameters (default: OFF) 204 !! nam_diaharm Harmonic analysis of tidal constituents (default: OFF)205 204 !! nam_diadct transports through some sections (default: OFF) 206 205 !! nam_diatmb Top Middle Bottom Output (default: OFF) -
NEMO/branches/2019/dev_r11943_MERGE_2019/tests/ICE_ADV2D/EXPREF/namelist_cfg
r11536 r12205 202 202 !! namdiu Cool skin and warm layer models (default: OFF) 203 203 !! namflo float parameters (default: OFF) 204 !! nam_diaharm Harmonic analysis of tidal constituents (default: OFF)205 204 !! nam_diadct transports through some sections (default: OFF) 206 205 !! nam_diatmb Top Middle Bottom Output (default: OFF) -
NEMO/branches/2019/dev_r11943_MERGE_2019/tests/ICE_AGRIF/EXPREF/1_namelist_cfg
r11536 r12205 202 202 !! namdiu Cool skin and warm layer models (default: OFF) 203 203 !! namflo float parameters (default: OFF) 204 !! nam_diaharm Harmonic analysis of tidal constituents (default: OFF)205 204 !! nam_diadct transports through some sections (default: OFF) 206 205 !! nam_diatmb Top Middle Bottom Output (default: OFF) -
NEMO/branches/2019/dev_r11943_MERGE_2019/tests/ICE_AGRIF/EXPREF/namelist_cfg
r11536 r12205 202 202 !! namdiu Cool skin and warm layer models (default: OFF) 203 203 !! namflo float parameters (default: OFF) 204 !! nam_diaharm Harmonic analysis of tidal constituents (default: OFF)205 204 !! nam_diadct transports through some sections (default: OFF) 206 205 !! nam_diatmb Top Middle Bottom Output (default: OFF) -
NEMO/branches/2019/dev_r11943_MERGE_2019/tests/ISOMIP/EXPREF/namelist_cfg
r12150 r12205 428 428 !! namdiu Cool skin and warm layer models (default: OFF) 429 429 !! namflo float parameters (default: OFF) 430 !! nam_diaharm Harmonic analysis of tidal constituents (default: OFF)431 430 !! nam_diadct transports through some sections (default: OFF) 432 431 !! nam_diatmb Top Middle Bottom Output (default: OFF) -
NEMO/branches/2019/dev_r11943_MERGE_2019/tests/LOCK_EXCHANGE/EXPREF/namelist_FCT2_flux_ubs_cfg
r12150 r12205 276 276 !! namdiu Cool skin and warm layer models (default: OFF) 277 277 !! namflo float parameters (default: OFF) 278 !! nam_diaharm Harmonic analysis of tidal constituents (default: OFF)279 278 !! nam_diadct transports through some sections (default: OFF) 280 279 !! nam_diatmb Top Middle Bottom Output (default: OFF) -
NEMO/branches/2019/dev_r11943_MERGE_2019/tests/OVERFLOW/EXPREF/namelist_zps_FCT4_flux_ubs_cfg
r11848 r12205 294 294 !! namdiu Cool skin and warm layer models (default: OFF) 295 295 !! namflo float parameters (default: OFF) 296 !! nam_diaharm Harmonic analysis of tidal constituents (default: OFF)297 296 !! nam_diadct transports through some sections (default: OFF) 298 297 !! nam_diatmb Top Middle Bottom Output (default: OFF) -
NEMO/branches/2019/dev_r11943_MERGE_2019/tests/VORTEX/EXPREF/1_namelist_cfg
r11536 r12205 273 273 !! namdiu Cool skin and warm layer models (default: OFF) 274 274 !! namflo float parameters (default: OFF) 275 !! nam_diaharm Harmonic analysis of tidal constituents (default: OFF)276 275 !! nam_diadct transports through some sections (default: OFF) 277 276 !! nam_diatmb Top Middle Bottom Output (default: OFF) -
NEMO/branches/2019/dev_r11943_MERGE_2019/tests/VORTEX/EXPREF/namelist_cfg
r11536 r12205 264 264 !! namdiu Cool skin and warm layer models (default: OFF) 265 265 !! namflo float parameters (default: OFF) 266 !! nam_diaharm Harmonic analysis of tidal constituents (default; OFF)267 266 !! nam_diadct transports through some sections (default: OFF) 268 267 !! nam_diatmb Top Middle Bottom Output (default: OFF) -
NEMO/branches/2019/dev_r11943_MERGE_2019/tests/WAD/EXPREF/namelist_cfg
r11949 r12205 422 422 !! namdiu Cool skin and warm layer models (default: OFF) 423 423 !! namflo float parameters (default: OFF) 424 !! nam_diaharm Harmonic analysis of tidal constituents (default: OFF)425 424 !! nam_diadct transports through some sections (default: OFF) 426 425 !! nam_diatmb Top Middle Bottom Output (default: OFF)
Note: See TracChangeset
for help on using the changeset viewer.