Homec4science
Diffusion SPClibs (feature/refactoring) Browse

spclibs/examples/mesh.f90
03540012fedefeature/refactoring

Last Change
Raw File

mesh.f90
View Options

!>
!> @file mesh.f90
!>
!> @brief
!>
!> @copyright
!> Copyright (©) 2021 EPFL (Ecole Polytechnique Fédérale de Lausanne)
!> SPC (Swiss Plasma Center)
!>
!> spclibs is free software: you can redistribute it and/or modify it under
!> the terms of the GNU Lesser General Public License as published by the Free
!> Software Foundation, either version 3 of the License, or (at your option)
!> any later version.
!>
!> spclibs 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.
!>
!> You should have received a copy of the GNU Lesser General Public License
!> along with this program. If not, see <https://www.gnu.org/licenses/>.
!>
!> @authors
!> (in alphabetical order)
!> @author Trach-Minh Tran <trach-minh.tran@epfl.ch>
!>
PROGRAM main
!
! Equidistant and non-equidistant mesh
!
USE bsplines
IMPLICIT NONE
INTEGER :: i, nx
DOUBLE PRECISION :: coefs(5), dev, dx
DOUBLE PRECISION, DIMENSION(:), ALLOCATABLE :: xgrid
!
INTERFACE
SUBROUTINE meshdist(coefs, x, nx)
DOUBLE PRECISION, INTENT(in) :: coefs(5)
INTEGER, INTENT(iN) :: nx
DOUBLE PRECISION, INTENT(inout) :: x(0:nx)
END SUBROUTINE meshdist
END INTERFACE
!
NAMELIST /newrun/ nx, coefs
!
nx = 8 ! Number oh intevals in x
coefs(1:5) = (/0.2d0, 1.d0, 0.d0, 0.d0, 1.d0/) ! Mesh point distribution function
! see function FDIST in MESHDIST
!
READ(*,newrun)
WRITE(*,newrun)
!
ALLOCATE(xgrid(0:nx))
xgrid(0) = 0.0d0
xgrid(nx) = 1.0d0
CALL meshdist(coefs, xgrid, nx)
PRINT*, 'Mesh is equidistant?', is_equid(xgrid, dev)
PRINT*, 'dev =', dev
!
dx = 1.0d0/REAL(nx,8)
xgrid = (/ (i*dx, i=0,nx) /)
PRINT*, 'Mesh is equidistant?', is_equid(xgrid, dev)
PRINT*, 'dev =', dev
!
DEALLOCATE(xgrid)
END PROGRAM main
c4science · Help