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.

dtadyn.F90 in trunk/NEMOGCM/NEMO/OFF_SRC – NEMO

Context Navigation

source: trunk/NEMOGCM/NEMO/OFF_SRC/dtadyn.F90 @ 2762

Last change on this file since 2762 was 2762, checked in by cetlod, 13 years ago
modify dtadyn.F90 routine to be used without key_ldfslp, see ticket #822
Property svn:keywords set to `Id`
File size: 34.8 KB

Rev	Line
[325]	1	MODULE dtadyn
	2	!!======================================================================
	3	!! * MODULE dtadyn *
[2528]	4	!! Off-line : interpolation of the physical fields
	5	!!======================================================================
	6	!! History : OPA ! 1992-01 (M. Imbard) Original code
	7	!! 8.0 ! 1998-04 (L.Bopp MA Foujols) slopes for isopyc.
	8	!! - ! 1998-05 (L. Bopp) read output of coupled run
	9	!! 8.2 ! 2001-01 (M. Levy et M. Benjelloul) add netcdf FORMAT
	10	!! NEMO 1.0 ! 2005-03 (O. Aumont and A. El Moussaoui) F90
	11	!! - ! 2005-12 (C. Ethe) Adapted for DEGINT
	12	!! 3.0 ! 2007-06 (C. Ethe) use of iom module
	13	!! - ! 2007-09 (C. Ethe) add swap_dyn_data
	14	!! 3.3 ! 2010-11 (C. Ethe) Full reorganization of the off-line: phasing with the on-line
	15	!!----------------------------------------------------------------------
[325]	16
	17	!!----------------------------------------------------------------------
	18	!! dta_dyn_init : initialization, namelist read, and parameters control
	19	!! dta_dyn : Interpolation of the fields
	20	!!----------------------------------------------------------------------
	21	USE oce ! ocean dynamics and tracers variables
[2528]	22	USE c1d ! 1D configuration: lk_c1d
	23	USE dom_oce ! ocean domain: variables
	24	USE zdf_oce ! ocean vertical physics: variables
	25	USE sbc_oce ! surface module: variables
[325]	26	USE phycst ! physical constants
[2528]	27	USE trabbl ! active tracer: bottom boundary layer
	28	USE ldfslp ! lateral diffusion: iso-neutral slopes
	29	USE ldfeiv ! eddy induced velocity coef.
[446]	30	USE ldftra_oce ! ocean tracer lateral physics
[2528]	31	USE zdfmxl ! vertical physics: mixed layer depth
	32	USE eosbn2 ! equation of state - Brunt Vaisala frequency
[325]	33	USE lbclnk ! ocean lateral boundary conditions (or mpp link)
[2528]	34	USE zpshde ! z-coord. with partial steps: horizontal derivatives
	35	USE in_out_manager ! I/O manager
	36	USE iom ! I/O library
[325]	37	USE lib_mpp ! distributed memory computing library
[2528]	38	USE prtctl ! print control
[325]	39
	40	IMPLICIT NONE
	41	PRIVATE
	42
[2528]	43	PUBLIC dta_dyn_init ! called by opa.F90
	44	PUBLIC dta_dyn ! called by step.F90
[325]	45
[2528]	46	LOGICAL, PUBLIC :: lperdyn = .TRUE. !: boolean for periodic fields or not
	47	LOGICAL, PUBLIC :: lfirdyn = .TRUE. !: boolean for the first call or not
[325]	48
[2528]	49	INTEGER, PUBLIC :: ndtadyn = 73 !: Number of dat in one year
	50	INTEGER, PUBLIC :: ndtatot = 73 !: Number of data in the input field
	51	INTEGER, PUBLIC :: nsptint = 1 !: type of spatial interpolation
[325]	52
[2528]	53	CHARACTER(len=45) :: cfile_grid_T = 'dyna_grid_T.nc' ! name of the grid_T file
	54	CHARACTER(len=45) :: cfile_grid_U = 'dyna_grid_U.nc' ! name of the grid_U file
	55	CHARACTER(len=45) :: cfile_grid_V = 'dyna_grid_V.nc' ! name of the grid_V file
	56	CHARACTER(len=45) :: cfile_grid_W = 'dyna_grid_W.nc' ! name of the grid_W file
[1735]	57
[2528]	58	REAL(wp) :: rnspdta ! number of time step per 2 consecutives data
	59	REAL(wp) :: rnspdta2 ! rnspdta * 0.5
[1735]	60
[2528]	61	INTEGER :: ndyn1, ndyn2 !
	62	INTEGER :: nlecoff = 0 ! switch for the first read
	63	INTEGER :: numfl_t, numfl_u, numfl_v, numfl_w
[325]	64
[2715]	65	REAL(wp), ALLOCATABLE, SAVE, DIMENSION(:,:,:,:) :: tdta ! temperature at two consecutive times
	66	REAL(wp), ALLOCATABLE, SAVE, DIMENSION(:,:,:,:) :: sdta ! salinity at two consecutive times
	67	REAL(wp), ALLOCATABLE, SAVE, DIMENSION(:,:,:,:) :: udta ! zonal velocity at two consecutive times
	68	REAL(wp), ALLOCATABLE, SAVE, DIMENSION(:,:,:,:) :: vdta ! meridional velocity at two consecutive times
	69	REAL(wp), ALLOCATABLE, SAVE, DIMENSION(:,:,:,:) :: wdta ! vertical velocity at two consecutive times
	70	REAL(wp), ALLOCATABLE, SAVE, DIMENSION(:,:,:,:) :: avtdta ! vertical diffusivity coefficient
[325]	71
[2715]	72	REAL(wp), ALLOCATABLE, SAVE, DIMENSION(:,:,:) :: hmlddta ! mixed layer depth at two consecutive times
	73	REAL(wp), ALLOCATABLE, SAVE, DIMENSION(:,:,:) :: wspddta ! wind speed at two consecutive times
	74	REAL(wp), ALLOCATABLE, SAVE, DIMENSION(:,:,:) :: frlddta ! sea-ice fraction at two consecutive times
	75	REAL(wp), ALLOCATABLE, SAVE, DIMENSION(:,:,:) :: empdta ! E-P at two consecutive times
	76	REAL(wp), ALLOCATABLE, SAVE, DIMENSION(:,:,:) :: qsrdta ! short wave heat flux at two consecutive times
	77	REAL(wp), ALLOCATABLE, SAVE, DIMENSION(:,:,:) :: bblxdta ! frequency of bbl in the x direction at 2 consecutive times
	78	REAL(wp), ALLOCATABLE, SAVE, DIMENSION(:,:,:) :: bblydta ! frequency of bbl in the y direction at 2 consecutive times
[2528]	79	LOGICAL :: l_offbbl
[2762]	80	#if defined key_ldfslp && ! defined key_c1d
[2715]	81	REAL(wp), ALLOCATABLE, SAVE, DIMENSION(:,:,:,:) :: uslpdta ! zonal isopycnal slopes
	82	REAL(wp), ALLOCATABLE, SAVE, DIMENSION(:,:,:,:) :: vslpdta ! meridional isopycnal slopes
	83	REAL(wp), ALLOCATABLE, SAVE, DIMENSION(:,:,:,:) :: wslpidta ! zonal diapycnal slopes
	84	REAL(wp), ALLOCATABLE, SAVE, DIMENSION(:,:,:,:) :: wslpjdta ! meridional diapycnal slopes
[325]	85	#endif
[2528]	86	#if ! defined key_degrad && defined key_traldf_c2d && defined key_traldf_eiv
[2715]	87	REAL(wp), ALLOCATABLE, SAVE, DIMENSION(:,:,:) :: aeiwdta ! G&M coefficient
[1501]	88	#endif
[2528]	89	#if defined key_degrad
[2715]	90	REAL(wp), ALLOCATABLE, SAVE, DIMENSION(:,:,:,:) :: ahtudta, ahtvdta, ahtwdta ! Lateral diffusivity
[2528]	91	# if defined key_traldf_eiv
[2715]	92	REAL(wp), ALLOCATABLE, SAVE, DIMENSION(:,:,:,:) :: aeiudta, aeivdta, aeiwdta ! G&M coefficient
[495]	93	# endif
[446]	94	#endif
	95
[325]	96	!! * Substitutions
	97	# include "domzgr_substitute.h90"
	98	# include "vectopt_loop_substitute.h90"
[343]	99	!!----------------------------------------------------------------------
[2528]	100	!! NEMO/OFF 3.3 , NEMO Consortium (2010)
	101	!! $Id$
	102	!! Software governed by the CeCILL licence (NEMOGCM/NEMO_CeCILL.txt)
[343]	103	!!----------------------------------------------------------------------
[325]	104	CONTAINS
	105
[1501]	106	SUBROUTINE dta_dyn( kt )
[325]	107	!!----------------------------------------------------------------------
	108	!! * ROUTINE dta_dyn *
	109	!!
[2528]	110	!! ** Purpose : Prepares dynamics and physics fields from an NEMO run
	111	!! for an off-line simulation of passive tracers
[325]	112	!!
	113	!! ** Method : calculates the position of DATA to read READ DATA
	114	!! (example month changement) computes slopes IF needed
	115	!! interpolates DATA IF needed
[2528]	116	!!----------------------------------------------------------------------
	117	INTEGER, INTENT(in) :: kt ! ocean time-step index
[325]	118	!!
[2528]	119	INTEGER :: iper, iperm1, iswap, izt ! local integers
	120	REAL(wp) :: zt, zweigh ! local scalars
[325]	121	!!----------------------------------------------------------------------
	122
[2528]	123	zt = ( REAL(kt,wp) + rnspdta2 ) / rnspdta
	124	izt = INT( zt )
	125	zweigh = zt - REAL( INT(zt), wp )
[325]	126
[2528]	127	IF( lperdyn ) THEN ; iperm1 = MOD( izt, ndtadyn )
	128	ELSE ; iperm1 = MOD( izt, ndtatot - 1 ) + 1
[431]	129	ENDIF
[325]	130
	131	iper = iperm1 + 1
[1291]	132	IF( iperm1 == 0 ) THEN
	133	IF( lperdyn ) THEN
[325]	134	iperm1 = ndtadyn
	135	ELSE
[1291]	136	IF( lfirdyn ) THEN
[2528]	137	IF(lwp) WRITE (numout,*) 'dta_dyn: dynamic file is not periodic with or without interpolation &
	138	& we take the first value for the previous period iperm1 = 0 '
[325]	139	END IF
	140	END IF
	141	END IF
	142
	143	iswap = 0
	144
	145	! 1. First call lfirdyn = true
	146	! ----------------------------
	147
[1501]	148	IF( lfirdyn ) THEN
[2528]	149	ndyn1 = iperm1 ! store the information of the period read
[1501]	150	ndyn2 = iper
	151
[2528]	152	IF(lwp) THEN
	153	WRITE (numout,*) ' dynamics data read for the period ndyn1 =', ndyn1, &
	154	& ' and for the period ndyn2 = ', ndyn2
[1501]	155	WRITE (numout,*) ' time step is : ', kt
[2528]	156	WRITE (numout,*) ' we have ndtadyn = ', ndtadyn, ' records in the dynamic file for one year'
[1501]	157	END IF
	158	!
[2528]	159	CALL dynrea( kt, MAX( 1, iperm1) ) ! data read for the iperm1 period
	160	!
	161	CALL swap_dyn_data ! swap from record 2 to 1
	162	!
[1501]	163	iswap = 1 ! indicates swap
	164	!
[2528]	165	CALL dynrea( kt, iper ) ! data read for the iper period
	166	!
	167	lfirdyn = .FALSE. ! trace the first call
[325]	168	ENDIF
	169	!
[1501]	170	! And now what we have to do at every time step
[325]	171	! check the validity of the period in memory
	172	!
[1501]	173	IF( iperm1 /= ndyn1 ) THEN
[2528]	174	!
	175	IF( iperm1 == 0 ) THEN
	176	IF(lwp) THEN
[1501]	177	WRITE (numout,*) ' dynamic file is not periodic with periodic interpolation'
	178	WRITE (numout,*) ' we take the last value for the last period '
	179	WRITE (numout,*) ' iperm1 = 12, iper = 13 '
	180	ENDIF
	181	iperm1 = 12
	182	iper = 13
	183	ENDIF
	184	!
[2528]	185	CALL swap_dyn_data ! We have to prepare a new read of data : swap from record 2 to 1
[1501]	186	!
[2528]	187	iswap = 1 ! indicates swap
	188	!
[1501]	189	CALL dynrea( kt, iper ) ! data read for the iper period
[2528]	190	!
	191	ndyn1 = ndyn2 ! store the information of the period read
[1501]	192	ndyn2 = iper
[2528]	193	!
	194	IF(lwp) THEN
	195	WRITE (numout,*) ' dynamics data read for the period ndyn1 =', ndyn1, &
	196	& ' and for the period ndyn2 = ', ndyn2
[1501]	197	WRITE (numout,*) ' time step is : ', kt
	198	END IF
	199	!
	200	END IF
[325]	201	!
[1501]	202	! Compute the data at the given time step
	203	!----------------------------------------
[495]	204
[2528]	205	IF( nsptint == 0 ) THEN ! No space interpolation, data are probably correct
	206	! ! We have to initialize data if we have changed the period
	207	CALL assign_dyn_data
	208	ELSEIF( nsptint == 1 ) THEN ! linear interpolation
[1501]	209	CALL linear_interp_dyn_data( zweigh )
[2528]	210	ELSE ! other interpolation
[1501]	211	WRITE (numout,*) ' this kind of interpolation do not exist at the moment : we stop'
	212	STOP 'dtadyn'
[325]	213	END IF
[2528]	214	!
	215	CALL eos( tsn, rhd, rhop ) ! In any case, we need rhop
	216	!
	217	#if ! defined key_degrad && defined key_traldf_c2d
	218	! ! In case of 2D varying coefficients, we need aeiv and aeiu
	219	IF( lk_traldf_eiv ) CALL dta_eiv( kt ) ! eddy induced velocity coefficient
[446]	220	#endif
[2528]	221	!
	222	IF( .NOT. l_offbbl ) THEN ! Compute bbl coefficients if needed
	223	tsb(:,:,:,:) = tsn(:,:,:,:)
	224	CALL bbl( kt, 'TRC')
	225	END IF
[2762]	226	!
[2528]	227	IF(ln_ctl) THEN
	228	CALL prt_ctl(tab3d_1=tsn(:,:,:,jp_tem), clinfo1=' tn - : ', mask1=tmask, ovlap=1, kdim=jpk )
	229	CALL prt_ctl(tab3d_1=tsn(:,:,:,jp_sal), clinfo1=' sn - : ', mask1=tmask, ovlap=1, kdim=jpk )
	230	CALL prt_ctl(tab3d_1=un , clinfo1=' un - : ', mask1=tmask, ovlap=1, kdim=jpk )
	231	CALL prt_ctl(tab3d_1=vn , clinfo1=' vn - : ', mask1=tmask, ovlap=1, kdim=jpk )
	232	CALL prt_ctl(tab3d_1=wn , clinfo1=' wn - : ', mask1=tmask, ovlap=1, kdim=jpk )
	233	CALL prt_ctl(tab3d_1=avt , clinfo1=' kz - : ', mask1=tmask, ovlap=1, kdim=jpk )
	234	CALL prt_ctl(tab2d_1=fr_i , clinfo1=' fr_i - : ', mask1=tmask, ovlap=1 )
	235	CALL prt_ctl(tab2d_1=hmld , clinfo1=' hmld - : ', mask1=tmask, ovlap=1 )
	236	CALL prt_ctl(tab2d_1=emps , clinfo1=' emps - : ', mask1=tmask, ovlap=1 )
	237	CALL prt_ctl(tab2d_1=wndm , clinfo1=' wspd - : ', mask1=tmask, ovlap=1 )
	238	CALL prt_ctl(tab2d_1=qsr , clinfo1=' qsr - : ', mask1=tmask, ovlap=1 )
	239	ENDIF
	240	!
[325]	241	END SUBROUTINE dta_dyn
	242
[2528]	243
[2715]	244	INTEGER FUNCTION dta_dyn_alloc()
	245	!!---------------------------------------------------------------------
	246	!! * ROUTINE dta_dyn_alloc *
	247	!!---------------------------------------------------------------------
	248
	249	ALLOCATE( tdta (jpi,jpj,jpk,2), sdta (jpi,jpj,jpk,2), &
	250	& udta (jpi,jpj,jpk,2), vdta (jpi,jpj,jpk,2), &
	251	& wdta (jpi,jpj,jpk,2), avtdta (jpi,jpj,jpk,2), &
[2762]	252	#if defined key_ldfslp && ! defined key_c1d
[2715]	253	& uslpdta (jpi,jpj,jpk,2), vslpdta (jpi,jpj,jpk,2), &
	254	& wslpidta(jpi,jpj,jpk,2), wslpjdta(jpi,jpj,jpk,2), &
	255	#endif
	256	#if defined key_degrad
	257	& ahtudta (jpi,jpj,jpk,2), ahtvdta (jpi,jpj,jpk,2), &
	258	& ahtwdta (jpi,jpj,jpk,2), &
	259	# if defined key_traldf_eiv
	260	& aeiudta (jpi,jpj,jpk,2), aeivdta (jpi,jpj,jpk,2), &
	261	& aeiwdta (jpi,jpj,jpk,2), &
	262	# endif
	263	#endif
	264	#if ! defined key_degrad && defined key_traldf_c2d && defined key_traldf_eiv
	265	& aeiwdta (jpi,jpj, 2), &
	266	#endif
	267	& hmlddta (jpi,jpj, 2), wspddta (jpi,jpj, 2), &
	268	& frlddta (jpi,jpj, 2), qsrdta (jpi,jpj, 2), &
	269	& empdta (jpi,jpj, 2), STAT=dta_dyn_alloc )
	270	!
	271	IF( dta_dyn_alloc /= 0 ) CALL ctl_warn('dta_dyn_alloc: failed to allocate facvol array')
	272	!
	273	END FUNCTION dta_dyn_alloc
	274
	275
[325]	276	SUBROUTINE dynrea( kt, kenr )
	277	!!----------------------------------------------------------------------
	278	!! * ROUTINE dynrea *
	279	!!
	280	!! ** Purpose : READ dynamics fiels from OPA9 netcdf output
	281	!!
[2528]	282	!! ** Method : READ the kenr records of DATA and store in udta(...,2), ....
[325]	283	!!----------------------------------------------------------------------
[2715]	284	USE wrk_nemo, ONLY: wrk_in_use, wrk_not_released
[2762]	285	USE wrk_nemo, ONLY: zu => wrk_3d_3 , zv => wrk_3d_4 , zw => wrk_3d_5
	286	USE wrk_nemo, ONLY: zt => wrk_3d_6 , zs => wrk_3d_7 , zavt => wrk_3d_8
	287	USE wrk_nemo, ONLY: zemp => wrk_2d_11 , zqsr => wrk_2d_12, zmld => wrk_2d_13
	288	USE wrk_nemo, ONLY: zice => wrk_2d_14 , zwspd => wrk_2d_15
	289	USE wrk_nemo, ONLY: ztaux => wrk_2d_16 , ztauy => wrk_2d_17
	290	USE wrk_nemo, ONLY: zbblx => wrk_2d_18 , zbbly => wrk_2d_19
[2715]	291	USE wrk_nemo, ONLY: zaeiw2d => wrk_2d_10
[2762]	292	USE wrk_nemo, ONLY: ztsn => wrk_4d_1
[2715]	293	!
[2528]	294	INTEGER, INTENT(in) :: kt, kenr ! time index
	295	!!
[1501]	296	INTEGER :: jkenr
[2762]	297	#if defined key_degrad
	298	REAL(wp), ALLOCATABLE, DIMENSION(:,:,:) :: zahtu, zahtv, zahtw ! Lateral diffusivity
	299	# if defined key_traldf_eiv
	300	REAL(wp), ALLOCATABLE, DIMENSION(:,:,:) :: zaeiu, zaeiv, zaeiw ! G&M coefficient
	301	# endif
	302	#endif
[2528]	303	!!----------------------------------------------------------------------
[2715]	304	!
[2762]	305	IF( wrk_in_use(3, 3,4,5,6,7,8) .OR. &
	306	wrk_in_use(4, 1) .OR. &
	307	wrk_in_use(2, 10,11,12,13,14,15,16,17,18,19) ) THEN
[2715]	308	CALL ctl_stop('domrea/dta_dyn: requested workspace arrays unavailable') ; RETURN
	309	ENDIF
[2762]	310
	311	#if defined key_degrad
	312	ALLOCATE( zahtu(jpi,jpj,jpk), zahtv(jpi,jpj,jpk), zahtw(jpi,jpj,jpk) )
	313	# if defined key_traldf_eiv
	314	ALLOCATE( zaeiu(jpi,jpj,jpk), zaeiv(jpi,jpj,jpk), zaeiw(jpi,jpj,jpk) )
	315	# endif
	316	#endif
[1323]	317
[325]	318	! cas d'un fichier non periodique : on utilise deux fois le premier et
	319	! le dernier champ temporel
	320
	321	jkenr = kenr
	322
	323	IF(lwp) THEN
	324	WRITE(numout,*)
[2528]	325	WRITE(numout,*) 'Dynrea : read dynamical fields, kenr = ', jkenr
	326	WRITE(numout,*) '~~~~~~~'
	327	#if defined key_degrad
[495]	328	WRITE(numout,*) ' Degraded fields'
	329	#endif
[325]	330	WRITE(numout,*)
	331	ENDIF
	332
	333
	334	IF( kt == nit000 .AND. nlecoff == 0 ) THEN
	335	nlecoff = 1
[1501]	336	CALL iom_open ( cfile_grid_T, numfl_t )
	337	CALL iom_open ( cfile_grid_U, numfl_u )
	338	CALL iom_open ( cfile_grid_V, numfl_v )
	339	CALL iom_open ( cfile_grid_W, numfl_w )
[325]	340	ENDIF
	341
[495]	342	! file grid-T
	343	!---------------
[1501]	344	CALL iom_get( numfl_t, jpdom_data, 'votemper', zt (:,:,:), jkenr )
	345	CALL iom_get( numfl_t, jpdom_data, 'vosaline', zs (:,:,:), jkenr )
	346	CALL iom_get( numfl_t, jpdom_data, 'somixhgt', zmld (:,: ), jkenr )
	347	CALL iom_get( numfl_t, jpdom_data, 'sowaflcd', zemp (:,: ), jkenr )
	348	CALL iom_get( numfl_t, jpdom_data, 'soshfldo', zqsr (:,: ), jkenr )
	349	CALL iom_get( numfl_t, jpdom_data, 'soicecov', zice (:,: ), jkenr )
	350	IF( iom_varid( numfl_t, 'sowindsp', ldstop = .FALSE. ) > 0 ) THEN
	351	CALL iom_get( numfl_t, jpdom_data, 'sowindsp', zwspd(:,:), jkenr )
	352	ELSE
	353	CALL iom_get( numfl_u, jpdom_data, 'sozotaux', ztaux(:,:), jkenr )
	354	CALL iom_get( numfl_v, jpdom_data, 'sometauy', ztauy(:,:), jkenr )
	355	CALL tau2wnd( ztaux, ztauy, zwspd )
	356	ENDIF
	357	! files grid-U / grid_V
	358	CALL iom_get( numfl_u, jpdom_data, 'vozocrtx', zu (:,:,:), jkenr )
	359	CALL iom_get( numfl_v, jpdom_data, 'vomecrty', zv (:,:,:), jkenr )
[2559]	360	#if defined key_trabbl
	361	IF( .NOT. lk_c1d .AND. nn_bbl_ldf == 1 ) THEN
[2528]	362	IF( iom_varid( numfl_u, 'sobblcox', ldstop = .FALSE. ) > 0 .AND. &
	363	& iom_varid( numfl_v, 'sobblcoy', ldstop = .FALSE. ) > 0 ) THEN
	364	CALL iom_get( numfl_u, jpdom_data, 'sobblcox', zbblx(:,:), jkenr )
	365	CALL iom_get( numfl_v, jpdom_data, 'sobblcoy', zbbly(:,:), jkenr )
	366	l_offbbl = .TRUE.
	367	ENDIF
[1501]	368	ENDIF
[2559]	369	#endif
[325]	370
[495]	371	! file grid-W
[2762]	372	! CALL iom_get ( numfl_w, jpdom_data, 'vovecrtz', zw (:,:,:), jkenr )
[1501]	373	! Computation of vertical velocity using horizontal divergence
[2762]	374	CALL wzv( zu, zv, zw )
[1501]	375
[2528]	376	IF( iom_varid( numfl_w, 'voddmavs', ldstop = .FALSE. ) > 0 ) THEN ! avs exist: it is used
	377	CALL iom_get( numfl_w, jpdom_data, 'voddmavs', zavt (:,:,:), jkenr )
	378	ELSE ! no avs: use avt
	379	CALL iom_get( numfl_w, jpdom_data, 'votkeavt', zavt (:,:,:), jkenr )
	380	ENDIF
[325]	381
[2528]	382	#if ! defined key_degrad && defined key_traldf_c2d && defined key_traldf_eiv
[2715]	383	CALL iom_get( numfl_w, jpdom_data, 'soleaeiw', zaeiw2d(:,: ), jkenr )
[1501]	384	#endif
[446]	385
[2528]	386	#if defined key_degrad
[1501]	387	CALL iom_get( numfl_u, jpdom_data, 'vozoahtu', zahtu(:,:,:), jkenr )
	388	CALL iom_get( numfl_v, jpdom_data, 'vomeahtv', zahtv(:,:,:), jkenr )
	389	CALL iom_get( numfl_w, jpdom_data, 'voveahtw', zahtw(:,:,:), jkenr )
[2528]	390	# if defined key_traldf_eiv
[1501]	391	CALL iom_get( numfl_u, jpdom_data, 'vozoaeiu', zaeiu(:,:,:), jkenr )
	392	CALL iom_get( numfl_v, jpdom_data, 'vomeaeiv', zaeiv(:,:,:), jkenr )
	393	CALL iom_get( numfl_w, jpdom_data, 'voveaeiw', zaeiw(:,:,:), jkenr )
[495]	394	# endif
[446]	395	#endif
	396
[2762]	397	udta (:,:,:,2) = zu (:,:,:) * umask(:,:,:)
	398	vdta (:,:,:,2) = zv (:,:,:) * vmask(:,:,:)
	399	wdta (:,:,:,2) = zw (:,:,:) * tmask(:,:,:)
	400	tdta (:,:,:,2) = zt (:,:,:) * tmask(:,:,:)
	401	sdta (:,:,:,2) = zs (:,:,:) * tmask(:,:,:)
[1501]	402	avtdta(:,:,:,2) = zavt(:,:,:) * tmask(:,:,:)
[612]	403
[2762]	404	#if defined key_ldfslp && ! defined key_c1d
	405	! Computes slopes (here tsn and avt are used as workspace)
	406	ztsn (:,:,:,jp_tem) = tdta (:,:,:,2)
	407	ztsn (:,:,:,jp_sal) = sdta (:,:,:,2)
	408	avt(:,:,:) = avtdta(:,:,:,2)
	409
	410	CALL eos( ztsn, rhd, rhop ) ! Time-filtered in situ density
	411	CALL bn2( ztsn, rn2 ) ! before Brunt-Vaisala frequency
	412	IF( ln_zps ) &
	413	& CALL zps_hde( kt, jpts, ztsn, gtsu, gtsv, & ! Partial steps: before Horizontal DErivative
	414	& rhd, gru , grv ) ! of t, s, rd at the bottom ocean level
	415	CALL zdf_mxl( kt ) ! mixed layer depth
	416	CALL ldf_slp( kt, rhd, rn2 )
	417
	418	uslpdta (:,:,:,2) = uslp (:,:,:)
	419	vslpdta (:,:,:,2) = vslp (:,:,:)
	420	wslpidta(:,:,:,2) = wslpi(:,:,:)
	421	wslpjdta(:,:,:,2) = wslpj(:,:,:)
	422	#endif
	423
[2528]	424	#if ! defined key_degrad && defined key_traldf_c2d && defined key_traldf_eiv
[2715]	425	aeiwdta(:,:,2) = zaeiw2d(:,:) * tmask(:,:,1)
[446]	426	#endif
[495]	427
[2528]	428	#if defined key_degrad
[495]	429	ahtudta(:,:,:,2) = zahtu(:,:,:) * umask(:,:,:)
	430	ahtvdta(:,:,:,2) = zahtv(:,:,:) * vmask(:,:,:)
	431	ahtwdta(:,:,:,2) = zahtw(:,:,:) * tmask(:,:,:)
[2528]	432	# if defined key_traldf_eiv
[495]	433	aeiudta(:,:,:,2) = zaeiu(:,:,:) * umask(:,:,:)
	434	aeivdta(:,:,:,2) = zaeiv(:,:,:) * vmask(:,:,:)
	435	aeiwdta(:,:,:,2) = zaeiw(:,:,:) * tmask(:,:,:)
	436	# endif
[325]	437	#endif
	438
[1501]	439	! fluxes
[325]	440	!
[1501]	441	wspddta(:,:,2) = zwspd(:,:) * tmask(:,:,1)
	442	frlddta(:,:,2) = MIN( 1., zice(:,:) ) * tmask(:,:,1)
	443	empdta (:,:,2) = zemp(:,:) * tmask(:,:,1)
	444	qsrdta (:,:,2) = zqsr(:,:) * tmask(:,:,1)
	445	hmlddta(:,:,2) = zmld(:,:) * tmask(:,:,1)
[2528]	446
[2559]	447	#if defined key_trabbl
[2528]	448	IF( l_offbbl ) THEN
	449	bblxdta(:,:,2) = MAX( 0., zbblx(:,:) )
	450	bblydta(:,:,2) = MAX( 0., zbbly(:,:) )
	451	WHERE( bblxdta(:,:,2) > 2. ) bblxdta(:,:,2) = 0.
	452	WHERE( bblydta(:,:,2) > 2. ) bblydta(:,:,2) = 0.
	453	ENDIF
[2559]	454	#endif
[495]	455
	456	IF( kt == nitend ) THEN
	457	CALL iom_close ( numfl_t )
	458	CALL iom_close ( numfl_u )
	459	CALL iom_close ( numfl_v )
	460	CALL iom_close ( numfl_w )
	461	ENDIF
[2528]	462	!
[2762]	463	IF( wrk_not_released(3, 3,4,5,6,7,8) .OR. &
	464	wrk_not_released(4, 1 ) .OR. &
	465	wrk_not_released(2, 10,11,12,13,14,15,16,17,18,19) ) THEN
[2715]	466	CALL ctl_stop('domrea/dta_dyn: failed to release workspace arrays')
	467	END IF
[2762]	468	#if defined key_degrad
	469	DEALLOCATE( zahtu ) ; DEALLOCATE( zahtv ) ; DEALLOCATE( zahtw )
	470	# if defined key_traldf_eiv
	471	DEALLOCATE( zaeiu ) ; DEALLOCATE( zaeiv ) ; DEALLOCATE( zaeiw )
	472	# endif
	473	#endif
[2715]	474	!
[325]	475	END SUBROUTINE dynrea
	476
[2528]	477
[1501]	478	SUBROUTINE dta_dyn_init
	479	!!----------------------------------------------------------------------
	480	!! * ROUTINE dta_dyn_init *
	481	!!
	482	!! ** Purpose : initializations of parameters for the interpolation
	483	!!
	484	!! ** Method :
	485	!!----------------------------------------------------------------------
[2715]	486	REAL(wp) :: znspyr !: number of time step per year
	487	!
[2528]	488	NAMELIST/namdyn/ ndtadyn, ndtatot, nsptint, lperdyn, &
[2715]	489	& cfile_grid_T, cfile_grid_U, cfile_grid_V, cfile_grid_W
[1501]	490	!!----------------------------------------------------------------------
[2715]	491	!
	492	IF( dta_dyn_alloc() /= 0 ) CALL ctl_stop( 'STOP', 'dta_dyn_alloc: unable to allocate standard ocean arrays' )
	493	!
[2528]	494	REWIND( numnam ) ! Read Namelist namdyn : Lateral physics on tracers
[1501]	495	READ ( numnam, namdyn )
[2715]	496	!
[2528]	497	IF(lwp) THEN ! control print
[1501]	498	WRITE(numout,*)
	499	WRITE(numout,*) 'namdyn : offline dynamical selection'
	500	WRITE(numout,*) '~~~~~~~'
	501	WRITE(numout,*) ' Namelist namdyn : set parameters for the lecture of the dynamical fields'
	502	WRITE(numout,*)
	503	WRITE(numout,*) ' number of elements in the FILE for a year ndtadyn = ' , ndtadyn
	504	WRITE(numout,*) ' total number of elements in the FILE ndtatot = ' , ndtatot
	505	WRITE(numout,*) ' type of interpolation nsptint = ' , nsptint
	506	WRITE(numout,*) ' loop on the same FILE lperdyn = ' , lperdyn
	507	WRITE(numout,*) ' '
	508	WRITE(numout,*) ' name of grid_T file cfile_grid_T = ', TRIM(cfile_grid_T)
	509	WRITE(numout,*) ' name of grid_U file cfile_grid_U = ', TRIM(cfile_grid_U)
	510	WRITE(numout,*) ' name of grid_V file cfile_grid_V = ', TRIM(cfile_grid_V)
	511	WRITE(numout,*) ' name of grid_W file cfile_grid_W = ', TRIM(cfile_grid_W)
	512	WRITE(numout,*) ' '
	513	ENDIF
[2528]	514	!
[1735]	515	znspyr = nyear_len(1) * rday / rdt
[2715]	516	rnspdta = znspyr / REAL( ndtadyn, wp )
[1735]	517	rnspdta2 = rnspdta * 0.5
[2528]	518	!
	519	CALL dta_dyn( nit000 )
	520	!
[1501]	521	END SUBROUTINE dta_dyn_init
	522
[2528]	523
[2762]	524	SUBROUTINE wzv( pu, pv, pw )
[1501]	525	!!----------------------------------------------------------------------
	526	!! * ROUTINE wzv *
	527	!!
	528	!! ** Purpose : Compute the now vertical velocity after the array swap
	529	!!
[2528]	530	!! ** Method : - compute the now divergence given by :
	531	!! * z-coordinate ONLY !!!!
[1501]	532	!! hdiv = 1/(e1t*e2t) [ di(e2u u) + dj(e1v v) ]
	533	!! - Using the incompressibility hypothesis, the vertical
	534	!! velocity is computed by integrating the horizontal divergence
	535	!! from the bottom to the surface.
[2528]	536	!! The boundary conditions are w=0 at the bottom (no flux).
	537	!!----------------------------------------------------------------------
	538	REAL(wp), DIMENSION(jpi,jpj,jpk), INTENT(in ) :: pu, pv !: horizontal velocities
	539	REAL(wp), DIMENSION(jpi,jpj,jpk), INTENT( out) :: pw !: verticla velocity
[1501]	540	!!
	541	INTEGER :: ji, jj, jk
	542	REAL(wp) :: zu, zu1, zv, zv1, zet
[2762]	543	REAL(wp), DIMENSION(jpi,jpj,jpk) :: zhdiv !: horizontal divergence
[2528]	544	!!----------------------------------------------------------------------
	545	!
[1501]	546	! Computation of vertical velocity using horizontal divergence
[2762]	547	zhdiv(:,:,:) = 0.
[1501]	548	DO jk = 1, jpkm1
	549	DO jj = 2, jpjm1
	550	DO ji = fs_2, fs_jpim1 ! vector opt.
	551	zu = pu(ji ,jj ,jk) * umask(ji ,jj ,jk) * e2u(ji ,jj ) * fse3u(ji ,jj ,jk)
	552	zu1 = pu(ji-1,jj ,jk) * umask(ji-1,jj ,jk) * e2u(ji-1,jj ) * fse3u(ji-1,jj ,jk)
	553	zv = pv(ji ,jj ,jk) * vmask(ji ,jj ,jk) * e1v(ji ,jj ) * fse3v(ji ,jj ,jk)
	554	zv1 = pv(ji ,jj-1,jk) * vmask(ji ,jj-1,jk) * e1v(ji ,jj-1) * fse3v(ji ,jj-1,jk)
	555	zet = 1. / ( e1t(ji,jj) * e2t(ji,jj) * fse3t(ji,jj,jk) )
[2762]	556	zhdiv(ji,jj,jk) = ( zu - zu1 + zv - zv1 ) * zet
[1501]	557	END DO
	558	END DO
[2528]	559	END DO
[2762]	560	CALL lbc_lnk( zhdiv, 'T', 1. ) ! Lateral boundary conditions on zhdiv
[2528]	561	!
[1501]	562	! computation of vertical velocity from the bottom
[2528]	563	pw(:,:,jpk) = 0._wp
[1501]	564	DO jk = jpkm1, 1, -1
[2762]	565	pw(:,:,jk) = pw(:,:,jk+1) - fse3t(:,:,jk) * zhdiv(:,:,jk)
[1501]	566	END DO
[2528]	567	!
[1501]	568	END SUBROUTINE wzv
	569
[2528]	570
	571	SUBROUTINE dta_eiv( kt )
	572	!!----------------------------------------------------------------------
	573	!! * ROUTINE dta_eiv *
	574	!!
	575	!! ** Purpose : Compute the eddy induced velocity coefficient from the
	576	!! growth rate of baroclinic instability.
	577	!!
	578	!! ** Method : Specific to the offline model. Computes the horizontal
	579	!! values from the vertical value
	580	!!----------------------------------------------------------------------
	581	INTEGER, INTENT( in ) :: kt ! ocean time-step inedx
	582	!!
	583	INTEGER :: ji, jj ! dummy loop indices
	584	!!----------------------------------------------------------------------
	585	!
	586	IF( kt == nit000 ) THEN
	587	IF(lwp) WRITE(numout,*)
	588	IF(lwp) WRITE(numout,*) 'dta_eiv : eddy induced velocity coefficients'
	589	IF(lwp) WRITE(numout,*) '~~~~~~~'
	590	ENDIF
	591	!
[2559]	592	#if defined key_ldfeiv
[2528]	593	! Average the diffusive coefficient at u- v- points
	594	DO jj = 2, jpjm1
	595	DO ji = fs_2, fs_jpim1 ! vector opt.
	596	aeiu(ji,jj) = .5 * ( aeiw(ji,jj) + aeiw(ji+1,jj ) )
	597	aeiv(ji,jj) = .5 * ( aeiw(ji,jj) + aeiw(ji ,jj+1) )
	598	END DO
	599	END DO
	600	CALL lbc_lnk( aeiu, 'U', 1. ) ; CALL lbc_lnk( aeiv, 'V', 1. ) ! lateral boundary condition
[2559]	601	#endif
[2528]	602	!
	603	END SUBROUTINE dta_eiv
	604
	605
[1501]	606	SUBROUTINE tau2wnd( ptaux, ptauy, pwspd )
	607	!!---------------------------------------------------------------------
	608	!! * ROUTINE sbc_tau2wnd *
	609	!!
	610	!! ** Purpose : Estimation of wind speed as a function of wind stress
	611	!!
	612	!! ** Method : \|tau\|=rhoaCd\|U\|^2
	613	!!---------------------------------------------------------------------
[2528]	614	REAL(wp), DIMENSION(jpi,jpj), INTENT(in ) :: ptaux, ptauy ! wind stress in i-j direction resp.
	615	REAL(wp), DIMENSION(jpi,jpj), INTENT( out) :: pwspd ! wind speed
	616	!!
	617	REAL(wp) :: zrhoa = 1.22_wp ! Air density kg/m3
	618	REAL(wp) :: zcdrag = 1.5e-3_wp ! drag coefficient
	619	REAL(wp) :: ztx, zty, ztau, zcoef ! temporary variables
	620	INTEGER :: ji, jj ! dummy indices
[1501]	621	!!---------------------------------------------------------------------
[1643]	622	zcoef = 1. / ( zrhoa * zcdrag )
[1501]	623	!CDIR NOVERRCHK
[1699]	624	DO jj = 2, jpjm1
[1501]	625	!CDIR NOVERRCHK
[1699]	626	DO ji = fs_2, fs_jpim1 ! vector opt.
	627	ztx = ptaux(ji,jj) * umask(ji,jj,1) + ptaux(ji-1,jj ) * umask(ji-1,jj ,1)
	628	zty = ptauy(ji,jj) * vmask(ji,jj,1) + ptauy(ji ,jj-1) * vmask(ji ,jj-1,1)
	629	ztau = 0.5 * SQRT( ztx * ztx + zty * zty )
[1501]	630	pwspd(ji,jj) = SQRT ( ztau * zcoef ) * tmask(ji,jj,1)
	631	END DO
	632	END DO
[1699]	633	CALL lbc_lnk( pwspd(:,:), 'T', 1. )
[2528]	634	!
[1501]	635	END SUBROUTINE tau2wnd
	636
	637
	638	SUBROUTINE swap_dyn_data
	639	!!----------------------------------------------------------------------
	640	!! * ROUTINE swap_dyn_data *
	641	!!
	642	!! ** Purpose : swap array data
	643	!!----------------------------------------------------------------------
[2528]	644	!
[1501]	645	! swap from record 2 to 1
	646	tdta (:,:,:,1) = tdta (:,:,:,2)
	647	sdta (:,:,:,1) = sdta (:,:,:,2)
	648	avtdta (:,:,:,1) = avtdta (:,:,:,2)
	649	udta (:,:,:,1) = udta (:,:,:,2)
	650	vdta (:,:,:,1) = vdta (:,:,:,2)
	651	wdta (:,:,:,1) = wdta (:,:,:,2)
[2528]	652	#if defined key_ldfslp && ! defined key_c1d
[1501]	653	uslpdta (:,:,:,1) = uslpdta (:,:,:,2)
	654	vslpdta (:,:,:,1) = vslpdta (:,:,:,2)
	655	wslpidta(:,:,:,1) = wslpidta(:,:,:,2)
	656	wslpjdta(:,:,:,1) = wslpjdta(:,:,:,2)
	657	#endif
	658	hmlddta(:,:,1) = hmlddta(:,:,2)
	659	wspddta(:,:,1) = wspddta(:,:,2)
	660	frlddta(:,:,1) = frlddta(:,:,2)
	661	empdta (:,:,1) = empdta (:,:,2)
	662	qsrdta (:,:,1) = qsrdta (:,:,2)
[2528]	663	IF( l_offbbl ) THEN
	664	bblxdta(:,:,1) = bblxdta(:,:,2)
	665	bblydta(:,:,1) = bblydta(:,:,2)
	666	ENDIF
[1501]	667
[2528]	668	#if ! defined key_degrad && defined key_traldf_c2d && defined key_traldf_eiv
[1501]	669	aeiwdta(:,:,1) = aeiwdta(:,:,2)
	670	#endif
	671
[2528]	672	#if defined key_degrad
[1501]	673	ahtudta(:,:,:,1) = ahtudta(:,:,:,2)
	674	ahtvdta(:,:,:,1) = ahtvdta(:,:,:,2)
	675	ahtwdta(:,:,:,1) = ahtwdta(:,:,:,2)
[2528]	676	# if defined key_traldf_eiv
[1501]	677	aeiudta(:,:,:,1) = aeiudta(:,:,:,2)
	678	aeivdta(:,:,:,1) = aeivdta(:,:,:,2)
	679	aeiwdta(:,:,:,1) = aeiwdta(:,:,:,2)
	680	# endif
	681	#endif
[2528]	682	!
	683	END SUBROUTINE swap_dyn_data
[1501]	684
	685
	686	SUBROUTINE assign_dyn_data
	687	!!----------------------------------------------------------------------
	688	!! * ROUTINE assign_dyn_data *
	689	!!
	690	!! ** Purpose : Assign dynamical data to the data that have been read
	691	!! without time interpolation
	692	!!
	693	!!----------------------------------------------------------------------
	694
[2528]	695	tsn(:,:,:,jp_tem) = tdta (:,:,:,2)
	696	tsn(:,:,:,jp_sal) = sdta (:,:,:,2)
	697	avt(:,:,:) = avtdta(:,:,:,2)
[1501]	698
	699	un (:,:,:) = udta (:,:,:,2)
	700	vn (:,:,:) = vdta (:,:,:,2)
	701	wn (:,:,:) = wdta (:,:,:,2)
	702
[2528]	703	#if defined key_ldfslp && ! defined key_c1d
[1501]	704	uslp (:,:,:) = uslpdta (:,:,:,2)
	705	vslp (:,:,:) = vslpdta (:,:,:,2)
	706	wslpi(:,:,:) = wslpidta(:,:,:,2)
	707	wslpj(:,:,:) = wslpjdta(:,:,:,2)
	708	#endif
	709
	710	hmld(:,:) = hmlddta(:,:,2)
	711	wndm(:,:) = wspddta(:,:,2)
	712	fr_i(:,:) = frlddta(:,:,2)
	713	emp (:,:) = empdta (:,:,2)
	714	emps(:,:) = emp(:,:)
	715	qsr (:,:) = qsrdta (:,:,2)
[2559]	716	#if defined key_trabbl
[2528]	717	IF( l_offbbl ) THEN
	718	ahu_bbl(:,:) = bblxdta(:,:,2)
	719	ahv_bbl(:,:) = bblydta(:,:,2)
	720	ENDIF
[2559]	721	#endif
[2528]	722	#if ! defined key_degrad && defined key_traldf_c2d && defined key_traldf_eiv
[1501]	723	aeiw(:,:) = aeiwdta(:,:,2)
	724	#endif
	725
[2528]	726	#if defined key_degrad
[1501]	727	ahtu(:,:,:) = ahtudta(:,:,:,2)
	728	ahtv(:,:,:) = ahtvdta(:,:,:,2)
	729	ahtw(:,:,:) = ahtwdta(:,:,:,2)
[2528]	730	# if defined key_traldf_eiv
[1501]	731	aeiu(:,:,:) = aeiudta(:,:,:,2)
	732	aeiv(:,:,:) = aeivdta(:,:,:,2)
	733	aeiw(:,:,:) = aeiwdta(:,:,:,2)
	734	# endif
	735	#endif
[2528]	736	!
[1501]	737	END SUBROUTINE assign_dyn_data
	738
[2528]	739
[1501]	740	SUBROUTINE linear_interp_dyn_data( pweigh )
	741	!!----------------------------------------------------------------------
[2528]	742	!! * ROUTINE linear_interp_dyn_data *
[1501]	743	!!
	744	!! ** Purpose : linear interpolation of data
[2528]	745	!!----------------------------------------------------------------------
	746	REAL(wp), INTENT(in) :: pweigh ! weigh
[1501]	747	!!
	748	REAL(wp) :: zweighm1
	749	!!----------------------------------------------------------------------
	750
	751	zweighm1 = 1. - pweigh
	752
[2528]	753	tsn(:,:,:,jp_tem) = zweighm1 * tdta (:,:,:,1) + pweigh * tdta (:,:,:,2)
	754	tsn(:,:,:,jp_sal) = zweighm1 * sdta (:,:,:,1) + pweigh * sdta (:,:,:,2)
	755	avt(:,:,:) = zweighm1 * avtdta(:,:,:,1) + pweigh * avtdta(:,:,:,2)
[1501]	756
	757	un (:,:,:) = zweighm1 * udta (:,:,:,1) + pweigh * udta (:,:,:,2)
	758	vn (:,:,:) = zweighm1 * vdta (:,:,:,1) + pweigh * vdta (:,:,:,2)
	759	wn (:,:,:) = zweighm1 * wdta (:,:,:,1) + pweigh * wdta (:,:,:,2)
	760
[2528]	761	#if defined key_ldfslp && ! defined key_c1d
[1501]	762	uslp (:,:,:) = zweighm1 * uslpdta (:,:,:,1) + pweigh * uslpdta (:,:,:,2)
	763	vslp (:,:,:) = zweighm1 * vslpdta (:,:,:,1) + pweigh * vslpdta (:,:,:,2)
	764	wslpi(:,:,:) = zweighm1 * wslpidta(:,:,:,1) + pweigh * wslpidta(:,:,:,2)
	765	wslpj(:,:,:) = zweighm1 * wslpjdta(:,:,:,1) + pweigh * wslpjdta(:,:,:,2)
	766	#endif
	767
	768	hmld(:,:) = zweighm1 * hmlddta(:,:,1) + pweigh * hmlddta(:,:,2)
	769	wndm(:,:) = zweighm1 * wspddta(:,:,1) + pweigh * wspddta(:,:,2)
	770	fr_i(:,:) = zweighm1 * frlddta(:,:,1) + pweigh * frlddta(:,:,2)
	771	emp (:,:) = zweighm1 * empdta (:,:,1) + pweigh * empdta (:,:,2)
	772	emps(:,:) = emp(:,:)
	773	qsr (:,:) = zweighm1 * qsrdta (:,:,1) + pweigh * qsrdta (:,:,2)
[2559]	774	#if defined key_trabbl
[2528]	775	IF( l_offbbl ) THEN
	776	ahu_bbl(:,:) = zweighm1 * bblxdta(:,:,1) + pweigh * bblxdta(:,:,2)
	777	ahv_bbl(:,:) = zweighm1 * bblydta(:,:,1) + pweigh * bblydta(:,:,2)
	778	ENDIF
[2559]	779	#endif
[1501]	780
[2528]	781	#if ! defined key_degrad && defined key_traldf_c2d && defined key_traldf_eiv
[1501]	782	aeiw(:,:) = zweighm1 * aeiwdta(:,:,1) + pweigh * aeiwdta(:,:,2)
	783	#endif
	784
[2528]	785	#if defined key_degrad
[1501]	786	ahtu(:,:,:) = zweighm1 * ahtudta(:,:,:,1) + pweigh * ahtudta(:,:,:,2)
	787	ahtv(:,:,:) = zweighm1 * ahtvdta(:,:,:,1) + pweigh * ahtvdta(:,:,:,2)
	788	ahtw(:,:,:) = zweighm1 * ahtwdta(:,:,:,1) + pweigh * ahtwdta(:,:,:,2)
[2528]	789	# if defined key_traldf_eiv
[1501]	790	aeiu(:,:,:) = zweighm1 * aeiudta(:,:,:,1) + pweigh * aeiudta(:,:,:,2)
	791	aeiv(:,:,:) = zweighm1 * aeivdta(:,:,:,1) + pweigh * aeivdta(:,:,:,2)
	792	aeiw(:,:,:) = zweighm1 * aeiwdta(:,:,:,1) + pweigh * aeiwdta(:,:,:,2)
	793	# endif
	794	#endif
[2528]	795	!
[1501]	796	END SUBROUTINE linear_interp_dyn_data
	797
[2528]	798	!!======================================================================
[325]	799	END MODULE dtadyn

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

Download in other formats: