Inline Harmonic Analysis

Update 6: The frequencies of 6 overtides/compound tides (MK3, M4,
MS4, S4, M6, M8) have been added.

Author: c2xu

src/parameterizations/lateral/MOM_tidal_forcing.F90

@@ -150,10 +150,10 @@ end function eq_phase
 !> Looks up angular frequencies for the main tidal constituents.
 !! Values used here are from previous versions of MOM.
 function tidal_frequency(constit)
-  character (len=2), intent(in) :: constit !> Constituent to look up
+  character (len=*), intent(in) :: constit !> Constituent to look up
   real :: tidal_frequency                  !> Angular frequency [rad s-1]
 
-  select case (constit)
+  select case (trim(constit))
     case ("M2")
       tidal_frequency = 1.4051890e-4
     case ("S2")
@@ -174,6 +174,18 @@ function tidal_frequency(constit)
       tidal_frequency = 0.053234e-4
     case ("MM")
       tidal_frequency = 0.026392e-4
+    case ("MK3")
+      tidal_frequency = 2.1344007e-4
+    case ("M4")
+      tidal_frequency = 2.8103780e-4
+    case ("MS4")
+      tidal_frequency = 2.8596300e-4
+    case ("S4")
+      tidal_frequency = 2.9088820e-4
+    case ("M6")
+      tidal_frequency = 4.2155670e-4
+    case ("M8")
+      tidal_frequency = 5.6207560e-4
     case default
       call MOM_error(FATAL, "tidal_frequency: unrecognized constituent")
   end select