source: branches/nemo_v3_3_beta/NEMOGCM/NEMO/C1D_SRC/dynnxt_c1d.F90 @ 2287

Last change on this file since 2287 was 2287, checked in by smasson, 10 years ago

update licence of all NEMO files…

  • Property svn:keywords set to Id
File size: 4.9 KB
Line 
1MODULE dynnxt_c1d
2   !!======================================================================
3   !!                       ***  MODULE  dynnxt_c1d  ***
4   !! Ocean dynamics: time stepping in 1D configuration
5   !!======================================================================
6   !! History :  2.0  !  2004-10  (C. Ethe)  Original code from dynnxt.F90
7   !!            3.0  !  2008-04  (G.madec)  Style only
8   !!----------------------------------------------------------------------
9#if defined key_c1d
10   !!----------------------------------------------------------------------
11   !!   'key_c1d'                                          1D Configuration
12   !!---------------------------------------------------------------------- 
13   !!----------------------------------------------------------------------
14   !!   dyn_nxt_c1d : update the horizontal velocity from the momentum trend
15   !!----------------------------------------------------------------------
16   USE oce             ! ocean dynamics and tracers
17   USE dom_oce         ! ocean space and time domain
18   USE in_out_manager  ! I/O manager
19   USE lbclnk          ! lateral boundary condition (or mpp link)
20   USE prtctl          ! Print control
21
22   IMPLICIT NONE
23   PRIVATE
24
25   PUBLIC dyn_nxt_c1d                ! routine called by step.F90
26   !!----------------------------------------------------------------------
27   !! NEMO/C1D 3.3 , NEMO Consortium (2010)
28   !! $Id$
29   !! Software governed by the CeCILL licence (NEMOGCM/NEMO_CeCILL.txt)
30   !!----------------------------------------------------------------------
31
32CONTAINS
33
34   SUBROUTINE dyn_nxt_c1d ( kt )
35      !!----------------------------------------------------------------------
36      !!                  ***  ROUTINE dyn_nxt_c1d  ***
37      !!                   
38      !! ** Purpose :   Compute the after horizontal velocity from the
39      !!      momentum trend.
40      !!
41      !! ** Method  :   Apply lateral boundary conditions on the trends (ua,va)
42      !!      through calls to routine lbc_lnk.
43      !!      After velocity is compute using a leap-frog scheme environment:
44      !!         (ua,va) = (ub,vb) + 2 rdt (ua,va)
45      !!      Time filter applied on now horizontal velocity to avoid the
46      !!      divergence of two consecutive time-steps and swap of dynamics
47      !!      arrays to start the next time step:
48      !!         (ub,vb) = (un,vn) + atfp [ (ub,vb) + (ua,va) - 2 (un,vn) ]
49      !!         (un,vn) = (ua,va)
50      !!
51      !! ** Action : - Update ub,vb arrays, the before horizontal velocity
52      !!             - Update un,vn arrays, the now horizontal velocity
53      !!----------------------------------------------------------------------
54      INTEGER, INTENT( in ) ::   kt      ! ocean time-step index
55      !!
56      INTEGER  ::   jk           ! dummy loop indices
57      REAL(wp) ::   z2dt         ! temporary scalar
58      !!----------------------------------------------------------------------
59
60      IF( kt == nit000 ) THEN
61         IF(lwp) WRITE(numout,*)
62         IF(lwp) WRITE(numout,*) 'dyn_nxt_c1d : time stepping on 1D configuation'
63         IF(lwp) WRITE(numout,*) '~~~~~~~'
64      ENDIF
65
66      ! Local constant initialization
67      z2dt = 2. * rdt
68      IF( neuler == 0 .AND. kt == nit000 )  z2dt = rdt
69
70      ! Lateral boundary conditions on ( ua, va )
71      CALL lbc_lnk( ua, 'U', -1. )
72      CALL lbc_lnk( va, 'V', -1. )
73
74      DO jk = 1, jpkm1      ! Next Velocity
75         ua(:,:,jk) = ( ub(:,:,jk) + z2dt * ua(:,:,jk) ) * umask(:,:,jk)
76         va(:,:,jk) = ( vb(:,:,jk) + z2dt * va(:,:,jk) ) * vmask(:,:,jk)
77      END DO
78 
79      DO jk = 1, jpkm1      ! Time filter and swap of dynamics arrays
80         IF( neuler == 0 .AND. kt == nit000 ) THEN      ! Euler (forward) time stepping
81             ub(:,:,jk) = un(:,:,jk)
82             vb(:,:,jk) = vn(:,:,jk)
83             un(:,:,jk) = ua(:,:,jk)
84             vn(:,:,jk) = va(:,:,jk)
85         ELSE                                           ! Leap-frog time stepping
86             ub(:,:,jk) = atfp * ( ub(:,:,jk) + ua(:,:,jk) ) + atfp1 * un(:,:,jk)
87             vb(:,:,jk) = atfp * ( vb(:,:,jk) + va(:,:,jk) ) + atfp1 * vn(:,:,jk)
88             un(:,:,jk) = ua(:,:,jk)
89             vn(:,:,jk) = va(:,:,jk)
90         ENDIF
91      END DO
92
93      IF(ln_ctl)   CALL prt_ctl( tab3d_1=un, clinfo1=' nxt_c1d  - Un: ', mask1=umask,   &
94         &                       tab3d_2=vn, clinfo2=' Vn: '           , mask2=vmask )
95      !
96   END SUBROUTINE dyn_nxt_c1d
97
98#else
99   !!----------------------------------------------------------------------
100   !!   Default key                                     NO 1D Config
101   !!----------------------------------------------------------------------
102CONTAINS
103   SUBROUTINE dyn_nxt_c1d ( kt )
104      WRITE(*,*) 'dyn_nxt_c1d: You should not have seen this print! error?', kt
105   END SUBROUTINE dyn_nxt_c1d
106#endif
107
108   !!======================================================================
109END MODULE dynnxt_c1d
Note: See TracBrowser for help on using the repository browser.