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.
dynzad.F90 in branches/2014/dev_CNRS1_10_TEOS10_Ediag/NEMOGCM/NEMO/OPA_SRC/DYN – NEMO

source: branches/2014/dev_CNRS1_10_TEOS10_Ediag/NEMOGCM/NEMO/OPA_SRC/DYN/dynzad.F90 @ 4915

Last change on this file since 4915 was 4619, checked in by gm, 10 years ago

#1294 : TEOS-10 and Ediag

  • Property svn:keywords set to Id
File size: 6.0 KB
RevLine 
[3]1MODULE dynzad
2   !!======================================================================
3   !!                       ***  MODULE  dynzad  ***
4   !! Ocean dynamics : vertical advection trend
5   !!======================================================================
[2715]6   !! History :  OPA  ! 1991-01  (G. Madec) Original code
7   !!            7.0  ! 1991-11  (G. Madec)
8   !!            7.5  ! 1996-01  (G. Madec) statement function for e3
9   !!   NEMO     0.5  ! 2002-07  (G. Madec) Free form, F90
[503]10   !!----------------------------------------------------------------------
[3]11   
12   !!----------------------------------------------------------------------
[503]13   !!   dyn_zad       : vertical advection momentum trend
[3]14   !!----------------------------------------------------------------------
[503]15   USE oce            ! ocean dynamics and tracers
16   USE dom_oce        ! ocean space and time domain
[888]17   USE sbc_oce        ! surface boundary condition: ocean
[4619]18   USE trd_oce        ! trends: ocean variables
19   USE trddyn         ! trend manager: dynamics
20   !
[719]21   USE in_out_manager ! I/O manager
[4619]22   USE lib_mpp        ! MPP library
[503]23   USE prtctl         ! Print control
[4619]24   USE wrk_nemo       ! Memory Allocation
25   USE timing         ! Timing
[3]26
27   IMPLICIT NONE
28   PRIVATE
29   
[503]30   PUBLIC   dyn_zad   ! routine called by step.F90
[3]31
32   !! * Substitutions
33#  include "domzgr_substitute.h90"
34#  include "vectopt_loop_substitute.h90"
35   !!----------------------------------------------------------------------
[2528]36   !! NEMO/OPA 3.3 , NEMO Consortium (2010)
[888]37   !! $Id$
[2715]38   !! Software governed by the CeCILL licence     (NEMOGCM/NEMO_CeCILL.txt)
[3]39   !!----------------------------------------------------------------------
40CONTAINS
41
42   SUBROUTINE dyn_zad ( kt )
43      !!----------------------------------------------------------------------
44      !!                  ***  ROUTINE dynzad  ***
45      !!
46      !! ** Purpose :   Compute the now vertical momentum advection trend and
47      !!      add it to the general trend of momentum equation.
48      !!
49      !! ** Method  :   The now vertical advection of momentum is given by:
50      !!         w dz(u) = ua + 1/(e1u*e2u*e3u) mk+1[ mi(e1t*e2t*wn) dk(un) ]
51      !!         w dz(v) = va + 1/(e1v*e2v*e3v) mk+1[ mj(e1t*e2t*wn) dk(vn) ]
52      !!      Add this trend to the general trend (ua,va):
53      !!         (ua,va) = (ua,va) + w dz(u,v)
54      !!
55      !! ** Action  : - Update (ua,va) with the vert. momentum adv. trends
[4619]56      !!              - Send the trends to trddyn for diagnostics (l_trddyn=T)
[3294]57      !!----------------------------------------------------------------------
[503]58      INTEGER, INTENT(in) ::   kt   ! ocean time-step inedx
[2715]59      !
[503]60      INTEGER  ::   ji, jj, jk      ! dummy loop indices
61      REAL(wp) ::   zua, zva        ! temporary scalars
[3294]62      REAL(wp), POINTER, DIMENSION(:,:,:) ::  zwuw , zwvw
63      REAL(wp), POINTER, DIMENSION(:,:  ) ::  zww
64      REAL(wp), POINTER, DIMENSION(:,:,:) ::  ztrdu, ztrdv
[3]65      !!----------------------------------------------------------------------
[3294]66      !
67      IF( nn_timing == 1 )  CALL timing_start('dyn_zad')
68      !
69      CALL wrk_alloc( jpi,jpj, zww ) 
70      CALL wrk_alloc( jpi,jpj,jpk, zwuw , zwvw ) 
71      !
[3]72      IF( kt == nit000 ) THEN
73         IF(lwp)WRITE(numout,*)
74         IF(lwp)WRITE(numout,*) 'dyn_zad : arakawa advection scheme'
75      ENDIF
[216]76
[503]77      IF( l_trddyn )   THEN         ! Save ua and va trends
[3294]78         CALL wrk_alloc( jpi, jpj, jpk, ztrdu, ztrdv ) 
[503]79         ztrdu(:,:,:) = ua(:,:,:) 
80         ztrdv(:,:,:) = va(:,:,:) 
[216]81      ENDIF
[3]82     
[503]83      DO jk = 2, jpkm1              ! Vertical momentum advection at level w and u- and v- vertical
84         DO jj = 2, jpj                   ! vertical fluxes
85            DO ji = fs_2, jpi             ! vector opt.
[3]86               zww(ji,jj) = 0.25 * e1t(ji,jj) * e2t(ji,jj) * wn(ji,jj,jk)
87            END DO
88         END DO
[503]89         DO jj = 2, jpjm1                 ! vertical momentum advection at w-point
90            DO ji = fs_2, fs_jpim1        ! vector opt.
[3]91               zwuw(ji,jj,jk) = ( zww(ji+1,jj  ) + zww(ji,jj) ) * ( un(ji,jj,jk-1)-un(ji,jj,jk) )
92               zwvw(ji,jj,jk) = ( zww(ji  ,jj+1) + zww(ji,jj) ) * ( vn(ji,jj,jk-1)-vn(ji,jj,jk) )
93            END DO 
94         END DO   
95      END DO
[503]96      DO jj = 2, jpjm1              ! Surface and bottom values set to zero
97         DO ji = fs_2, fs_jpim1           ! vector opt.
[3]98            zwuw(ji,jj, 1 ) = 0.e0
99            zwvw(ji,jj, 1 ) = 0.e0
100            zwuw(ji,jj,jpk) = 0.e0
101            zwvw(ji,jj,jpk) = 0.e0
102         END DO 
103      END DO
104
[503]105      DO jk = 1, jpkm1              ! Vertical momentum advection at u- and v-points
[3]106         DO jj = 2, jpjm1
[503]107            DO ji = fs_2, fs_jpim1       ! vector opt.
108               !                         ! vertical momentum advective trends
[3]109               zua = - ( zwuw(ji,jj,jk) + zwuw(ji,jj,jk+1) ) / ( e1u(ji,jj) * e2u(ji,jj) * fse3u(ji,jj,jk) )
110               zva = - ( zwvw(ji,jj,jk) + zwvw(ji,jj,jk+1) ) / ( e1v(ji,jj) * e2v(ji,jj) * fse3v(ji,jj,jk) )
[503]111               !                         ! add the trends to the general momentum trends
[3]112               ua(ji,jj,jk) = ua(ji,jj,jk) + zua
113               va(ji,jj,jk) = va(ji,jj,jk) + zva
114            END DO 
115         END DO 
116      END DO
117
[503]118      IF( l_trddyn ) THEN           ! save the vertical advection trends for diagnostic
119         ztrdu(:,:,:) = ua(:,:,:) - ztrdu(:,:,:)
120         ztrdv(:,:,:) = va(:,:,:) - ztrdv(:,:,:)
[4619]121         CALL trd_dyn( ztrdu, ztrdv, jpdyn_zad, kt )
[3294]122         CALL wrk_dealloc( jpi, jpj, jpk, ztrdu, ztrdv ) 
[216]123      ENDIF
[503]124      !                             ! Control print
125      IF(ln_ctl)   CALL prt_ctl( tab3d_1=ua, clinfo1=' zad  - Ua: ', mask1=umask,   &
126         &                       tab3d_2=va, clinfo2=       ' Va: ', mask2=vmask, clinfo3='dyn' )
127      !
[3294]128      CALL wrk_dealloc( jpi,jpj, zww ) 
129      CALL wrk_dealloc( jpi,jpj,jpk, zwuw , zwvw ) 
[2715]130      !
[3294]131      IF( nn_timing == 1 )  CALL timing_stop('dyn_zad')
132      !
[3]133   END SUBROUTINE dyn_zad
134
[503]135   !!======================================================================
[3]136END MODULE dynzad
Note: See TracBrowser for help on using the repository browser.