construct_grid Subroutine

public subroutine construct_grid(grid_type, n, g, rmin, rmax, l, mapping)

general function to generate grids grid_type= LGL, LG, or FGH mapping_type= linear or algebraic(reciprocal) This function does 1) set parameters of N and grid_type 2) allocate memory for r, rp, D1, D2, and so on 2) make computational grids 3) make real grids by means of a mapping function 4) construct D1 and D2

Arguments

Type IntentOptional Attributes Name
character(len=*) :: grid_type
integer, intent(in) :: n
type(grid), intent(out) :: g
real(kind=long), intent(in), optional :: rmin
real(kind=long), intent(in), optional :: rmax
real(kind=long), intent(in), optional :: l
character(len=*), intent(in), optional :: mapping