source: branches/UKMO/r6232_tracer_advection/NEMOGCM/NEMO/OPA_SRC/SBC/sbctide.F90 @ 9295

Last change on this file since 9295 was 9295, checked in by jcastill, 3 years ago

Remove svn keywords

File size: 5.3 KB
Line 
1MODULE sbctide
2   !!======================================================================
3   !!                       ***  MODULE  sbctide  ***
4   !! Initialization of tidal forcing
5   !!======================================================================
6   !! History :  9.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
11   USE daymod
12   USE dynspg_oce
13   USE tideini
14   !
15   USE iom
16   USE in_out_manager  ! I/O units
17   USE ioipsl          ! NetCDF IPSL library
18   USE lbclnk          ! ocean lateral boundary conditions (or mpp link)
19
20   IMPLICIT NONE
21   PUBLIC
22
23   REAL(wp), PUBLIC, ALLOCATABLE, DIMENSION(:,:) ::   pot_astro   !
24
25#if defined key_tide
26   !!----------------------------------------------------------------------
27   !!   'key_tide' :                                        tidal potential
28   !!----------------------------------------------------------------------
29   !!   sbc_tide            :
30   !!   tide_init_potential :
31   !!----------------------------------------------------------------------
32
33   LOGICAL, PUBLIC, PARAMETER ::   lk_tide  = .TRUE.
34   REAL(wp), ALLOCATABLE, DIMENSION(:,:,:) ::   amp_pot, phi_pot
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 sbc_tide( kt )
44      !!----------------------------------------------------------------------
45      !!                 ***  ROUTINE sbc_tide  ***
46      !!----------------------------------------------------------------------     
47      INTEGER, INTENT( in ) ::   kt     ! ocean time-step
48      INTEGER               ::   jk     ! dummy loop index
49      !!----------------------------------------------------------------------
50
51      IF( nsec_day == NINT(0.5_wp * rdttra(1)) ) THEN      ! start a new day
52         !
53         IF( kt == nit000 ) THEN
54            ALLOCATE( amp_pot(jpi,jpj,nb_harmo),                      &
55               &      phi_pot(jpi,jpj,nb_harmo), pot_astro(jpi,jpj)   )
56         ENDIF
57         !
58         amp_pot(:,:,:) = 0._wp
59         phi_pot(:,:,:) = 0._wp
60         pot_astro(:,:) = 0._wp
61         !
62         CALL tide_harmo( omega_tide, v0tide, utide, ftide, ntide, nb_harmo )
63         !
64         kt_tide = kt
65         !
66         IF(lwp) THEN
67            WRITE(numout,*)
68            WRITE(numout,*) 'sbc_tide : Update of the components and (re)Init. the potential at kt=', kt
69            WRITE(numout,*) '~~~~~~~~ '
70            DO jk = 1, nb_harmo
71               WRITE(numout,*) Wave(ntide(jk))%cname_tide, utide(jk), ftide(jk), v0tide(jk), omega_tide(jk)
72            END DO
73         ENDIF
74         !
75         IF( ln_tide_pot )   CALL tide_init_potential
76         !
77      ENDIF
78      !
79   END SUBROUTINE sbc_tide
80
81
82   SUBROUTINE tide_init_potential
83      !!----------------------------------------------------------------------
84      !!                 ***  ROUTINE tide_init_potential  ***
85      !!----------------------------------------------------------------------
86      INTEGER  ::   ji, jj, jk   ! dummy loop indices
87      REAL(wp) ::   zcons, ztmp1, ztmp2, zlat, zlon, ztmp, zamp, zcs   ! local scalar
88      !!----------------------------------------------------------------------
89
90      DO jk = 1, nb_harmo
91         zcons = 0.7_wp * Wave(ntide(jk))%equitide * ftide(jk)
92         DO ji = 1, jpi
93            DO jj = 1, jpj
94               ztmp1 =  amp_pot(ji,jj,jk) * COS( phi_pot(ji,jj,jk) )
95               ztmp2 = -amp_pot(ji,jj,jk) * SIN( phi_pot(ji,jj,jk) )
96               zlat = gphit(ji,jj)*rad !! latitude en radian
97               zlon = glamt(ji,jj)*rad !! longitude en radian
98               ztmp = v0tide(jk) + utide(jk) + Wave(ntide(jk))%nutide * zlon
99               ! le potentiel est composé des effets des astres:
100               IF    ( Wave(ntide(jk))%nutide == 1 )  THEN  ;  zcs = zcons * SIN( 2._wp*zlat )
101               ELSEIF( Wave(ntide(jk))%nutide == 2 )  THEN  ;  zcs = zcons * COS( zlat )**2
102               ELSE                                         ;  zcs = 0._wp
103               ENDIF
104               ztmp1 = ztmp1 + zcs * COS( ztmp )
105               ztmp2 = ztmp2 - zcs * SIN( ztmp )
106               zamp = SQRT( ztmp1*ztmp1 + ztmp2*ztmp2 )
107               amp_pot(ji,jj,jk) = zamp
108               phi_pot(ji,jj,jk) = ATAN2( -ztmp2 / MAX( 1.e-10_wp , zamp ) ,   &
109                  &                        ztmp1 / MAX( 1.e-10_wp,  zamp )   )
110            END DO
111         END DO
112      END DO
113      !
114   END SUBROUTINE tide_init_potential
115
116#else
117  !!----------------------------------------------------------------------
118  !!   Default case :   Empty module
119  !!----------------------------------------------------------------------
120  LOGICAL, PUBLIC, PARAMETER ::   lk_tide = .FALSE.
121CONTAINS
122  SUBROUTINE sbc_tide( kt )      ! Empty routine
123    INTEGER         , INTENT(in) ::   kt         ! ocean time-step
124    WRITE(*,*) 'sbc_tide: You should not have seen this print! error?', kt
125  END SUBROUTINE sbc_tide
126#endif
127
128  !!======================================================================
129END MODULE sbctide
Note: See TracBrowser for help on using the repository browser.