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.

step.F90 in NEMO/trunk/src/OCE – NEMO

Context Navigation

source: NEMO/trunk/src/OCE/step.F90 @ 14201

Last change on this file since 14201 was 14201, checked in by smasson, 3 years ago
trunk: phase step.F90 and stpmlf.F90...
Property svn:keywords set to `Id`
File size: 24.2 KB

Line
1	MODULE step
2	!!======================================================================
3	!! * MODULE step *
4	!! Time-stepping : manager of the ocean, tracer and ice time stepping
5	!!======================================================================
6	!! History : OPA ! 1991-03 (G. Madec) Original code
7	!! - ! 1991-11 (G. Madec)
8	!! - ! 1992-06 (M. Imbard) add a first output record
9	!! - ! 1996-04 (G. Madec) introduction of dynspg
10	!! - ! 1996-04 (M.A. Foujols) introduction of passive tracer
11	!! 8.0 ! 1997-06 (G. Madec) new architecture of call
12	!! 8.2 ! 1997-06 (G. Madec, M. Imbard, G. Roullet) free surface
13	!! - ! 1999-02 (G. Madec, N. Grima) hpg implicit
14	!! - ! 2000-07 (J-M Molines, M. Imbard) Open Bondary Conditions
15	!! NEMO 1.0 ! 2002-06 (G. Madec) free form, suppress macro-tasking
16	!! - ! 2004-08 (C. Talandier) New trends organization
17	!! - ! 2005-01 (C. Ethe) Add the KPP closure scheme
18	!! - ! 2005-11 (G. Madec) Reorganisation of tra and dyn calls
19	!! - ! 2006-01 (L. Debreu, C. Mazauric) Agrif implementation
20	!! - ! 2006-07 (S. Masson) restart using iom
21	!! 3.2 ! 2009-02 (G. Madec, R. Benshila) reintroduicing z*-coordinate
22	!! - ! 2009-06 (S. Masson, G. Madec) TKE restart compatible with key_cpl
23	!! 3.3 ! 2010-05 (K. Mogensen, A. Weaver, M. Martin, D. Lea) Assimilation interface
24	!! - ! 2010-10 (C. Ethe, G. Madec) reorganisation of initialisation phase + merge TRC-TRA
25	!! 3.4 ! 2011-04 (G. Madec, C. Ethe) Merge of dtatem and dtasal
26	!! 3.6 ! 2012-07 (J. Simeon, G. Madec. C. Ethe) Online coarsening of outputs
27	!! 3.6 ! 2014-04 (F. Roquet, G. Madec) New equations of state
28	!! 3.6 ! 2014-10 (E. Clementi, P. Oddo) Add Qiao vertical mixing in case of waves
29	!! 3.7 ! 2014-10 (G. Madec) LDF simplication
30	!! - ! 2014-12 (G. Madec) remove KPP scheme
31	!! - ! 2015-11 (J. Chanut) free surface simplification (remove filtered free surface)
32	!! 4.0 ! 2017-05 (G. Madec) introduction of the vertical physics manager (zdfphy)
33	!! 4.1 ! 2019-08 (A. Coward, D. Storkey) rewrite in preparation for new timestepping scheme
34	!!----------------------------------------------------------------------
35	#if defined key_qco \|\| defined key_linssh
36	!!----------------------------------------------------------------------
37	!! 'key_qco' EMPTY MODULE Quasi-Eulerian vertical coordinate
38	!! OR
39	!! 'key_linssh EMPTY MODULE Fixed in time vertical coordinate
40	!!----------------------------------------------------------------------
41	#else
42	!!----------------------------------------------------------------------
43	!! stp : OPA system time-stepping
44	!!----------------------------------------------------------------------
45	USE step_oce ! time stepping definition modules
46
47	IMPLICIT NONE
48	PRIVATE
49
50	PUBLIC stp ! called by nemogcm.F90
51
52	! ! time level indices !
53	INTEGER, PUBLIC :: Nbb, Nnn, Naa, Nrhs !: used by nemo_init
54
55	!! * Substitutions
56	# include "do_loop_substitute.h90"
57	!!----------------------------------------------------------------------
58	!! NEMO/OCE 4.0 , NEMO Consortium (2018)
59	!! $Id$
60	!! Software governed by the CeCILL license (see ./LICENSE)
61	!!----------------------------------------------------------------------
62	CONTAINS
63
64	#if defined key_agrif
65	RECURSIVE SUBROUTINE stp( )
66	INTEGER :: kstp ! ocean time-step index
67	#else
68	SUBROUTINE stp( kstp )
69	INTEGER, INTENT(in) :: kstp ! ocean time-step index
70	#endif
71	!!----------------------------------------------------------------------
72	!! * ROUTINE stp *
73	!!
74	!! ** Purpose : - Time stepping of OPA (momentum and active tracer eqs.)
75	!! - Time stepping of SI3 (dynamic and thermodynamic eqs.)
76	!! - Time stepping of TRC (passive tracer eqs.)
77	!!
78	!! ** Method : -1- Update forcings and data
79	!! -2- Update ocean physics
80	!! -3- Compute the t and s trends
81	!! -4- Update t and s
82	!! -5- Compute the momentum trends
83	!! -6- Update the horizontal velocity
84	!! -7- Compute the diagnostics variables (rd,N2, hdiv,w)
85	!! -8- Outputs and diagnostics
86	!!----------------------------------------------------------------------
87	INTEGER :: ji, jj, jk, jtile ! dummy loop indice
88	!! ---------------------------------------------------------------------
89	#if defined key_agrif
90	IF( nstop > 0 ) RETURN ! avoid to go further if an error was detected during previous time step (child grid)
91	kstp = nit000 + Agrif_Nb_Step()
92	Kbb_a = Nbb; Kmm_a = Nnn; Krhs_a = Nrhs ! agrif_oce module copies of time level indices
93	IF( lk_agrif_debug ) THEN
94	IF( Agrif_Root() .and. lwp) WRITE(,) '---'
95	IF(lwp) WRITE(,) 'Grid Number', Agrif_Fixed(),' time step ', kstp, 'int tstep', Agrif_NbStepint()
96	ENDIF
97	IF( kstp == nit000 + 1 ) lk_agrif_fstep = .FALSE.
98	# if defined key_iomput
99	IF( Agrif_Nbstepint() == 0 ) CALL iom_swap( cxios_context )
100	# endif
101	#endif
102	!
103	IF( ln_timing ) CALL timing_start('stp')
104	!
105	!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
106	! model timestep
107	!<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<
108	!
109	IF( l_1st_euler ) THEN ! start or restart with Euler 1st time-step
110	rDt = rn_Dt
111	r1_Dt = 1._wp / rDt
112	ENDIF
113	!
114	!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
115	! update I/O and calendar
116	!<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<
117	!
118	IF( kstp == nit000 ) THEN ! initialize IOM context (must be done after nemo_init for AGRIF+XIOS+OASIS)
119	CALL iom_init( cxios_context, ld_closedef=.FALSE. ) ! for model grid (including possible AGRIF zoom)
120	IF( lk_diamlr ) CALL dia_mlr_iom_init ! with additional setup for multiple-linear-regression analysis
121	CALL iom_init_closedef
122	IF( ln_crs ) CALL iom_init( TRIM(cxios_context)//"_crs" ) ! for coarse grid
123	ENDIF
124	IF( kstp == nitrst .AND. lwxios ) THEN
125	CALL iom_swap( cw_ocerst_cxt )
126	CALL iom_init_closedef( cw_ocerst_cxt )
127	CALL iom_setkt( kstp - nit000 + 1, cw_ocerst_cxt )
128	#if defined key_top
129	CALL iom_swap( cw_toprst_cxt )
130	CALL iom_init_closedef( cw_toprst_cxt )
131	CALL iom_setkt( kstp - nit000 + 1, cw_toprst_cxt )
132	#endif
133	ENDIF
134	#if defined key_si3
135	IF( kstp + nn_fsbc - 1 == nitrst .AND. lwxios ) THEN
136	CALL iom_swap( cw_icerst_cxt )
137	CALL iom_init_closedef( cw_icerst_cxt )
138	CALL iom_setkt( kstp - nit000 + 1, cw_icerst_cxt )
139	ENDIF
140	#endif
141	IF( kstp /= nit000 ) CALL day( kstp ) ! Calendar (day was already called at nit000 in day_init)
142	CALL iom_setkt( kstp - nit000 + 1, cxios_context ) ! tell IOM we are at time step kstp
143	IF( ln_crs ) CALL iom_setkt( kstp - nit000 + 1, TRIM(cxios_context)//"_crs" ) ! tell IOM we are at time step kstp
144
145	!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
146	! Update external forcing (tides, open boundaries, ice shelf interaction and surface boundary condition (including sea-ice)
147	!<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<
148	IF( ln_tide ) CALL tide_update( kstp ) ! update tide potential
149	IF( ln_apr_dyn ) CALL sbc_apr ( kstp ) ! atmospheric pressure (NB: call before bdy_dta which needs ssh_ib)
150	IF( ln_bdy ) CALL bdy_dta ( kstp, Nnn ) ! update dynamic & tracer data at open boundaries
151	IF( ln_isf ) CALL isf_stp ( kstp, Nnn )
152	CALL sbc ( kstp, Nbb, Nnn ) ! Sea Boundary Condition (including sea-ice)
153
154	!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
155	! Update stochastic parameters and random T/S fluctuations
156	!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
157	IF( ln_sto_eos ) CALL sto_par( kstp ) ! Stochastic parameters
158	IF( ln_sto_eos ) CALL sto_pts( ts(:,:,:,:,Nnn) ) ! Random T/S fluctuations
159
160	!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
161	! Ocean physics update
162	!<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<
163	! THERMODYNAMICS
164	CALL eos_rab( ts(:,:,:,:,Nbb), rab_b, Nnn ) ! before local thermal/haline expension ratio at T-points
165	CALL eos_rab( ts(:,:,:,:,Nnn), rab_n, Nnn ) ! now local thermal/haline expension ratio at T-points
166	CALL bn2 ( ts(:,:,:,:,Nbb), rab_b, rn2b, Nnn ) ! before Brunt-Vaisala frequency
167	CALL bn2 ( ts(:,:,:,:,Nnn), rab_n, rn2, Nnn ) ! now Brunt-Vaisala frequency
168
169	! VERTICAL PHYSICS
170	CALL zdf_phy( kstp, Nbb, Nnn, Nrhs ) ! vertical physics update (top/bot drag, avt, avs, avm + MLD)
171
172	! LATERAL PHYSICS
173	!
174	IF( l_ldfslp ) THEN ! slope of lateral mixing
175	CALL eos( ts(:,:,:,:,Nbb), rhd, gdept_0(:,:,:) ) ! before in situ density
176
177	IF( ln_zps .AND. .NOT. ln_isfcav) &
178	& CALL zps_hde ( kstp, Nnn, jpts, ts(:,:,:,:,Nbb), gtsu, gtsv, & ! Partial steps: before horizontal gradient
179	& rhd, gru , grv ) ! of t, s, rd at the last ocean level
180
181	IF( ln_zps .AND. ln_isfcav) &
182	& CALL zps_hde_isf( kstp, Nnn, jpts, ts(:,:,:,:,Nbb), gtsu, gtsv, gtui, gtvi, & ! Partial steps for top cell (ISF)
183	& rhd, gru , grv , grui, grvi ) ! of t, s, rd at the first ocean level
184	IF( ln_traldf_triad ) THEN
185	CALL ldf_slp_triad( kstp, Nbb, Nnn ) ! before slope for triad operator
186	ELSE
187	CALL ldf_slp ( kstp, rhd, rn2b, Nbb, Nnn ) ! before slope for standard operator
188	ENDIF
189	ENDIF
190	! ! eddy diffusivity coeff.
191	IF( l_ldftra_time .OR. l_ldfeiv_time ) CALL ldf_tra( kstp, Nbb, Nnn ) ! and/or eiv coeff.
192	IF( l_ldfdyn_time ) CALL ldf_dyn( kstp, Nbb ) ! eddy viscosity coeff.
193
194	!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
195	! Ocean dynamics : hdiv, ssh, e3, u, v, w
196	!<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<
197
198	CALL ssh_nxt ( kstp, Nbb, Nnn, ssh, Naa ) ! after ssh (includes call to div_hor)
199	IF( .NOT.ln_linssh ) &
200	& CALL dom_vvl_sf_nxt( kstp, Nbb, Nnn, Naa ) ! after vertical scale factors
201	CALL wzv ( kstp, Nbb, Nnn, Naa, ww ) ! now cross-level velocity
202	IF( ln_zad_Aimp ) CALL wAimp ( kstp, Nnn ) ! Adaptive-implicit vertical advection partitioning
203	CALL eos ( ts(:,:,:,:,Nnn), rhd, rhop, gdept(:,:,:,Nnn) ) ! now in situ density for hpg computation
204
205
206	uu(:,:,:,Nrhs) = 0._wp ! set dynamics trends to zero
207	vv(:,:,:,Nrhs) = 0._wp
208
209	IF( lk_asminc .AND. ln_asmiau .AND. ln_dyninc ) &
210	& CALL dyn_asm_inc ( kstp, Nbb, Nnn, uu, vv, Nrhs ) ! apply dynamics assimilation increment
211	IF( ln_bdy ) CALL bdy_dyn3d_dmp ( kstp, Nbb, uu, vv, Nrhs ) ! bdy damping trends
212	#if defined key_agrif
213	IF(.NOT. Agrif_Root()) &
214	& CALL Agrif_Sponge_dyn ! momentum sponge
215	#endif
216	CALL dyn_adv( kstp, Nbb, Nnn , uu, vv, Nrhs ) ! advection (VF or FF) ==> RHS
217	CALL dyn_vor( kstp, Nnn , uu, vv, Nrhs ) ! vorticity ==> RHS
218	CALL dyn_ldf( kstp, Nbb, Nnn , uu, vv, Nrhs ) ! lateral mixing
219	IF( ln_zdfosm ) CALL dyn_osm( kstp, Nnn , uu, vv, Nrhs ) ! OSMOSIS non-local velocity fluxes ==> RHS
220	CALL dyn_hpg( kstp, Nnn , uu, vv, Nrhs ) ! horizontal gradient of Hydrostatic pressure
221	CALL dyn_spg( kstp, Nbb, Nnn, Nrhs, uu, vv, ssh, uu_b, vv_b, Naa ) ! surface pressure gradient
222
223	! With split-explicit free surface, since now transports have been updated and ssh(:,:,Nrhs) as well
224	IF( ln_dynspg_ts ) THEN ! vertical scale factors and vertical velocity need to be updated
225	CALL div_hor ( kstp, Nbb, Nnn ) ! Horizontal divergence (2nd call in time-split case)
226	IF(.NOT.ln_linssh) CALL dom_vvl_sf_nxt( kstp, Nbb, Nnn, Naa, kcall=2 ) ! after vertical scale factors (update depth average component)
227	ENDIF
228	CALL dyn_zdf ( kstp, Nbb, Nnn, Nrhs, uu, vv, Naa ) ! vertical diffusion
229	IF( ln_dynspg_ts ) THEN ! vertical scale factors and vertical velocity need to be updated
230	CALL wzv ( kstp, Nbb, Nnn, Naa, ww ) ! Nnn cross-level velocity
231	IF( ln_zad_Aimp ) CALL wAimp ( kstp, Nnn ) ! Adaptive-implicit vertical advection partitioning
232	ENDIF
233
234
235	!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
236	! cool skin
237	!<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<
238	IF ( ln_diurnal ) CALL diurnal_layers( kstp )
239
240	!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
241	! diagnostics and outputs
242	!<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<
243	IF( ln_floats ) CALL flo_stp ( kstp, Nbb, Nnn ) ! drifting Floats
244	IF( ln_diacfl ) CALL dia_cfl ( kstp, Nnn ) ! Courant number diagnostics
245	CALL dia_hth ( kstp, Nnn ) ! Thermocline depth (20 degres isotherm depth)
246	IF( ln_diadct ) CALL dia_dct ( kstp, Nnn ) ! Transports
247	CALL dia_ar5 ( kstp, Nnn ) ! ar5 diag
248	CALL dia_ptr ( kstp, Nnn ) ! Poleward adv/ldf TRansports diagnostics
249	CALL dia_wri ( kstp, Nnn ) ! ocean model: outputs
250	IF( ln_crs ) CALL crs_fld ( kstp, Nnn ) ! ocean model: online field coarsening & output
251	IF( lk_diadetide ) CALL dia_detide( kstp ) ! Weights computation for daily detiding of model diagnostics
252	IF( lk_diamlr ) CALL dia_mlr ! Update time used in multiple-linear-regression analysis
253
254	#if defined key_top
255	!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
256	! Passive Tracer Model
257	!<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<
258	CALL trc_stp ( kstp, Nbb, Nnn, Nrhs, Naa ) ! time-stepping
259	#endif
260
261	!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
262	! Active tracers
263	!<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<
264	! Loop over tile domains
265	DO jtile = 1, nijtile
266	IF( ln_tile ) CALL dom_tile( ntsi, ntsj, ntei, ntej, ktile = jtile )
267
268	DO_3D( 0, 0, 0, 0, 1, jpk )
269	ts(ji,jj,jk,:,Nrhs) = 0._wp ! set tracer trends to zero
270	END_3D
271
272	IF( lk_asminc .AND. ln_asmiau .AND. &
273	& ln_trainc ) CALL tra_asm_inc( kstp, Nbb, Nnn, ts, Nrhs ) ! apply tracer assimilation increment
274	CALL tra_sbc ( kstp, Nnn, ts, Nrhs ) ! surface boundary condition
275	IF( ln_traqsr ) CALL tra_qsr ( kstp, Nnn, ts, Nrhs ) ! penetrative solar radiation qsr
276	IF( ln_isf ) CALL tra_isf ( kstp, Nnn, ts, Nrhs ) ! ice shelf heat flux
277	IF( ln_trabbc ) CALL tra_bbc ( kstp, Nnn, ts, Nrhs ) ! bottom heat flux
278	IF( ln_trabbl ) CALL tra_bbl ( kstp, Nbb, Nnn, ts, Nrhs ) ! advective (and/or diffusive) bottom boundary layer scheme
279	IF( ln_tradmp ) CALL tra_dmp ( kstp, Nbb, Nnn, ts, Nrhs ) ! internal damping trends
280	IF( ln_bdy ) CALL bdy_tra_dmp( kstp, Nbb, ts, Nrhs ) ! bdy damping trends
281	END DO
282
283	#if defined key_agrif
284	IF(.NOT. Agrif_Root()) THEN
285	IF( ln_tile ) CALL dom_tile( ntsi, ntsj, ntei, ntej, ktile = 0 )
286	CALL Agrif_Sponge_tra ! tracers sponge
287	ENDIF
288	#endif
289
290	! TEMP: [tiling] Separate loop over tile domains (due to tra_adv workarounds for tiling)
291	DO jtile = 1, nijtile
292	IF( ln_tile ) CALL dom_tile( ntsi, ntsj, ntei, ntej, ktile = jtile )
293
294	CALL tra_adv ( kstp, Nbb, Nnn, ts, Nrhs ) ! hor. + vert. advection ==> RHS
295	IF( ln_zdfmfc ) CALL tra_mfc ( kstp, Nbb, ts, Nrhs ) ! Mass Flux Convection
296	IF( ln_zdfosm ) CALL tra_osm ( kstp, Nnn, ts, Nrhs ) ! OSMOSIS non-local tracer fluxes ==> RHS
297	IF( lrst_oce .AND. ln_zdfosm ) &
298	& CALL osm_rst ( kstp, Nnn, 'WRITE' ) ! write OSMOSIS outputs + ww (so must do here) to restarts
299	CALL tra_ldf ( kstp, Nbb, Nnn, ts, Nrhs ) ! lateral mixing
300
301	CALL tra_zdf ( kstp, Nbb, Nnn, Nrhs, ts, Naa ) ! vertical mixing and after tracer fields
302	IF( ln_zdfnpc ) CALL tra_npc ( kstp, Nnn, Nrhs, ts, Naa ) ! update after fields by non-penetrative convection
303	END DO
304
305	IF( ln_tile ) CALL dom_tile( ntsi, ntsj, ntei, ntej, ktile = 0 ) ! Revert to tile over full domain
306	!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
307	! Set boundary conditions, time filter and swap time levels
308	!<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<
309	!!jc1: For agrif, it would be much better to finalize tracers/momentum here (e.g. bdy conditions) and move the swap
310	!! (and time filtering) after Agrif update. Then restart would be done after and would contain updated fields.
311	!! If so:
312	!! (i) no need to call agrif update at initialization time
313	!! (ii) no need to update "before" fields
314	!!
315	!! Apart from creating new tra_swp/dyn_swp routines, this however:
316	!! (i) makes boundary conditions at initialization time computed from updated fields which is not the case between
317	!! two restarts => restartability issue. One can circumvent this, maybe, by assuming "interface separation",
318	!! e.g. a shift of the feedback interface inside child domain.
319	!! (ii) requires that all restart outputs of updated variables by agrif (e.g. passive tracers/tke/barotropic arrays) are done at the same
320	!! place.
321	!!
322	!!jc2: dynnxt must be the latest call. e3t(:,:,:,Nbb) are indeed updated in that routine
323	CALL tra_atf ( kstp, Nbb, Nnn, Naa, ts ) ! time filtering of "now" tracer arrays
324	CALL dyn_atf ( kstp, Nbb, Nnn, Naa, uu, vv, e3t, e3u, e3v ) ! time filtering of "now" velocities and scale factors
325	CALL ssh_atf ( kstp, Nbb, Nnn, Naa, ssh ) ! time filtering of "now" sea surface height
326	!
327	! Swap time levels
328	Nrhs = Nbb
329	Nbb = Nnn
330	Nnn = Naa
331	Naa = Nrhs
332	!
333	IF(.NOT.ln_linssh) CALL dom_vvl_sf_update( kstp, Nbb, Nnn, Naa ) ! recompute vertical scale factors
334	!
335	IF( ln_diahsb ) CALL dia_hsb ( kstp, Nbb, Nnn ) ! - ML - global conservation diagnostics
336
337	!!gm : This does not only concern the dynamics ==>>> add a new title
338	!!gm2: why ouput restart before AGRIF update?
339	!!
340	!!jc: That would be better, but see comment above
341	!!
342	IF( lrst_oce ) CALL rst_write ( kstp, Nbb, Nnn ) ! write output ocean restart file
343	IF( ln_sto_eos ) CALL sto_rst_write( kstp ) ! write restart file for stochastic parameters
344
345	#if defined key_agrif
346	!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
347	! AGRIF recursive integration
348	!<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<
349	Kbb_a = Nbb; Kmm_a = Nnn; Krhs_a = Nrhs ! agrif_oce module copies of time level indices
350	CALL Agrif_Integrate_ChildGrids( stp ) ! allows to finish all the Child Grids before updating
351
352	#endif
353	!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
354	! Control
355	!<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<
356	CALL stp_ctl ( kstp, Nnn )
357
358	#if defined key_agrif
359	!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
360	! AGRIF update
361	!<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<
362	IF( Agrif_NbStepint() == 0 .AND. nstop == 0 ) &
363	& CALL Agrif_update_all( ) ! Update all components
364
365	#endif
366	IF( ln_diaobs .AND. nstop == 0 ) &
367	& CALL dia_obs( kstp, Nnn ) ! obs-minus-model (assimilation) diags (after dynamics update)
368
369	!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
370	! File manipulation at the end of the first time step
371	!<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<
372	IF( kstp == nit000 ) THEN ! 1st time step only
373	CALL iom_close( numror ) ! close input ocean restart file
374	IF( lrxios ) CALL iom_context_finalize( cr_ocerst_cxt )
375	IF(lwm) CALL FLUSH ( numond ) ! flush output namelist oce
376	IF(lwm .AND. numoni /= -1 ) CALL FLUSH ( numoni ) ! flush output namelist ice (if exist)
377	ENDIF
378
379	!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
380	! Coupled mode
381	!<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<
382	IF( lk_oasis .AND. nstop == 0 ) CALL sbc_cpl_snd( kstp, Nbb, Nnn ) ! coupled mode : field exchanges
383	!
384	#if defined key_iomput
385	!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
386	! Finalize contextes if end of simulation or error detected
387	!<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<
388	IF( kstp == nitend .OR. nstop > 0 ) THEN
389	CALL iom_context_finalize( cxios_context ) ! needed for XIOS+AGRIF
390	IF( ln_crs ) CALL iom_context_finalize( trim(cxios_context)//"_crs" ) !
391	ENDIF
392	#endif
393	!
394	IF( l_1st_euler ) THEN ! recover Leap-frog timestep
395	rDt = 2._wp * rn_Dt
396	r1_Dt = 1._wp / rDt
397	l_1st_euler = .FALSE.
398	ENDIF
399	!
400	IF( ln_timing ) CALL timing_stop('stp')
401	!
402	END SUBROUTINE stp
403	!
404	#endif
405	!!======================================================================
406	END MODULE step

Note: See TracBrowser for help on using the repository browser.

Download in other formats: