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.

tranxt.F90 in branches/2017/dev_merge_2017/NEMOGCM/NEMO/OPA_SRC/TRA – NEMO

Context Navigation

source: branches/2017/dev_merge_2017/NEMOGCM/NEMO/OPA_SRC/TRA/tranxt.F90 @ 9093

Last change on this file since 9093 was 9023, checked in by timgraham, 7 years ago
Merged METO_MERCATOR branch and resolved all conflicts in OPA_SRC
Property svn:keywords set to `Id`
File size: 19.2 KB

Rev	Line
[3]	1	MODULE tranxt
	2	!!======================================================================
	3	!! * MODULE tranxt *
	4	!! Ocean active tracers: time stepping on temperature and salinity
	5	!!======================================================================
[1110]	6	!! History : OPA ! 1991-11 (G. Madec) Original code
	7	!! 7.0 ! 1993-03 (M. Guyon) symetrical conditions
	8	!! 8.0 ! 1996-02 (G. Madec & M. Imbard) opa release 8.0
	9	!! - ! 1996-04 (A. Weaver) Euler forward step
	10	!! 8.2 ! 1999-02 (G. Madec, N. Grima) semi-implicit pressure grad.
	11	!! NEMO 1.0 ! 2002-08 (G. Madec) F90: Free form and module
	12	!! - ! 2002-11 (C. Talandier, A-M Treguier) Open boundaries
	13	!! - ! 2005-04 (C. Deltel) Add Asselin trend in the ML budget
	14	!! 2.0 ! 2006-02 (L. Debreu, C. Mazauric) Agrif implementation
	15	!! 3.0 ! 2008-06 (G. Madec) time stepping always done in trazdf
[1438]	16	!! 3.1 ! 2009-02 (G. Madec, R. Benshila) re-introduce the vvl option
[2528]	17	!! 3.3 ! 2010-04 (M. Leclair, G. Madec) semi-implicit hpg with asselin filter + modified LF-RA
	18	!! - ! 2010-05 (C. Ethe, G. Madec) merge TRC-TRA
[3]	19	!!----------------------------------------------------------------------
[503]	20
	21	!!----------------------------------------------------------------------
[2528]	22	!! tra_nxt : time stepping on tracers
	23	!! tra_nxt_fix : time stepping on tracers : fixed volume case
	24	!! tra_nxt_vvl : time stepping on tracers : variable volume case
[3]	25	!!----------------------------------------------------------------------
	26	USE oce ! ocean dynamics and tracers variables
	27	USE dom_oce ! ocean space and time domain variables
[2528]	28	USE sbc_oce ! surface boundary condition: ocean
[5467]	29	USE sbcrnf ! river runoffs
[6140]	30	USE sbcisf ! ice shelf melting
[4990]	31	USE zdf_oce ! ocean vertical mixing
[1438]	32	USE domvvl ! variable volume
[4990]	33	USE trd_oce ! trends: ocean variables
	34	USE trdtra ! trends manager: tracers
	35	USE traqsr ! penetrative solar radiation (needed for nksr)
	36	USE phycst ! physical constant
[9019]	37	USE ldftra ! lateral physics : tracers
	38	USE ldfslp ! lateral physics : slopes
	39	USE bdy_oce , ONLY : ln_bdy
[3294]	40	USE bdytra ! open boundary condition (bdy_tra routine)
[4990]	41	!
[3]	42	USE in_out_manager ! I/O manager
	43	USE lbclnk ! ocean lateral boundary conditions (or mpp link)
[258]	44	USE prtctl ! Print control
[4990]	45	USE timing ! Timing
[2528]	46	#if defined key_agrif
[389]	47	USE agrif_opa_interp
[2528]	48	#endif
[3]	49
	50	IMPLICIT NONE
	51	PRIVATE
	52
[2528]	53	PUBLIC tra_nxt ! routine called by step.F90
	54	PUBLIC tra_nxt_fix ! to be used in trcnxt
	55	PUBLIC tra_nxt_vvl ! to be used in trcnxt
[592]	56
	57	!! * Substitutions
[6140]	58	# include "vectopt_loop_substitute.h90"
[3]	59	!!----------------------------------------------------------------------
[2528]	60	!! NEMO/OPA 3.3 , NEMO-Consortium (2010)
[2715]	61	!! $Id$
	62	!! Software governed by the CeCILL licence (NEMOGCM/NEMO_CeCILL.txt)
[3]	63	!!----------------------------------------------------------------------
	64	CONTAINS
	65
	66	SUBROUTINE tra_nxt( kt )
	67	!!----------------------------------------------------------------------
	68	!! * ROUTINE tranxt *
	69	!!
[1110]	70	!! ** Purpose : Apply the boundary condition on the after temperature
	71	!! and salinity fields, achieved the time stepping by adding
	72	!! the Asselin filter on now fields and swapping the fields.
[3]	73	!!
[1110]	74	!! ** Method : At this stage of the computation, ta and sa are the
	75	!! after temperature and salinity as the time stepping has
	76	!! been performed in trazdf_imp or trazdf_exp module.
[3]	77	!!
[1110]	78	!! - Apply lateral boundary conditions on (ta,sa)
	79	!! at the local domain boundaries through lbc_lnk call,
[7646]	80	!! at the one-way open boundaries (ln_bdy=T),
[4990]	81	!! at the AGRIF zoom boundaries (lk_agrif=T)
[1110]	82	!!
[1438]	83	!! - Update lateral boundary conditions on AGRIF children
	84	!! domains (lk_agrif=T)
[1110]	85	!!
[6140]	86	!! ** Action : - tsb & tsn ready for the next time step
[503]	87	!!----------------------------------------------------------------------
	88	INTEGER, INTENT(in) :: kt ! ocean time-step index
	89	!!
[6140]	90	INTEGER :: ji, jj, jk, jn ! dummy loop indices
	91	REAL(wp) :: zfact ! local scalars
[9019]	92	REAL(wp), ALLOCATABLE, DIMENSION(:,:,:) :: ztrdt, ztrds
[3]	93	!!----------------------------------------------------------------------
[3294]	94	!
[9019]	95	IF( ln_timing ) CALL timing_start( 'tra_nxt')
[3294]	96	!
[1110]	97	IF( kt == nit000 ) THEN
	98	IF(lwp) WRITE(numout,*)
	99	IF(lwp) WRITE(numout,*) 'tra_nxt : achieve the time stepping by Asselin filter and array swap'
	100	IF(lwp) WRITE(numout,*) '~~~~~~~'
[592]	101	ENDIF
	102
[1110]	103	! Update after tracer on domain lateral boundaries
	104	!
[5656]	105	#if defined key_agrif
	106	CALL Agrif_tra ! AGRIF zoom boundaries
	107	#endif
	108	!
[4230]	109	CALL lbc_lnk( tsa(:,:,:,jp_tem), 'T', 1._wp ) ! local domain boundaries (T-point, unchanged sign)
	110	CALL lbc_lnk( tsa(:,:,:,jp_sal), 'T', 1._wp )
[1110]	111	!
[7646]	112	IF( ln_bdy ) CALL bdy_tra( kt ) ! BDY open boundaries
[1438]	113
	114	! set time step size (Euler/Leapfrog)
[9019]	115	IF( neuler == 0 .AND. kt == nit000 ) THEN ; r2dt = rdt ! at nit000 (Euler)
[6140]	116	ELSEIF( kt <= nit000 + 1 ) THEN ; r2dt = 2._wp* rdt ! at nit000 or nit000+1 (Leapfrog)
[1438]	117	ENDIF
[3]	118
[1110]	119	! trends computation initialisation
[7646]	120	IF( l_trdtra ) THEN
[9019]	121	ALLOCATE( ztrdt(jpi,jpj,jpk) , ztrds(jpi,jpj,jpk) )
[8698]	122	ztrdt(:,:,jpk) = 0._wp
	123	ztrds(:,:,jpk) = 0._wp
[4990]	124	IF( ln_traldf_iso ) THEN ! diagnose the "pure" Kz diffusive trend
	125	CALL trd_tra( kt, 'TRA', jp_tem, jptra_zdfp, ztrdt )
	126	CALL trd_tra( kt, 'TRA', jp_sal, jptra_zdfp, ztrds )
	127	ENDIF
[7646]	128	! total trend for the non-time-filtered variables.
[8698]	129	zfact = 1.0 / rdt
	130	! G Nurser 23 Mar 2017. Recalculate trend as Delta(e3t*T)/e3tn; e3tn cancel from tsn terms
[7646]	131	DO jk = 1, jpkm1
[8698]	132	ztrdt(:,:,jk) = ( tsa(:,:,jk,jp_tem)e3t_a(:,:,jk) / e3t_n(:,:,jk) - tsn(:,:,jk,jp_tem)) zfact
	133	ztrds(:,:,jk) = ( tsa(:,:,jk,jp_sal)e3t_a(:,:,jk) / e3t_n(:,:,jk) - tsn(:,:,jk,jp_sal)) zfact
[7646]	134	END DO
	135	CALL trd_tra( kt, 'TRA', jp_tem, jptra_tot, ztrdt )
	136	CALL trd_tra( kt, 'TRA', jp_sal, jptra_tot, ztrds )
[9019]	137	IF( ln_linssh ) THEN ! linear sea surface height only
[8698]	138	! Store now fields before applying the Asselin filter
	139	! in order to calculate Asselin filter trend later.
	140	ztrdt(:,:,:) = tsn(:,:,:,jp_tem)
	141	ztrds(:,:,:) = tsn(:,:,:,jp_sal)
	142	ENDIF
[1110]	143	ENDIF
	144
[2528]	145	IF( neuler == 0 .AND. kt == nit000 ) THEN ! Euler time-stepping at first time-step (only swap)
	146	DO jn = 1, jpts
	147	DO jk = 1, jpkm1
[7753]	148	tsn(:,:,jk,jn) = tsa(:,:,jk,jn)
[2528]	149	END DO
	150	END DO
[9019]	151	IF (l_trdtra .AND. .NOT. ln_linssh ) THEN ! Zero Asselin filter contribution must be explicitly written out since for vvl
	152	! ! Asselin filter is output by tra_nxt_vvl that is not called on this time step
[8698]	153	ztrdt(:,:,:) = 0._wp
	154	ztrds(:,:,:) = 0._wp
	155	CALL trd_tra( kt, 'TRA', jp_tem, jptra_atf, ztrdt )
	156	CALL trd_tra( kt, 'TRA', jp_sal, jptra_atf, ztrds )
	157	END IF
[6140]	158	!
[2528]	159	ELSE ! Leap-Frog + Asselin filter time stepping
	160	!
[6140]	161	IF( ln_linssh ) THEN ; CALL tra_nxt_fix( kt, nit000, 'TRA', tsb, tsn, tsa, jpts ) ! linear free surface
	162	ELSE ; CALL tra_nxt_vvl( kt, nit000, rdt, 'TRA', tsb, tsn, tsa, &
	163	& sbc_tsc, sbc_tsc_b, jpts ) ! non-linear free surface
[2528]	164	ENDIF
[6140]	165	!
	166	DO jn = 1, jpts
	167	CALL lbc_lnk( tsb(:,:,:,jn), 'T', 1._wp )
	168	CALL lbc_lnk( tsn(:,:,:,jn), 'T', 1._wp )
	169	CALL lbc_lnk( tsa(:,:,:,jn), 'T', 1._wp )
	170	END DO
[5656]	171	ENDIF
[2715]	172	!
[8698]	173	IF( l_trdtra .AND. ln_linssh ) THEN ! trend of the Asselin filter (tb filtered - tb)/dt
	174	zfact = 1._wp / r2dt
[1110]	175	DO jk = 1, jpkm1
[7753]	176	ztrdt(:,:,jk) = ( tsb(:,:,jk,jp_tem) - ztrdt(:,:,jk) ) * zfact
	177	ztrds(:,:,jk) = ( tsb(:,:,jk,jp_sal) - ztrds(:,:,jk) ) * zfact
[1110]	178	END DO
[4990]	179	CALL trd_tra( kt, 'TRA', jp_tem, jptra_atf, ztrdt )
	180	CALL trd_tra( kt, 'TRA', jp_sal, jptra_atf, ztrds )
[1438]	181	END IF
[9019]	182	IF( l_trdtra ) DEALLOCATE( ztrdt , ztrds )
[2715]	183	!
[1438]	184	! ! control print
[2528]	185	IF(ln_ctl) CALL prt_ctl( tab3d_1=tsn(:,:,:,jp_tem), clinfo1=' nxt - Tn: ', mask1=tmask, &
	186	& tab3d_2=tsn(:,:,:,jp_sal), clinfo2= ' Sn: ', mask2=tmask )
[1438]	187	!
[9019]	188	IF( ln_timing ) CALL timing_stop('tra_nxt')
[3294]	189	!
[1438]	190	END SUBROUTINE tra_nxt
	191
	192
[3294]	193	SUBROUTINE tra_nxt_fix( kt, kit000, cdtype, ptb, ptn, pta, kjpt )
[1438]	194	!!----------------------------------------------------------------------
	195	!! * ROUTINE tra_nxt_fix *
	196	!!
	197	!! ** Purpose : fixed volume: apply the Asselin time filter and
	198	!! swap the tracer fields.
	199	!!
	200	!! ** Method : - Apply a Asselin time filter on now fields.
	201	!! - swap tracer fields to prepare the next time_step.
	202	!!
[6140]	203	!! ** Action : - tsb & tsn ready for the next time step
[1438]	204	!!----------------------------------------------------------------------
[6140]	205	INTEGER , INTENT(in ) :: kt ! ocean time-step index
	206	INTEGER , INTENT(in ) :: kit000 ! first time step index
	207	CHARACTER(len=3) , INTENT(in ) :: cdtype ! =TRA or TRC (tracer indicator)
	208	INTEGER , INTENT(in ) :: kjpt ! number of tracers
	209	REAL(wp), DIMENSION(jpi,jpj,jpk,kjpt), INTENT(inout) :: ptb ! before tracer fields
	210	REAL(wp), DIMENSION(jpi,jpj,jpk,kjpt), INTENT(inout) :: ptn ! now tracer fields
	211	REAL(wp), DIMENSION(jpi,jpj,jpk,kjpt), INTENT(inout) :: pta ! tracer trend
[2715]	212	!
[2528]	213	INTEGER :: ji, jj, jk, jn ! dummy loop indices
	214	REAL(wp) :: ztn, ztd ! local scalars
[1438]	215	!!----------------------------------------------------------------------
[6140]	216	!
[3294]	217	IF( kt == kit000 ) THEN
[1438]	218	IF(lwp) WRITE(numout,*)
[3294]	219	IF(lwp) WRITE(numout,*) 'tra_nxt_fix : time stepping', cdtype
[1438]	220	IF(lwp) WRITE(numout,*) '~~~~~~~~~~~'
	221	ENDIF
	222	!
[2528]	223	DO jn = 1, kjpt
[1438]	224	!
[2528]	225	DO jk = 1, jpkm1
[6140]	226	DO jj = 2, jpjm1
	227	DO ji = fs_2, fs_jpim1
[2528]	228	ztn = ptn(ji,jj,jk,jn)
[6140]	229	ztd = pta(ji,jj,jk,jn) - 2._wp * ztn + ptb(ji,jj,jk,jn) ! time laplacian on tracers
[2528]	230	!
[6140]	231	ptb(ji,jj,jk,jn) = ztn + atfp * ztd ! ptb <-- filtered ptn
	232	ptn(ji,jj,jk,jn) = pta(ji,jj,jk,jn) ! ptn <-- pta
[3]	233	END DO
[2528]	234	END DO
	235	END DO
[1110]	236	!
[2528]	237	END DO
[1438]	238	!
	239	END SUBROUTINE tra_nxt_fix
[3]	240
[1110]	241
[5385]	242	SUBROUTINE tra_nxt_vvl( kt, kit000, p2dt, cdtype, ptb, ptn, pta, psbc_tc, psbc_tc_b, kjpt )
[1438]	243	!!----------------------------------------------------------------------
	244	!! * ROUTINE tra_nxt_vvl *
	245	!!
	246	!! ** Purpose : Time varying volume: apply the Asselin time filter
	247	!! and swap the tracer fields.
	248	!!
	249	!! ** Method : - Apply a thickness weighted Asselin time filter on now fields.
	250	!! - swap tracer fields to prepare the next time_step.
[6140]	251	!! tb = ( e3t_ntn + atfp[ e3t_btb - 2 e3t_ntn + e3t_a*ta ] )
	252	!! /( e3t_n + atfp*[ e3t_b - 2 e3t_n + e3t_a ] )
	253	!! tn = ta
[1438]	254	!!
[6140]	255	!! ** Action : - tsb & tsn ready for the next time step
[1438]	256	!!----------------------------------------------------------------------
[6140]	257	INTEGER , INTENT(in ) :: kt ! ocean time-step index
	258	INTEGER , INTENT(in ) :: kit000 ! first time step index
	259	REAL(wp) , INTENT(in ) :: p2dt ! time-step
	260	CHARACTER(len=3) , INTENT(in ) :: cdtype ! =TRA or TRC (tracer indicator)
	261	INTEGER , INTENT(in ) :: kjpt ! number of tracers
	262	REAL(wp), DIMENSION(jpi,jpj,jpk,kjpt), INTENT(inout) :: ptb ! before tracer fields
	263	REAL(wp), DIMENSION(jpi,jpj,jpk,kjpt), INTENT(inout) :: ptn ! now tracer fields
	264	REAL(wp), DIMENSION(jpi,jpj,jpk,kjpt), INTENT(inout) :: pta ! tracer trend
	265	REAL(wp), DIMENSION(jpi,jpj ,kjpt), INTENT(in ) :: psbc_tc ! surface tracer content
	266	REAL(wp), DIMENSION(jpi,jpj ,kjpt), INTENT(in ) :: psbc_tc_b ! before surface tracer content
	267	!
[5930]	268	LOGICAL :: ll_traqsr, ll_rnf, ll_isf ! local logical
[2528]	269	INTEGER :: ji, jj, jk, jn ! dummy loop indices
[8698]	270	REAL(wp) :: zfact, zfact1, ztc_a , ztc_n , ztc_b , ztc_f , ztc_d ! local scalar
[2715]	271	REAL(wp) :: zfact2, ze3t_b, ze3t_n, ze3t_a, ze3t_f, ze3t_d ! - -
[9019]	272	REAL(wp), ALLOCATABLE, DIMENSION(:,:,:,:) :: ztrd_atf
[1438]	273	!!----------------------------------------------------------------------
[3294]	274	!
	275	IF( kt == kit000 ) THEN
[1438]	276	IF(lwp) WRITE(numout,*)
[3294]	277	IF(lwp) WRITE(numout,*) 'tra_nxt_vvl : time stepping', cdtype
[1438]	278	IF(lwp) WRITE(numout,*) '~~~~~~~~~~~'
	279	ENDIF
[2528]	280	!
	281	IF( cdtype == 'TRA' ) THEN
	282	ll_traqsr = ln_traqsr ! active tracers case and solar penetration
[5467]	283	ll_rnf = ln_rnf ! active tracers case and river runoffs
[6140]	284	ll_isf = ln_isf ! active tracers case and ice shelf melting
	285	ELSE ! passive tracers case
	286	ll_traqsr = .FALSE. ! NO solar penetration
	287	ll_rnf = .FALSE. ! NO river runoffs ???? !!gm BUG ?
	288	ll_isf = .FALSE. ! NO ice shelf melting/freezing !!gm BUG ??
[2528]	289	ENDIF
	290	!
[9019]	291	IF( ( l_trdtra .AND. cdtype == 'TRA' ) .OR. ( l_trdtrc .AND. cdtype == 'TRC' ) ) THEN
	292	ALLOCATE( ztrd_atf(jpi,jpj,jpk,kjpt) )
[8698]	293	ztrd_atf(:,:,:,:) = 0.0_wp
	294	ENDIF
	295	zfact = 1._wp / r2dt
[9019]	296	zfact1 = atfp * p2dt
	297	zfact2 = zfact1 * r1_rau0
[2528]	298	DO jn = 1, kjpt
	299	DO jk = 1, jpkm1
[6140]	300	DO jj = 2, jpjm1
	301	DO ji = fs_2, fs_jpim1
	302	ze3t_b = e3t_b(ji,jj,jk)
	303	ze3t_n = e3t_n(ji,jj,jk)
	304	ze3t_a = e3t_a(ji,jj,jk)
[2528]	305	! ! tracer content at Before, now and after
	306	ztc_b = ptb(ji,jj,jk,jn) * ze3t_b
	307	ztc_n = ptn(ji,jj,jk,jn) * ze3t_n
	308	ztc_a = pta(ji,jj,jk,jn) * ze3t_a
	309	!
	310	ze3t_d = ze3t_a - 2. * ze3t_n + ze3t_b
	311	ztc_d = ztc_a - 2. * ztc_n + ztc_b
	312	!
	313	ze3t_f = ze3t_n + atfp * ze3t_d
	314	ztc_f = ztc_n + atfp * ztc_d
	315	!
[5643]	316	IF( jk == mikt(ji,jj) ) THEN ! first level
	317	ze3t_f = ze3t_f - zfact2 * ( (emp_b(ji,jj) - emp(ji,jj) ) &
	318	& + (fwfisf_b(ji,jj) - fwfisf(ji,jj)) )
[5385]	319	ztc_f = ztc_f - zfact1 * ( psbc_tc(ji,jj,jn) - psbc_tc_b(ji,jj,jn) )
[2528]	320	ENDIF
[9023]	321	IF( ln_rnf_depth ) THEN
	322	! Rivers are not just at the surface must go down to nk_rnf(ji,jj)
	323	IF( mikt(ji,jj) <=jk .and. jk <= nk_rnf(ji,jj) ) THEN
	324	ze3t_f = ze3t_f - zfact2 * ( - (rnf_b(ji,jj) - rnf(ji,jj) ) ) &
	325	& * ( e3t_n(ji,jj,jk) / h_rnf(ji,jj) )
	326	ENDIF
	327	ELSE
	328	IF( jk == mikt(ji,jj) ) THEN ! first level
	329	ze3t_f = ze3t_f - zfact2 * ( - (rnf_b(ji,jj) - rnf(ji,jj) ) )
	330	ENDIF
	331	ENDIF
	332
[6140]	333	!
[5643]	334	! solar penetration (temperature only)
	335	IF( ll_traqsr .AND. jn == jp_tem .AND. jk <= nksr ) &
[2528]	336	& ztc_f = ztc_f - zfact1 * ( qsr_hc(ji,jj,jk) - qsr_hc_b(ji,jj,jk) )
[6140]	337	!
[5643]	338	! river runoff
	339	IF( ll_rnf .AND. jk <= nk_rnf(ji,jj) ) &
[5467]	340	& ztc_f = ztc_f - zfact1 * ( rnf_tsc(ji,jj,jn) - rnf_tsc_b(ji,jj,jn) ) &
[6140]	341	& * e3t_n(ji,jj,jk) / h_rnf(ji,jj)
	342	!
[5643]	343	! ice shelf
	344	IF( ll_isf ) THEN
	345	! level fully include in the Losch_2008 ice shelf boundary layer
	346	IF ( jk >= misfkt(ji,jj) .AND. jk < misfkb(ji,jj) ) &
	347	ztc_f = ztc_f - zfact1 * ( risf_tsc(ji,jj,jn) - risf_tsc_b(ji,jj,jn) ) &
[6140]	348	& * e3t_n(ji,jj,jk) * r1_hisf_tbl (ji,jj)
[5643]	349	! level partially include in Losch_2008 ice shelf boundary layer
	350	IF ( jk == misfkb(ji,jj) ) &
	351	ztc_f = ztc_f - zfact1 * ( risf_tsc(ji,jj,jn) - risf_tsc_b(ji,jj,jn) ) &
[6140]	352	& * e3t_n(ji,jj,jk) * r1_hisf_tbl (ji,jj) * ralpha(ji,jj)
[5643]	353	END IF
[6140]	354	!
[5467]	355	ze3t_f = 1.e0 / ze3t_f
	356	ptb(ji,jj,jk,jn) = ztc_f * ze3t_f ! ptb <-- ptn filtered
	357	ptn(ji,jj,jk,jn) = pta(ji,jj,jk,jn) ! ptn <-- pta
	358	!
[8698]	359	IF( ( l_trdtra .and. cdtype == 'TRA' ) .OR. ( l_trdtrc .and. cdtype == 'TRC' ) ) THEN
	360	ztrd_atf(ji,jj,jk,jn) = (ztc_f - ztc_n) * zfact/ze3t_n
	361	ENDIF
	362	!
[1438]	363	END DO
	364	END DO
[2528]	365	END DO
	366	!
	367	END DO
[503]	368	!
[9019]	369	IF( ( l_trdtra .AND. cdtype == 'TRA' ) .OR. ( l_trdtrc .AND. cdtype == 'TRC' ) ) THEN
	370	IF( l_trdtra .AND. cdtype == 'TRA' ) THEN
	371	CALL trd_tra( kt, cdtype, jp_tem, jptra_atf, ztrd_atf(:,:,:,jp_tem) )
	372	CALL trd_tra( kt, cdtype, jp_sal, jptra_atf, ztrd_atf(:,:,:,jp_sal) )
	373	ENDIF
	374	IF( l_trdtrc .AND. cdtype == 'TRC' ) THEN
	375	DO jn = 1, kjpt
	376	CALL trd_tra( kt, cdtype, jn, jptra_atf, ztrd_atf(:,:,:,jn) )
	377	END DO
	378	ENDIF
	379	DEALLOCATE( ztrd_atf )
[8698]	380	ENDIF
	381	!
[1438]	382	END SUBROUTINE tra_nxt_vvl
[3]	383
	384	!!======================================================================
	385	END MODULE tranxt

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

Download in other formats: