The concurrent simulation of ocean circulations and ocean tides is carried out with the Max-Planck Institute Ocean Model (MPIOM/TP6M L40 mpiom-1.6.3.) forced by the full luni-solar tidal potential as an additional body force and by the surface fluxes of momentum, heat and freshwater derived the NCEP/NCAR reanalysis for the period 1981-2012. A tripolar grid with a horizontal resolution of about 0.1 degrees is used. There are total 40 vertical levels in z-coordinates. More details can be find in Li and von Storch (2020). The file name of the data sets is composed as follows. STORMTIDE2_TP6ML40_<variable_acronym>_3d_1hr_<date>_<run>.nc There are five 3-dimensional hourly variables: sea_water_potential_temperature (tho), sea_water_salinity (sao), sea_water_x_velocity (u), sea_water_y_velocity (v) and upward_sea_water_velocity (w), for January, April, July and October of 2012. 2012 is the last year of the simulation. Storing four months hourly was doable at the time when the simulation is produced. One run is provided.

For the development of the earth-system model ICON, carried out at the Max-Planck Institute for Meteorology in cooperation with the German meteorological service (DWD), the ocean compartment plays a major role. It is now used to push the frontiers of simulating the climate based on first principles by resolving its major energetic motions, including those at small scales such as eddies and waves in the ocean. This new development, however, still requires substantial efforts. One such a process is internal tides that are generated as (barotropic) tides flow over topographic features such as underwater ridges and sea mounts. To ensure a realistic representation of internal tides, one needs to make sure that the model used is capable to realistically simulate the tides. A set of experiments has been carried to assess the ability of ICON-O to simulate open ocean tides. Tides are tested in the global ocean model ICON-O in various horizontal resolutions: R2B6 (40km), R2B8 (10km), and BCT (BaseCampTelescope: a telescoping -inhomogeneous- grid ranging from 80km to 8km with a focus region in the South Atlantic). All runs are performed with a relatively high vertical resolution (128 vertical levels) and additionally the vertical coordinate z* is used and different heights of top levels are applied. To improve the simulated tides two new parameterizations are introduced in ICON-O. Self-Attraction and Loading (SAL) and a parameterization for Tidal-Bottom-Drag (TBD) are introduced individually and together in all 3 horizontal resolutions. 100 years of spin-up simulations with OMIP forcing and without tides have been performed for all set-ups. After the climate has spun-up reasonably, 3 more years of simulation now with tides have been performed. For the tides to spin-up, only the last year is evaluated. The raw model output (hourly sea surface elevation) and the spin-up runs are also available (see references). One year (the last) of the ICON-O simulation with tides is analyzed by the least-square-fitting method of Foreman et al. (2009). From this harmonic analysis the harmonic constants of the eight major diurnal and semi-diurnal tidal constituents (M2, S2, N2, K2, K1, O1, P1, and Q1) are derived. Phases and amplitudes on the original ICON-O grid (XX_harmonics.nc) and a regridded regular 0.1-degree grid (XX_harmonics_r3600x1800.nc) are archived.