-
Notifications
You must be signed in to change notification settings - Fork 0
/
SwanPrepComp.ftn90
96 lines (96 loc) · 3.06 KB
/
SwanPrepComp.ftn90
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
subroutine SwanPrepComp ( cross )
!
! --|-----------------------------------------------------------|--
! | Delft University of Technology |
! | Faculty of Civil Engineering and Geosciences |
! | Environmental Fluid Mechanics Section |
! | P.O. Box 5048, 2600 GA Delft, The Netherlands |
! | |
! | Programmer: Marcel Zijlema |
! --|-----------------------------------------------------------|--
!
!
! SWAN (Simulating WAves Nearshore); a third generation wave model
! Copyright (C) 1993-2016 Delft University of Technology
!
! This program is free software; you can redistribute it and/or
! modify it under the terms of the GNU General Public License as
! published by the Free Software Foundation; either version 2 of
! the License, or (at your option) any later version.
!
! This program is distributed in the hope that it will be useful,
! but WITHOUT ANY WARRANTY; without even the implied warranty of
! MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
! GNU General Public License for more details.
!
! A copy of the GNU General Public License is available at
! http://www.gnu.org/copyleft/gpl.html#SEC3
! or by writing to the Free Software Foundation, Inc.,
! 59 Temple Place, Suite 330, Boston, MA 02111-1307 USA
!
!
! Authors
!
! 40.80: Marcel Zijlema
!
! Updates
!
! 40.80, July 2007: New subroutine
!
! Purpose
!
! Does some preparations before computation is started
!
! Modules used
!
use ocpcomm4
use swcomm3
use SwanGriddata
!PUN use SwanGridobjects
!
implicit none
!
! Argument variables
!
integer, dimension(nfaces), intent(out) :: cross ! contains sequence number of obstacles for each face
! where they crossing or zero if no crossing
!
! Local variables
!
integer, save :: ient = 0 ! number of entries in this subroutine
!PUN integer :: ivert ! loop counter over vertices
!
!PUN type(verttype), dimension(:), pointer :: vert ! datastructure for vertices with their attributes
!
! Structure
!
! Description of the pseudo code
!
! Source text
!
if (ltrace) call strace (ient,'SwanPrepComp')
!
!PUN ! point to vertex object
!PUN !
!PUN vert => gridobject%vert_grid
!PUN !
! deallocate arrays kvertc and kvertf (we don't use them anymore!)
!
if (allocated(kvertc)) deallocate(kvertc)
if (allocated(kvertf)) deallocate(kvertf)
!
!PUN ! ghost vertices are regarded as vertices with boundary condition
!PUN !
!PUN do ivert = 1, nverts
!PUN if ( vmark(ivert) == 999 ) vert(ivert)%atti(VBC) = 1
!PUN enddo
!PUN !
! setup a vertex list
!
call SwanVertlist
!
! find obstacles in computational grid, if present
!
if ( NUMOBS > 0 ) call SwanFindObstacles ( cross )
!
end subroutine SwanPrepComp