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.

sshwzv.F90 in NEMO/trunk/src/OCE/DYN – NEMO

Context Navigation

source: NEMO/trunk/src/OCE/DYN/sshwzv.F90 @ 12489

Last change on this file since 12489 was 12489, checked in by davestorkey, 4 years ago

Preparation for new timestepping scheme #2390.
Main changes:

Initial euler timestep now handled in stp and not in TRA/DYN routines.
Renaming of all timestep parameters. In summary, the namelist parameter is now rn_Dt and the current timestep is rDt (and rDt_ice, rDt_trc etc).
Renaming of a few miscellaneous parameters, eg. atfp -> rn_atfp (namelist parameter used everywhere) and rau0 -> rho0.

This version gives bit-comparable results to the previous version of the trunk.

Property svn:keywords set to Id

File size: 18.8 KB

Rev	Line
[1565]	1	MODULE sshwzv
[3]	2	!!==============================================================================
[1438]	3	!! * MODULE sshwzv *
	4	!! Ocean dynamics : sea surface height and vertical velocity
[3]	5	!!==============================================================================
[1438]	6	!! History : 3.1 ! 2009-02 (G. Madec, M. Leclair) Original code
[2528]	7	!! 3.3 ! 2010-04 (M. Leclair, G. Madec) modified LF-RA
	8	!! - ! 2010-05 (K. Mogensen, A. Weaver, M. Martin, D. Lea) Assimilation interface
	9	!! - ! 2010-09 (D.Storkey and E.O'Dea) bug fixes for BDY module
[4292]	10	!! 3.3 ! 2011-10 (M. Leclair) split former ssh_wzv routine and remove all vvl related work
[11414]	11	!! 4.0 ! 2018-12 (A. Coward) add mixed implicit/explicit advection
[12377]	12	!! 4.1 ! 2019-08 (A. Coward, D. Storkey) Rename ssh_nxt -> ssh_atf. Now only does time filtering.
[3]	13	!!----------------------------------------------------------------------
[1438]	14
[3]	15	!!----------------------------------------------------------------------
[6140]	16	!! ssh_nxt : after ssh
[12377]	17	!! ssh_atf : time filter the ssh arrays
[6140]	18	!! wzv : compute now vertical velocity
[1438]	19	!!----------------------------------------------------------------------
[6140]	20	USE oce ! ocean dynamics and tracers variables
[12377]	21	USE isf_oce ! ice shelf
[6140]	22	USE dom_oce ! ocean space and time domain variables
	23	USE sbc_oce ! surface boundary condition: ocean
	24	USE domvvl ! Variable volume
	25	USE divhor ! horizontal divergence
	26	USE phycst ! physical constants
[9019]	27	USE bdy_oce , ONLY : ln_bdy, bdytmask ! Open BounDarY
[6140]	28	USE bdydyn2d ! bdy_ssh routine
[2528]	29	#if defined key_agrif
[9570]	30	USE agrif_oce_interp
[2528]	31	#endif
[6140]	32	!
[10364]	33	USE iom
[6140]	34	USE in_out_manager ! I/O manager
	35	USE restart ! only for lrst_oce
	36	USE prtctl ! Print control
	37	USE lbclnk ! ocean lateral boundary condition (or mpp link)
	38	USE lib_mpp ! MPP library
	39	USE timing ! Timing
[9023]	40	USE wet_dry ! Wetting/Drying flux limiting
[592]	41
[3]	42	IMPLICIT NONE
	43	PRIVATE
	44
[1438]	45	PUBLIC ssh_nxt ! called by step.F90
[4292]	46	PUBLIC wzv ! called by step.F90
[10364]	47	PUBLIC wAimp ! called by step.F90
[12377]	48	PUBLIC ssh_atf ! called by step.F90
[3]	49
	50	!! * Substitutions
[12377]	51	# include "do_loop_substitute.h90"
[3]	52	!!----------------------------------------------------------------------
[9598]	53	!! NEMO/OCE 4.0 , NEMO Consortium (2018)
[888]	54	!! $Id$
[10068]	55	!! Software governed by the CeCILL license (see ./LICENSE)
[592]	56	!!----------------------------------------------------------------------
[3]	57	CONTAINS
	58
[12377]	59	SUBROUTINE ssh_nxt( kt, Kbb, Kmm, pssh, Kaa )
[3]	60	!!----------------------------------------------------------------------
[4292]	61	!! * ROUTINE ssh_nxt *
[1438]	62	!!
[12377]	63	!! ** Purpose : compute the after ssh (ssh(Kaa))
[3]	64	!!
[4292]	65	!! ** Method : - Using the incompressibility hypothesis, the ssh increment
	66	!! is computed by integrating the horizontal divergence and multiply by
	67	!! by the time step.
[3]	68	!!
[12377]	69	!! ** action : ssh(:,:,Kaa), after sea surface height
[2528]	70	!!
	71	!! Reference : Leclair, M., and G. Madec, 2009, Ocean Modelling.
[3]	72	!!----------------------------------------------------------------------
[12377]	73	INTEGER , INTENT(in ) :: kt ! time step
	74	INTEGER , INTENT(in ) :: Kbb, Kmm, Kaa ! time level index
	75	REAL(wp), DIMENSION(jpi,jpj,jpt), INTENT(inout) :: pssh ! sea-surface height
[4292]	76	!
[12489]	77	INTEGER :: jk ! dummy loop index
	78	REAL(wp) :: zcoef ! local scalar
[9019]	79	REAL(wp), DIMENSION(jpi,jpj) :: zhdiv ! 2D workspace
[3]	80	!!----------------------------------------------------------------------
[3294]	81	!
[9019]	82	IF( ln_timing ) CALL timing_start('ssh_nxt')
[3294]	83	!
[3]	84	IF( kt == nit000 ) THEN
	85	IF(lwp) WRITE(numout,*)
[4292]	86	IF(lwp) WRITE(numout,*) 'ssh_nxt : after sea surface height'
[1438]	87	IF(lwp) WRITE(numout,*) '~~~~~~~ '
[3]	88	ENDIF
[2528]	89	!
[12489]	90	zcoef = 0.5_wp * r1_rho0
[3]	91
[1438]	92	! !------------------------------!
	93	! ! After Sea Surface Height !
	94	! !------------------------------!
[9023]	95	IF(ln_wd_il) THEN
[12489]	96	CALL wad_lmt(pssh(:,:,Kbb), zcoef * (emp_b(:,:) + emp(:,:)), rDt, Kmm, uu, vv )
[7753]	97	ENDIF
[7646]	98
[12377]	99	CALL div_hor( kt, Kbb, Kmm ) ! Horizontal divergence
[7646]	100	!
[7753]	101	zhdiv(:,:) = 0._wp
[1438]	102	DO jk = 1, jpkm1 ! Horizontal divergence of barotropic transports
[12377]	103	zhdiv(:,:) = zhdiv(:,:) + e3t(:,:,jk,Kmm) * hdiv(:,:,jk)
[1438]	104	END DO
	105	! ! Sea surface elevation time stepping
[4338]	106	! In time-split case we need a first guess of the ssh after (using the baroclinic timestep) in order to
	107	! compute the vertical velocity which can be used to compute the non-linear terms of the momentum equations.
	108	!
[12489]	109	pssh(:,:,Kaa) = ( pssh(:,:,Kbb) - rDt * ( zcoef * ( emp_b(:,:) + emp(:,:) ) + zhdiv(:,:) ) ) * ssmask(:,:)
[9023]	110	!
	111	#if defined key_agrif
[12377]	112	Kbb_a = Kbb; Kmm_a = Kmm; Krhs_a = Kaa; CALL agrif_ssh( kt )
[9023]	113	#endif
	114	!
[5930]	115	IF ( .NOT.ln_dynspg_ts ) THEN
[7646]	116	IF( ln_bdy ) THEN
[12377]	117	CALL lbc_lnk( 'sshwzv', pssh(:,:,Kaa), 'T', 1. ) ! Not sure that's necessary
	118	CALL bdy_ssh( pssh(:,:,Kaa) ) ! Duplicate sea level across open boundaries
[5930]	119	ENDIF
[4292]	120	ENDIF
	121	! !------------------------------!
	122	! ! outputs !
	123	! !------------------------------!
	124	!
[12377]	125	IF(sn_cfctl%l_prtctl) CALL prt_ctl( tab2d_1=pssh(:,:,Kaa), clinfo1=' pssh(:,:,Kaa) - : ', mask1=tmask )
[4292]	126	!
[9019]	127	IF( ln_timing ) CALL timing_stop('ssh_nxt')
[4292]	128	!
	129	END SUBROUTINE ssh_nxt
	130
	131
[12377]	132	SUBROUTINE wzv( kt, Kbb, Kmm, pww, Kaa )
[4292]	133	!!----------------------------------------------------------------------
	134	!! * ROUTINE wzv *
	135	!!
	136	!! ** Purpose : compute the now vertical velocity
	137	!!
	138	!! ** Method : - Using the incompressibility hypothesis, the vertical
	139	!! velocity is computed by integrating the horizontal divergence
	140	!! from the bottom to the surface minus the scale factor evolution.
	141	!! The boundary conditions are w=0 at the bottom (no flux) and.
	142	!!
[12377]	143	!! ** action : pww : now vertical velocity
[4292]	144	!!
	145	!! Reference : Leclair, M., and G. Madec, 2009, Ocean Modelling.
	146	!!----------------------------------------------------------------------
[12377]	147	INTEGER , INTENT(in) :: kt ! time step
	148	INTEGER , INTENT(in) :: Kbb, Kmm, Kaa ! time level indices
	149	REAL(wp), DIMENSION(jpi,jpj,jpk), INTENT(inout) :: pww ! now vertical velocity
[4292]	150	!
[5836]	151	INTEGER :: ji, jj, jk ! dummy loop indices
[9019]	152	REAL(wp), ALLOCATABLE, DIMENSION(:,:,:) :: zhdiv
[4292]	153	!!----------------------------------------------------------------------
	154	!
[9019]	155	IF( ln_timing ) CALL timing_start('wzv')
[5836]	156	!
[4292]	157	IF( kt == nit000 ) THEN
	158	IF(lwp) WRITE(numout,*)
	159	IF(lwp) WRITE(numout,*) 'wzv : now vertical velocity '
	160	IF(lwp) WRITE(numout,*) '~~~~~ '
	161	!
[12377]	162	pww(:,:,jpk) = 0._wp ! bottom boundary condition: w=0 (set once for all)
[4292]	163	ENDIF
	164	! !------------------------------!
	165	! ! Now Vertical Velocity !
	166	! !------------------------------!
	167	!
	168	IF( ln_vvl_ztilde .OR. ln_vvl_layer ) THEN ! z_tilde and layer cases
[9019]	169	ALLOCATE( zhdiv(jpi,jpj,jpk) )
[4292]	170	!
	171	DO jk = 1, jpkm1
	172	! horizontal divergence of thickness diffusion transport ( velocity multiplied by e3t)
[4338]	173	! - ML - note: computation already done in dom_vvl_sf_nxt. Could be optimized (not critical and clearer this way)
[12377]	174	DO_2D_00_00
	175	zhdiv(ji,jj,jk) = r1_e1e2t(ji,jj) * ( un_td(ji,jj,jk) - un_td(ji-1,jj,jk) + vn_td(ji,jj,jk) - vn_td(ji,jj-1,jk) )
	176	END_2D
[592]	177	END DO
[10425]	178	CALL lbc_lnk('sshwzv', zhdiv, 'T', 1.) ! - ML - Perhaps not necessary: not used for horizontal "connexions"
[4292]	179	! ! Is it problematic to have a wrong vertical velocity in boundary cells?
[12377]	180	! ! Same question holds for hdiv. Perhaps just for security
[4292]	181	DO jk = jpkm1, 1, -1 ! integrate from the bottom the hor. divergence
	182	! computation of w
[12377]	183	pww(:,:,jk) = pww(:,:,jk+1) - ( e3t(:,:,jk,Kmm) * hdiv(:,:,jk) + zhdiv(:,:,jk) &
[12489]	184	& + r1_Dt * ( e3t(:,:,jk,Kaa) - e3t(:,:,jk,Kbb) ) ) * tmask(:,:,jk)
[4292]	185	END DO
[12377]	186	! IF( ln_vvl_layer ) pww(:,:,:) = 0.e0
[9019]	187	DEALLOCATE( zhdiv )
[4292]	188	ELSE ! z_star and linear free surface cases
	189	DO jk = jpkm1, 1, -1 ! integrate from the bottom the hor. divergence
[7753]	190	! computation of w
[12377]	191	pww(:,:,jk) = pww(:,:,jk+1) - ( e3t(:,:,jk,Kmm) * hdiv(:,:,jk) &
[12489]	192	& + r1_Dt * ( e3t(:,:,jk,Kaa) - e3t(:,:,jk,Kbb) ) ) * tmask(:,:,jk)
[4292]	193	END DO
[1438]	194	ENDIF
[592]	195
[7646]	196	IF( ln_bdy ) THEN
[4327]	197	DO jk = 1, jpkm1
[12377]	198	pww(:,:,jk) = pww(:,:,jk) * bdytmask(:,:)
[4327]	199	END DO
	200	ENDIF
[4292]	201	!
[9023]	202	#if defined key_agrif
	203	IF( .NOT. AGRIF_Root() ) THEN
[12377]	204	IF ((nbondi == 1).OR.(nbondi == 2)) pww(nlci-1 , : ,:) = 0.e0 ! east
	205	IF ((nbondi == -1).OR.(nbondi == 2)) pww(2 , : ,:) = 0.e0 ! west
	206	IF ((nbondj == 1).OR.(nbondj == 2)) pww(: ,nlcj-1 ,:) = 0.e0 ! north
	207	IF ((nbondj == -1).OR.(nbondj == 2)) pww(: ,2 ,:) = 0.e0 ! south
[9023]	208	ENDIF
	209	#endif
[5836]	210	!
[9124]	211	IF( ln_timing ) CALL timing_stop('wzv')
[9023]	212	!
[5836]	213	END SUBROUTINE wzv
[592]	214
	215
[12377]	216	SUBROUTINE ssh_atf( kt, Kbb, Kmm, Kaa, pssh )
[1438]	217	!!----------------------------------------------------------------------
[12377]	218	!! * ROUTINE ssh_atf *
[1438]	219	!!
[12377]	220	!! ** Purpose : Apply Asselin time filter to now SSH.
[1438]	221	!!
[2528]	222	!! ** Method : - apply Asselin time fiter to now ssh (excluding the forcing
	223	!! from the filter, see Leclair and Madec 2010) and swap :
[12489]	224	!! pssh(:,:,Kmm) = pssh(:,:,Kaa) + rn_atfp * ( pssh(:,:,Kbb) -2 pssh(:,:,Kmm) + pssh(:,:,Kaa) )
	225	!! - rn_atfp * rn_Dt * ( emp_b - emp ) / rho0
[1438]	226	!!
[12377]	227	!! ** action : - pssh(:,:,Kmm) time filtered
[2528]	228	!!
	229	!! Reference : Leclair, M., and G. Madec, 2009, Ocean Modelling.
[1438]	230	!!----------------------------------------------------------------------
[12377]	231	INTEGER , INTENT(in ) :: kt ! ocean time-step index
	232	INTEGER , INTENT(in ) :: Kbb, Kmm, Kaa ! ocean time level indices
	233	REAL(wp), DIMENSION(jpi,jpj,jpt), INTENT(inout) :: pssh ! SSH field
[6140]	234	!
	235	REAL(wp) :: zcoef ! local scalar
[1438]	236	!!----------------------------------------------------------------------
[3294]	237	!
[12377]	238	IF( ln_timing ) CALL timing_start('ssh_atf')
[3294]	239	!
[1438]	240	IF( kt == nit000 ) THEN
	241	IF(lwp) WRITE(numout,*)
[12377]	242	IF(lwp) WRITE(numout,*) 'ssh_atf : Asselin time filter of sea surface height'
[1438]	243	IF(lwp) WRITE(numout,*) '~~~~~~~ '
	244	ENDIF
[6140]	245	! !== Euler time-stepping: no filter, just swap ==!
[12489]	246	IF ( .NOT.( l_1st_euler ) ) THEN ! Only do time filtering for leapfrog timesteps
[12377]	247	! ! filtered "now" field
[12489]	248	pssh(:,:,Kmm) = pssh(:,:,Kmm) + rn_atfp * ( pssh(:,:,Kbb) - 2 * pssh(:,:,Kmm) + pssh(:,:,Kaa) )
[12377]	249	IF( .NOT.ln_linssh ) THEN ! "now" <-- with forcing removed
[12489]	250	zcoef = rn_atfp * rn_Dt * r1_rho0
[12377]	251	pssh(:,:,Kmm) = pssh(:,:,Kmm) - zcoef * ( emp_b(:,:) - emp (:,:) &
	252	& - rnf_b(:,:) + rnf (:,:) &
	253	& + fwfisf_cav_b(:,:) - fwfisf_cav(:,:) &
	254	& + fwfisf_par_b(:,:) - fwfisf_par(:,:) ) * ssmask(:,:)
	255
	256	! ice sheet coupling
[12489]	257	IF ( ln_isf .AND. ln_isfcpl .AND. kt == nit000+1) pssh(:,:,Kbb) = pssh(:,:,Kbb) - rn_atfp * rn_Dt * ( risfcpl_ssh(:,:) - 0.0 ) * ssmask(:,:)
[12377]	258
[6140]	259	ENDIF
[1438]	260	ENDIF
	261	!
[12377]	262	IF(sn_cfctl%l_prtctl) CALL prt_ctl( tab2d_1=pssh(:,:,Kmm), clinfo1=' pssh(:,:,Kmm) - : ', mask1=tmask )
[2528]	263	!
[12377]	264	IF( ln_timing ) CALL timing_stop('ssh_atf')
[3294]	265	!
[12377]	266	END SUBROUTINE ssh_atf
[3]	267
[12377]	268	SUBROUTINE wAimp( kt, Kmm )
[10364]	269	!!----------------------------------------------------------------------
	270	!! * ROUTINE wAimp *
	271	!!
	272	!! ** Purpose : compute the Courant number and partition vertical velocity
	273	!! if a proportion needs to be treated implicitly
	274	!!
	275	!! ** Method : -
	276	!!
[12377]	277	!! ** action : ww : now vertical velocity (to be handled explicitly)
[10364]	278	!! : wi : now vertical velocity (for implicit treatment)
	279	!!
[11414]	280	!! Reference : Shchepetkin, A. F. (2015): An adaptive, Courant-number-dependent
	281	!! implicit scheme for vertical advection in oceanic modeling.
	282	!! Ocean Modelling, 91, 38-69.
[10364]	283	!!----------------------------------------------------------------------
	284	INTEGER, INTENT(in) :: kt ! time step
[12377]	285	INTEGER, INTENT(in) :: Kmm ! time level index
[10364]	286	!
	287	INTEGER :: ji, jj, jk ! dummy loop indices
[11293]	288	REAL(wp) :: zCu, zcff, z1_e3t ! local scalars
[10364]	289	REAL(wp) , PARAMETER :: Cu_min = 0.15_wp ! local parameters
[11407]	290	REAL(wp) , PARAMETER :: Cu_max = 0.30_wp ! local parameters
[10364]	291	REAL(wp) , PARAMETER :: Cu_cut = 2._wp*Cu_max - Cu_min ! local parameters
	292	REAL(wp) , PARAMETER :: Fcu = 4._wpCu_max(Cu_max-Cu_min) ! local parameters
	293	!!----------------------------------------------------------------------
	294	!
	295	IF( ln_timing ) CALL timing_start('wAimp')
	296	!
	297	IF( kt == nit000 ) THEN
	298	IF(lwp) WRITE(numout,*)
	299	IF(lwp) WRITE(numout,*) 'wAimp : Courant number-based partitioning of now vertical velocity '
	300	IF(lwp) WRITE(numout,*) '~~~~~ '
[11293]	301	wi(:,:,:) = 0._wp
[10364]	302	ENDIF
	303	!
[11414]	304	! Calculate Courant numbers
	305	IF( ln_vvl_ztilde .OR. ln_vvl_layer ) THEN
[12377]	306	DO_3D_00_00( 1, jpkm1 )
	307	z1_e3t = 1._wp / e3t(ji,jj,jk,Kmm)
[12489]	308	! 2*rn_Dt and not rDt (for restartability)
	309	Cu_adv(ji,jj,jk) = 2._wp * rn_Dt * ( ( MAX( ww(ji,jj,jk) , 0._wp ) - MIN( ww(ji,jj,jk+1) , 0._wp ) ) &
[12377]	310	& + ( MAX( e2u(ji ,jj)e3u(ji ,jj,jk,Kmm)uu(ji ,jj,jk,Kmm) + un_td(ji ,jj,jk), 0._wp ) - &
	311	& MIN( e2u(ji-1,jj)e3u(ji-1,jj,jk,Kmm)uu(ji-1,jj,jk,Kmm) + un_td(ji-1,jj,jk), 0._wp ) ) &
	312	& * r1_e1e2t(ji,jj) &
	313	& + ( MAX( e1v(ji,jj )e3v(ji,jj ,jk,Kmm)vv(ji,jj ,jk,Kmm) + vn_td(ji,jj ,jk), 0._wp ) - &
	314	& MIN( e1v(ji,jj-1)e3v(ji,jj-1,jk,Kmm)vv(ji,jj-1,jk,Kmm) + vn_td(ji,jj-1,jk), 0._wp ) ) &
	315	& * r1_e1e2t(ji,jj) &
	316	& ) * z1_e3t
	317	END_3D
[11414]	318	ELSE
[12377]	319	DO_3D_00_00( 1, jpkm1 )
	320	z1_e3t = 1._wp / e3t(ji,jj,jk,Kmm)
[12489]	321	! 2*rn_Dt and not rDt (for restartability)
	322	Cu_adv(ji,jj,jk) = 2._wp * rn_Dt * ( ( MAX( ww(ji,jj,jk) , 0._wp ) - MIN( ww(ji,jj,jk+1) , 0._wp ) ) &
[12377]	323	& + ( MAX( e2u(ji ,jj)e3u(ji ,jj,jk,Kmm)uu(ji ,jj,jk,Kmm), 0._wp ) - &
	324	& MIN( e2u(ji-1,jj)e3u(ji-1,jj,jk,Kmm)uu(ji-1,jj,jk,Kmm), 0._wp ) ) &
	325	& * r1_e1e2t(ji,jj) &
	326	& + ( MAX( e1v(ji,jj )e3v(ji,jj ,jk,Kmm)vv(ji,jj ,jk,Kmm), 0._wp ) - &
	327	& MIN( e1v(ji,jj-1)e3v(ji,jj-1,jk,Kmm)vv(ji,jj-1,jk,Kmm), 0._wp ) ) &
	328	& * r1_e1e2t(ji,jj) &
	329	& ) * z1_e3t
	330	END_3D
[11414]	331	ENDIF
[10907]	332	CALL lbc_lnk( 'sshwzv', Cu_adv, 'T', 1. )
[10364]	333	!
	334	CALL iom_put("Courant",Cu_adv)
	335	!
	336	IF( MAXVAL( Cu_adv(:,:,:) ) > Cu_min ) THEN ! Quick check if any breaches anywhere
[12377]	337	DO_3DS_11_11( jpkm1, 2, -1 )
	338	!
	339	zCu = MAX( Cu_adv(ji,jj,jk) , Cu_adv(ji,jj,jk-1) )
[11293]	340	! alt:
[12377]	341	! IF ( ww(ji,jj,jk) > 0._wp ) THEN
[11293]	342	! zCu = Cu_adv(ji,jj,jk)
	343	! ELSE
	344	! zCu = Cu_adv(ji,jj,jk-1)
	345	! ENDIF
[12377]	346	!
	347	IF( zCu <= Cu_min ) THEN !<-- Fully explicit
	348	zcff = 0._wp
	349	ELSEIF( zCu < Cu_cut ) THEN !<-- Mixed explicit
	350	zcff = ( zCu - Cu_min )**2
	351	zcff = zcff / ( Fcu + zcff )
	352	ELSE !<-- Mostly implicit
	353	zcff = ( zCu - Cu_max )/ zCu
	354	ENDIF
	355	zcff = MIN(1._wp, zcff)
	356	!
	357	wi(ji,jj,jk) = zcff * ww(ji,jj,jk)
	358	ww(ji,jj,jk) = ( 1._wp - zcff ) * ww(ji,jj,jk)
	359	!
	360	Cu_adv(ji,jj,jk) = zcff ! Reuse array to output coefficient below and in stp_ctl
	361	END_3D
[11293]	362	Cu_adv(:,:,1) = 0._wp
[10364]	363	ELSE
	364	! Fully explicit everywhere
[11407]	365	Cu_adv(:,:,:) = 0._wp ! Reuse array to output coefficient below and in stp_ctl
[10907]	366	wi (:,:,:) = 0._wp
[10364]	367	ENDIF
	368	CALL iom_put("wimp",wi)
	369	CALL iom_put("wi_cff",Cu_adv)
[12377]	370	CALL iom_put("wexp",ww)
[10364]	371	!
	372	IF( ln_timing ) CALL timing_stop('wAimp')
	373	!
	374	END SUBROUTINE wAimp
[3]	375	!!======================================================================
[1565]	376	END MODULE sshwzv

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

Download in other formats: