Context Navigation

step.F90 @ 10764

Last change on this file since 10764 was 10764, checked in by jcastill, 5 years ago
Minimal set of changes for coupling order
File size: 20.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
[5930]	31	!! - ! 2015-11 (J. Chanut) free surface simplification
[3]	32	!!----------------------------------------------------------------------
[888]	33
	34	!!----------------------------------------------------------------------
[2528]	35	!! stp : OPA system time-stepping
[3]	36	!!----------------------------------------------------------------------
[3764]	37	USE step_oce ! time stepping definition modules
[6140]	38	!
	39	USE iom ! xIOs server
[389]	40
[3]	41	IMPLICIT NONE
	42	PRIVATE
	43
[5836]	44	PUBLIC stp ! called by nemogcm.F90
[3]	45
	46	!!----------------------------------------------------------------------
[5836]	47	!! NEMO/OPA 3.7 , NEMO Consortium (2015)
[10763]	48	!! $Id$
[2528]	49	!! Software governed by the CeCILL licence (NEMOGCM/NEMO_CeCILL.txt)
[3]	50	!!----------------------------------------------------------------------
	51	CONTAINS
	52
[508]	53	#if defined key_agrif
[5656]	54	RECURSIVE SUBROUTINE stp( )
[1438]	55	INTEGER :: kstp ! ocean time-step index
[508]	56	#else
[467]	57	SUBROUTINE stp( kstp )
[1438]	58	INTEGER, INTENT(in) :: kstp ! ocean time-step index
[467]	59	#endif
	60	!!----------------------------------------------------------------------
[3]	61	!! * ROUTINE stp *
[3764]	62	!!
[3]	63	!! ** Purpose : - Time stepping of OPA (momentum and active tracer eqs.)
	64	!! - Time stepping of LIM (dynamic and thermodynamic eqs.)
	65	!! - Tme stepping of TRC (passive tracer eqs.)
[3764]	66	!!
	67	!! ** Method : -1- Update forcings and data
	68	!! -2- Update ocean physics
	69	!! -3- Compute the t and s trends
	70	!! -4- Update t and s
[3]	71	!! -5- Compute the momentum trends
	72	!! -6- Update the horizontal velocity
[5836]	73	!! -7- Compute the diagnostics variables (rd,N2, hdiv,w)
[3]	74	!! -8- Outputs and diagnostics
	75	!!----------------------------------------------------------------------
[7753]	76	INTEGER :: ji,jj,jk ! dummy loop indice
[3]	77	INTEGER :: indic ! error indicator if < 0
[4338]	78	INTEGER :: kcall ! optional integer argument (dom_vvl_sf_nxt)
[3]	79	!! ---------------------------------------------------------------------
[392]	80	#if defined key_agrif
[389]	81	kstp = nit000 + Agrif_Nb_Step()
[5836]	82	IF( lk_agrif_debug ) THEN
	83	IF( Agrif_Root() .and. lwp) WRITE(,) '---'
	84	IF(lwp) WRITE(,) 'Grid Number', Agrif_Fixed(),' time step ', kstp, 'int tstep', Agrif_NbStepint()
[5656]	85	ENDIF
[5836]	86	IF( kstp == nit000 + 1 ) lk_agrif_fstep = .FALSE.
[1793]	87	# if defined key_iomput
[5407]	88	IF( Agrif_Nbstepint() == 0 ) CALL iom_swap( cxios_context )
[3764]	89	# endif
	90	#endif
[5836]	91	!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
	92	! update I/O and calendar
	93	!<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<
	94	indic = 0 ! reset to no error condition
	95
	96	IF( kstp == nit000 ) THEN ! initialize IOM context (must be done after nemo_init for AGRIF+XIOS+OASIS)
	97	CALL iom_init( cxios_context ) ! for model grid (including passible AGRIF zoom)
	98	IF( ln_crs ) CALL iom_init( TRIM(cxios_context)//"_crs" ) ! for coarse grid
[4152]	99	ENDIF
[1725]	100	IF( kstp /= nit000 ) CALL day( kstp ) ! Calendar (day was already called at nit000 in day_init)
[5836]	101	CALL iom_setkt( kstp - nit000 + 1, cxios_context ) ! tell IOM we are at time step kstp
	102	IF( ln_crs ) CALL iom_setkt( kstp - nit000 + 1, TRIM(cxios_context)//"_crs" ) ! tell IOM we are at time step kstp
[3]	103
	104	!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
[5836]	105	! Update external forcing (tides, open boundaries, and surface boundary condition (including sea-ice)
[3]	106	!<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<
[7646]	107	IF( ln_tide ) CALL sbc_tide( kstp ) ! update tide potential
[5836]	108	IF( ln_apr_dyn ) CALL sbc_apr ( kstp ) ! atmospheric pressure (NB: call before bdy_dta which needs ssh_ib)
[7646]	109	IF( ln_bdy ) CALL bdy_dta ( kstp, time_offset=+1 ) ! update dynamic & tracer data at open boundaries
[5836]	110	CALL sbc ( kstp ) ! Sea Boundary Condition (including sea-ice)
	111
[3]	112	!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
[5329]	113	! Update stochastic parameters and random T/S fluctuations
	114	!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
[6464]	115	IF( ln_sto_eos ) CALL sto_par( kstp ) ! Stochastic parameters
	116	IF( ln_sto_eos ) CALL sto_pts( tsn ) ! Random T/S fluctuations
[5329]	117
	118	!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
[3294]	119	! Ocean physics update (ua, va, tsa used as workspace)
[3]	120	!<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<
[4990]	121	! THERMODYNAMICS
	122	CALL eos_rab( tsb, rab_b ) ! before local thermal/haline expension ratio at T-points
	123	CALL eos_rab( tsn, rab_n ) ! now local thermal/haline expension ratio at T-points
	124	CALL bn2 ( tsb, rab_b, rn2b ) ! before Brunt-Vaisala frequency
	125	CALL bn2 ( tsn, rab_n, rn2 ) ! now Brunt-Vaisala frequency
[5836]	126
[1492]	127	!
[3764]	128	! VERTICAL PHYSICS
[4369]	129	CALL zdf_bfr( kstp ) ! bottom friction (if quadratic)
[1492]	130	! ! Vertical eddy viscosity and diffusivity coefficients
[7646]	131	IF( lk_zdfric ) CALL zdf_ric ( kstp ) ! Richardson number dependent Kz
	132	IF( lk_zdftke ) CALL zdf_tke ( kstp ) ! TKE closure scheme for Kz
	133	IF( lk_zdfgls ) CALL zdf_gls ( kstp ) ! GLS closure scheme for Kz
	134	IF( ln_zdfqiao ) CALL zdf_qiao( kstp ) ! Qiao vertical mixing
	135	!
	136	IF( lk_zdfcst ) THEN ! Constant Kz (reset avt, avm[uv] to the background value)
[7753]	137	avt (:,:,:) = rn_avt0 * wmask (:,:,:)
	138	avmu(:,:,:) = rn_avm0 * wumask(:,:,:)
	139	avmv(:,:,:) = rn_avm0 * wvmask(:,:,:)
[677]	140	ENDIF
[5836]	141
[1492]	142	IF( ln_rnf_mouth ) THEN ! increase diffusivity at rivers mouths
[7753]	143	DO jk = 2, nkrnf ; avt(:,:,jk) = avt(:,:,jk) + 2._wp * rn_avt_rnf * rnfmsk(:,:) * tmask(:,:,jk) ; END DO
[3]	144	ENDIF
[1492]	145	IF( ln_zdfevd ) CALL zdf_evd( kstp ) ! enhanced vertical eddy diffusivity
[3]	146
[1492]	147	IF( lk_zdftmx ) CALL zdf_tmx( kstp ) ! tidal vertical mixing
[3]	148
[5836]	149	IF( lk_zdfddm ) CALL zdf_ddm( kstp ) ! double diffusive mixing
[3764]	150
[1492]	151	CALL zdf_mxl( kstp ) ! mixed layer depth
[3]	152
[2528]	153	! write TKE or GLS information in the restart file
[1531]	154	IF( lrst_oce .AND. lk_zdftke ) CALL tke_rst( kstp, 'WRITE' )
[2528]	155	IF( lrst_oce .AND. lk_zdfgls ) CALL gls_rst( kstp, 'WRITE' )
[1492]	156	!
[3764]	157	! LATERAL PHYSICS
[1492]	158	!
[5836]	159	IF( l_ldfslp ) THEN ! slope of lateral mixing
[5120]	160	CALL eos( tsb, rhd, gdept_0(:,:,:) ) ! before in situ density
[5836]	161
[5120]	162	IF( ln_zps .AND. .NOT. ln_isfcav) &
	163	& CALL zps_hde ( kstp, jpts, tsb, gtsu, gtsv, & ! Partial steps: before horizontal gradient
	164	& rhd, gru , grv ) ! of t, s, rd at the last ocean level
[5836]	165
[5120]	166	IF( ln_zps .AND. ln_isfcav) &
[6140]	167	& CALL zps_hde_isf( kstp, jpts, tsb, gtsu, gtsv, gtui, gtvi, & ! Partial steps for top cell (ISF)
	168	& rhd, gru , grv , grui, grvi ) ! of t, s, rd at the first ocean level
[5836]	169	IF( ln_traldf_triad ) THEN
	170	CALL ldf_slp_triad( kstp ) ! before slope for triad operator
	171	ELSE
	172	CALL ldf_slp ( kstp, rhd, rn2b ) ! before slope for standard operator
[2528]	173	ENDIF
[1481]	174	ENDIF
[6381]	175	! ! eddy diffusivity coeff.
	176	IF( l_ldftra_time .OR. l_ldfeiv_time ) CALL ldf_tra( kstp ) ! and/or eiv coeff.
	177	IF( l_ldfdyn_time ) CALL ldf_dyn( kstp ) ! eddy viscosity coeff.
[3]	178
[4369]	179	!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
[5930]	180	! Ocean dynamics : hdiv, ssh, e3, u, v, w
[4369]	181	!<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<
[5930]	182
[6140]	183	CALL ssh_nxt ( kstp ) ! after ssh (includes call to div_hor)
	184	IF(.NOT.ln_linssh ) CALL dom_vvl_sf_nxt( kstp ) ! after vertical scale factors
	185	CALL wzv ( kstp ) ! now cross-level velocity
	186	CALL eos ( tsn, rhd, rhop, gdept_n(:,:,:) ) ! now in situ density for hpg computation
	187
[5930]	188	!!jc: fs simplification
[6140]	189	!!jc: lines below are useless if ln_linssh=F. Keep them here (which maintains a bug if ln_linssh=T and ln_zps=T, cf ticket #1636)
[5930]	190	!! but ensures reproductible results
	191	!! with previous versions using split-explicit free surface
[6140]	192	IF( ln_zps .AND. .NOT. ln_isfcav ) &
	193	& CALL zps_hde ( kstp, jpts, tsn, gtsu, gtsv, & ! Partial steps: before horizontal gradient
	194	& rhd, gru , grv ) ! of t, s, rd at the last ocean level
	195	IF( ln_zps .AND. ln_isfcav ) &
	196	& CALL zps_hde_isf( kstp, jpts, tsn, 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
[5930]	198	!!jc: fs simplification
[7753]	199
	200	ua(:,:,:) = 0._wp ! set dynamics trends to zero
	201	va(:,:,:) = 0._wp
[4369]	202
[5930]	203	IF( lk_asminc .AND. ln_asmiau .AND. ln_dyninc ) &
[6381]	204	& CALL dyn_asm_inc ( kstp ) ! apply dynamics assimilation increment
[7646]	205	IF( ln_bdy ) CALL bdy_dyn3d_dmp ( kstp ) ! bdy damping trends
[4369]	206	#if defined key_agrif
[5930]	207	IF(.NOT. Agrif_Root()) &
	208	& CALL Agrif_Sponge_dyn ! momentum sponge
[4369]	209	#endif
[5930]	210	CALL dyn_adv ( kstp ) ! advection (vector or flux form)
	211	CALL dyn_vor ( kstp ) ! vorticity term including Coriolis
	212	CALL dyn_ldf ( kstp ) ! lateral mixing
	213	CALL dyn_hpg ( kstp ) ! horizontal gradient of Hydrostatic pressure
	214	CALL dyn_spg ( kstp ) ! surface pressure gradient
[4369]	215
[5930]	216	! With split-explicit free surface, since now transports have been updated and ssha as well
	217	IF( ln_dynspg_ts ) THEN ! vertical scale factors and vertical velocity need to be updated
	218	CALL div_hor ( kstp ) ! Horizontal divergence (2nd call in time-split case)
[6140]	219	IF(.NOT.ln_linssh) CALL dom_vvl_sf_nxt( kstp, kcall=2 ) ! after vertical scale factors (update depth average component)
[5930]	220	CALL wzv ( kstp ) ! now cross-level velocity
[4369]	221	ENDIF
	222
[5930]	223	CALL dyn_bfr ( kstp ) ! bottom friction
	224	CALL dyn_zdf ( kstp ) ! vertical diffusion
	225
[1482]	226	!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
[6140]	227	! cool skin
	228	!<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<
	229	IF ( ln_diurnal ) CALL stp_diurnal( kstp )
	230
	231	!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
	232	! diagnostics and outputs (ua, va, tsa used as workspace)
[1482]	233	!<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<
[6140]	234	IF( lk_floats ) CALL flo_stp( kstp ) ! drifting Floats
	235	IF( nn_diacfl == 1 ) CALL dia_cfl( kstp ) ! Courant number diagnostics
	236	IF( lk_diahth ) CALL dia_hth( kstp ) ! Thermocline depth (20 degres isotherm depth)
	237	IF( lk_diadct ) CALL dia_dct( kstp ) ! Transports
[7646]	238	CALL dia_ar5( kstp ) ! ar5 diag
[6140]	239	IF( lk_diaharm ) CALL dia_harm( kstp ) ! Tidal harmonic analysis
	240	CALL dia_wri( kstp ) ! ocean model: outputs
[4152]	241	!
[5930]	242	IF( ln_crs ) CALL crs_fld ( kstp ) ! ocean model: online field coarsening & output
[6140]	243
[923]	244	#if defined key_top
[3]	245	!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
	246	! Passive Tracer Model
	247	!<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<
[5930]	248	CALL trc_stp ( kstp ) ! time-stepping
[3]	249	#endif
	250
	251	!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
[5930]	252	! Active tracers
[3]	253	!<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<
[7753]	254	tsa(:,:,:,:) = 0._wp ! set tracer trends to zero
[3]	255
[5656]	256	IF( lk_asminc .AND. ln_asmiau .AND. &
[5930]	257	& ln_trainc ) CALL tra_asm_inc ( kstp ) ! apply tracer assimilation increment
	258	CALL tra_sbc ( kstp ) ! surface boundary condition
	259	IF( ln_traqsr ) CALL tra_qsr ( kstp ) ! penetrative solar radiation qsr
	260	IF( ln_trabbc ) CALL tra_bbc ( kstp ) ! bottom heat flux
	261	IF( lk_trabbl ) CALL tra_bbl ( kstp ) ! advective (and/or diffusive) bottom boundary layer scheme
	262	IF( ln_tradmp ) CALL tra_dmp ( kstp ) ! internal damping trends
[7646]	263	IF( ln_bdy ) CALL bdy_tra_dmp ( kstp ) ! bdy damping trends
[392]	264	#if defined key_agrif
[5930]	265	IF(.NOT. Agrif_Root()) &
	266	& CALL Agrif_Sponge_tra ! tracers sponge
[389]	267	#endif
[5930]	268	CALL tra_adv ( kstp ) ! horizontal & vertical advection
	269	CALL tra_ldf ( kstp ) ! lateral mixing
[1239]	270
[5930]	271	!!gm : why CALL to dia_ptr has been moved here??? (use trends info?)
	272	IF( ln_diaptr ) CALL dia_ptr ! Poleward adv/ldf TRansports diagnostics
[5836]	273	!!gm
[5930]	274	CALL tra_zdf ( kstp ) ! vertical mixing and after tracer fields
	275	IF( ln_zdfnpc ) CALL tra_npc ( kstp ) ! update after fields by non-penetrative convection
[1481]	276
[3]	277	!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
[5930]	278	! Set boundary conditions and Swap
[3]	279	!<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<
[5930]	280	!!jc1: For agrif, it would be much better to finalize tracers/momentum here (e.g. bdy conditions) and move the swap
	281	!! (and time filtering) after Agrif update. Then restart would be done after and would contain updated fields.
	282	!! If so:
	283	!! (i) no need to call agrif update at initialization time
	284	!! (ii) no need to update "before" fields
	285	!!
	286	!! Apart from creating new tra_swp/dyn_swp routines, this however:
	287	!! (i) makes boundary conditions at initialization time computed from updated fields which is not the case between
	288	!! two restarts => restartability issue. One can circumvent this, maybe, by assuming "interface separation",
	289	!! e.g. a shift of the feedback interface inside child domain.
	290	!! (ii) requires that all restart outputs of updated variables by agrif (e.g. passive tracers/tke/barotropic arrays) are done at the same
	291	!! place.
	292	!!
[6140]	293	!!jc2: dynnxt must be the latest call. e3t_b are indeed updated in that routine
[5930]	294	CALL tra_nxt ( kstp ) ! finalize (bcs) tracer fields at next time step and swap
	295	CALL dyn_nxt ( kstp ) ! finalize (bcs) velocities at next time step and swap
	296	CALL ssh_swp ( kstp ) ! swap of sea surface height
[6140]	297	IF(.NOT.ln_linssh) CALL dom_vvl_sf_swp( kstp ) ! swap of vertical scale factors
[5656]	298	!
[7646]	299	IF( ln_diahsb ) CALL dia_hsb( kstp ) ! - ML - global conservation diagnostics
[5836]	300
	301	!!gm : This does not only concern the dynamics ==>>> add a new title
	302	!!gm2: why ouput restart before AGRIF update?
[5930]	303	!!
	304	!!jc: That would be better, but see comment above
	305	!!
[6464]	306	IF( lrst_oce ) CALL rst_write ( kstp ) ! write output ocean restart file
	307	IF( ln_sto_eos ) CALL sto_rst_write( kstp ) ! write restart file for stochastic parameters
[593]	308
[5656]	309	#if defined key_agrif
	310	!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
	311	! AGRIF
	312	!<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<
[5930]	313	CALL Agrif_Integrate_ChildGrids( stp )
[5656]	314
[6140]	315	IF( Agrif_NbStepint() == 0 ) THEN ! AGRIF Update
	316	!!jc in fact update is useless at last time step, but do it for global diagnostics
[5930]	317	CALL Agrif_Update_Tra() ! Update active tracers
	318	CALL Agrif_Update_Dyn() ! Update momentum
[5656]	319	ENDIF
	320	#endif
[6140]	321	IF( ln_diaobs ) CALL dia_obs( kstp ) ! obs-minus-model (assimilation) diagnostics (call after dynamics update)
[2528]	322
[3]	323	!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
[5656]	324	! Control
[3]	325	!<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<
[5930]	326	CALL stp_ctl ( kstp, indic )
[6140]	327	IF( indic < 0 ) THEN
[5930]	328	CALL ctl_stop( 'step: indic < 0' )
	329	CALL dia_wri_state( 'output.abort', kstp )
[10764]	330	CALL ctl_stop('STOP','NEMO failure in stp')
[1482]	331	ENDIF
[6140]	332	IF( kstp == nit000 ) THEN
[4624]	333	CALL iom_close( numror ) ! close input ocean restart file
	334	IF(lwm) CALL FLUSH ( numond ) ! flush output namelist oce
[5836]	335	IF(lwm.AND.numoni /= -1 ) &
	336	& CALL FLUSH ( numoni ) ! flush output namelist ice (if exist)
[4147]	337	ENDIF
[3]	338
[508]	339	!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
[75]	340	! Coupled mode
	341	!<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<
[10764]	342	!IF( lk_oasis ) CALL sbc_cpl_snd( kstp ) ! coupled mode : field exchanges
[508]	343	!
[4153]	344	#if defined key_iomput
	345	IF( kstp == nitend .OR. indic < 0 ) THEN
[5407]	346	CALL iom_context_finalize( cxios_context ) ! needed for XIOS+AGRIF
	347	IF( ln_crs ) CALL iom_context_finalize( trim(cxios_context)//"_crs" ) !
[4152]	348	ENDIF
[4153]	349	#endif
[3769]	350	!
[3294]	351	IF( nn_timing == 1 .AND. kstp == nit000 ) CALL timing_reset
[1481]	352	!
[3]	353	END SUBROUTINE stp
[6140]	354
[3]	355	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/branches/UKMO/r8395_coupling_sequence/NEMOGCM/NEMO/OPA_SRC/step.F90 @ 10764

Download in other formats: