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/branches/2020/dev_r12377_KERNEL-06_techene_e3/src/OCE – NEMO

Context Navigation

source: NEMO/branches/2020/dev_r12377_KERNEL-06_techene_e3/src/OCE/step.F90 @ 12724

Last change on this file since 12724 was 12724, checked in by techene, 4 years ago
branch KERNEL-06 : merge with trunk@12698 #2385 - in duplcated files : changes to comply to the new trunk variables and some loop bug fixes
Property svn:keywords set to `Id`
File size: 22.9 KB

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

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

Download in other formats: