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_r10721_KERNEL-02_Storkey_Coward_IMMERSE_first_steps/src/OCE – NEMO

Context Navigation

source: NEMO/branches/2019/dev_r10721_KERNEL-02_Storkey_Coward_IMMERSE_first_steps/src/OCE/step.F90 @ 10806

Last change on this file since 10806 was 10806, checked in by davestorkey, 5 years ago
2019/dev_r10721_KERNEL-02_Storkey_Coward_IMMERSE_first_steps branch: Latest updates. Make sure all time-dependent 3D variables are converted in previously modified modules.
Property svn:keywords set to `Id`
File size: 21.8 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
[10743]	46	PUBLIC update_pointers ! called by nemo_init
[3]	47
	48	!!----------------------------------------------------------------------
[9598]	49	!! NEMO/OCE 4.0 , NEMO Consortium (2018)
[888]	50	!! $Id$
[10068]	51	!! Software governed by the CeCILL license (see ./LICENSE)
[3]	52	!!----------------------------------------------------------------------
	53	CONTAINS
	54
[508]	55	#if defined key_agrif
[5656]	56	RECURSIVE SUBROUTINE stp( )
[1438]	57	INTEGER :: kstp ! ocean time-step index
[508]	58	#else
[467]	59	SUBROUTINE stp( kstp )
[1438]	60	INTEGER, INTENT(in) :: kstp ! ocean time-step index
[467]	61	#endif
	62	!!----------------------------------------------------------------------
[3]	63	!! * ROUTINE stp *
[3764]	64	!!
[9019]	65	!! ** Purpose : - Time stepping of OPA (momentum and active tracer eqs.)
[9570]	66	!! - Time stepping of SI3 (dynamic and thermodynamic eqs.)
[9019]	67	!! - Time stepping of TRC (passive tracer eqs.)
[3764]	68	!!
	69	!! ** Method : -1- Update forcings and data
	70	!! -2- Update ocean physics
	71	!! -3- Compute the t and s trends
	72	!! -4- Update t and s
[3]	73	!! -5- Compute the momentum trends
	74	!! -6- Update the horizontal velocity
[5836]	75	!! -7- Compute the diagnostics variables (rd,N2, hdiv,w)
[3]	76	!! -8- Outputs and diagnostics
	77	!!----------------------------------------------------------------------
[9019]	78	INTEGER :: ji, jj, jk ! dummy loop indice
	79	INTEGER :: indic ! error indicator if < 0
	80	!!gm kcall can be removed, I guess
	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
[9124]	94	!
	95	IF( ln_timing ) CALL timing_start('stp')
	96	!
[5836]	97	!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
	98	! update I/O and calendar
	99	!<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<
	100	indic = 0 ! reset to no error condition
	101
	102	IF( kstp == nit000 ) THEN ! initialize IOM context (must be done after nemo_init for AGRIF+XIOS+OASIS)
	103	CALL iom_init( cxios_context ) ! for model grid (including passible AGRIF zoom)
	104	IF( ln_crs ) CALL iom_init( TRIM(cxios_context)//"_crs" ) ! for coarse grid
[4152]	105	ENDIF
[1725]	106	IF( kstp /= nit000 ) CALL day( kstp ) ! Calendar (day was already called at nit000 in day_init)
[5836]	107	CALL iom_setkt( kstp - nit000 + 1, cxios_context ) ! tell IOM we are at time step kstp
	108	IF( ln_crs ) CALL iom_setkt( kstp - nit000 + 1, TRIM(cxios_context)//"_crs" ) ! tell IOM we are at time step kstp
[3]	109
	110	!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
[5836]	111	! Update external forcing (tides, open boundaries, and surface boundary condition (including sea-ice)
[3]	112	!<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<
[7646]	113	IF( ln_tide ) CALL sbc_tide( kstp ) ! update tide potential
[5836]	114	IF( ln_apr_dyn ) CALL sbc_apr ( kstp ) ! atmospheric pressure (NB: call before bdy_dta which needs ssh_ib)
[7646]	115	IF( ln_bdy ) CALL bdy_dta ( kstp, time_offset=+1 ) ! update dynamic & tracer data at open boundaries
[5836]	116	CALL sbc ( kstp ) ! Sea Boundary Condition (including sea-ice)
	117
[3]	118	!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
[5329]	119	! Update stochastic parameters and random T/S fluctuations
	120	!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
[6464]	121	IF( ln_sto_eos ) CALL sto_par( kstp ) ! Stochastic parameters
	122	IF( ln_sto_eos ) CALL sto_pts( tsn ) ! Random T/S fluctuations
[5329]	123
	124	!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
[9490]	125	! Ocean physics update
[3]	126	!<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<
[4990]	127	! THERMODYNAMICS
	128	CALL eos_rab( tsb, rab_b ) ! before local thermal/haline expension ratio at T-points
	129	CALL eos_rab( tsn, rab_n ) ! now local thermal/haline expension ratio at T-points
	130	CALL bn2 ( tsb, rab_b, rn2b ) ! before Brunt-Vaisala frequency
	131	CALL bn2 ( tsn, rab_n, rn2 ) ! now Brunt-Vaisala frequency
[5836]	132
[3764]	133	! VERTICAL PHYSICS
[9019]	134	CALL zdf_phy( kstp ) ! vertical physics update (top/bot drag, avt, avs, avm + MLD)
[5836]	135
[3764]	136	! LATERAL PHYSICS
[1492]	137	!
[5836]	138	IF( l_ldfslp ) THEN ! slope of lateral mixing
[5120]	139	CALL eos( tsb, rhd, gdept_0(:,:,:) ) ! before in situ density
[5836]	140
[5120]	141	IF( ln_zps .AND. .NOT. ln_isfcav) &
	142	& CALL zps_hde ( kstp, jpts, tsb, gtsu, gtsv, & ! Partial steps: before horizontal gradient
	143	& rhd, gru , grv ) ! of t, s, rd at the last ocean level
[5836]	144
[5120]	145	IF( ln_zps .AND. ln_isfcav) &
[6140]	146	& CALL zps_hde_isf( kstp, jpts, tsb, gtsu, gtsv, gtui, gtvi, & ! Partial steps for top cell (ISF)
	147	& rhd, gru , grv , grui, grvi ) ! of t, s, rd at the first ocean level
[5836]	148	IF( ln_traldf_triad ) THEN
	149	CALL ldf_slp_triad( kstp ) ! before slope for triad operator
	150	ELSE
	151	CALL ldf_slp ( kstp, rhd, rn2b ) ! before slope for standard operator
[2528]	152	ENDIF
[1481]	153	ENDIF
[6381]	154	! ! eddy diffusivity coeff.
	155	IF( l_ldftra_time .OR. l_ldfeiv_time ) CALL ldf_tra( kstp ) ! and/or eiv coeff.
	156	IF( l_ldfdyn_time ) CALL ldf_dyn( kstp ) ! eddy viscosity coeff.
[3]	157
[4369]	158	!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
[5930]	159	! Ocean dynamics : hdiv, ssh, e3, u, v, w
[4369]	160	!<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<
[5930]	161
[6140]	162	CALL ssh_nxt ( kstp ) ! after ssh (includes call to div_hor)
[10364]	163	IF( .NOT.ln_linssh ) CALL dom_vvl_sf_nxt( kstp ) ! after vertical scale factors
[6140]	164	CALL wzv ( kstp ) ! now cross-level velocity
[10364]	165	IF( ln_zad_Aimp ) CALL wAimp ( kstp ) ! Adaptive-implicit vertical advection partitioning
[6140]	166	CALL eos ( tsn, rhd, rhop, gdept_n(:,:,:) ) ! now in situ density for hpg computation
	167
[5930]	168	!!jc: fs simplification
[6140]	169	!!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]	170	!! but ensures reproductible results
	171	!! with previous versions using split-explicit free surface
[6140]	172	IF( ln_zps .AND. .NOT. ln_isfcav ) &
	173	& CALL zps_hde ( kstp, jpts, tsn, gtsu, gtsv, & ! Partial steps: before horizontal gradient
	174	& rhd, gru , grv ) ! of t, s, rd at the last ocean level
	175	IF( ln_zps .AND. ln_isfcav ) &
	176	& CALL zps_hde_isf( kstp, jpts, tsn, gtsu, gtsv, gtui, gtvi, & ! Partial steps for top cell (ISF)
	177	& rhd, gru , grv , grui, grvi ) ! of t, s, rd at the first ocean level
[5930]	178	!!jc: fs simplification
[7753]	179
[10806]	180	uu(:,:,:,Nrhs) = 0._wp ! set dynamics trends to zero
	181	vv(:,:,:,Nrhs) = 0._wp
[4369]	182
[5930]	183	IF( lk_asminc .AND. ln_asmiau .AND. ln_dyninc ) &
[6381]	184	& CALL dyn_asm_inc ( kstp ) ! apply dynamics assimilation increment
[7646]	185	IF( ln_bdy ) CALL bdy_dyn3d_dmp ( kstp ) ! bdy damping trends
[4369]	186	#if defined key_agrif
[5930]	187	IF(.NOT. Agrif_Root()) &
	188	& CALL Agrif_Sponge_dyn ! momentum sponge
[4369]	189	#endif
[10789]	190	CALL dyn_adv ( kstp, Nm1, Nnn, uu(:,:,:,Nrhs), vv(:,:,:,Nrhs) ) ! advection (vector or flux form)
	191	CALL dyn_vor ( kstp, Nnn, uu(:,:,:,Nrhs), vv(:,:,:,Nrhs) ) ! vorticity term including Coriolis
[10806]	192	CALL dyn_ldf ( kstp, Nm1, Nnn, uu(:,:,:,Nrhs), vv(:,:,:,Nrhs) ) ! lateral mixing
[9019]	193	IF( ln_zdfosm ) CALL dyn_osm ( kstp ) ! OSMOSIS non-local velocity fluxes
[5930]	194	CALL dyn_hpg ( kstp ) ! horizontal gradient of Hydrostatic pressure
	195	CALL dyn_spg ( kstp ) ! surface pressure gradient
[4369]	196
[5930]	197	! With split-explicit free surface, since now transports have been updated and ssha as well
	198	IF( ln_dynspg_ts ) THEN ! vertical scale factors and vertical velocity need to be updated
	199	CALL div_hor ( kstp ) ! Horizontal divergence (2nd call in time-split case)
[6140]	200	IF(.NOT.ln_linssh) CALL dom_vvl_sf_nxt( kstp, kcall=2 ) ! after vertical scale factors (update depth average component)
[5930]	201	CALL wzv ( kstp ) ! now cross-level velocity
[10364]	202	IF( ln_zad_Aimp ) CALL wAimp ( kstp ) ! Adaptive-implicit vertical advection partitioning
[4369]	203	ENDIF
[9019]	204
[5930]	205	CALL dyn_zdf ( kstp ) ! vertical diffusion
	206
[1482]	207	!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
[6140]	208	! cool skin
	209	!<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<
	210	IF ( ln_diurnal ) CALL stp_diurnal( kstp )
	211
	212	!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
[9490]	213	! diagnostics and outputs
[1482]	214	!<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<
[9490]	215	IF( lk_floats ) CALL flo_stp ( kstp ) ! drifting Floats
	216	IF( ln_diacfl ) CALL dia_cfl ( kstp ) ! Courant number diagnostics
	217	IF( lk_diahth ) CALL dia_hth ( kstp ) ! Thermocline depth (20 degres isotherm depth)
	218	IF( lk_diadct ) CALL dia_dct ( kstp ) ! Transports
	219	CALL dia_ar5 ( kstp ) ! ar5 diag
[6140]	220	IF( lk_diaharm ) CALL dia_harm( kstp ) ! Tidal harmonic analysis
[9490]	221	CALL dia_wri ( kstp ) ! ocean model: outputs
[4152]	222	!
[5930]	223	IF( ln_crs ) CALL crs_fld ( kstp ) ! ocean model: online field coarsening & output
[6140]	224
[923]	225	#if defined key_top
[3]	226	!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
	227	! Passive Tracer Model
	228	!<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<
[5930]	229	CALL trc_stp ( kstp ) ! time-stepping
[3]	230	#endif
	231
	232	!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
[5930]	233	! Active tracers
[3]	234	!<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<
[10806]	235	ts(:,:,:,:,Nrhs) = 0._wp ! set tracer trends to zero
[3]	236
[5656]	237	IF( lk_asminc .AND. ln_asmiau .AND. &
[5930]	238	& ln_trainc ) CALL tra_asm_inc ( kstp ) ! apply tracer assimilation increment
[10806]	239	CALL tra_sbc ( kstp, Nnn, ts(:,:,:,:,Nrhs) ) ! surface boundary condition
	240	IF( ln_traqsr ) CALL tra_qsr ( kstp, Nnn, Nnn_2lev, ts(:,:,:,:,Nrhs) ) ! penetrative solar radiation qsr
	241	IF( ln_trabbc ) CALL tra_bbc ( kstp, Nnn, ts(:,:,:,:,Nrhs) ) ! bottom heat flux
	242	IF( ln_trabbl ) CALL tra_bbl ( kstp, Nm1, Nnn, Nnn_2lev, ts(:,:,:,:,Nrhs) ) ! advective (and/or diffusive) bottom boundary layer scheme
	243	IF( ln_tradmp ) CALL tra_dmp ( kstp, Nm1, Nnn_2lev, ts(:,:,:,:,Nrhs) ) ! internal damping trends
[7646]	244	IF( ln_bdy ) CALL bdy_tra_dmp ( kstp ) ! bdy damping trends
[392]	245	#if defined key_agrif
[5930]	246	IF(.NOT. Agrif_Root()) &
	247	& CALL Agrif_Sponge_tra ! tracers sponge
[389]	248	#endif
[10806]	249	CALL tra_adv ( kstp, Nm1, Nnn, Np1, Nnn_2lev, ts(:,:,:,:,Nrhs) ) ! horizontal & vertical advection
[9019]	250	IF( ln_zdfosm ) CALL tra_osm ( kstp ) ! OSMOSIS non-local tracer fluxes
	251	IF( lrst_oce .AND. ln_zdfosm ) &
	252	& CALL osm_rst( kstp, 'WRITE' )! write OSMOSIS outputs + wn (so must do here) to restarts
[10806]	253	CALL tra_ldf ( kstp, Nm1, Nnn, Nnn_2lev, ts(:,:,:,:,Nrhs) ) ! lateral mixing
[1239]	254
[5930]	255	!!gm : why CALL to dia_ptr has been moved here??? (use trends info?)
	256	IF( ln_diaptr ) CALL dia_ptr ! Poleward adv/ldf TRansports diagnostics
[5836]	257	!!gm
[5930]	258	CALL tra_zdf ( kstp ) ! vertical mixing and after tracer fields
	259	IF( ln_zdfnpc ) CALL tra_npc ( kstp ) ! update after fields by non-penetrative convection
[1481]	260
[3]	261	!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
[5930]	262	! Set boundary conditions and Swap
[3]	263	!<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<
[5930]	264	!!jc1: For agrif, it would be much better to finalize tracers/momentum here (e.g. bdy conditions) and move the swap
	265	!! (and time filtering) after Agrif update. Then restart would be done after and would contain updated fields.
	266	!! If so:
	267	!! (i) no need to call agrif update at initialization time
	268	!! (ii) no need to update "before" fields
	269	!!
	270	!! Apart from creating new tra_swp/dyn_swp routines, this however:
	271	!! (i) makes boundary conditions at initialization time computed from updated fields which is not the case between
	272	!! two restarts => restartability issue. One can circumvent this, maybe, by assuming "interface separation",
	273	!! e.g. a shift of the feedback interface inside child domain.
	274	!! (ii) requires that all restart outputs of updated variables by agrif (e.g. passive tracers/tke/barotropic arrays) are done at the same
	275	!! place.
	276	!!
[6140]	277	!!jc2: dynnxt must be the latest call. e3t_b are indeed updated in that routine
[5930]	278	CALL tra_nxt ( kstp ) ! finalize (bcs) tracer fields at next time step and swap
[10799]	279	CALL dyn_nxt ( kstp ) !
[5930]	280	CALL ssh_swp ( kstp ) ! swap of sea surface height
[10799]	281	IF(.NOT.ln_linssh) CALL dom_vvl_sf_swp( kstp ) ! interpolate vertical scale factors for Nnn time level
[5656]	282	!
[9210]	283	IF( ln_diahsb ) CALL dia_hsb ( kstp ) ! - ML - global conservation diagnostics
[5836]	284
	285	!!gm : This does not only concern the dynamics ==>>> add a new title
	286	!!gm2: why ouput restart before AGRIF update?
[5930]	287	!!
	288	!!jc: That would be better, but see comment above
	289	!!
[9210]	290	IF( lrst_oce ) CALL rst_write ( kstp ) ! write output ocean restart file
	291	IF( ln_sto_eos ) CALL sto_rst_write( kstp ) ! write restart file for stochastic parameters
[593]	292
[5656]	293	#if defined key_agrif
	294	!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
	295	! AGRIF
	296	!<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<
[9485]	297	CALL Agrif_Integrate_ChildGrids( stp ) ! allows to finish all the Child Grids before updating
[5656]	298
[9780]	299	IF( Agrif_NbStepint() == 0 ) CALL Agrif_update_all( ) ! Update all components
[9485]	300	#endif
[9210]	301	IF( ln_diaobs ) CALL dia_obs ( kstp ) ! obs-minus-model (assimilation) diagnostics (call after dynamics update)
[2528]	302
[3]	303	!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
[5656]	304	! Control
[3]	305	!<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<
[9210]	306	CALL stp_ctl ( kstp, indic )
	307
	308	IF( kstp == nit000 ) THEN ! 1st time step only
	309	CALL iom_close( numror ) ! close input ocean restart file
	310	IF(lwm) CALL FLUSH ( numond ) ! flush output namelist oce
	311	IF(lwm .AND. numoni /= -1 ) CALL FLUSH ( numoni ) ! flush output namelist ice (if exist)
[1482]	312	ENDIF
[3]	313
[508]	314	!>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
[75]	315	! Coupled mode
	316	!<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<
[5836]	317	!!gm why lk_oasis and not lk_cpl ????
[5930]	318	IF( lk_oasis ) CALL sbc_cpl_snd( kstp ) ! coupled mode : field exchanges
[508]	319	!
[4153]	320	#if defined key_iomput
	321	IF( kstp == nitend .OR. indic < 0 ) THEN
[5407]	322	CALL iom_context_finalize( cxios_context ) ! needed for XIOS+AGRIF
[9367]	323	IF(lrxios) CALL iom_context_finalize( crxios_context )
[5407]	324	IF( ln_crs ) CALL iom_context_finalize( trim(cxios_context)//"_crs" ) !
[4152]	325	ENDIF
[4153]	326	#endif
[3769]	327	!
[9210]	328	IF( ln_timing ) CALL timing_stop('stp')
[1481]	329	!
[3]	330	END SUBROUTINE stp
[6140]	331
[10743]	332	SUBROUTINE update_pointers
	333	!!----------------------------------------------------------------------
	334	!! * ROUTINE update_pointers *
	335	!!
	336	!! ** Purpose : Associate temporary pointer arrays.
	337	!! For IMMERSE development phase only - to be deleted
	338	!!
	339	!! ** Method :
	340	!!----------------------------------------------------------------------
	341
	342	ub => uu(:,:,:,Nm1); un => uu(:,:,:,Nnn); ua => uu(:,:,:,Np1)
	343	vb => vv(:,:,:,Nm1); vn => vv(:,:,:,Nnn); va => vv(:,:,:,Np1)
[10789]	344	wn => ww(:,:,:)
[10806]	345	hdivn => hdiv(:,:,:)
	346	rab_b => r_ab (:,:,:,:,Nm1_2lev); rab_n => r_ab (:,:,:,:,Nnn_2lev)
	347	rn2b => r_n2 (:,:,:,Nm1_2lev) ; rn2 => r_n2 (:,:,:,Nnn_2lev)
[10743]	348
[10802]	349	tsb => ts(:,:,:,:,Nm1); tsn => ts(:,:,:,:,Nnn); tsa => ts(:,:,:,:,Np1)
	350
[10743]	351	e3t_b => e3t(:,:,:,Nm1); e3t_n => e3t(:,:,:,Nnn); e3t_a => e3t(:,:,:,Np1)
	352	e3u_b => e3u(:,:,:,Nm1); e3u_n => e3u(:,:,:,Nnn); e3u_a => e3u(:,:,:,Np1)
	353	e3v_b => e3v(:,:,:,Nm1); e3v_n => e3v(:,:,:,Nnn); e3v_a => e3v(:,:,:,Np1)
	354
	355	e3f_n => e3f(:,:,:)
	356
	357	e3w_b => e3w (:,:,:,Nm1_2lev); e3w_n => e3w (:,:,:,Nnn_2lev)
	358	e3uw_b => e3uw(:,:,:,Nm1_2lev); e3uw_n => e3uw(:,:,:,Nnn_2lev)
	359	e3vw_b => e3vw(:,:,:,Nm1_2lev); e3vw_n => e3vw(:,:,:,Nnn_2lev)
	360
[10787]	361	gdept_b => gdept(:,:,:,Nm1_2lev); gdept_n => gdept(:,:,:,Nnn_2lev)
	362	gdepw_b => gdepw(:,:,:,Nm1_2lev); gdepw_n => gdepw(:,:,:,Nnn_2lev)
	363	gde3w_n => gde3w(:,:,:)
	364
[10743]	365	END SUBROUTINE update_pointers
	366
[9019]	367	!!======================================================================
[3]	368	END MODULE step

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

Download in other formats: