Changeset 15422 for NEMO/branches/UKMO/r14075_India_uncoupled/src/OCE/SBC/sbctide.F90

Timestamp:

2021-10-21T11:19:25+02:00 (3 years ago)

Author:

jcastill

Message:

Changes tested so that they can merged with the CO9 Met Office branch - jpmax_harmo should be 34 with FES14 tides, but the last components are not used anyway

File:

• NEMO/branches/UKMO/r14075_India_uncoupled/src/OCE/SBC/sbctide.F90 (modified) (3 diffs)

NEMO/branches/UKMO/r14075_India_uncoupled/src/OCE/SBC/sbctide.F90

@@ -16 +16 @@
    USE ioipsl         ! NetCDF IPSL library
    USE lbclnk         ! ocean lateral boundary conditions (or mpp link)
+   USE bdytides
 
    IMPLICIT NONE
@@ -108 +109 @@
 
       DO jk = 1, nb_harmo
-         zcons = 0.7_wp * Wave(ntide(jk))%equitide * ftide(jk)
+         ! love number now provides in tide namelist  
+         zcons = dn_love_number * Wave(ntide(jk))%equitide * ftide(jk)
          DO ji = 1, jpi
             DO jj = 1, jpj
@@ -119 +121 @@
                IF    ( Wave(ntide(jk))%nutide == 1 )  THEN  ;  zcs = zcons * SIN( 2._wp*zlat )
                ELSEIF( Wave(ntide(jk))%nutide == 2 )  THEN  ;  zcs = zcons * COS( zlat )**2
+               ! Add tide potential for long period tides  
+               ELSEIF( Wave(ntide(jk))%nutide == 0 )  THEN  ;  zcs = zcons * (0.5_wp-1.5_wp*SIN(zlat)**2._wp)
                ELSE                                         ;  zcs = 0._wp
                ENDIF