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 tags/nemo_v3_2/nemo_v3_2/NEMO/OFF_SRC/DOM – NEMO

source: tags/nemo_v3_2/nemo_v3_2/NEMO/OFF_SRC/DOM/domstp.F90 @ 1878

Last change on this file since 1878 was 1878, checked in by flavoni, 14 years ago

initial test for nemogcm
File size: 4.9 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 :
11   !!        !  90-10  (O. Marti)  Original code
12   !!        !  96-01  (G. Madec)  terrain following coordinates
13   !!   8.5  !  02-08  (G. Madec)  F90: Free form and module
14   !!----------------------------------------------------------------------
15   !! * Modules used
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
20   IMPLICIT NONE
21   PRIVATE
22
23   !! * routine accessibility
24   PUBLIC dom_stp        ! routine called by inidom.F90
25
26   !! * Substitutions
27#  include "domzgr_substitute.h90"
28   !!----------------------------------------------------------------------
29   !!   OPA 9.0 , LOCEAN-IPSL (2005)
30   !! $Id: domstp.F90 1152 2008-06-26 14:11:13Z rblod $
31   !! Software governed by the CeCILL licence (modipsl/doc/NEMO_CeCILL.txt)
32   !!----------------------------------------------------------------------
33
34CONTAINS
35
36   SUBROUTINE dom_stp
37      !!----------------------------------------------------------------------
38      !!                    ***  ROUTINE dom_stp  ***
39      !!         
40      !! ** Purpose :   Intialize ocean time step for the run
41      !!
42      !! ** Method  : - Initialization of a coef. use in the Asselin time
43      !!      filter:  atfp1 = 1 - 2 * atfp  where atfp is the Asselin time
44      !!      filter parameter read in namelist
45      !!              - Model time step:
46      !!      nacc = 0 : synchronous time intergration.
47      !!      There is one time step only, defined by: rdt, rdttra(k)=rdt
48      !!      nacc = 1 : accelerating the convergence. There is 2 different
49      !!      time steps for dynamics and tracers:
50      !!        rdt      : dynamical part
51      !!        rdttra(k): temperature and salinity
52      !!      The tracer time step is a function of vertical level. the model
53      !!      reference time step ( i.e. for wind stress, surface heat and
54      !!      salt fluxes) is the surface tracer time step is rdttra(1).
55      !!         N.B. depth dependent acceleration of convergence is not im-
56      !!      plemented for s-coordinate.
57      !!
58      !! ** Action  : - rdttra   : vertical profile of tracer time step
59      !!              - atfp1    : = 1 - 2*atfp
60      !!
61      !! References :
62      !!      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 )   &
91                 & CALL ctl_stop ( ' depth dependent acceleration of convergence not implemented in s-coordinates' )
92            IF(lwp) WRITE(numout,*)'         tracers   time step :  dt (hours)  level'
93
94            DO jk = 1, jpk
95               IF( gdept_0(jk) <= rdth ) rdttra(jk) = rdtmin
96               IF( gdept_0(jk) >  rdth ) THEN
97                  rdttra(jk) = rdtmin + ( rdtmax - rdtmin )   &
98                                      * ( EXP( ( gdept_0(jk ) - rdth ) / rdth ) - 1. )   &
99                                      / ( EXP( ( gdept_0(jpk) - rdth ) / rdth ) - 1. )
100               ENDIF
101               IF(lwp) WRITE(numout,"(36x,f5.2,5x,i3)") rdttra(jk)/3600., jk
102            END DO 
103
104         CASE DEFAULT              ! E R R O R
105
106            WRITE(ctmp1,*) ' nacc value e r r o r, nacc= ',nacc
107            CALL ctl_stop( ctmp1 )
108
109      END SELECT
110
111   END SUBROUTINE dom_stp
112
113   !!======================================================================
114END MODULE domstp
Note: See TracBrowser for help on using the repository browser.