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/2019/dev_r12072_MERGE_OPTION2_2019/src/OCE – NEMO

Context Navigation

source: NEMO/branches/2019/dev_r12072_MERGE_OPTION2_2019/src/OCE/step.F90 @ 12256

Last change on this file since 12256 was 12256, checked in by smasson, 4 years ago
rev12232_dev_r12072_MERGE_OPTION2_2019: bugfix in fldread offset in bdy, see #2263. all sette tests ok
Property svn:keywords set to `Id`
File size: 19.5 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)
[3]	33	!!----------------------------------------------------------------------
[888]	34
	35	!!----------------------------------------------------------------------
[2528]	36	!! stp : OPA system time-stepping
[3]	37	!!----------------------------------------------------------------------
[3764]	38	USE step_oce ! time stepping definition modules
[6140]	39	!
	40	USE iom ! xIOs server
[389]	41
[3]	42	IMPLICIT NONE
	43	PRIVATE
	44
[5836]	45	PUBLIC stp ! called by nemogcm.F90
[3]	46
	47	!!----------------------------------------------------------------------
[9598]	48	!! NEMO/OCE 4.0 , NEMO Consortium (2018)
[888]	49	!! $Id$
[10068]	50	!! Software governed by the CeCILL license (see ./LICENSE)
[3]	51	!!----------------------------------------------------------------------
	52	CONTAINS
	53
[508]	54	#if defined key_agrif
[5656]	55	RECURSIVE SUBROUTINE stp( )
[1438]	56	INTEGER :: kstp ! ocean time-step index
[508]	57	#else
[467]	58	SUBROUTINE stp( kstp )
[1438]	59	INTEGER, INTENT(in) :: kstp ! ocean time-step index
[467]	60	#endif
	61	!!----------------------------------------------------------------------
[3]	62	!! * ROUTINE stp *
[3764]	63	!!
[9019]	64	!! ** Purpose : - Time stepping of OPA (momentum and active tracer eqs.)
[9570]	65	!! - Time stepping of SI3 (dynamic and thermodynamic eqs.)
[9019]	66	!! - Time stepping of TRC (passive tracer eqs.)
[3764]	67	!!
	68	!! ** Method : -1- Update forcings and data
	69	!! -2- Update ocean physics
	70	!! -3- Compute the t and s trends
	71	!! -4- Update t and s
[3]	72	!! -5- Compute the momentum trends
	73	!! -6- Update the horizontal velocity
[5836]	74	!! -7- Compute the diagnostics variables (rd,N2, hdiv,w)
[3]	75	!! -8- Outputs and diagnostics
	76	!!----------------------------------------------------------------------
[9019]	77	INTEGER :: ji, jj, jk ! dummy loop indice
	78	INTEGER :: indic ! error indicator if < 0
	79	!!gm kcall can be removed, I guess
	80	INTEGER :: kcall ! optional integer argument (dom_vvl_sf_nxt)
[3]	81	!! ---------------------------------------------------------------------
[392]	82	#if defined key_agrif
[389]	83	kstp = nit000 + Agrif_Nb_Step()
[5836]	84	IF( lk_agrif_debug ) THEN
	85	IF( Agrif_Root() .and. lwp) WRITE(,) '---'
	86	IF(lwp) WRITE(,) 'Grid Number', Agrif_Fixed(),' time step ', kstp, 'int tstep', Agrif_NbStepint()
[5656]	87	ENDIF
[5836]	88	IF( kstp == nit000 + 1 ) lk_agrif_fstep = .FALSE.
[1793]	89	# if defined key_iomput
[5407]	90	IF( Agrif_Nbstepint() == 0 ) CALL iom_swap( cxios_context )
[3764]	91	# endif
	92	#endif
[9124]	93	!
	94	IF( ln_timing ) CALL timing_start('stp')
	95	!
[5836]	96	!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
	97	! update I/O and calendar
	98	!<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<
	99	indic = 0 ! reset to no error condition
	100
	101	IF( kstp == nit000 ) THEN ! initialize IOM context (must be done after nemo_init for AGRIF+XIOS+OASIS)
[12184]	102	CALL iom_init( cxios_context, ld_closedef=.FALSE. ) ! for model grid (including passible AGRIF zoom)
	103	IF( lk_diamlr ) CALL dia_mlr_iom_init ! with additional setup for multiple-linear-regression analysis
	104	CALL iom_init_closedef
[5836]	105	IF( ln_crs ) CALL iom_init( TRIM(cxios_context)//"_crs" ) ! for coarse grid
[4152]	106	ENDIF
[1725]	107	IF( kstp /= nit000 ) CALL day( kstp ) ! Calendar (day was already called at nit000 in day_init)
[5836]	108	CALL iom_setkt( kstp - nit000 + 1, cxios_context ) ! tell IOM we are at time step kstp
	109	IF( ln_crs ) CALL iom_setkt( kstp - nit000 + 1, TRIM(cxios_context)//"_crs" ) ! tell IOM we are at time step kstp
[3]	110
	111	!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
[5836]	112	! Update external forcing (tides, open boundaries, and surface boundary condition (including sea-ice)
[3]	113	!<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<
[12184]	114	IF( ln_tide ) CALL tide_update( kstp ) ! update tide potential
[5836]	115	IF( ln_apr_dyn ) CALL sbc_apr ( kstp ) ! atmospheric pressure (NB: call before bdy_dta which needs ssh_ib)
[12256]	116	IF( ln_bdy ) CALL bdy_dta ( kstp ) ! update dynamic & tracer data at open boundaries
[12166]	117	IF( ln_isf ) CALL isf_stp ( kstp ) ! ice shelf/ocean boundary condition
[5836]	118	CALL sbc ( kstp ) ! Sea Boundary Condition (including sea-ice)
	119
[3]	120	!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
[5329]	121	! Update stochastic parameters and random T/S fluctuations
	122	!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
[6464]	123	IF( ln_sto_eos ) CALL sto_par( kstp ) ! Stochastic parameters
	124	IF( ln_sto_eos ) CALL sto_pts( tsn ) ! Random T/S fluctuations
[5329]	125
	126	!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
[9490]	127	! Ocean physics update
[3]	128	!<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<
[4990]	129	! THERMODYNAMICS
	130	CALL eos_rab( tsb, rab_b ) ! before local thermal/haline expension ratio at T-points
	131	CALL eos_rab( tsn, rab_n ) ! now local thermal/haline expension ratio at T-points
	132	CALL bn2 ( tsb, rab_b, rn2b ) ! before Brunt-Vaisala frequency
	133	CALL bn2 ( tsn, rab_n, rn2 ) ! now Brunt-Vaisala frequency
[5836]	134
[3764]	135	! VERTICAL PHYSICS
[9019]	136	CALL zdf_phy( kstp ) ! vertical physics update (top/bot drag, avt, avs, avm + MLD)
[5836]	137
[3764]	138	! LATERAL PHYSICS
[1492]	139	!
[5836]	140	IF( l_ldfslp ) THEN ! slope of lateral mixing
[5120]	141	CALL eos( tsb, rhd, gdept_0(:,:,:) ) ! before in situ density
[5836]	142
[5120]	143	IF( ln_zps .AND. .NOT. ln_isfcav) &
	144	& CALL zps_hde ( kstp, jpts, tsb, gtsu, gtsv, & ! Partial steps: before horizontal gradient
	145	& rhd, gru , grv ) ! of t, s, rd at the last ocean level
[5836]	146
[5120]	147	IF( ln_zps .AND. ln_isfcav) &
[6140]	148	& CALL zps_hde_isf( kstp, jpts, tsb, gtsu, gtsv, gtui, gtvi, & ! Partial steps for top cell (ISF)
	149	& rhd, gru , grv , grui, grvi ) ! of t, s, rd at the first ocean level
[5836]	150	IF( ln_traldf_triad ) THEN
	151	CALL ldf_slp_triad( kstp ) ! before slope for triad operator
	152	ELSE
	153	CALL ldf_slp ( kstp, rhd, rn2b ) ! before slope for standard operator
[2528]	154	ENDIF
[1481]	155	ENDIF
[6381]	156	! ! eddy diffusivity coeff.
	157	IF( l_ldftra_time .OR. l_ldfeiv_time ) CALL ldf_tra( kstp ) ! and/or eiv coeff.
	158	IF( l_ldfdyn_time ) CALL ldf_dyn( kstp ) ! eddy viscosity coeff.
[4369]	159	!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
[5930]	160	! Ocean dynamics : hdiv, ssh, e3, u, v, w
[4369]	161	!<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<
[5930]	162
[6140]	163	CALL ssh_nxt ( kstp ) ! after ssh (includes call to div_hor)
[10364]	164	IF( .NOT.ln_linssh ) CALL dom_vvl_sf_nxt( kstp ) ! after vertical scale factors
[6140]	165	CALL wzv ( kstp ) ! now cross-level velocity
[10364]	166	IF( ln_zad_Aimp ) CALL wAimp ( kstp ) ! Adaptive-implicit vertical advection partitioning
[6140]	167	CALL eos ( tsn, rhd, rhop, gdept_n(:,:,:) ) ! now in situ density for hpg computation
	168
[7753]	169
	170	ua(:,:,:) = 0._wp ! set dynamics trends to zero
	171	va(:,:,:) = 0._wp
[4369]	172
[5930]	173	IF( lk_asminc .AND. ln_asmiau .AND. ln_dyninc ) &
[6381]	174	& CALL dyn_asm_inc ( kstp ) ! apply dynamics assimilation increment
[7646]	175	IF( ln_bdy ) CALL bdy_dyn3d_dmp ( kstp ) ! bdy damping trends
[4369]	176	#if defined key_agrif
[5930]	177	IF(.NOT. Agrif_Root()) &
	178	& CALL Agrif_Sponge_dyn ! momentum sponge
[4369]	179	#endif
[5930]	180	CALL dyn_adv ( kstp ) ! advection (vector or flux form)
	181	CALL dyn_vor ( kstp ) ! vorticity term including Coriolis
	182	CALL dyn_ldf ( kstp ) ! lateral mixing
[9019]	183	IF( ln_zdfosm ) CALL dyn_osm ( kstp ) ! OSMOSIS non-local velocity fluxes
[5930]	184	CALL dyn_hpg ( kstp ) ! horizontal gradient of Hydrostatic pressure
	185	CALL dyn_spg ( kstp ) ! surface pressure gradient
[4369]	186
[5930]	187	! With split-explicit free surface, since now transports have been updated and ssha as well
	188	IF( ln_dynspg_ts ) THEN ! vertical scale factors and vertical velocity need to be updated
	189	CALL div_hor ( kstp ) ! Horizontal divergence (2nd call in time-split case)
[6140]	190	IF(.NOT.ln_linssh) CALL dom_vvl_sf_nxt( kstp, kcall=2 ) ! after vertical scale factors (update depth average component)
[11287]	191	ENDIF
[12109]	192	CALL dyn_zdf ( kstp ) ! vertical diffusion
	193	IF( ln_dynspg_ts ) THEN ! vertical scale factors and vertical velocity need to be updated
[5930]	194	CALL wzv ( kstp ) ! now cross-level velocity
[10364]	195	IF( ln_zad_Aimp ) CALL wAimp ( kstp ) ! Adaptive-implicit vertical advection partitioning
[4369]	196	ENDIF
[12109]	197
[5930]	198
[1482]	199	!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
[6140]	200	! cool skin
	201	!<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<
	202	IF ( ln_diurnal ) CALL stp_diurnal( kstp )
	203
	204	!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
[9490]	205	! diagnostics and outputs
[1482]	206	!<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<
[11536]	207	IF( ln_floats ) CALL flo_stp ( kstp ) ! drifting Floats
[9490]	208	IF( ln_diacfl ) CALL dia_cfl ( kstp ) ! Courant number diagnostics
[12109]	209	CALL dia_hth ( kstp ) ! Thermocline depth (20 degres isotherm depth)
[11536]	210	IF( ln_diadct ) CALL dia_dct ( kstp ) ! Transports
[9490]	211	CALL dia_ar5 ( kstp ) ! ar5 diag
[12109]	212	CALL dia_ptr ( kstp ) ! Poleward adv/ldf TRansports diagnostics
[9490]	213	CALL dia_wri ( kstp ) ! ocean model: outputs
[4152]	214	!
[12184]	215	IF( lk_diadetide ) CALL dia_detide( kstp ) ! Weights computation for daily detiding of model diagnostics
	216	IF( lk_diamlr ) CALL dia_mlr ! Update time used in multiple-linear-regression analysis
[5930]	217	IF( ln_crs ) CALL crs_fld ( kstp ) ! ocean model: online field coarsening & output
[6140]	218
[923]	219	#if defined key_top
[3]	220	!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
	221	! Passive Tracer Model
	222	!<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<
[5930]	223	CALL trc_stp ( kstp ) ! time-stepping
[3]	224	#endif
	225
	226	!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
[5930]	227	! Active tracers
[3]	228	!<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<
[7753]	229	tsa(:,:,:,:) = 0._wp ! set tracer trends to zero
[3]	230
[5656]	231	IF( lk_asminc .AND. ln_asmiau .AND. &
[5930]	232	& ln_trainc ) CALL tra_asm_inc ( kstp ) ! apply tracer assimilation increment
	233	CALL tra_sbc ( kstp ) ! surface boundary condition
	234	IF( ln_traqsr ) CALL tra_qsr ( kstp ) ! penetrative solar radiation qsr
[12166]	235	IF( ln_isf ) CALL tra_isf ( kstp ) ! ice shelf heat flux
[5930]	236	IF( ln_trabbc ) CALL tra_bbc ( kstp ) ! bottom heat flux
[9019]	237	IF( ln_trabbl ) CALL tra_bbl ( kstp ) ! advective (and/or diffusive) bottom boundary layer scheme
[5930]	238	IF( ln_tradmp ) CALL tra_dmp ( kstp ) ! internal damping trends
[7646]	239	IF( ln_bdy ) CALL bdy_tra_dmp ( kstp ) ! bdy damping trends
[392]	240	#if defined key_agrif
[5930]	241	IF(.NOT. Agrif_Root()) &
	242	& CALL Agrif_Sponge_tra ! tracers sponge
[389]	243	#endif
[5930]	244	CALL tra_adv ( kstp ) ! horizontal & vertical advection
[9019]	245	IF( ln_zdfosm ) CALL tra_osm ( kstp ) ! OSMOSIS non-local tracer fluxes
	246	IF( lrst_oce .AND. ln_zdfosm ) &
	247	& CALL osm_rst( kstp, 'WRITE' )! write OSMOSIS outputs + wn (so must do here) to restarts
[5930]	248	CALL tra_ldf ( kstp ) ! lateral mixing
[1239]	249
[5930]	250	CALL tra_zdf ( kstp ) ! vertical mixing and after tracer fields
	251	IF( ln_zdfnpc ) CALL tra_npc ( kstp ) ! update after fields by non-penetrative convection
[1481]	252
[3]	253	!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
[5930]	254	! Set boundary conditions and Swap
[3]	255	!<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<
[5930]	256	!!jc1: For agrif, it would be much better to finalize tracers/momentum here (e.g. bdy conditions) and move the swap
	257	!! (and time filtering) after Agrif update. Then restart would be done after and would contain updated fields.
	258	!! If so:
	259	!! (i) no need to call agrif update at initialization time
	260	!! (ii) no need to update "before" fields
	261	!!
	262	!! Apart from creating new tra_swp/dyn_swp routines, this however:
	263	!! (i) makes boundary conditions at initialization time computed from updated fields which is not the case between
	264	!! two restarts => restartability issue. One can circumvent this, maybe, by assuming "interface separation",
	265	!! e.g. a shift of the feedback interface inside child domain.
	266	!! (ii) requires that all restart outputs of updated variables by agrif (e.g. passive tracers/tke/barotropic arrays) are done at the same
	267	!! place.
	268	!!
[6140]	269	!!jc2: dynnxt must be the latest call. e3t_b are indeed updated in that routine
[5930]	270	CALL tra_nxt ( kstp ) ! finalize (bcs) tracer fields at next time step and swap
[9250]	271	CALL dyn_nxt ( kstp ) ! finalize (bcs) velocities at next time step and swap (always called after tra_nxt)
[5930]	272	CALL ssh_swp ( kstp ) ! swap of sea surface height
[6140]	273	IF(.NOT.ln_linssh) CALL dom_vvl_sf_swp( kstp ) ! swap of vertical scale factors
[5656]	274	!
[9210]	275	IF( ln_diahsb ) CALL dia_hsb ( kstp ) ! - ML - global conservation diagnostics
[5836]	276
	277	!!gm : This does not only concern the dynamics ==>>> add a new title
	278	!!gm2: why ouput restart before AGRIF update?
[5930]	279	!!
	280	!!jc: That would be better, but see comment above
	281	!!
[9210]	282	IF( lrst_oce ) CALL rst_write ( kstp ) ! write output ocean restart file
	283	IF( ln_sto_eos ) CALL sto_rst_write( kstp ) ! write restart file for stochastic parameters
[593]	284
[5656]	285	#if defined key_agrif
	286	!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
	287	! AGRIF
	288	!<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<
[9485]	289	CALL Agrif_Integrate_ChildGrids( stp ) ! allows to finish all the Child Grids before updating
[5656]	290
[9780]	291	IF( Agrif_NbStepint() == 0 ) CALL Agrif_update_all( ) ! Update all components
[9485]	292	#endif
[9210]	293	IF( ln_diaobs ) CALL dia_obs ( kstp ) ! obs-minus-model (assimilation) diagnostics (call after dynamics update)
[2528]	294
[3]	295	!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
[5656]	296	! Control
[3]	297	!<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<
[9210]	298	CALL stp_ctl ( kstp, indic )
	299
	300	IF( kstp == nit000 ) THEN ! 1st time step only
	301	CALL iom_close( numror ) ! close input ocean restart file
	302	IF(lwm) CALL FLUSH ( numond ) ! flush output namelist oce
	303	IF(lwm .AND. numoni /= -1 ) CALL FLUSH ( numoni ) ! flush output namelist ice (if exist)
[1482]	304	ENDIF
[3]	305
[508]	306	!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
[75]	307	! Coupled mode
	308	!<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<
[5836]	309	!!gm why lk_oasis and not lk_cpl ????
[5930]	310	IF( lk_oasis ) CALL sbc_cpl_snd( kstp ) ! coupled mode : field exchanges
[508]	311	!
[4153]	312	#if defined key_iomput
	313	IF( kstp == nitend .OR. indic < 0 ) THEN
[5407]	314	CALL iom_context_finalize( cxios_context ) ! needed for XIOS+AGRIF
[9367]	315	IF(lrxios) CALL iom_context_finalize( crxios_context )
[5407]	316	IF( ln_crs ) CALL iom_context_finalize( trim(cxios_context)//"_crs" ) !
[4152]	317	ENDIF
[4153]	318	#endif
[3769]	319	!
[9210]	320	IF( ln_timing ) CALL timing_stop('stp')
[1481]	321	!
[3]	322	END SUBROUTINE stp
[6140]	323
[9019]	324	!!======================================================================
[3]	325	END MODULE step

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

Download in other formats: