source: flexpart.git/src/interpol_vdep_nests.f90

10.4.1_peseiGFS_025bugfixes+enhancementsdevrelease-10release-10.4.1scaling-bug
Last change on this file was 92fab65, checked in by Ignacio Pisso <ip@…>, 4 years ago

add SPDX-License-Identifier to all .f90 files

  • Property mode set to 100644
File size: 2.4 KB
Line 
1! SPDX-FileCopyrightText: FLEXPART 1998-2019, see flexpart_license.txt
2! SPDX-License-Identifier: GPL-3.0-or-later
3
4subroutine interpol_vdep_nests(level,vdepo)
5  !                                 i     o
6  !****************************************************************************
7  !                                                                           *
8  !  Interpolation of the deposition velocity on 2-d model layer.             *
9  !  In horizontal direction bilinear interpolation interpolation is used.    *
10  !  Temporally a linear interpolation is used.                               *
11  !                                                                           *
12  !  1 first time                                                             *
13  !  2 second time                                                            *
14  !                                                                           *
15  !                                                                           *
16  !     Author: A. Stohl                                                      *
17  !                                                                           *
18  !     30 May 1994                                                           *
19  !                                                                           *
20  !****************************************************************************
21  !                                                                           *
22  ! Variables:                                                                *
23  !                                                                           *
24  ! level                number of species for which interpolation is done    *
25  !                                                                           *
26  !****************************************************************************
27
28  use par_mod
29  use com_mod
30  use interpol_mod
31
32  implicit none
33
34  integer :: level,indexh,m
35  real :: y(2),vdepo
36
37  ! a) Bilinear horizontal interpolation
38
39  do m=1,2
40    indexh=memind(m)
41
42    y(m)=p1*vdepn(ix ,jy ,level,indexh,ngrid) &
43         +p2*vdepn(ixp,jy ,level,indexh,ngrid) &
44         +p3*vdepn(ix ,jyp,level,indexh,ngrid) &
45         +p4*vdepn(ixp,jyp,level,indexh,ngrid)
46  end do
47
48
49  ! b) Temporal interpolation
50
51  vdepo=(y(1)*dt2+y(2)*dt1)*dtt
52
53  depoindicator(level)=.false.
54
55
56end subroutine interpol_vdep_nests
Note: See TracBrowser for help on using the repository browser.
hosted by ZAMG