source: NEMO/branches/2019/ENHANCE-02_ISF_nemo/src/OCE/DYN/divhor.F90 @ 11987

Last change on this file since 11987 was 11987, checked in by mathiot, 11 months ago

ENHANCE-02_ISF_nemo: changes needed after Dave's review

  • Property svn:keywords set to Id
File size: 5.3 KB
RevLine 
[5777]1MODULE divhor
[3]2   !!==============================================================================
[5777]3   !!                       ***  MODULE  divhor  ***
4   !! Ocean diagnostic variable : now horizontal divergence
[3]5   !!==============================================================================
[5777]6   !! History :  1.0  ! 2002-09  (G. Madec, E. Durand)  Free form, F90
[2528]7   !!             -   ! 2005-01  (J. Chanut) Unstructured open boundaries
8   !!             -   ! 2003-08  (G. Madec)  merged of cur and div, free form, F90
9   !!             -   ! 2005-01  (J. Chanut, A. Sellar) unstructured open boundaries
10   !!            3.3  ! 2010-09  (D.Storkey and E.O'Dea) bug fixes for BDY module
11   !!             -   ! 2010-10  (R. Furner, G. Madec) runoff and cla added directly here
[5777]12   !!            3.7  ! 2014-01  (G. Madec) suppression of velocity curl from in-core memory
13   !!             -   ! 2014-12  (G. Madec) suppression of cross land advection option
14   !!             -   ! 2015-10  (G. Madec) add velocity and rnf flag in argument of div_hor
[2528]15   !!----------------------------------------------------------------------
[3]16
17   !!----------------------------------------------------------------------
[5777]18   !!   div_hor    : Compute the horizontal divergence field
[3]19   !!----------------------------------------------------------------------
20   USE oce             ! ocean dynamics and tracers
21   USE dom_oce         ! ocean space and time domain
[11987]22   USE sbc_oce, ONLY : ln_rnf      ! river runoff
23   USE sbcrnf , ONLY : sbc_rnf_div ! river runoff
24   USE isf_oce, ONLY : ln_isf      ! ice shelf
25   USE isfhdiv, ONLY : isf_hdiv    ! ice shelf
[9023]26#if defined key_asminc   
27   USE asminc          ! Assimilation increment
28#endif
[5777]29   !
[3]30   USE in_out_manager  ! I/O manager
31   USE lbclnk          ! ocean lateral boundary conditions (or mpp link)
[2715]32   USE lib_mpp         ! MPP library
[3294]33   USE timing          ! Timing
[3]34
35   IMPLICIT NONE
36   PRIVATE
37
[5777]38   PUBLIC   div_hor    ! routine called by step.F90 and istate.F90
[3]39
40   !! * Substitutions
41#  include "vectopt_loop_substitute.h90"
42   !!----------------------------------------------------------------------
[9598]43   !! NEMO/OCE 4.0 , NEMO Consortium (2018)
[1152]44   !! $Id$
[10068]45   !! Software governed by the CeCILL license (see ./LICENSE)
[3]46   !!----------------------------------------------------------------------
47CONTAINS
48
[5777]49   SUBROUTINE div_hor( kt )
[3]50      !!----------------------------------------------------------------------
[5777]51      !!                  ***  ROUTINE div_hor  ***
[3]52      !!                   
[5777]53      !! ** Purpose :   compute the horizontal divergence at now time-step
[3]54      !!
[5777]55      !! ** Method  :   the now divergence is computed as :
56      !!         hdivn = 1/(e1e2t*e3t) ( di[e2u*e3u un] + dj[e1v*e3v vn] )
57      !!      and correct with runoff inflow (div_rnf) and cross land flow (div_cla)
[3]58      !!
[5777]59      !! ** Action  : - update hdivn, the now horizontal divergence
[3]60      !!----------------------------------------------------------------------
[2715]61      INTEGER, INTENT(in) ::   kt   ! ocean time-step index
[2528]62      !
[2715]63      INTEGER  ::   ji, jj, jk    ! dummy loop indices
64      REAL(wp) ::   zraur, zdep   ! local scalars
[11541]65      REAL(wp), DIMENSION(jpi,jpj) :: ztmp
[3]66      !!----------------------------------------------------------------------
[3294]67      !
[9019]68      IF( ln_timing )   CALL timing_start('div_hor')
[3294]69      !
[3]70      IF( kt == nit000 ) THEN
71         IF(lwp) WRITE(numout,*)
[5777]72         IF(lwp) WRITE(numout,*) 'div_hor : horizontal velocity divergence '
73         IF(lwp) WRITE(numout,*) '~~~~~~~   '
[3]74      ENDIF
[5777]75      !
76      DO jk = 1, jpkm1                                      !==  Horizontal divergence  ==!
[3]77         DO jj = 2, jpjm1
78            DO ji = fs_2, fs_jpim1   ! vector opt.
[9019]79               hdivn(ji,jj,jk) = (  e2u(ji  ,jj) * e3u_n(ji  ,jj,jk) * un(ji  ,jj,jk)      &
80                  &               - e2u(ji-1,jj) * e3u_n(ji-1,jj,jk) * un(ji-1,jj,jk)      &
81                  &               + e1v(ji,jj  ) * e3v_n(ji,jj  ,jk) * vn(ji,jj  ,jk)      &
82                  &               - e1v(ji,jj-1) * e3v_n(ji,jj-1,jk) * vn(ji,jj-1,jk)  )   &
83                  &            * r1_e1e2t(ji,jj) / e3t_n(ji,jj,jk)
[3]84            END DO 
85         END DO 
[5777]86      END DO
[11823]87      !
[9019]88#if defined key_agrif
89      IF( .NOT. Agrif_Root() ) THEN
[9116]90         IF( nbondi == -1 .OR. nbondi == 2 )   hdivn(   2   ,  :   ,:) = 0._wp      ! west
91         IF( nbondi ==  1 .OR. nbondi == 2 )   hdivn( nlci-1,  :   ,:) = 0._wp      ! east
92         IF( nbondj == -1 .OR. nbondj == 2 )   hdivn(   :   ,  2   ,:) = 0._wp      ! south
93         IF( nbondj ==  1 .OR. nbondj == 2 )   hdivn(   :   ,nlcj-1,:) = 0._wp      ! north
[9019]94      ENDIF
95#endif
[2715]96      !
[11852]97      IF( ln_rnf )                      CALL sbc_rnf_div( hdivn )       !==  runoffs           ==!   (update hdivn field)
[2528]98      !
[9023]99#if defined key_asminc 
100      IF( ln_sshinc .AND. ln_asmiau )   CALL ssh_asm_div( kt, hdivn )   !==  SSH assimilation  ==!   (update hdivn field)
101      !
102#endif
[3294]103      !
[11852]104      IF( ln_isf )                      CALL isf_hdiv( kt, hdivn )      !==  ice shelf         ==!   (update hdivn field)
[11529]105      !
[10425]106      CALL lbc_lnk( 'divhor', hdivn, 'T', 1. )   !   (no sign change)
[6140]107      !
[9019]108      IF( ln_timing )   CALL timing_stop('div_hor')
[5777]109      !
110   END SUBROUTINE div_hor
[3]111   
112   !!======================================================================
[5777]113END MODULE divhor
Note: See TracBrowser for help on using the repository browser.