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

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

Changeset 12205


Ignore:
Timestamp:
2019-12-12T11:52:50+01:00 (4 years ago)
Author:
acc
Message:

2019/dev_r11943_MERGE_2019: Merge in dev_r11879_ENHANCE-05_SimonM-Harmonic_Analysis. SETTE tested

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  
    381381!!   namdiu       Cool skin and warm layer models                       (default: OFF) 
    382382!!   namflo       float parameters                                      (default: OFF) 
    383 !!   nam_diaharm  Harmonic analysis of tidal constituents               (default: OFF) 
    384383!!   nam_diadct   transports through some sections                      (default: OFF) 
    385384!!   nam_diatmb   Top Middle Bottom Output                              (default: OFF) 
  • NEMO/branches/2019/dev_r11943_MERGE_2019/cfgs/AGRIF_DEMO/EXPREF/2_namelist_cfg

    r12182 r12205  
    330330!!   namdiu       Cool skin and warm layer models                       (default: OFF) 
    331331!!   namflo       float parameters                                      (default: OFF) 
    332 !!   nam_diaharm  Harmonic analysis of tidal constituents               (default: OFF) 
    333332!!   nam_diadct   transports through some sections                      (default: OFF) 
    334333!!   nam_diatmb   Top Middle Bottom Output                              (default: OFF) 
  • NEMO/branches/2019/dev_r11943_MERGE_2019/cfgs/AGRIF_DEMO/EXPREF/3_namelist_cfg

    r12182 r12205  
    312312!!   namdiu       Cool skin and warm layer models                       (default: OFF) 
    313313!!   namflo       float parameters                                      (default: OFF) 
    314 !!   nam_diaharm  Harmonic analysis of tidal constituents               (default: OFF) 
    315314!!   nam_diadct   transports through some sections                      (default: OFF) 
    316315!!   nam_diatmb   Top Middle Bottom Output                              (default: OFF) 
  • NEMO/branches/2019/dev_r11943_MERGE_2019/cfgs/AGRIF_DEMO/EXPREF/namelist_cfg

    r12182 r12205  
    381381!!   namdiu       Cool skin and warm layer models                       (default: OFF) 
    382382!!   namflo       float parameters                                      (default: OFF) 
    383 !!   nam_diaharm  Harmonic analysis of tidal constituents               (default: OFF) 
    384383!!   nam_diadct   transports through some sections                      (default: OFF) 
    385384!!   nam_diatmb   Top Middle Bottom Output                              (default: OFF) 
  • NEMO/branches/2019/dev_r11943_MERGE_2019/cfgs/AMM12/EXPREF/context_nemo.xml

    r12193 r12205  
    2222    <field_definition src="./field_def_nemo-oce.xml"/>   <!--  NEMO ocean dynamics                     --> 
    2323 
     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 
    2477<!-- Files definition --> 
    2578    <file_definition src="./file_def_nemo-oce.xml"/>     <!--  NEMO ocean dynamics                     --> 
    2679 
     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 
    27100<!-- Axis definition --> 
    28101    <axis_definition src="./axis_def_nemo.xml"/> 
    29   
     102 
    30103<!-- Domain definition --> 
    31104    <domain_definition src="./domain_def_nemo.xml"/> 
  • NEMO/branches/2019/dev_r11943_MERGE_2019/cfgs/AMM12/EXPREF/file_def_nemo-oce.xml

    r10641 r12205  
    4949 
    5050      <file_group id="25h_mean" output_freq="1d" output_level="10" enabled=".TRUE."> 
    51  
    5251 
    5352     <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  
    154154&nam_tide      !   tide parameters                                      (default: OFF) 
    155155!----------------------------------------------------------------------- 
    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' 
    174173/ 
    175174!----------------------------------------------------------------------- 
     
    351350!!   namdiu       Cool skin and warm layer models                       (default: OFF) 
    352351!!   namflo       float parameters                                      (default: OFF) 
    353 !!   nam_diaharm  Harmonic analysis of tidal constituents               (default: OFF) 
    354352!!   nam_diadct   transports through some sections                      (default: OFF) 
    355353!!   nam_diatmb   Top Middle Bottom Output                              (default: OFF) 
  • NEMO/branches/2019/dev_r11943_MERGE_2019/cfgs/C1D_PAPA/EXPREF/namelist_cfg

    r12182 r12205  
    420420!!   namdiu       Cool skin and warm layer models                       (default: OFF) 
    421421!!   namflo       float parameters                                      (default: OFF) 
    422 !!   nam_diaharm  Harmonic analysis of tidal constituents               (default: OFF) 
    423422!!   nam_diadct   transports through some sections                      (default: OFF) 
    424423!!   nam_diatmb   Top Middle Bottom Output                              (default: OFF) 
     
    448447/ 
    449448!----------------------------------------------------------------------- 
    450 &nam_diaharm   !   Harmonic analysis of tidal constituents              (default: OFF) 
    451 !----------------------------------------------------------------------- 
    452 / 
    453 !----------------------------------------------------------------------- 
    454449&nam_diadct    ! transports through some sections                       (default: OFF) 
    455450!----------------------------------------------------------------------- 
  • NEMO/branches/2019/dev_r11943_MERGE_2019/cfgs/GYRE_BFM/EXPREF/namelist_cfg

    r11536 r12205  
    228228!!   namdiu       Cool skin and warm layer models                       (default: OFF) 
    229229!!   namflo       float parameters                                      (default: OFF) 
    230 !!   nam_diaharm  Harmonic analysis of tidal constituents               (default: OFF) 
    231230!!   nam_diadct   transports through some sections                      (default: OFF) 
    232231!!   nam_diatmb   Top Middle Bottom Output                              (default: OFF) 
  • NEMO/branches/2019/dev_r11943_MERGE_2019/cfgs/GYRE_PISCES/EXPREF/namelist_cfg

    r11536 r12205  
    222222!!   namdiu       Cool skin and warm layer models                       (default: OFF) 
    223223!!   namflo       float parameters                                      (default: OFF) 
    224 !!   nam_diaharm  Harmonic analysis of tidal constituents               (default: OFF) 
    225224!!   nam_diadct   transports through some sections                      (default: OFF) 
    226225!!   nam_diatmb   Top Middle Bottom Output                              (default: OFF) 
  • NEMO/branches/2019/dev_r11943_MERGE_2019/cfgs/ORCA2_ICE_PISCES/EXPREF/namelist_cfg

    r12182 r12205  
    391391!!   namdiu       Cool skin and warm layer models                       (default: OFF) 
    392392!!   namflo       float parameters                                      (default: OFF) 
    393 !!   nam_diaharm  Harmonic analysis of tidal constituents               (default: OFF) 
    394393!!   nam_diadct   transports through some sections                      (default: OFF) 
    395394!!   nam_diatmb   Top Middle Bottom Output                              (default: OFF) 
  • NEMO/branches/2019/dev_r11943_MERGE_2019/cfgs/ORCA2_OFF_PISCES/EXPREF/namelist_cfg

    r11536 r12205  
    377377!!   namdiu       Cool skin and warm layer models                       (default: OFF) 
    378378!!   namflo       float parameters                                      (default: OFF) 
    379 !!   nam_diaharm  Harmonic analysis of tidal constituents               (default: OFF) 
    380379!!   nam_diadct   transports through some sections                      (default: OFF) 
    381380!!   nam_diatmb   Top Middle Bottom Output                              (default: OFF) 
     
    405404/ 
    406405!----------------------------------------------------------------------- 
    407 &nam_diaharm   !   Harmonic analysis of tidal constituents              (default: OFF) 
    408 !----------------------------------------------------------------------- 
    409 / 
    410 !----------------------------------------------------------------------- 
    411406&nam_diadct    !   transports through some sections                     (default: OFF) 
    412407!----------------------------------------------------------------------- 
  • NEMO/branches/2019/dev_r11943_MERGE_2019/cfgs/ORCA2_OFF_TRC/EXPREF/namelist_cfg

    r11536 r12205  
    375375!!   namdiu       Cool skin and warm layer models                       (default: OFF) 
    376376!!   namflo       float parameters                                      (default: OFF) 
    377 !!   nam_diaharm  Harmonic analysis of tidal constituents               ("key_diaharm") 
    378377!!   namdct       transports through some sections                      ("key_diadct") 
    379378!!   nam_diatmb   Top Middle Bottom Output                              (default: OFF) 
     
    403402/ 
    404403!----------------------------------------------------------------------- 
    405 &nam_diaharm   !   Harmonic analysis of tidal constituents              ("key_diaharm") 
    406 !----------------------------------------------------------------------- 
    407 / 
    408 !----------------------------------------------------------------------- 
    409404&namdct        ! transports through some sections                       ("key_diadct") 
    410405!----------------------------------------------------------------------- 
  • NEMO/branches/2019/dev_r11943_MERGE_2019/cfgs/ORCA2_SAS_ICE/EXPREF/namelist_cfg

    r12182 r12205  
    196196!!   namdiu       Cool skin and warm layer models                       (default: OFF) 
    197197!!   namflo       float parameters                                      (default: OFF) 
    198 !!   nam_diaharm  Harmonic analysis of tidal constituents               (default: OFF) 
    199198!!   nam_diadct   transports through some sections                      (default: OFF) 
    200199!!   nam_diatmb   Top Middle Bottom Output                              (default: OFF) 
  • NEMO/branches/2019/dev_r11943_MERGE_2019/cfgs/SHARED/field_def_nemo-oce.xml

    r12193 r12205  
    11<?xml version="1.0"?>  
    22    <!-- $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> 
    375     
    476    <!--  
     
    880============================================================================================================ 
    981    --> 
    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> 
    11102 
    12103    <!--  
     
    138229 
    139230      <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" /> 
    213252      </field_group> 
    214253 
  • NEMO/branches/2019/dev_r11943_MERGE_2019/cfgs/SHARED/grid_def_nemo.xml

    r12193 r12205  
     1<?xml version="1.0"?> 
    12<!-- 
    23============================================================================================================ 
     
    45============================================================================================================ 
    56    --> 
    6   
     7 
    78    <grid_definition>     
    89 
     
    5960        <!--  --> 
    6061       <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        <!--  --> 
    7866       <grid id="grid_EqT" > 
    7967         <domain domain_ref="EqT" /> 
     
    237225      </grid> 
    238226 
    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"> 
    244236        <axis axis_ref="section_ice"/> 
    245      </grid> 
    246  
    247      <grid id="grid_4strait_ice"> 
    248        <axis axis_ref="section_ice"/> 
    249      </grid> 
     237      </grid> 
    250238  
    251239      <!-- scalars --> 
     
    253241        <scalar/> 
    254242      </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  
    642642!----------------------------------------------------------------------- 
    643643   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 
    644647      ln_tide_pot   = .false.               !  use tidal potential forcing 
     648         rn_tide_gamma = 0.7                   ! Tidal tilt factor 
    645649         ln_scal_load  = .false.               ! Use scalar approximation for 
    646650            rn_scal_load = 0.094               !     load potential 
     
    649653            ! 
    650654      ln_tide_ramp  = .false.               !  Use linear ramp for tides at startup 
    651          rdttideramp   =    0.                 !  ramp duration in days 
    652       clname(1)    = 'DUMMY'               !  name of constituent - all tidal components must be set in namelist_cfg 
     655         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 
    653657/ 
    654658!----------------------------------------------------------------------- 
     
    726730   filtide          = 'bdydta/amm12_bdytide_'   !  file name root of tidal forcing files 
    727731   ln_bdytide_2ddta = .false.                   ! 
    728    ln_bdytide_conj  = .false.                   ! 
    729732/ 
    730733 
     
    11921195!!   namdiu       Cool skin and warm layer models                       (default: OFF) 
    11931196!!   namflo       float parameters                                      (default: OFF) 
    1194 !!   nam_diaharm  Harmonic analysis of tidal constituents               (default: OFF) 
    11951197!!   nam_diadct   transports through some sections                      (default: OFF) 
    11961198!!   nam_diatmb   Top Middle Bottom Output                              (default: OFF) 
     
    12451247      ln_ariane   = .true.    !    Input with Ariane tool convention(T) 
    12461248      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 not 
    1252        nit000_han = 1      !    First time step used for harmonic analysis 
    1253        nitend_han = 75     !    Last time step used for harmonic analysis 
    1254        nstep_han  = 15     !    Time step frequency for harmonic analysis 
    1255        tname(1)   = 'M2'   !    Name of tidal constituents 
    1256        tname(2)   = 'K1'   !              --- 
    12571249/ 
    12581250!----------------------------------------------------------------------- 
  • NEMO/branches/2019/dev_r11943_MERGE_2019/cfgs/SPITZ12/EXPREF/namelist_cfg

    r12182 r12205  
    138138&nam_tide      !   tide parameters                                      (default: OFF) 
    139139!----------------------------------------------------------------------- 
    140    ln_tide     = .true.      ! Activate tides 
    141       ln_tide_pot   = .false.                !  use tidal potential forcing 
    142       clname(1)     =   'M2'    !  name of constituent 
    143       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' 
    150150/ 
    151151!----------------------------------------------------------------------- 
     
    348348!!   namdiu       Cool skin and warm layer models                       (default: OFF) 
    349349!!   namflo       float parameters                                      (default: OFF) 
    350 !!   nam_diaharm  Harmonic analysis of tidal constituents               (default: OFF) 
    351350!!   nam_diadct   transports through some sections                      (default: OFF) 
    352351!!   nam_diatmb   Top Middle Bottom Output                              (default: OFF) 
  • NEMO/branches/2019/dev_r11943_MERGE_2019/doc/latex/NEMO/subfiles/chap_DIA.tex

    r11693 r12205  
    15801580 
    15811581%% ================================================================================================= 
    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 analysis 
    1597  
    1598  - \np{nitend_han}{nitend\_han} is the  last time step used for harmonic analysis 
    1599  
    1600  - \np{nstep_han}{nstep\_han}  is the  time step frequency for harmonic analysis 
    1601  
    1602 % - \np{nb_ana}{nb\_ana}     is the number of harmonics to analyse 
    1603  
    1604  - \np{tname}{tname}       is an array with names of tidal constituents to analyse 
    1605  
    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 %% ================================================================================================= 
    16281582\section[Transports across sections (\texttt{\textbf{key\_diadct}})]{Transports across sections (\protect\key{diadct})} 
    16291583\label{sec:DIA_diag_dct} 
  • NEMO/branches/2019/dev_r11943_MERGE_2019/src/OCE/BDY/bdydta.F90

    r12145 r12205  
    2323   USE phycst         ! physical constants 
    2424   USE sbcapr         ! atmospheric pressure forcing 
    25    USE sbctide        ! Tidal forcing or not 
     25   USE tide_mod, ONLY: ln_tide ! tidal forcing 
    2626   USE bdy_oce        ! ocean open boundary conditions   
    2727   USE bdytides       ! tidal forcing at boundaries 
     
    7575CONTAINS 
    7676 
    77    SUBROUTINE bdy_dta( kt, Kmm, kit, kt_offset ) 
     77   SUBROUTINE bdy_dta( kt, Kmm, kt_offset ) 
    7878      !!---------------------------------------------------------------------- 
    7979      !!                   ***  SUBROUTINE bdy_dta  *** 
     
    8686      INTEGER, INTENT(in)           ::   kt           ! ocean time-step index  
    8787      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 
    9089      !                                               ! is present then units = subcycle timesteps. 
    9190      !                                               ! kt_offset = 0 => get data at "now" time level 
     
    106105      ! Initialise data arrays once for all from initial conditions where required 
    107106      !--------------------------------------------------------------------------- 
    108       IF( kt == nit000 .AND. .NOT.PRESENT(kit) ) THEN 
     107      IF( kt == nit000 ) THEN 
    109108 
    110109         ! Calculate depth-mean currents 
     
    217216         ! read/update all bdy data 
    218217         ! ------------------------ 
    219          CALL fld_read( kt, 1, bf_alias, kit = kit, kt_offset = kt_offset ) 
     218         CALL fld_read( kt, 1, bf_alias, kt_offset = kt_offset ) 
    220219         ! apply some corrections in some specific cases... 
    221220         ! -------------------------------------------------- 
     
    275274            END DO 
    276275         ENDIF   ! ltotvel 
    277  
    278          ! update tidal harmonic forcing 
    279          IF( PRESENT(kit) .AND. nn_dyn2d_dta(jbdy) .GE. 2 ) THEN 
    280             CALL bdytide_update( kt = kt, idx = idx_bdy(jbdy), dta = dta_alias, td = tides(jbdy),   &  
    281                &                 kit = kit, kt_offset = kt_offset ) 
    282          ENDIF 
    283276 
    284277         !  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  
    2222   USE bdydta         ! open boundary cond. setting   (bdy_dta_init routine) 
    2323   USE bdytides       ! open boundary cond. setting   (bdytide_init routine) 
    24    USE sbctide        ! Tidal forcing or not 
     24   USE tide_mod, ONLY: ln_tide ! tidal forcing 
    2525   USE phycst   , ONLY: rday 
    2626   ! 
  • NEMO/branches/2019/dev_r11943_MERGE_2019/src/OCE/BDY/bdytides.F90

    r11960 r12205  
    1818   USE phycst         ! physical constants 
    1919   USE bdy_oce        ! ocean open boundary conditions 
    20    USE tideini        !  
     20   USE tide_mod       !  
    2121   USE daymod         ! calendar 
    2222   ! 
     
    3030 
    3131   PUBLIC   bdytide_init     ! routine called in bdy_init 
    32    PUBLIC   bdytide_update   ! routine called in bdy_dta 
    3332   PUBLIC   bdy_dta_tides    ! routine called in dyn_spg_ts 
    3433 
     
    4544   TYPE(OBC_DATA)  , PUBLIC, DIMENSION(jp_bdy) :: dta_bdy_s  !: bdy external data (slow component) 
    4645 
     46   INTEGER ::   kt_tide 
     47 
    4748   !!---------------------------------------------------------------------- 
    4849   !! NEMO/OCE 4.0 , NEMO Consortium (2018) 
     
    6465      CHARACTER(len=80)                         ::   filtide             !: Filename root for tidal input files 
    6566      LOGICAL                                   ::   ln_bdytide_2ddta    !: If true, read 2d harmonic data 
    66       LOGICAL                                   ::   ln_bdytide_conj     !: If true, assume complex conjugate tidal data 
    6767      !! 
    6868      INTEGER                                   ::   ib_bdy, itide, ib   !: dummy loop indices 
     
    7979      TYPE(TIDES_DATA),  POINTER                ::   td                  !: local short cut    
    8080      !! 
    81       NAMELIST/nambdy_tide/filtide, ln_bdytide_2ddta, ln_bdytide_conj 
     81      NAMELIST/nambdy_tide/filtide, ln_bdytide_2ddta 
    8282      !!---------------------------------------------------------------------- 
    8383      ! 
     
    119119            IF(lwp) WRITE(numout,*) '          Namelist nambdy_tide : tidal harmonic forcing at open boundaries' 
    120120            IF(lwp) WRITE(numout,*) '             read tidal data in 2d files: ', ln_bdytide_2ddta 
    121             IF(lwp) WRITE(numout,*) '             assume complex conjugate   : ', ln_bdytide_conj 
    122121            IF(lwp) WRITE(numout,*) '             Number of tidal components to read: ', nb_harmo 
    123122            IF(lwp) THEN  
    124123                    WRITE(numout,*) '             Tidal components: '  
    125124               DO itide = 1, nb_harmo 
    126                   WRITE(numout,*)  '                 ', Wave(ntide(itide))%cname_tide  
     125                  WRITE(numout,*)  '                 ', tide_harmonics(itide)%cname_tide  
    127126               END DO 
    128127            ENDIF  
     
    165164               igrd = 1                       ! Everything is at T-points here 
    166165               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(:,:) )  
    169168                  DO ib = 1, ilen0(igrd) 
    170169                     ii = idx_bdy(ib_bdy)%nbi(ib,igrd) 
     
    182181               igrd = 2                       ! Everything is at U-points here 
    183182               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(:,:) ) 
    186185                  DO ib = 1, ilen0(igrd) 
    187186                     ii = idx_bdy(ib_bdy)%nbi(ib,igrd) 
     
    199198               igrd = 3                       ! Everything is at V-points here 
    200199               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(:,:) ) 
    203202                  DO ib = 1, ilen0(igrd) 
    204203                     ii = idx_bdy(ib_bdy)%nbi(ib,igrd) 
     
    224223               DO itide = 1, nb_harmo 
    225224                  !                                                              ! 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' 
    227226                  CALL iom_open( clfile, inum ) 
    228227                  CALL fld_map( inum, 'z1' , dta_read(1:ilen0(1),1:1,1:1) , 1, idx_bdy(ib_bdy)%nbmap(:,1) ) 
     
    232231                  CALL iom_close( inum ) 
    233232                  !                                                              ! 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' 
    235234                  CALL iom_open( clfile, inum ) 
    236235                  CALL fld_map( inum, 'u1' , dta_read(1:ilen0(2),1:1,1:1) , 1, idx_bdy(ib_bdy)%nbmap(:,2) ) 
     
    240239                  CALL iom_close( inum ) 
    241240                  !                                                              ! 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' 
    243242                  CALL iom_open( clfile, inum ) 
    244243                  CALL fld_map( inum, 'v1' , dta_read(1:ilen0(3),1:1,1:1) , 1, idx_bdy(ib_bdy)%nbmap(:,3) ) 
     
    254253            ENDIF ! ln_bdytide_2ddta=.true. 
    255254            ! 
    256             IF( ln_bdytide_conj ) THEN    ! assume complex conjugate in data files 
    257                td%ssh0(:,:,2) = - td%ssh0(:,:,2) 
    258                td%u0  (:,:,2) = - td%u0  (:,:,2) 
    259                td%v0  (:,:,2) = - td%v0  (:,:,2) 
    260             ENDIF 
    261             ! 
    262255            ! Allocate slow varying data in the case of time splitting: 
    263256            ! Do it anyway because at this stage knowledge of free surface scheme is unknown 
     
    274267      ! 
    275268   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 counter 
    286       TYPE(OBC_INDEX)  , INTENT(in   ) ::   idx         ! OBC indices 
    287       TYPE(OBC_DATA)   , INTENT(inout) ::   dta         ! OBC external data 
    288       TYPE(TIDES_DATA) , INTENT(inout) ::   td          ! tidal harmonics data 
    289       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 kit 
    291       !                                                 ! is present then units = subcycle timesteps. 
    292       !                                                 ! kt_offset = 0  => get data at "now"    time level 
    293       !                                                 ! kt_offset = -1 => get data at "before" time level 
    294       !                                                 ! kt_offset = +1 => get data at "after"  time level 
    295       !                                                 ! etc. 
    296       ! 
    297       INTEGER  ::   itide, igrd, ib       ! dummy loop indices 
    298       INTEGER  ::   time_add              ! time offset in units of timesteps 
    299       INTEGER, DIMENSION(3) ::   ilen0    ! length of boundary data (from OBC arrays) 
    300       REAL(wp) ::   z_arg, z_sarg, zflag, zramp   ! local scalars     
    301       REAL(wp), DIMENSION(jpmax_harmo) :: z_sist, z_cost 
    302       !!---------------------------------------------------------------------- 
    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=1 
    309       IF ( PRESENT(kit) ) THEN 
    310         IF ( kit /= 1 ) zflag=0 
    311       ENDIF 
    312  
    313       IF ( (nsec_day == NINT(0.5_wp * rdt) .OR. kt==nit000) .AND. zflag==1 ) THEN 
    314         ! 
    315         kt_tide = kt - (nsec_day - 0.5_wp * rdt)/rdt 
    316         ! 
    317         IF(lwp) THEN 
    318            WRITE(numout,*) 
    319            WRITE(numout,*) 'bdytide_update : (re)Initialization of the tidal bdy forcing at kt=',kt 
    320            WRITE(numout,*) '~~~~~~~~~~~~~~ ' 
    321         ENDIF 
    322         ! 
    323         CALL tide_init_elevation ( idx, td ) 
    324         CALL tide_init_velocities( idx, td ) 
    325         ! 
    326       ENDIF  
    327  
    328       time_add = 0 
    329       IF( PRESENT(kt_offset) ) THEN 
    330          time_add = kt_offset 
    331       ENDIF 
    332           
    333       IF( PRESENT(kit) ) THEN   
    334          z_arg = ((kt-kt_tide) * rdt + (kit+0.5_wp*(time_add-1)) * rdt / REAL(nn_baro,wp) ) 
    335       ELSE                               
    336          z_arg = ((kt-kt_tide)+time_add) * rdt 
    337       ENDIF 
    338  
    339       ! Linear ramp on tidal component at open boundaries  
    340       zramp = 1._wp 
    341       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_harmo 
    344          z_sarg = z_arg * omega_tide(itide) 
    345          z_cost(itide) = COS( z_sarg ) 
    346          z_sist(itide) = SIN( z_sarg ) 
    347       END DO 
    348  
    349       DO itide = 1, nb_harmo 
    350          igrd=1                              ! SSH on tracer grid 
    351          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 DO 
    354          igrd=2                              ! U grid 
    355          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 DO 
    358          igrd=3                              ! V grid 
    359          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 DO 
    362       END DO 
    363       ! 
    364    END SUBROUTINE bdytide_update 
    365269 
    366270 
     
    406310      ! Linear ramp on tidal component at open boundaries  
    407311      zramp = 1. 
    408       IF (ln_tide_ramp) zramp = MIN(MAX( (z_arg - nit000*rdt)/(rdttideramp*rday),0.),1.) 
     312      IF (ln_tide_ramp) zramp = MIN(MAX( (z_arg - nit000*rdt)/(rn_tide_ramp_dt*rday),0.),1.) 
    409313 
    410314      DO ib_bdy = 1,nb_bdy 
     
    447351            DO itide = 1, nb_harmo 
    448352               ! 
    449                z_sarg = (z_arg + zoff) * omega_tide(itide) 
     353               z_sarg = (z_arg + zoff) * tide_harmonics(itide)%omega 
    450354               z_cost = zramp * COS( z_sarg ) 
    451355               z_sist = zramp * SIN( z_sarg ) 
     
    505409         END DO 
    506410         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 
    509413         ENDDO 
    510414         DO ib = 1 , ilen0(igrd) 
     
    544448         END DO 
    545449         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 
    548452         ENDDO 
    549453         DO ib = 1, ilen0(igrd) 
     
    565469         END DO 
    566470         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 
    569473         ENDDO 
    570474         DO ib = 1, ilen0(igrd) 
  • NEMO/branches/2019/dev_r11943_MERGE_2019/src/OCE/DYN/dynspg.F90

    r11960 r12205  
    2121   USE dynspg_exp     ! surface pressure gradient     (dyn_spg_exp routine) 
    2222   USE dynspg_ts      ! surface pressure gradient     (dyn_spg_ts  routine) 
    23    USE sbctide        !  
    24    USE updtide        !  
     23   USE tide_mod       ! 
    2524   USE trd_oce        ! trends: ocean variables 
    2625   USE trddyn         ! trend manager: dynamics 
     
    4342   INTEGER, PARAMETER ::   np_EXP = 0   !       explicit time stepping 
    4443   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 
    4546 
    4647   !! * Substitutions 
     
    116117         IF( .NOT.ln_dynspg_ts .AND. ( ln_tide_pot .AND. ln_tide )  ) THEN   ! N.B. added directly at sub-time-step in ts-case 
    117118            ! 
    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) 
    119122            ! 
    120123            DO jj = 2, jpjm1                         ! add tide potential forcing 
  • NEMO/branches/2019/dev_r11943_MERGE_2019/src/OCE/DYN/dynspg_ts.F90

    r12150 r12205  
    4444   USE bdytides        ! open boundary condition data 
    4545   USE bdydyn2d        ! open boundary conditions on barotropic variables 
    46    USE sbctide         ! tides 
    47    USE updtide         ! tide potential 
     46   USE tide_mod        ! 
    4847   USE sbcwave         ! surface wave 
    4948   USE diatmb          ! Top,middle,bottom output 
     
    173172      REAL(wp), ALLOCATABLE, DIMENSION(:,:) :: ztwdmask, zuwdmask, zvwdmask ! ROMS wetting and drying masks at t,u,v points 
    174173      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 
    175175      !!---------------------------------------------------------------------- 
    176176      ! 
     
    444444         ! 
    445445         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 
    447451         ! 
    448452         !                    !==  extrapolation at mid-step  ==!   (jn+1/2) 
  • NEMO/branches/2019/dev_r11943_MERGE_2019/src/OCE/IOM/iom.F90

    r12193 r12205  
    5555   LOGICAL, PUBLIC, PARAMETER ::   lk_iomput = .FALSE.       !: iom_put flag 
    5656#endif 
    57    PUBLIC iom_init, iom_swap, iom_open, iom_close, iom_setkt, iom_varid, iom_get, iom_get_var 
     57   PUBLIC iom_init, iom_init_closedef, iom_swap, iom_open, iom_close, iom_setkt, iom_varid, iom_get, iom_get_var 
    5858   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_val 
     59   PUBLIC iom_use, iom_context_finalize, iom_update_file_name, iom_miss_val 
    6060 
    6161   PRIVATE iom_rp0d, iom_rp1d, iom_rp2d, iom_rp3d 
     
    6464#if defined key_iomput 
    6565   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_sdate 
     66   PRIVATE set_grid, set_grid_bounds, set_scalar, set_xmlatt, set_mooring, iom_sdate 
    6767   PRIVATE iom_set_rst_context, iom_set_rstw_active, iom_set_rstr_active 
    6868# endif 
     
    9292CONTAINS 
    9393 
    94    SUBROUTINE iom_init( cdname, fname, ld_tmppatch )  
     94   SUBROUTINE iom_init( cdname, fname, ld_tmppatch, ld_closedef )  
    9595      !!---------------------------------------------------------------------- 
    9696      !!                     ***  ROUTINE   *** 
     
    102102      CHARACTER(len=*), OPTIONAL, INTENT(in)  :: fname 
    103103      LOGICAL         , OPTIONAL, INTENT(in)  :: ld_tmppatch 
     104      LOGICAL         , OPTIONAL, INTENT(in)  :: ld_closedef 
    104105#if defined key_iomput 
    105106      ! 
     
    116117      INTEGER ::   nldi_save, nlei_save    !:      and close boundaries in output files 
    117118      INTEGER ::   nldj_save, nlej_save    !: 
     119      LOGICAL ::   ll_closedef = .TRUE. 
    118120      !!---------------------------------------------------------------------- 
    119121      ! 
     
    130132         IF( njmpp + jpj - 1 == jpjglo ) nlej = jpj 
    131133      ENDIF 
     134      IF ( PRESENT(ld_closedef) ) ll_closedef = ld_closedef 
    132135      ! 
    133136      ALLOCATE( zt_bnds(2,jpk), zw_bnds(2,jpk) ) 
     
    268271      ENDIF 
    269272      ! 
    270       ! end file definition 
     273      ! set time step length 
    271274      dtime%second = rdt 
    272275      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 
    275279      ! 
    276280      DEALLOCATE( zt_bnds, zw_bnds ) 
     
    283287      ! 
    284288   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 
    285303 
    286304   SUBROUTINE iom_set_rstw_var_active(field) 
  • NEMO/branches/2019/dev_r11943_MERGE_2019/src/OCE/TDE/tide_mod.F90

    r12065 r12205  
    1717   !!---------------------------------------------------------------------- 
    1818 
    19    USE oce, ONLY : sshn       ! sea-surface height 
     19   USE oce, ONLY : ssh        ! sea-surface height 
    2020   USE par_oce                ! ocean parameters 
    2121   USE phycst, ONLY : rpi     ! pi 
     
    132132      sn_tide_cnames(:) = '' 
    133133      ! Read Namelist nam_tide 
    134       REWIND( numnam_ref )              ! Namelist nam_tide in reference namelist : Tides 
    135134      READ  ( numnam_ref, nam_tide, IOSTAT = ios, ERR = 901) 
    136135901   IF( ios /= 0 )   CALL ctl_nam ( ios , 'nam_tide in reference namelist' ) 
    137136      ! 
    138       REWIND( numnam_cfg )              ! Namelist nam_tide in configuration namelist : Tides 
    139137      READ  ( numnam_cfg, nam_tide, IOSTAT = ios, ERR = 902 ) 
    140138902   IF( ios >  0 )   CALL ctl_nam ( ios , 'nam_tide in configuration namelist' ) 
     
    715713 
    716714 
    717    SUBROUTINE upd_tide(pdelta) 
     715   SUBROUTINE upd_tide(pdelta, Kmm) 
    718716      !!---------------------------------------------------------------------- 
    719717      !!                 ***  ROUTINE upd_tide  *** 
     
    726724      !!----------------------------------------------------------------------       
    727725      REAL, INTENT(in)              ::   pdelta      ! Temporal offset in seconds 
     726      INTEGER, INTENT(IN)           ::   Kmm         ! Time level index 
    728727      INTEGER                       ::   jk          ! Dummy loop index 
    729728      REAL(wp)                      ::   zt, zramp   ! Local scalars 
     
    755754      ! 
    756755      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 * sshn(:,:) ) 
     756         IF ( iom_use( "tide_pot" ) ) CALL iom_put( "tide_pot", pot_astro(:,:) + rn_scal_load * ssh(:,:,Kmm) ) 
    758757      END IF 
    759758      ! 
  • NEMO/branches/2019/dev_r11943_MERGE_2019/src/OCE/nemogcm.F90

    r12193 r12205  
    4646   USE closea         ! treatment of closed seas (for ln_closea) 
    4747   USE usrdef_nam     ! user defined configuration 
    48    USE tideini        ! tidal components initialization   (tide_ini routine) 
     48   USE tide_mod, ONLY : tide_init ! tidal components initialization   (tide_init routine) 
    4949   USE bdy_oce,  ONLY : ln_bdy 
    5050   USE bdyini         ! open boundary cond. setting       (bdy_init routine) 
     
    5959   USE diaobs         ! Observation diagnostics       (dia_obs_init routine) 
    6060   USE diacfl         ! CFL diagnostics               (dia_cfl_init routine) 
    61    USE diaharm        ! tidal harmonics diagnostics  (dia_harm_init routine) 
     61   USE diamlr         ! IOM context management for multiple-linear-regression analysis 
    6262   USE step           ! NEMO time-stepping                 (stp     routine) 
    6363   USE isfstp         ! ice shelf                     (isf_stp_init routine) 
     
    7575   USE diatmb         ! Top,middle,bottom output 
    7676   USE dia25h         ! 25h mean output 
     77   USE diadetide      ! Weights computation for daily detiding of model diagnostics 
    7778   USE sbc_oce , ONLY : lk_oasis 
    7879   USE wet_dry        ! Wetting and drying setting   (wad_init routine) 
     
    492493                           CALL dia_tmb_init    ! TMB outputs 
    493494                           CALL dia_25h_init( Nbb )    ! 25h mean  outputs 
    494                            CALL dia_harm_init   ! tidal harmonics outputs 
     495                           CALL dia_detide_init ! Weights computation for daily detiding of model diagnostics 
    495496      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 
    496498 
    497499      !                                      ! Assimilation increments 
  • NEMO/branches/2019/dev_r11943_MERGE_2019/src/OCE/step.F90

    r12193 r12205  
    107107                              
    108108      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 
    110112         IF( ln_crs      )   CALL iom_init( TRIM(cxios_context)//"_crs" )  ! for coarse grid 
    111113      ENDIF 
     
    117119      ! Update external forcing (tides, open boundaries, ice shelf interaction and surface boundary condition (including sea-ice) 
    118120      !<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<< 
    119       IF( ln_tide    )   CALL sbc_tide( kstp )                        ! update tide potential 
     121      IF( ln_tide    )   CALL tide_update( kstp )                     ! update tide potential 
    120122      IF( ln_apr_dyn )   CALL sbc_apr ( kstp )                        ! atmospheric pressure (NB: call before bdy_dta which needs ssh_ib)  
    121123      IF( ln_bdy     )   CALL bdy_dta ( kstp, Nnn, kt_offset = +1 )   ! update dynamic & tracer data at open boundaries 
     
    211213      ! diagnostics and outputs 
    212214      !<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<< 
    213       IF( ln_floats  )   CALL flo_stp ( kstp, Nbb, Nnn )        ! drifting Floats 
    214       IF( ln_diacfl  )   CALL dia_cfl ( kstp,      Nnn )        ! Courant number diagnostics 
    215                          CALL dia_hth ( kstp,      Nnn )        ! Thermocline depth (20 degres isotherm depth) 
    216       IF( ln_diadct  )   CALL dia_dct ( kstp,      Nnn )        ! Transports 
    217                          CALL dia_ar5 ( kstp,      Nnn )        ! ar5 diag 
    218                          CALL dia_ptr ( kstp,      Nnn )        ! Poleward adv/ldf TRansports diagnostics 
    219       IF( ln_diaharm )   CALL dia_harm( kstp,      Nnn )        ! Tidal harmonic analysis 
    220                          CALL dia_wri ( kstp,      Nnn )        ! ocean model: outputs 
    221       ! 
    222       IF( ln_crs     )   CALL crs_fld ( kstp,      Nnn )   ! ocean model: online field coarsening & output 
     215      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 
    223225       
    224226#if defined key_top 
  • NEMO/branches/2019/dev_r11943_MERGE_2019/src/OCE/step_oce.F90

    r12150 r12205  
    1919   USE sbccpl          ! surface boundary condition: coupled formulation (call send at end of step) 
    2020   USE sbcapr          ! surface boundary condition: atmospheric pressure 
    21    USE sbctide         ! Tide initialisation 
     21   USE tide_mod, ONLY : ln_tide, tide_update 
    2222   USE sbcwave         ! Wave intialisation 
    2323 
     
    8181   USE diahth          ! thermocline depth                (dia_hth routine) 
    8282   USE diahsb          ! heat, salt and volume budgets    (dia_hsb routine) 
    83    USE diaharm 
    8483   USE diacfl 
    8584   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 
    8687   USE flo_oce         ! floats variables 
    8788   USE floats          ! floats computation               (flo_stp routine) 
  • NEMO/branches/2019/dev_r11943_MERGE_2019/tests/CANAL/EXPREF/field_def_nemo-oce.xml

    r9572 r12205  
    110110         <field id="tpt_dep"      long_name="T-point depth"                  standard_name="depth_below_geoid"   unit="m"   grid_ref="grid_T_3D" /> 
    111111         <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"      />   
    190112      </field_group> 
    191113 
  • NEMO/branches/2019/dev_r11943_MERGE_2019/tests/CANAL/EXPREF/namelist_cfg

    r11536 r12205  
    276276!!   namdiu       Cool skin and warm layer models                       (default: OFF) 
    277277!!   namflo       float parameters                                      (default: OFF) 
    278 !!   nam_diaharm  Harmonic analysis of tidal constituents               (default: OFF) 
    279278!!   nam_diadct   transports through some sections                      (default: OFF) 
    280279!!   nam_diatmb   Top Middle Bottom Output                              (default: OFF) 
  • NEMO/branches/2019/dev_r11943_MERGE_2019/tests/ICE_ADV1D/EXPREF/namelist_cfg

    r11536 r12205  
    202202!!   namdiu       Cool skin and warm layer models                       (default: OFF) 
    203203!!   namflo       float parameters                                      (default: OFF) 
    204 !!   nam_diaharm  Harmonic analysis of tidal constituents               (default: OFF) 
    205204!!   nam_diadct   transports through some sections                      (default: OFF) 
    206205!!   nam_diatmb   Top Middle Bottom Output                              (default: OFF) 
  • NEMO/branches/2019/dev_r11943_MERGE_2019/tests/ICE_ADV1D/EXPREF/namelist_cfg_120pts

    r11536 r12205  
    202202!!   namdiu       Cool skin and warm layer models                       (default: OFF) 
    203203!!   namflo       float parameters                                      (default: OFF) 
    204 !!   nam_diaharm  Harmonic analysis of tidal constituents               (default: OFF) 
    205204!!   nam_diadct   transports through some sections                      (default: OFF) 
    206205!!   nam_diatmb   Top Middle Bottom Output                              (default: OFF) 
  • NEMO/branches/2019/dev_r11943_MERGE_2019/tests/ICE_ADV1D/EXPREF/namelist_cfg_240pts

    r11536 r12205  
    202202!!   namdiu       Cool skin and warm layer models                       (default: OFF) 
    203203!!   namflo       float parameters                                      (default: OFF) 
    204 !!   nam_diaharm  Harmonic analysis of tidal constituents               (default: OFF) 
    205204!!   nam_diadct   transports through some sections                      (default: OFF) 
    206205!!   nam_diatmb   Top Middle Bottom Output                              (default: OFF) 
  • NEMO/branches/2019/dev_r11943_MERGE_2019/tests/ICE_ADV1D/EXPREF/namelist_cfg_60pts

    r11536 r12205  
    202202!!   namdiu       Cool skin and warm layer models                       (default: OFF) 
    203203!!   namflo       float parameters                                      (default: OFF) 
    204 !!   nam_diaharm  Harmonic analysis of tidal constituents               (default: OFF) 
    205204!!   nam_diadct   transports through some sections                      (default: OFF) 
    206205!!   nam_diatmb   Top Middle Bottom Output                              (default: OFF) 
  • NEMO/branches/2019/dev_r11943_MERGE_2019/tests/ICE_ADV2D/EXPREF/namelist_cfg

    r11536 r12205  
    202202!!   namdiu       Cool skin and warm layer models                       (default: OFF) 
    203203!!   namflo       float parameters                                      (default: OFF) 
    204 !!   nam_diaharm  Harmonic analysis of tidal constituents               (default: OFF) 
    205204!!   nam_diadct   transports through some sections                      (default: OFF) 
    206205!!   nam_diatmb   Top Middle Bottom Output                              (default: OFF) 
  • NEMO/branches/2019/dev_r11943_MERGE_2019/tests/ICE_AGRIF/EXPREF/1_namelist_cfg

    r11536 r12205  
    202202!!   namdiu       Cool skin and warm layer models                       (default: OFF) 
    203203!!   namflo       float parameters                                      (default: OFF) 
    204 !!   nam_diaharm  Harmonic analysis of tidal constituents               (default: OFF) 
    205204!!   nam_diadct   transports through some sections                      (default: OFF) 
    206205!!   nam_diatmb   Top Middle Bottom Output                              (default: OFF) 
  • NEMO/branches/2019/dev_r11943_MERGE_2019/tests/ICE_AGRIF/EXPREF/namelist_cfg

    r11536 r12205  
    202202!!   namdiu       Cool skin and warm layer models                       (default: OFF) 
    203203!!   namflo       float parameters                                      (default: OFF) 
    204 !!   nam_diaharm  Harmonic analysis of tidal constituents               (default: OFF) 
    205204!!   nam_diadct   transports through some sections                      (default: OFF) 
    206205!!   nam_diatmb   Top Middle Bottom Output                              (default: OFF) 
  • NEMO/branches/2019/dev_r11943_MERGE_2019/tests/ISOMIP/EXPREF/namelist_cfg

    r12150 r12205  
    428428!!   namdiu       Cool skin and warm layer models                       (default: OFF) 
    429429!!   namflo       float parameters                                      (default: OFF) 
    430 !!   nam_diaharm  Harmonic analysis of tidal constituents               (default: OFF) 
    431430!!   nam_diadct   transports through some sections                      (default: OFF) 
    432431!!   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  
    276276!!   namdiu       Cool skin and warm layer models                       (default: OFF) 
    277277!!   namflo       float parameters                                      (default: OFF) 
    278 !!   nam_diaharm  Harmonic analysis of tidal constituents               (default: OFF) 
    279278!!   nam_diadct   transports through some sections                      (default: OFF) 
    280279!!   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  
    294294!!   namdiu       Cool skin and warm layer models                       (default: OFF) 
    295295!!   namflo       float parameters                                      (default: OFF) 
    296 !!   nam_diaharm  Harmonic analysis of tidal constituents               (default: OFF) 
    297296!!   nam_diadct   transports through some sections                      (default: OFF) 
    298297!!   nam_diatmb   Top Middle Bottom Output                              (default: OFF) 
  • NEMO/branches/2019/dev_r11943_MERGE_2019/tests/VORTEX/EXPREF/1_namelist_cfg

    r11536 r12205  
    273273!!   namdiu       Cool skin and warm layer models                       (default: OFF) 
    274274!!   namflo       float parameters                                      (default: OFF) 
    275 !!   nam_diaharm  Harmonic analysis of tidal constituents               (default: OFF) 
    276275!!   nam_diadct   transports through some sections                      (default: OFF) 
    277276!!   nam_diatmb   Top Middle Bottom Output                              (default: OFF) 
  • NEMO/branches/2019/dev_r11943_MERGE_2019/tests/VORTEX/EXPREF/namelist_cfg

    r11536 r12205  
    264264!!   namdiu       Cool skin and warm layer models                       (default: OFF) 
    265265!!   namflo       float parameters                                      (default: OFF) 
    266 !!   nam_diaharm  Harmonic analysis of tidal constituents               (default; OFF) 
    267266!!   nam_diadct   transports through some sections                      (default: OFF) 
    268267!!   nam_diatmb   Top Middle Bottom Output                              (default: OFF) 
  • NEMO/branches/2019/dev_r11943_MERGE_2019/tests/WAD/EXPREF/namelist_cfg

    r11949 r12205  
    422422!!   namdiu       Cool skin and warm layer models                       (default: OFF) 
    423423!!   namflo       float parameters                                      (default: OFF) 
    424 !!   nam_diaharm  Harmonic analysis of tidal constituents               (default: OFF) 
    425424!!   nam_diadct   transports through some sections                      (default: OFF) 
    426425!!   nam_diatmb   Top Middle Bottom Output                              (default: OFF) 
Note: See TracChangeset for help on using the changeset viewer.