source: NEMO/branches/2020/dev_r12563_ASINTER-06_ABL_improvement/src/SAS/step.F90 @ 13159

Last change on this file since 13159 was 13159, checked in by gsamson, 3 months ago

merge trunk@r13136 into ASINTER-06 branch; pass all SETTE tests; results identical to trunk@r13136; ticket #2419

  • Property svn:keywords set to Id
File size: 7.2 KB
RevLine 
[3331]1MODULE step
2   !!======================================================================
3   !!                       ***  MODULE step  ***
4   !! Time-stepping    : manager of the ocean, tracer and ice time stepping
5   !!                    version for standalone surface scheme
6   !!======================================================================
7   !! History :  OPA  !  1991-03  (G. Madec)  Original code
[3362]8   !!             .   !    .                                                     
9   !!             .   !    .                                                     
10   !!   NEMO     3.5  !  2012-03  (S. Alderson)
[3331]11   !!----------------------------------------------------------------------
12
13   !!----------------------------------------------------------------------
14   !!   stp             : OPA system time-stepping
15   !!----------------------------------------------------------------------
16   USE oce              ! ocean dynamics and tracers variables
17   USE dom_oce          ! ocean space and time domain variables
18   USE daymod           ! calendar                         (day     routine)
[6140]19   USE sbc_oce          ! surface boundary condition: fields
[3331]20   USE sbcmod           ! surface boundary condition       (sbc     routine)
21   USE sbcrnf           ! surface boundary condition: runoff variables
[6140]22   USE sbccpl           ! surface boundary condition: coupled interface
[3331]23   USE eosbn2           ! equation of state                (eos_bn2 routine)
24   USE diawri           ! Standard run outputs             (dia_wri routine)
[7646]25   USE bdy_oce   , ONLY: ln_bdy
[9019]26   USE bdydta           ! mandatory for sea-ice
[6140]27   USE stpctl           ! time stepping control            (stp_ctl routine)
28   !
29   USE in_out_manager   ! I/O manager
30   USE prtctl           ! Print control                    (prt_ctl routine)
31   USE iom              !
32   USE lbclnk           !
33   USE timing           ! Timing           
34#if defined key_iomput
35   USE xios
36#endif
[5510]37
[7646]38#if defined key_agrif
[9019]39   USE agrif_oce, ONLY: lk_agrif_debug
[9570]40#if defined key_si3
[9596]41   USE agrif_ice_update
[7646]42#endif
[9510]43#endif
[7646]44   
[3331]45   IMPLICIT NONE
46   PRIVATE
47
[6140]48   PUBLIC   stp   ! called by nemogcm.F90
[3331]49
50   !!----------------------------------------------------------------------
[12377]51   !! time level indices
52   !!----------------------------------------------------------------------
53   INTEGER, PUBLIC :: Nbb, Nnn, Naa, Nrhs          !! used by nemo_init
54   !!----------------------------------------------------------------------
[10068]55   !! NEMO/SAS 4.0 , NEMO Consortium (2018)
[5215]56   !! $Id$
[10068]57   !! Software governed by the CeCILL license (see ./LICENSE)
[3331]58   !!----------------------------------------------------------------------
59CONTAINS
60
61#if defined key_agrif
[7761]62   RECURSIVE SUBROUTINE stp( )
[3331]63      INTEGER             ::   kstp   ! ocean time-step index
64#else
65   SUBROUTINE stp( kstp )
66      INTEGER, INTENT(in) ::   kstp   ! ocean time-step index
67#endif
68      !!----------------------------------------------------------------------
69      !!                     ***  ROUTINE stp  ***
70      !!                     
[3358]71      !! ** Purpose : - Time stepping of SBC (surface boundary)
[3331]72      !!
73      !! ** Method  : -1- Update forcings and data 
[3358]74      !!              -2- Outputs and diagnostics
[3331]75      !!----------------------------------------------------------------------
76
77#if defined key_agrif
[13159]78      IF( nstop > 0 ) RETURN   ! avoid to go further if an error was detected during previous time step (child grid)
[3331]79      kstp = nit000 + Agrif_Nb_Step()
[12377]80      Kbb_a = Nbb; Kmm_a = Nnn; Krhs_a = Nrhs   ! agrif_oce module copies of time level indices
[13159]81      IF( lk_agrif_debug ) THEN
82         IF( Agrif_Root() .and. lwp)   WRITE(*,*) '---'
83         IF(lwp)   WRITE(*,*) 'Grid Number', Agrif_Fixed(),' time step ', kstp, 'int tstep', Agrif_NbStepint()
[7646]84      ENDIF
[13159]85      IF( kstp == nit000 + 1 )   lk_agrif_fstep = .FALSE.
[3331]86# if defined key_iomput
[5407]87      IF( Agrif_Nbstepint() == 0 )   CALL iom_swap( cxios_context )
[3331]88# endif   
89#endif   
[5407]90      IF( kstp == nit000 )   CALL iom_init( cxios_context ) ! iom_put initialization (must be done after nemo_init for AGRIF+XIOS+OASIS)
[3331]91      IF( kstp /= nit000 )   CALL day( kstp )             ! Calendar (day was already called at nit000 in day_init)
[5407]92                             CALL iom_setkt( kstp - nit000 + 1, cxios_context )   ! tell iom we are at time step kstp
[3331]93
[9019]94      ! ==> clem: open boundaries is mandatory for sea-ice because ice BDY is not decoupled from 
[5510]95      !           the environment of ocean BDY. Therefore bdy is called in both OPA and SAS modules.
96      !           From SAS: ocean bdy data are wrong  (but we do not care) and ice bdy data are OK. 
97      !           This is not clean and should be changed in the future.
98      ! ==>
[12377]99      IF( ln_bdy     )       CALL bdy_dta( kstp,      Nnn )                   ! update dynamic & tracer data at open boundaries
100                             CALL sbc    ( kstp, Nbb, Nnn )                   ! Sea Boundary Condition (including sea-ice)
[3331]101
[12377]102                             CALL dia_wri( kstp,      Nnn )                   ! ocean model: outputs
[3331]103
[7646]104#if defined key_agrif
105      !>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
[13159]106      ! AGRIF recursive integration
[7646]107      !<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<     
[13159]108                             CALL Agrif_Integrate_ChildGrids( stp )
109                             
110#endif                             
111      !>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
112      ! Control
113      !<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<
114                             CALL stp_ctl( kstp, Nnn )
[9510]115
[13159]116#if defined key_agrif
117      !>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
118      ! AGRIF update
119      !<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<     
120      IF( Agrif_NbStepint() == 0 .AND. nstop == 0 ) THEN                       ! AGRIF Update from zoom N to zoom 1 then to Parent
[9570]121#if defined key_si3
[9781]122                             CALL Agrif_Update_ice( )   ! update sea-ice
[7646]123#endif
[9510]124      ENDIF
[13159]125
[9510]126#endif
[7646]127      !>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
[13159]128      ! File manipulation at the end of the first time step
129      !<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<                         
130      IF( kstp == nit000   ) CALL iom_close( numror )                          ! close input  ocean restart file
[7646]131     
[5407]132      !>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
133      ! Coupled mode
134      !<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<
[13159]135      IF( lk_oasis .AND. nstop == 0 ) CALL sbc_cpl_snd( kstp, Nbb, Nnn )       ! coupled mode : field exchanges if OASIS-coupled ice
[5407]136
[3769]137#if defined key_iomput
[9367]138      IF( kstp == nitrst ) THEN
139         IF(.NOT.lwxios) THEN
140            CALL iom_close( numrow )     
141         ELSE
142            CALL iom_context_finalize( cwxios_context )
143         ENDIF
144         lrst_oce = .FALSE.
145      ENDIF
[13159]146      IF( kstp == nitend .OR. nstop > 0 ) THEN
147         CALL iom_context_finalize( cxios_context ) ! needed for XIOS+AGRIF
[5407]148      ENDIF
[3769]149#endif
[3331]150      !
[9019]151      IF( ln_timing .AND.  kstp == nit000  )   CALL timing_reset
[3331]152      !
153   END SUBROUTINE stp
154
155   !!======================================================================
156END MODULE step
Note: See TracBrowser for help on using the repository browser.