Changes between Version 14 and Version 15 of 2019WP/ENHANCE-05_SimonM-Harmonic_Analysis
- Timestamp:
- 2019-11-06T14:59:49+01:00 (4 years ago)
Legend:
- Unmodified
- Added
- Removed
- Modified
-
2019WP/ENHANCE-05_SimonM-Harmonic_Analysis
v14 v15 23 23 === Current implementation of tidal harmonic-analysis diagnostics and available enhancements 24 24 25 Tidal harmonic analysis of model fields based on multiple linear least-squares regression is available in the current trunk version of NEMO in module `diaharm` (source:/NEMO/trunk/src/OCE/DIA/diaharm.F90@10835). It is restricted to the analysis of sea-surface height and barotropic velocity fields, to a hard-coded maximum number of time slices within uninterrupted model runs, and its memory footprint during the whole analysis period can be very large. The current implementation has aspects that are deprecated (use of a pre-processor key, non-conventional namelist variable names, and statical allocation of potentially oversized arrays) and appears to be computationally inefficient. Further, it lacks sought-after features, such as the analysis of a wider range of model fields (including three-dimensional fields) and the possibility to span the analysis time period across restarted model runs, which have become available in more efficient alternative implementations of the same analysis method by N. Bruneau (pers. comm.) and E. O'Dea (https://forge.ipsl.jussieu.fr/nemo/browser/branches/UKMO/dev_5518_tid _analysis_restart) in previous NEMO versions. In addition, the harmonic analysis diagnostics of module `diaharm` remains deactivated, and thus its compilability untested, during the standard tests of the SETTE test suite.25 Tidal harmonic analysis of model fields based on multiple linear least-squares regression is available in the current trunk version of NEMO in module `diaharm` (source:/NEMO/trunk/src/OCE/DIA/diaharm.F90@10835). It is restricted to the analysis of sea-surface height and barotropic velocity fields, to a hard-coded maximum number of time slices within uninterrupted model runs, and its memory footprint during the whole analysis period can be very large. The current implementation has aspects that are deprecated (use of a pre-processor key, non-conventional namelist variable names, and statical allocation of potentially oversized arrays) and appears to be computationally inefficient. Further, it lacks sought-after features, such as the analysis of a wider range of model fields (including three-dimensional fields) and the possibility to span the analysis time period across restarted model runs, which have become available in more efficient alternative implementations of the same analysis method by N. Bruneau (pers. comm.) and E. O'Dea (https://forge.ipsl.jussieu.fr/nemo/browser/branches/UKMO/dev_5518_tide_analysis_restart) in previous NEMO versions. In addition, the harmonic analysis diagnostics of module `diaharm` remains deactivated, and thus its compilability untested, during the standard tests of the SETTE test suite. 26 26 27 27 === Proposed replacement of the tidal harmonic-analysis diagnostics by a generic formulation for multiple linear regression analysis