source: branches/2013/dev_r3858_CNRS3_Ediag/NEMOGCM/NEMO/OPA_SRC/DYN/dynspg_exp.F90 @ 3876

Last change on this file since 3876 was 3876, checked in by gm, 8 years ago

dev_r3858_CNRS3_Ediag: #927 phasing with 2011/dev_r3309_LOCEAN12_Ediag branche + mxl diag update

  • Property svn:keywords set to Id
File size: 5.5 KB
Line 
1MODULE dynspg_exp
2   !!======================================================================
3   !!                   ***  MODULE  dynspg_exp  ***
4   !! Ocean dynamics:  surface pressure gradient trend
5   !!======================================================================
6   !! History :  2.0  !  2005-11  (V. Garnier, G. Madec, L. Bessieres) Original code
7   !!            3.2  !  2009-06  (G. Madec, M. Leclair, R. Benshila) introduce sshwzv module
8   !!----------------------------------------------------------------------
9#if defined key_dynspg_exp   ||   defined key_esopa
10   !!----------------------------------------------------------------------
11   !!   'key_dynspg_exp'                              explicit free surface
12   !!----------------------------------------------------------------------
13   !!   dyn_spg_exp  : update the momentum trend with the surface
14   !!                      pressure gradient in the free surface constant 
15   !!                      volume case with vector optimization
16   !!----------------------------------------------------------------------
17   USE oce             ! ocean dynamics and tracers
18   USE dom_oce         ! ocean space and time domain
19   USE sbc_oce         ! surface boundary condition: ocean
20   USE obc_oce         ! Lateral open boundary condition
21   USE phycst          ! physical constants
22   USE obc_par         ! open boundary condition parameters
23   USE obcdta          ! open boundary condition data     (bdy_dta_bt routine)
24   !
25   USE in_out_manager  ! I/O manager
26   USE lib_mpp         ! distributed memory computing library
27   USE lbclnk          ! ocean lateral boundary conditions (or mpp link)
28   USE prtctl          ! Print control
29   USE iom             ! I/O library
30   USE timing          ! Timing
31
32
33   IMPLICIT NONE
34   PRIVATE
35
36   PUBLIC   dyn_spg_exp   ! routine called by step.F90
37
38   !! * Substitutions
39#  include "domzgr_substitute.h90"
40#  include "vectopt_loop_substitute.h90"
41   !!----------------------------------------------------------------------
42   !! NEMO/OPA 3.3 , NEMO Consortium (2010)
43   !! $Id$
44   !! Software governed by the CeCILL licence     (NEMOGCM/NEMO_CeCILL.txt)
45   !!----------------------------------------------------------------------
46CONTAINS
47
48   SUBROUTINE dyn_spg_exp( kt )
49      !!----------------------------------------------------------------------
50      !!                  ***  routine dyn_spg_exp  ***
51      !!
52      !! ** Purpose :   Compute the now trend due to the surface pressure
53      !!              gradient in case of explicit free surface formulation and
54      !!              add it to the general trend of momentum equation.
55      !!
56      !! ** Method  :   Explicit free surface formulation. Add to the general
57      !!              momentum trend the surface pressure gradient :
58      !!                      (ua,va) = (ua,va) + (spgu,spgv)
59      !!              where spgu = -1/rau0 d/dx(ps) = -g/e1u di( sshn )
60      !!                    spgv = -1/rau0 d/dy(ps) = -g/e2v dj( sshn )
61      !!
62      !! ** Action :   (ua,va)   trend of horizontal velocity increased by
63      !!                         the surf. pressure gradient trend
64      !!---------------------------------------------------------------------
65      INTEGER, INTENT(in)  ::   kt   ! ocean time-step index
66      !!
67      INTEGER ::   ji, jj, jk   ! dummy loop indices
68      !!----------------------------------------------------------------------
69      !
70      IF( nn_timing == 1 )  CALL timing_start('dyn_spg_exp')
71      !
72      IF( kt == nit000 ) THEN
73         IF(lwp) WRITE(numout,*)
74         IF(lwp) WRITE(numout,*) 'dyn_spg_exp : surface pressure gradient trend'
75         IF(lwp) WRITE(numout,*) '~~~~~~~~~~~   (explicit free surface)'
76         !
77         spgu(:,:) = 0._wp   ;   spgv(:,:) = 0._wp
78         !
79         IF( lk_vvl .AND. lwp ) WRITE(numout,*) '              lk_vvl=T : spg is included in dynhpg'
80      ENDIF
81
82
83!!gm bug ??  Rachid we have to discuss of the call below. I don't understand why it is here and not in ssh_wzv
84      IF( lk_obc )   CALL obc_dta_bt( kt, 0 )      ! OBC: read or estimate ssh and vertically integrated velocities
85!!gm
86
87      IF( .NOT. lk_vvl ) THEN          !* fixed volume : add the surface pressure gradient trend
88         !
89         DO jj = 2, jpjm1                    ! now surface pressure gradient
90            DO ji = fs_2, fs_jpim1   ! vector opt.
91               spgu(ji,jj) = - grav * ( sshn(ji+1,jj) - sshn(ji,jj) ) / e1u(ji,jj)
92               spgv(ji,jj) = - grav * ( sshn(ji,jj+1) - sshn(ji,jj) ) / e2v(ji,jj)
93            END DO
94         END DO
95         DO jk = 1, jpkm1                    ! Add it to the general trend
96            DO jj = 2, jpjm1
97               DO ji = fs_2, fs_jpim1   ! vector opt.
98                  ua(ji,jj,jk) = ua(ji,jj,jk) + spgu(ji,jj)
99                  va(ji,jj,jk) = va(ji,jj,jk) + spgv(ji,jj)
100               END DO
101            END DO
102         END DO
103         !
104      ENDIF
105      !
106      IF( nn_timing == 1 )  CALL timing_stop('dyn_spg_exp')
107      !
108   END SUBROUTINE dyn_spg_exp
109
110#else
111   !!----------------------------------------------------------------------
112   !!   Default case :   Empty module   No standart explicit free surface
113   !!----------------------------------------------------------------------
114CONTAINS
115   SUBROUTINE dyn_spg_exp( kt )       ! Empty routine
116      WRITE(*,*) 'dyn_spg_exp: You should not have seen this print! error?', kt
117   END SUBROUTINE dyn_spg_exp
118#endif
119
120   !!======================================================================
121END MODULE dynspg_exp
Note: See TracBrowser for help on using the repository browser.