New URL for NEMO forge!   http://forge.nemo-ocean.eu

Since March 2022 along with NEMO 4.2 release, the code development moved to a self-hosted GitLab.
This present forge is now archived and remained online for history.
tideini.F90 in branches/UKMO/r8395_India_uncoupled/NEMOGCM/NEMO/OPA_SRC/SBC – NEMO

source: branches/UKMO/r8395_India_uncoupled/NEMOGCM/NEMO/OPA_SRC/SBC/tideini.F90 @ 10685

Last change on this file since 10685 was 10685, checked in by jcastill, 5 years ago

Changes as Ash's files
File size: 5.0 KB
Line 
1MODULE tideini
2   !!======================================================================
3   !!                       ***  MODULE  tideini  ***
4   !! Initialization of tidal forcing
5   !!======================================================================
6   !! History :  1.0  !  2007  (O. Le Galloudec)  Original code
7   !!----------------------------------------------------------------------
8   USE oce            ! ocean dynamics and tracers variables
9   USE dom_oce        ! ocean space and time domain
10   USE phycst         ! physical constant
11   USE daymod         ! calandar
12   USE tide_mod       !
13   !
14   USE in_out_manager ! I/O units
15   USE iom            ! xIOs server
16   USE ioipsl         ! NetCDF IPSL library
17   USE lbclnk         ! ocean lateral boundary conditions (or mpp link)
18
19   IMPLICIT NONE
20   PUBLIC
21
22   REAL(wp), PUBLIC, ALLOCATABLE, DIMENSION(:) ::   omega_tide   !:
23   REAL(wp), PUBLIC, ALLOCATABLE, DIMENSION(:) ::   v0tide       !:
24   REAL(wp), PUBLIC, ALLOCATABLE, DIMENSION(:) ::   utide        !:
25   REAL(wp), PUBLIC, ALLOCATABLE, DIMENSION(:) ::   ftide        !:
26
27   LOGICAL , PUBLIC ::   ln_tide         !:
28   LOGICAL , PUBLIC ::   ln_tide_pot     !:
29   LOGICAL , PUBLIC ::   ln_tide_ramp    !:
30   INTEGER , PUBLIC ::   nb_harmo        !:
31   INTEGER , PUBLIC ::   kt_tide         !:
32   REAL(wp), PUBLIC ::   rdttideramp     !:
33   REAL(wp), PUBLIC ::   dn_love_number  !:
34   INTEGER , PUBLIC, ALLOCATABLE, DIMENSION(:) ::   ntide   !:
35
36   !!----------------------------------------------------------------------
37   !! NEMO/OPA 3.5 , NEMO Consortium (2013)
38   !! $Id$
39   !! Software governed by the CeCILL licence     (NEMOGCM/NEMO_CeCILL.txt)
40   !!----------------------------------------------------------------------
41CONTAINS
42   
43   SUBROUTINE tide_init
44      !!----------------------------------------------------------------------
45      !!                 ***  ROUTINE tide_init  ***
46      !!----------------------------------------------------------------------     
47      INTEGER  :: ji, jk
48      CHARACTER(LEN=4), DIMENSION(jpmax_harmo) :: clname
49      INTEGER  ::   ios                 ! Local integer output status for namelist read
50      !
51      NAMELIST/nam_tide/ln_tide, ln_tide_pot, ln_tide_ramp, rdttideramp, dn_love_number, clname
52      !!----------------------------------------------------------------------
53      !
54      ! Read Namelist nam_tide
55      REWIND( numnam_ref )              ! Namelist nam_tide in reference namelist : Tides
56      READ  ( numnam_ref, nam_tide, IOSTAT = ios, ERR = 901)
57901   IF( ios /= 0 )   CALL ctl_nam ( ios , 'nam_tide in reference namelist', lwp )
58      !
59      REWIND( numnam_cfg )              ! Namelist nam_tide in configuration namelist : Tides
60      READ  ( numnam_cfg, nam_tide, IOSTAT = ios, ERR = 902 )
61902   IF( ios /= 0 )   CALL ctl_nam ( ios , 'nam_tide in configuration namelist', lwp )
62      IF(lwm) WRITE ( numond, nam_tide )
63      !
64      IF (ln_tide) THEN
65         IF (lwp) THEN
66            WRITE(numout,*)
67            WRITE(numout,*) 'tide_init : Initialization of the tidal components'
68            WRITE(numout,*) '~~~~~~~~~ '
69            WRITE(numout,*) '   Namelist nam_tide'
70            WRITE(numout,*) '              Use tidal components : ln_tide      = ', ln_tide
71            WRITE(numout,*) '      Apply astronomical potential : ln_tide_pot  = ', ln_tide_pot
72            WRITE(numout,*) '                                     ln_tide_ramp = ', ln_tide_ramp
73            WRITE(numout,*) '                                     dn_love_number = ', dn_love_number
74         ENDIF
75      ELSE
76         IF(lwp) WRITE(numout,*)
77         IF(lwp) WRITE(numout,*) 'tide_init : tidal components not used (ln_tide = F)'
78         IF(lwp) WRITE(numout,*) '~~~~~~~~~ '
79         RETURN
80      ENDIF
81      !
82      CALL tide_init_Wave
83      !
84      nb_harmo=0
85      DO jk = 1, jpmax_harmo
86         DO ji = 1,jpmax_harmo
87            IF( TRIM(clname(jk)) == Wave(ji)%cname_tide )   nb_harmo = nb_harmo + 1
88         END DO
89      END DO
90      IF (ln_tide .and.lwp) WRITE(numout,*) '                                     nb_harmo     = ', nb_harmo
91
92      ! Ensure that tidal components have been set in namelist_cfg
93      IF( nb_harmo == 0 )   CALL ctl_stop( 'tide_init : No tidal components set in nam_tide' )
94      !
95      IF( ln_tide_ramp.AND.((nitend-nit000+1)*rdt/rday < rdttideramp) )   &
96         &   CALL ctl_stop('rdttideramp must be lower than run duration')
97      IF( ln_tide_ramp.AND.(rdttideramp<0.) ) &
98         &   CALL ctl_stop('rdttideramp must be positive')
99      !
100      ALLOCATE( ntide(nb_harmo) )
101      DO jk = 1, nb_harmo
102         DO ji = 1, jpmax_harmo
103            IF( TRIM(clname(jk)) == Wave(ji)%cname_tide ) THEN
104               ntide(jk) = ji
105               EXIT
106            ENDIF
107         END DO
108      END DO
109      !
110      ALLOCATE( omega_tide(nb_harmo), v0tide    (nb_harmo),   &
111         &      utide     (nb_harmo), ftide     (nb_harmo)  )
112      kt_tide = nit000
113      !
114   END SUBROUTINE tide_init
115     
116   !!======================================================================
117END MODULE tideini
Note: See TracBrowser for help on using the repository browser.