Context Navigation

← Previous Revision
Latest Revision
Next Revision →
Normal
Revision Log

step.F90 @ 8520

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

Download in other formats: