Skip to content

seg fault when executing on fusiont5 #3

@batchelordb

Description

@batchelordb

I’m trying to get genray going on fusiont5. Wael figured out how to build genray on fusiont5. He modified one of the archive makefiles. I will attach his version of the makefile and the diffs between his version and yours. I replicated what he did and got it to build with no errors reported. However when I run genray I get a seg fault. I will paste the log.genray output below. It seems to have trouble reading the namelists. However this is the genray.in file from one of your regression tests and worked in the past on edison and cori. Has the input interface changed recently?

DBB

MPI: qsize= 1

           =================================================
           GENRAY VERSION: genray_v10.12_180912
           =================================================

*********** code parameteres ************
number of rays in common/cone/ EC-cone: nraymax = 240
maximum number of EC cones: nconea = 20
maximum number of radial bins at the launching disk : n_mesh_disk_radial_bin_a = 5
nxeqda = 257
nyeqda = 257
nlimit = 101
nx4a=nxeqda+4 = 261
ny4a=nyeqda+4 = 261
nrya=ny4a = 261
nlim4=nlimit+4 = 105
It should be nrya=max(nxeqda,nyeqda)+4
nves = 62
number of points along each contour in the poloidal direction: nteta = 100
number of contours: npsi = 50
accuracy for the determination of contourspoins coordinates(zpsi,rpsi): epspsi = 1.0000000000000000E-004
nteta1=nteta+1 = 101
number of contours for the ploting of tokamakcross section: NL = 5
number of points along each contours in poloidal direction for plotting the tokamak cross section: NP = 100
for common/grill/
ngrilla = 10
nnkprmax = 60
nraymaxl = 240
,nthinmax = 160
for N_toroidal N_poloidal initial condition at i_n_poloidal=4
max{i=1,ngril1}nnktor(i): nnktormax = 60
max{i=1,ngril1}nnkpol(i): nnkpolmax = 60
max(nbulk): nbulka = 8
nbulkma=nbulka-1 = 7
maximal number of bin boundaries in the small radius direction
for power and CD radial profilesPower and current are tabulated at (NR-1) bin centers :NRA = 201
npsi4=npsi+4 = 54
nteta14=nteta1+4 = 105
nzy=nteta1+4 = 105
max number of points in arrays with plasma density, temperature, zeff. tpop, and vflow: ndensa = 201
ndens4a=ndensa+4 = 205
max number of output points along each ray: nrelta = 10000
n_relt_harma = 11
n_relt_harm1a = -5
n_relt_harm2a = 5
namber of points in rhoscat(nscat_n) and scatd(0:nscat_n): nscat_n = 10
to read distribution function f(iya,jx,lrz,ngen) written by CQL3D code
number of the points in the pitch angle mesh: iya = 120
iya1=iya+1 = 121
nrelta4=nrelta+4 = 10004
n_plot_dispa = 100
n_plot_disp_colda = 100
m_r_nperp_a = 20
m_i_nperp_a = 20
n_contour_plot_disp_a = 20
n_hot_roots_a = 4
n_power_switch_resonance_a = 3
gzonemax = 20
maximal numbes of wall points: n_wall_a = 200
maximal numbes of limiter points: n_limiter_a = 200
maximal number of limiters: max_limiters_a = 1
maximal number of wall points having the given poloidal angle theta: n_rho_wall_a = 10
to read normalized exponential density falls at poloidal mesh for edge_prof_nml,i
maximal number of mesh points: n_pol_edge_dens_a = 100
to create fine meshes of wall coordinates with additional points: n_wall_add_a = 10001
in subroutine read_all_namelists 1
in read_all_namelists after open file=genray.dat kode 0
======== USING genray.dat =========
in read_all_namelists after read genr kode= 0
print data from /genr/ obtained from input
prepare_genray_input after read tokamak kode= 0
Tokamak data
deltripl,nloop,i_ripple 0.0000000000000000 24 1
in read_write_genray_input.f after read(numercl)
&NUMERCL
IRKMETH= 2,
NDIM1= 6,
ISOLV= 1,
IDIF= 1,
NRELT= 1000,
PRMT1= 0.0000000000000000 ,
PRMT2= 999900.00000000000 ,
PRMT3= 1.0000000000000000E-004,
PRMT4= 1.0000000000000001E-005,
PRMT6= 3.7350000000000000E-003,
PRMT9= 1000000000000000.0 ,
ICORRECT= 1,
IOUT3D="enable ",
MAXSTEPS_RK= 10000,
I_OUTPUT= 1,
I_UH_SWITCH= 0,
UH_SWITCH= 1.5000000000000000 ,
PRMT6_UH_SWITCH= 1.0000000000000001E-005,
TOLL_HAMILT= 1.0000000000000000E-003,
I_POWER_SWITCH_RESONANCE= 0,
PRMT6_POWER_SWITCH_RESONANCE= 2.0000000000000000E-003,
N_POWER_SWITCH_RESONANCE= 1,
Y_POWER_SWITCH_RESONANCE= 0.50000000000000000 , 20.0000000000000000 ,
DEL_Y_POWER_SWITCH_RESONANCE= 1.0000000000000000E-002,
I_RESONANCE_CURVE_INTEGRATION_METHOD= 4,
EPSI= 1.0000000000000001E-005,
EPS_DELTA_POW= 10000000000.000000 ,
/
in read_all prmt= 0.0000000000000000 999900.00000000000 1.0000000000000000E-004 1.0000000000000001E-005 0.0000000000000000 3.7350000000000000E-003 0.0000000000000000 0.0000000000000000 1000000000000000.0
&OUTPUT
IWCNTR= 0,
IWOPEN= 1,
IWJ= 1,
ITOOLS= 0,
I_PLOT_B= 0,
I_PLOT_D= 0,
N_PLOT_DISP= 0,
R_PLOT_DISP= 100
0.0000000000000000 ,
ID_PLOT_DISP= 4,
Z_PLOT_DISP= 1000.0000000000000000 ,
N_PARALLEL_PLOT_DISP= 100
0.0000000000000000 ,
MAX_R_NPERP_PLOT_DISP= 1000.0000000000000000 ,
MIN_R_NPERP_PLOT_DISP= 100
0.0000000000000000 ,
MAX_I_NPERP_PLOT_DISP= 1000.0000000000000000 ,
MIN_I_NPERP_PLOT_DISP= 100
0.0000000000000000 ,
S_POLOID_PLOT_DISP= 1000.0000000000000000 ,
POINT_PLOT_DISP="poloidl_dist",
I_PLOT_DISP_COLD= 0,
N_PLOT_DISP_COLD= 0,
S_POLOID_PLOT_DISP_COLD= 100
0.0000000000000000 ,
R_PLOT_DISP_COLD= 1000.0000000000000000 ,
POINT_PLOT_DISP_COLD="poloidl_dist",
I_PLOT_WAVE_NORMAL_COLD= 0,
NUMBER_MAP_POINTS_REAL_NPERP= 10,
NUMBER_MAP_POINTS_IMAGE_NPERP= 10,
RATIO_MIN_R_NPERP= 0.50000000000000000 ,
RATIO_MAX_R_NPERP= 1.5000000000000000 ,
RATIO_MIN_I_NPERP= 0.0000000000000000 ,
RATIO_MAX_I_NPERP= 2.5000000000000000 ,
N_CONTOUR_PLOT_DISP= 20,
R_FREQ= 1.4900000000000000 ,
Z_FREQ= 0.0000000000000000 ,
ALPHA_FREQ= 180.00000000000000 ,
BETA_FREQ= 90.000000000000000 ,
DIST_FREQ= 1.2800000000000000 ,
MAX_PLOT_FREQ= 200.00000000000000 ,
NSTEPS_FREQ= 780,
N_EC_HARMONICS_FREQ= 6,
NPAR_FREQ= 0.0000000000000000 ,
/
nbulk1= 1
prepare_genray_input: partner=disabled
read_all_namelists (dens reading) kode= 0
at temperature reading kode= 0
at tpop reading kode= 0
at vflow reading kode= 0
&ZEFTAB
ZEFF1= 1.7100000000000000 , 1.7110000000000001 , 1.7100000000000000 , 1.7070000000000001 , 1.7010000000000001 ,
1.6940000000000000 , 3
1.6899999999999999 , 1.6850000000000001 , 1.6790000000000000 , 1.6739999999999999 , 1.6699999999999999 ,
1.6639999999999999 , 1.6599999999999999 , 1.6590000000000000 , 1.6599999999999999 , 1.6610000000000000 , 1.6470000000000000 ,
1.6100000000000001 , 1.5600000000000001 , 180*1.0000000000000000 ,
/
after read zeff1 1.7100000000000000 1.7110000000000001 1.7100000000000000 1.7070000000000001 1.7010000000000001 1.6940000000000000 1.6899999999999999 1.6899999999999999 1.6899999999999999 1.6850000000000001 1.6790000000000000 1.6739999999999999 1.6699999999999999 1.6639999999999999 1.6599999999999999 1.6590000000000000 1.6599999999999999 1.6610000000000000 1.6470000000000000 1.6100000000000001 1.5600000000000001 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000
at zeff reading kode= 0
in read_all_namelists after
read_uniform_column_profile
or read_nonuniform_line_profile
nj_tab_zeff1 201 201 201 201 201 201 201 201
before read_diskf
after read_diskf
before emission
after emission


check_read has found the negative value of iostat
The end of input in the given namelist emission was detected.
It can be that namelist emission is absent in genray.in file
Check data in emission in the file genray.in


before ox
after ox


check_read has found the negative value of iostat
The end of input in the given namelist ox was detected.
It can be that namelist ox is absent in genray.in file
Check data in ox in the file genray.in


before read adj_nml
after read adj_nml


check_read has found the negative value of iostat
The end of input in the given namelist adj_nml was detected.
It can be that namelist adj_nml is absent in genray.in file
Check data in adj_nml in the file genray.in


in read_write_genray_input.f i_adj= 0
Resetting i_adj to 0 for ieffic= 3
before edge_prof_nml


check_read has found the negative value of iostat
The end of input in the given namelist edge_prof_nml was detected.
It can be that namelist edge_prof_nml is absent in genray.in file
Check data in edge_prof_nml in the file genray.in


after edge_prof_nml
&EDGE_PROF_NML
I_EDGE_DENS_ANAL= 0,
I_EDGE_DENS_RZ_MESH= 0,
N_POL_EDGE_DENS= 100,
SIGMEDGN= 2.0000000000000000E-002,
SIGMEDGT= 2.0000000000000000E-002,
THETA_POL_EDGE_DENS_AR_DEGREE= 0.0000000000000000 , 3.6363636363636362 , 7.2727272727272725 , 10.909090909090908 , 14.545454545454545 ,
18.181818181818180 , 21.818181818181817 , 25.454545454545453 , 29.090909090909090 , 32.727272727272727 , 36.363636363636360 ,
40.000000000000000 , 43.636363636363633 , 47.272727272727273 , 50.909090909090907 , 54.545454545454547 , 58.181818181818180 ,
61.818181818181813 , 65.454545454545453 , 69.090909090909093 , 72.727272727272720 , 76.363636363636360 , 80.000000000000000 ,
83.636363636363640 , 87.272727272727266 , 90.909090909090907 , 94.545454545454547 , 98.181818181818173 , 101.81818181818181 ,
105.45454545454545 , 109.09090909090909 , 112.72727272727272 , 116.36363636363636 , 120.00000000000000 , 123.63636363636363 ,
127.27272727272727 , 130.90909090909091 , 134.54545454545453 , 138.18181818181819 , 141.81818181818181 , 145.45454545454544 ,
149.09090909090909 , 152.72727272727272 , 156.36363636363635 , 160.00000000000000 , 163.63636363636363 , 167.27272727272728 ,
170.90909090909091 , 174.54545454545453 , 178.18181818181819 , 181.81818181818181 , 185.45454545454544 , 189.09090909090909 ,
192.72727272727272 , 196.36363636363635 , 200.00000000000000 , 203.63636363636363 , 207.27272727272725 , 210.90909090909091 ,
214.54545454545453 , 218.18181818181819 , 221.81818181818181 , 225.45454545454544 , 229.09090909090909 , 232.72727272727272 ,
236.36363636363635 , 240.00000000000000 , 243.63636363636363 , 247.27272727272725 , 250.90909090909091 , 254.54545454545453 ,
258.18181818181819 , 261.81818181818181 , 265.45454545454544 , 269.09090909090907 , 272.72727272727269 , 276.36363636363637 ,
280.00000000000000 , 283.63636363636363 , 287.27272727272725 , 290.90909090909088 , 294.54545454545456 , 298.18181818181819 ,
301.81818181818181 , 305.45454545454544 , 309.09090909090907 , 312.72727272727269 , 316.36363636363637 , 320.00000000000000 ,
323.63636363636363 , 327.27272727272725 , 330.90909090909088 , 334.54545454545456 , 338.18181818181819 , 341.81818181818181 ,
345.45454545454544 , 349.09090909090907 , 352.72727272727269 , 356.36363636363637 , 360.00000000000000 ,
SIGMEDGN_AR= 100*2.0000000000000000E-002 ,
DENS_MIN_EDGE= 9.9999999999999995E-007,
TEMP_MIN_EDGE= 1.0000000000000000E-003,
THETA_POL_EDGE_1_DEGREE= 90.000000000000000 ,
THETA_POL_EDGE_2_DEGREE= 270.00000000000000 ,
SIGMA_THETA_POL_EDGE_1_DEGREE= 90.000000000000000 ,
SIGMA_THETA_POL_EDGE_2_DEGREE= 90.000000000000000 ,
SIGMA_EDGEN_0= 2.0000000000000000E-002,
SIGMA_EDGEN_1= 2.0000000000000000E-002,
SIGMA_EDGEN_2= 2.0000000000000000E-002,
SIGMA_WALL_N= 1.0000000000000000E-003,
SIGMA_LIM_TOROIDAL_DEGREE= 1.0000000000000000E-003,
NXEQD_ADD= 1000,
NYEQD_ADD= 1000,
/

Program received signal SIGSEGV: Segmentation fault - invalid memory reference.

Backtrace for this error:
#0 0x7fbd568a12da in ???
#1 0x7fbd568a0503 in ???
#2 0x7fbd55f1df1f in ???
#3 0x55796da84fbf in ???
#4 0x55796da8660e in ???
#5 0x7fbd55f00b96 in ???
#6 0x55796d7f52e9 in ???
#7 0xffffffffffffffff in ???

diff_make_files_wael_ComX.txt
makefile_mpi.gfortran64.txt

Metadata

Metadata

Assignees

No one assigned

    Labels

    No labels
    No labels

    Projects

    No projects

    Milestone

    No milestone

    Relationships

    None yet

    Development

    No branches or pull requests

    Issue actions