I'm stumbling trying to use the new CREATE_MIDPOINT grid option in GYRE 3.0. An example inlist an model is attached, and here's the terminal output:
Code: Select all
gyre_nad [3.0]
==============
Compiler : GCC version 4.9.0 20131223 (experimental)
Compiler options : -I /home/astro84/wball/mesa/sdk/include -march=sandybridge -mmmx -mno-3dnow -msse -msse2 -msse3 -mssse3 -mno-sse4a -mcx16 -msahf -mno-movbe -maes -mpclmul -mpopcnt -mno-abm -mno-lwp -mno-fma -mno-fma4 -mno-xop -mno-bmi -mno-tbm -mno-avx2 -msse4.2 -msse4.1 -mno-lzcnt -mno-rtm -mno-hle -mno-rdrnd -mno-f16c -mno-fsgsbase -mno-rdseed -mno-prfchw -mno-adx -mfxsr -mxsave -mxsaveopt -mno-avx512f -mno-avx512er -mno-avx512cd -mno-avx512pf --param l1-cache-size=32 --param l1-cache-line-size=64 --param l2-cache-size=10240 -mtune=sandybridge -mno-avx -mno-bmi2 -O2 -std=f2008 -ffpe-summary=none -finit-real=snan -ffpe-trap=invalid,overflow,zero -fbacktrace -fopenmp
OpenMP Threads : 4
MPI Processors : 1
Initialization
==============
Reading from MESA file best.gyre
Detected new-variant file
Adding central point
Building omega grid
ABORT at line 195 <gyre_grid:grid_range>:
Invalid op_type (the first op_type must be CREATE_*)
STOP Program aborted
Code: Select all
subroutine grid_range (gp, ml, op, x_in, x_i, x_o)
type(gridpar_t), intent(in) :: gp(:)
class(model_t), intent(in) :: ml
type(oscpar_t), intent(in) :: op
real(WP), allocatable, intent(in) :: x_in(:)
real(WP), intent(out) :: x_i
real(WP), intent(out) :: x_o
if(.NOT. (SIZE(gp) >= 1)) then
write(UNIT=ERROR_UNIT, FMT=*) 'ASSERT ''SIZE(gp) >= 1'' failed at line 173 <gyre_grid:grid_range>:'
write(UNIT=ERROR_UNIT, FMT=*) 'Empty gridpars'
stop
endif
! Determine the range spanned by the grid
select case (gp(1)%op_type)
case ('CREATE_CLONE')
if(.NOT. (ALLOCATED(x_in))) then
write(UNIT=ERROR_UNIT, FMT=*) 'ASSERT ''ALLOCATED(x_in)'' failed at line 179 <gyre_grid:grid_range>:'
write(UNIT=ERROR_UNIT, FMT=*) 'No grid to clone'
stop
endif
x_i = x_in(1)
x_o = x_in(SIZE(x_in))
case ('CREATE_UNIFORM')
x_i = 0._WP
x_o = 1._WP
case ('CREATE_GEOM')
x_i = 0._WP
x_o = 1._WP
case ('CREATE_LOG')
x_i = 0._WP
x_o = 1._WP
case ('CREATE_FROM_FILE')
x_i = 0._WP
x_o = 1._WP
case default
write(UNIT=ERROR_UNIT, FMT=*) 'ABORT at line 195 <gyre_grid:grid_range>:'
write(UNIT=ERROR_UNIT, FMT=*) 'Invalid op_type (the first op_type must be CREATE_*)'
stop 'Program aborted'
end select
! Finish
return
end subroutine grid_range
Code: Select all
case ('CREATE_MIDPOINT')
x_i = 0._WP
x_o = 1._WP
Scratch the last comment. Adding the extra case compiles but throws this error instead:
Code: Select all
gyre_nad [3.0]
==============
Compiler : GCC version 4.9.0 20131223 (experimental)
Compiler options : -I /home/astro84/wball/mesa/sdk/include -march=sandybridge -mmmx -mno-3dnow -msse -msse2 -msse3 -mssse3 -mno-sse4a -mcx16 -msahf -mno-movbe -maes -mpclmul -mpopcnt -mno-abm -mno-lwp -mno-fma -mno-fma4 -mno-xop -mno-bmi -mno-tbm -mno-avx2 -msse4.2 -msse4.1 -mno-lzcnt -mno-rtm -mno-hle -mno-rdrnd -mno-f16c -mno-fsgsbase -mno-rdseed -mno-prfchw -mno-adx -mfxsr -mxsave -mxsaveopt -mno-avx512f -mno-avx512er -mno-avx512cd -mno-avx512pf --param l1-cache-size=32 --param l1-cache-line-size=64 --param l2-cache-size=10240 -mtune=sandybridge -mno-avx -mno-bmi2 -O2 -std=f2008 -ffpe-summary=none -finit-real=snan -ffpe-trap=invalid,overflow,zero -fbacktrace -fopenmp
OpenMP Threads : 4
MPI Processors : 1
Initialization
==============
Reading from MESA file best.gyre
Detected new-variant file
Adding central point
Building omega grid
omega points : 500
omega range : 0.5001159574944368E+01 -> 0.4501043617449931E+02
ASSERT 'i > 0 .AND. i < this%n' failed at line 414 <core_spline:interp_v_>:
Out-of-bounds interpolation
Warrick
