1 | include(dom.inc)
2 | ! ===============================================================
3 | ! Copyright (c) CERFACS (all rights reserved)
4 | ! ===============================================================
5 |
6 |
7 | ! ================================================================
8 | ! ================================================================
9 | module avbp_coor
10 |
11 | implicit none
12 |
13 | character*80 coorfile
14 |
15 | DOM_REAL,dimension(:,:),allocatable:: coor
16 |
17 | DOM_INT coor_neq, &
18 | & coor_ndim, &
19 | & coor_nnode, &
20 | & coor_nblock, &
21 | & coor_nnbegin, &
22 | & coor_nnlength, &
23 | & coor_mynodes
24 |
25 | logical :: coor_itexists = .false.
26 | logical :: coor_show = .true.
27 |
28 | character*1,dimension(1:3),parameter :: &
29 | & coor_vname = (/'x','y','z'/)
30 |
31 | contains
32 |
33 | ! ================================================================
34 | subroutine destroycoormemory
35 |
36 | if(allocated(coor)) then
37 | deallocate(coor)
38 | endif
39 |
40 | end subroutine destroycoormemory
41 | ! ================================================================
42 |
43 |
44 | ! ================================================================
45 | subroutine getcoorinfo
46 |
47 | implicit none
48 | DOM_INT n
49 |
50 | inquire(file=coorfile,exist=coor_itexists)
51 |
52 | if(coor_itexists) then
53 | open( unit=1, file=coorfile, form='unformatted' )
54 | read(1) coor_ndim, coor_neq
55 | read(1) coor_nblock, coor_nnode
56 | close(1)
57 | else
58 | n = len_trim(coorfile)
59 | write(*,*) ' could not find the file ',coorfile(1:n)
60 | endif
61 |
62 | end subroutine getcoorinfo
63 | ! ================================================================
64 |
65 |
66 | ! ================================================================
67 | subroutine readcoor
68 |
69 | implicit none
70 | DOM_INT k,n
71 |
72 | n = len_trim(coorfile)
73 | if (coor_show) then
74 | write(*,*) ' >>>> Reading mesh file ',coorfile(1:n)
75 | end if
76 |
77 | inquire(file=coorfile,exist=coor_itexists)
78 |
79 |
80 | if(coor_itexists) then
81 | open( unit=1, file=coorfile, form='unformatted' )
82 |
83 | read(1) coor_ndim, coor_neq
84 | read(1) coor_nblock, coor_nnode
85 |
86 | if(allocated(coor)) deallocate(coor)
87 |
88 | allocate(coor(1:coor_ndim,1:coor_nnode))
89 |
90 | read(1) coor_nnbegin,coor_nnlength,coor_mynodes
91 | read(1) ( (coor(k,n), n=1,coor_nnode),k=1,coor_ndim )
92 |
93 | close(1)
94 | else
95 | n = len_trim(coorfile)
96 | write(*,*) ' could not find the file ',coorfile(1:n)
97 | endif
98 |
99 | end subroutine readcoor
100 | ! ================================================================
101 |
102 |
103 | ! ================================================================
104 | subroutine writecoor
105 |
106 | implicit none
107 | DOM_INT k,n
108 |
109 | n = len_trim(coorfile)
110 | write(*,*) ' >>>> Writing mesh file ',coorfile(1:n)
111 |
112 | open( unit=1, file=coorfile, form='unformatted' )
113 |
114 | write(1) coor_ndim, coor_neq
115 | write(1) coor_nblock, coor_nnode
116 |
117 | write(1) coor_nnbegin,coor_nnlength,coor_mynodes
118 | write(1) ( (coor(k,n), n=1,coor_nnode),k=1,coor_ndim )
119 |
120 | close(1)
121 |
122 | end subroutine writecoor
123 | ! ================================================================
124 |
125 |
126 | end module avbp_coor
127 | ! ================================================================
128 | ! ================================================================
129 |