Changeset 14166 for NEMO/branches/NERC/NEMO_4.0.2_CO9_package_tides/src/OCE/SBC/sbctide.F90

Timestamp:

2020-12-13T20:36:10+01:00 (4 years ago)

Author:

dbyrne

Message:

Long period tide forcing, variable love number, new nodal equation 20 added

File:

• NEMO/branches/NERC/NEMO_4.0.2_CO9_package_tides/src/OCE/SBC/sbctide.F90 (modified) (3 diffs)

NEMO/branches/NERC/NEMO_4.0.2_CO9_package_tides/src/OCE/SBC/sbctide.F90

@@ -16 +16 @@
    USE ioipsl         ! NetCDF IPSL library
    USE lbclnk         ! ocean lateral boundary conditions (or mpp link)
+
+   USE bdytides ! davbyr - Access to love number
 
    IMPLICIT NONE
@@ -108 +110 @@
 
       DO jk = 1, nb_harmo
-         zcons = 0.7_wp * Wave(ntide(jk))%equitide * ftide(jk)
+         ! davbyr - Insert variable Love number where once was 0.7
+         zcons = dn_love_number * Wave(ntide(jk))%equitide * ftide(jk)
+         ! END davbyr
          DO ji = 1, jpi
             DO jj = 1, jpj
@@ -119 +123 @@
                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
+               ! davbyr - Include long period tidal forcing
+               ELSEIF( Wave(ntide(jk))%nutide == 0 )  THEN  ;  zcs = zcons * (0.5_wp-1.5_wp*SIN(zlat)**2._wp)
+               ! END - davbyr
                ELSE                                         ;  zcs = 0._wp
                ENDIF