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.
domstp.F90 in branches/2011/DEV_r2739_STFC_dCSE/NEMOGCM/NEMO/OPA_SRC/DOM – NEMO

source: branches/2011/DEV_r2739_STFC_dCSE/NEMOGCM/NEMO/OPA_SRC/DOM/domstp.F90 @ 4400

Last change on this file since 4400 was 3211, checked in by spickles2, 12 years ago

Stephen Pickles, 11 Dec 2011

Commit to bring the rest of the DCSE NEMO development branch
in line with the latest development version. This includes
array index re-ordering of all OPA_SRC/.
  • Property svn:keywords set to Id
File size: 5.0 KB
Line 
1MODULE domstp
2   !!==============================================================================
3   !!                       ***  MODULE domstp   ***
4   !! Ocean initialization : time domain
5   !!==============================================================================
6
7   !!----------------------------------------------------------------------
8   !!   dom_stp        : ocean time domain initialization
9   !!----------------------------------------------------------------------
10   !! History :  OPA  ! 1990-10  (O. Marti)  Original code
11   !!                 ! 1996-01  (G. Madec)  terrain following coordinates
12   !!   NEMO     1.0  ! 2002-08  (G. Madec)  F90: Free form and module
13   !!----------------------------------------------------------------------
14   USE oce            ! ocean dynamics and tracers
15   USE dom_oce        ! ocean space and time domain
16   USE in_out_manager ! I/O manager
17   USE lib_mpp        ! MPP library
18
19   IMPLICIT NONE
20   PRIVATE
21
22   PUBLIC   dom_stp   ! routine called by inidom.F90
23
24   !! * Control permutation of array indices
25#  include "oce_ftrans.h90"
26#  include "dom_oce_ftrans.h90"
27
28   !! * Substitutions
29#  include "domzgr_substitute.h90"
30   !!----------------------------------------------------------------------
31   !! NEMO/OPA 3.3 , NEMO Consortium (2010)
32   !! $Id$
33   !! Software governed by the CeCILL licence     (NEMOGCM/NEMO_CeCILL.txt)
34   !!----------------------------------------------------------------------
35CONTAINS
36
37   SUBROUTINE dom_stp
38      !!----------------------------------------------------------------------
39      !!                    ***  ROUTINE dom_stp  ***
40      !!         
41      !! ** Purpose :   Intialize ocean time step for the run
42      !!
43      !! ** Method  : - Initialization of a coef. use in the Asselin time
44      !!      filter:  atfp1 = 1 - 2 * atfp  where atfp is the Asselin time
45      !!      filter parameter read in namelist
46      !!              - Model time step:
47      !!      nacc = 0 : synchronous time intergration.
48      !!      There is one time step only, defined by: rdt, rdttra(k)=rdt
49      !!      nacc = 1 : accelerating the convergence. There is 2 different
50      !!      time steps for dynamics and tracers:
51      !!        rdt      : dynamical part
52      !!        rdttra(k): temperature and salinity
53      !!      The tracer time step is a function of vertical level. the model
54      !!      reference time step ( i.e. for wind stress, surface heat and
55      !!      salt fluxes) is the surface tracer time step is rdttra(1).
56      !!         N.B. depth dependent acceleration of convergence is not im-
57      !!      plemented for s-coordinate.
58      !!
59      !! ** Action  : - rdttra   : vertical profile of tracer time step
60      !!              - atfp1    : = 1 - 2*atfp
61      !!
62      !! References :   Bryan, K., 1984, J. Phys. Oceanogr., 14, 666-673.
63      !!----------------------------------------------------------------------
64      INTEGER ::   jk              ! dummy loop indice
65      !!----------------------------------------------------------------------
66
67      IF(lwp) THEN
68         WRITE(numout,*)
69         WRITE(numout,*) 'dom_stp : time stepping setting'
70         WRITE(numout,*) '~~~~~~~'
71      ENDIF
72
73      ! 0. Asselin Time filter
74      ! ----------------------
75     
76      atfp1 = 1. - 2. * atfp
77
78      SELECT CASE ( nacc )
79
80         CASE ( 0 )                ! Synchronous time stepping
81            IF(lwp) WRITE(numout,*)'               synchronous time stepping'
82            IF(lwp) WRITE(numout,*)'               dynamics and tracer time step = ', rdt/3600., ' hours'
83
84            rdttra(:) = rdt
85
86         CASE ( 1 )                ! Accelerating the convergence
87            IF(lwp) WRITE(numout,*) '              no tracer damping in the turbocline'
88            IF(lwp) WRITE(numout,*)'               accelerating the convergence'
89            IF(lwp) WRITE(numout,*)'               dynamics time step = ', rdt/3600., ' hours'
90            IF( ln_sco .AND. rdtmin /= rdtmax .AND. lk_vvl )   &
91                 & CALL ctl_stop ( ' depth dependent acceleration of convergence not implemented in s-coordinates &
92                 &                   nor in variable volume' )
93            IF(lwp) WRITE(numout,*)'         tracers   time step :  dt (hours)  level'
94
95            DO jk = 1, jpk
96               IF( gdept_0(jk) <= rdth ) rdttra(jk) = rdtmin
97               IF( gdept_0(jk) >  rdth ) THEN
98                  rdttra(jk) = rdtmin + ( rdtmax - rdtmin )   &
99                                      * ( EXP( ( gdept_0(jk ) - rdth ) / rdth ) - 1. )   &
100                                      / ( EXP( ( gdept_0(jpk) - rdth ) / rdth ) - 1. )
101               ENDIF
102               IF(lwp) WRITE(numout,"(36x,f5.2,5x,i3)") rdttra(jk)/3600., jk
103            END DO 
104
105         CASE DEFAULT              ! E R R O R
106
107            WRITE(ctmp1,*) ' nacc value e r r o r, nacc= ',nacc
108            CALL ctl_stop( ctmp1 )
109
110      END SELECT
111
112   END SUBROUTINE dom_stp
113
114   !!======================================================================
115END MODULE domstp
Note: See TracBrowser for help on using the repository browser.