Changeset 12453 for NEMO/branches/UKMO/r12083_India_uncoupled/src/OCE/SBC

Timestamp:

2020-02-25T16:29:34+01:00 (4 years ago)

Author:

jcastill

Message:

First implementation of the branch - compiling after merge

Location:

NEMO/branches/UKMO/r12083_India_uncoupled/src/OCE/SBC

Files:

• sbctide.F90 (modified) (3 diffs)
• tide_FES14.h90 (added)
• tide_mod.F90 (modified) (3 diffs)
• tideini.F90 (modified) (4 diffs)

NEMO/branches/UKMO/r12083_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

NEMO/branches/UKMO/r12083_India_uncoupled/src/OCE/SBC/tide_mod.F90

@@ -16 +16 @@
    PUBLIC   tide_init_Wave   ! called by tideini and diaharm modules
 
+# if defined key_FES14_tides 
+   INTEGER, PUBLIC, PARAMETER ::   jpmax_harmo = 34   !: maximum number of harmonic 
+# else
    INTEGER, PUBLIC, PARAMETER ::   jpmax_harmo = 19   !: maximum number of harmonic
+# endif
 
    TYPE, PUBLIC ::    tide
@@ -41 +45 @@
 
    SUBROUTINE tide_init_Wave
+# if defined key_FES14_tides 
+#     include "tide_FES14.h90" 
+# else
 #     include "tide.h90"
+# endif
    END SUBROUTINE tide_init_Wave
 
@@ -331 +339 @@
          zf = zf * zf1 * zf1
          !
+      CASE( 20 )                 !==  formule 20,  compound waves ( 78 x 78 x 78 x 78 ) 
+         zf1 = nodal_factort(78) 
+         zf  = zf1 * zf1 * zf1 * zf1 
+         !
       CASE( 73 )                 !==  formule 73
          zs = sin(sh_I)

NEMO/branches/UKMO/r12083_India_uncoupled/src/OCE/SBC/tideini.F90

@@ -33 +33 @@
    INTEGER , PUBLIC ::   kt_tide         !:
    REAL(wp), PUBLIC ::   rdttideramp     !:
+   REAL(wp), PUBLIC ::   dn_love_number  !:
    REAL(wp), PUBLIC ::   rn_scal_load    !:
    CHARACTER(lc), PUBLIC ::   cn_tide_load   !: 
@@ -54 +55 @@
       !
       NAMELIST/nam_tide/ln_tide, ln_tide_pot, ln_scal_load, ln_read_load, cn_tide_load, &
-                  &     ln_tide_ramp, rn_scal_load, rdttideramp, clname
+                  &     ln_tide_ramp, rn_scal_load, rdttideramp, dn_love_number, clname
       !!----------------------------------------------------------------------
       !
@@ -78 +79 @@
             WRITE(numout,*) '         Read load potential from file           ln_read_load = ', ln_read_load
             WRITE(numout,*) '         Apply ramp on tides at startup          ln_tide_ramp = ', ln_tide_ramp
+            WRITE(numout,*) '                                              dn_love_number  = ', dn_love_number
             WRITE(numout,*) '         Fraction of SSH used in scal. approx.   rn_scal_load = ', rn_scal_load
             WRITE(numout,*) '         Duration (days) of ramp                 rdttideramp  = ', rdttideramp
@@ -99 +101 @@
       END DO
       !       
+      IF (ln_tide .and.lwp) WRITE(numout,*) '                                     nb_harmo     = ', nb_harmo
+
       ! Ensure that tidal components have been set in namelist_cfg
       IF( nb_harmo == 0 )   CALL ctl_stop( 'tide_init : No tidal components set in nam_tide' )