Context Navigation

← Previous Revision
Next Revision →
Normal
Revision Log

step.F90

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

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

New URL for NEMO forge! http://forge.nemo-ocean.eu

Context Navigation

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

Download in other formats: