Neko  0.9.99
A portable framework for high-order spectral element flow simulations
mesh_field.f90
Go to the documentation of this file.
1 ! Copyright (c) 2019-2023, The Neko Authors
2 ! All rights reserved.
3 !
4 ! Redistribution and use in source and binary forms, with or without
5 ! modification, are permitted provided that the following conditions
6 ! are met:
7 !
8 ! * Redistributions of source code must retain the above copyright
9 ! notice, this list of conditions and the following disclaimer.
10 !
11 ! * Redistributions in binary form must reproduce the above
12 ! copyright notice, this list of conditions and the following
13 ! disclaimer in the documentation and/or other materials provided
14 ! with the distribution.
15 !
16 ! * Neither the name of the authors nor the names of its
17 ! contributors may be used to endorse or promote products derived
18 ! from this software without specific prior written permission.
19 !
20 ! THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS
21 ! "AS IS" AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT
22 ! LIMITED TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS
23 ! FOR A PARTICULAR PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL THE
24 ! COPYRIGHT OWNER OR CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT,
25 ! INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING,
26 ! BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES;
27 ! LOSS OF USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION) HOWEVER
28 ! CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT
29 ! LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN
30 ! ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE
31 ! POSSIBILITY OF SUCH DAMAGE.
32 !
35 module mesh_field
36  use mesh, only : mesh_t
37  implicit none
38  private
39 
41  type, public :: mesh_fld_t
42  integer, allocatable :: data(:)
43  type(mesh_t), pointer :: msh
44  character(len=80) :: name
45  end type mesh_fld_t
46 
48 
49 contains
50 
51  subroutine mesh_field_init(fld, msh, fld_name)
52  type(mesh_fld_t), intent(inout) :: fld
53  type(mesh_t), target, intent(in) :: msh
54  character(len=*), optional :: fld_name
55 
56  call mesh_field_free(fld)
57 
58  fld%msh => msh
59  if (.not. allocated(fld%data)) then
60  allocate(fld%data(msh%nelv))
61  end if
62 
63  if (present(fld_name)) then
64  fld%name = fld_name
65  else
66  fld%name = 'MeshField'
67  end if
68 
69  fld%data = 0
70  end subroutine mesh_field_init
71 
72  subroutine mesh_field_free(fld)
73  type(mesh_fld_t), intent(inout) :: fld
74 
75  if (allocated(fld%data)) then
76  deallocate(fld%data)
77  end if
78 
79  nullify(fld%msh)
80  end subroutine mesh_field_free
81 
82 end module mesh_field
Defines a mesh field.
Definition: mesh_field.f90:35
subroutine, public mesh_field_free(fld)
Definition: mesh_field.f90:73
subroutine, public mesh_field_init(fld, msh, fld_name)
Definition: mesh_field.f90:52
Defines a mesh.
Definition: mesh.f90:34