search

PrincetonUniversity / STELLOPT

This is the GitHub repository for STELLOPT, the state-of-the-art stellarator optimization code.
https://princetonuniversity.github.io/STELLOPT/
MIT License
69 stars 19 forks source link

FIDASIM output is incorrect #201

Closed lazersos closed 1 year ago

lazersos commented 1 year ago

The -fidasim flag produces output which cannot be read by FIDASIM.

Others may also be missing.

FIDASIM Version:

   ____ ____ ___   ___    ____ ____ __  ___
  / __//  _// _ \ / _ |  / __//  _//  |/  /
 / _/ _/ / / // // __ | _\ \ _/ / / /|_/ /
/_/  /___//____//_/ |_|/___//___//_/  /_/

Version: v3.0.0-dev-32-g63c15db8-dirty
lazersos commented 1 year ago

OK after switching to branch v2.0.0 and recreating my atomic_tables.h5 file I now get a bizzare crash, but I don't think this is a BEAMS3D issue.

/Users/lazerson/Sims_Work/FIDASIM/fidasim fidasim_W7X_KJM_n10_e30_i16_OPTEH2_slow_inputs.dat 1
   ____ ____ ___   ___    ____ ____ __  ___
  / __//  _// _ \ / _ |  / __//  _//  |/  /
 / _/ _/ / / // // __ | _\ \ _/ / / /|_/ /
/_/  /___//____//_/ |_|/___//___//_/  /_/

Version: v2.0.0-dirty

FIDASIM is released as open source code under the MIT Licence.
For more information visit http://d3denergetic.github.io/FIDASIM/

---- Shot settings ----
 Shot:  3271890
 Time:    0 [ms]
 Runid: none

---- Input files ----
 Tables file: /Users/lazerson/Sims_Work/FIDASIM/atomic_tables.h5
 Geometry file: /Users/lazerson/Sims_Work/W7X/W7X_NBI/OPTEMIST/FIDASIM/fidasim_W7X_KJM_n10_e30_i16_OPTEH2_birth_geometry.h5
 Equilibrium file: /Users/lazerson/Sims_Work/W7X/W7X_NBI/OPTEMIST/FIDASIM/fidasim_W7X_KJM_n10_e30_i16_OPTEH2_birth_equilibrium.h5
 Distribution file: /Users/lazerson/Sims_Work/W7X/W7X_NBI/OPTEMIST/FIDASIM/fidasim_W7X_KJM_n10_e30_i16_OPTEH2_slow_distribution.h5

---- OpenMP settings ----
 Number of threads:  1

---- Random Number Generator settings ----
 RNG Seed: 1751101603

---- Atomic tables settings ----
 Maximum n/m:  6
 Beam/Fast-ion mass:  1.007 [amu]
 Thermal/Bulk-ion mass:  1.007 [amu]
 Impurity mass: 10.810 [amu]

---- Interpolation grid settings ----
 Nr: 128
 Nz: 128
 Nphi:  32
 dA:  3.43 [cm^2]

zsh: segmentation fault  /Users/lazerson/Sims_Work/FIDASIM/fidasim  1
lazersos commented 1 year ago

OK ran in debug mode and learned a bit more:

   ____ ____ ___   ___    ____ ____ __  ___
  / __//  _// _ \ / _ |  / __//  _//  |/  /
 / _/ _/ / / // // __ | _\ \ _/ / / /|_/ /
/_/  /___//____//_/ |_|/___//___//_/  /_/

Version: v2.0.0-dirty

FIDASIM is released as open source code under the MIT Licence.
For more information visit http://d3denergetic.github.io/FIDASIM/

########################### ATTENTION ###########################
# Running in debug mode. All optimizations have been turned off #
#################################################################

---- Shot settings ----
 Shot:  3271890
 Time:    0 [ms]
 Runid: none

---- Input files ----
 Tables file: /Users/lazerson/Sims_Work/FIDASIM/atomic_tables.h5
 Geometry file: /Users/lazerson/Sims_Work/W7X/W7X_NBI/OPTEMIST/FIDASIM/fidasim_W7X_KJM_n10_e30_i16_OPTEH2_birth_geometry.h5
 Equilibrium file: /Users/lazerson/Sims_Work/W7X/W7X_NBI/OPTEMIST/FIDASIM/fidasim_W7X_KJM_n10_e30_i16_OPTEH2_birth_equilibrium.h5
 Distribution file: /Users/lazerson/Sims_Work/W7X/W7X_NBI/OPTEMIST/FIDASIM/fidasim_W7X_KJM_n10_e30_i16_OPTEH2_slow_distribution.h5

---- OpenMP settings ----
 Number of threads:  1

---- Random Number Generator settings ----
 RNG Seed: 2098470796

---- Atomic tables settings ----
 Maximum n/m:  6
 Beam/Fast-ion mass:  1.007 [amu]
 Thermal/Bulk-ion mass:  1.007 [amu]
 Impurity mass: 10.810 [amu]

---- Interpolation grid settings ----
 Nr: 128
 Nz: 128
 Nphi:  32
 dA:  3.43 [cm^2]

At line 3042 of file fidasim.f90
Fortran runtime warning: An array temporary was created
At line 3043 of file fidasim.f90
Fortran runtime warning: An array temporary was created
At line 3044 of file fidasim.f90
Fortran runtime warning: An array temporary was created
At line 3045 of file fidasim.f90
Fortran runtime warning: An array temporary was created
At line 3046 of file fidasim.f90
Fortran runtime warning: An array temporary was created
At line 3047 of file fidasim.f90
Fortran runtime warning: An array temporary was created
At line 3048 of file fidasim.f90
Fortran runtime warning: An array temporary was created
At line 3123 of file fidasim.f90
Fortran runtime warning: An array temporary was created
At line 3124 of file fidasim.f90
Fortran runtime warning: An array temporary was created
At line 3125 of file fidasim.f90
Fortran runtime warning: An array temporary was created
At line 3126 of file fidasim.f90
Fortran runtime warning: An array temporary was created
At line 3127 of file fidasim.f90
Fortran runtime warning: An array temporary was created
At line 3128 of file fidasim.f90
Fortran runtime warning: An array temporary was created
At line 3132 of file fidasim.f90
Fortran runtime warning: An array temporary was created
At line 3133 of file fidasim.f90
Fortran runtime warning: An array temporary was created
At line 3133 of file fidasim.f90
Fortran runtime warning: An array temporary was created
At line 3133 of file fidasim.f90
Fortran runtime warning: An array temporary was created
At line 3134 of file fidasim.f90
Fortran runtime warning: An array temporary was created
At line 3135 of file fidasim.f90
Fortran runtime warning: An array temporary was created
At line 3135 of file fidasim.f90
Fortran runtime warning: An array temporary was created
At line 3135 of file fidasim.f90
Fortran runtime warning: An array temporary was created
At line 3136 of file fidasim.f90
Fortran runtime warning: An array temporary was created
At line 3137 of file fidasim.f90
Fortran runtime warning: An array temporary was created
At line 3137 of file fidasim.f90
Fortran runtime warning: An array temporary was created
At line 3137 of file fidasim.f90
Fortran runtime warning: An array temporary was created
---- Beam grid settings ----
 Nx: 170
 Ny: 150
 Nz:  50
 dV:  3.00 [cm^3]
 alpha: -2.07 [rad]
 beta:   0.00 [rad]
 gamma:  0.00 [rad]
 origin: [   0.00,   0.00,   0.00] [cm]
 Number of cells in plasma:  1275000

---- Neutral beam settings ----
 Beam: OPTEMIST-Diag
 Power:    1.00 [MW]
 Voltage: 150.00 [keV]

---- Fast-ion distribution settings ----
HDF5-DIAG: Error detected in HDF5 (1.14.1-2) thread 0:
  #000: H5D.c line 402 in H5Dopen2(): unable to synchronously open dataset
    major: Dataset
    minor: Can't open object
  #001: H5D.c line 363 in H5D__open_api_common(): unable to open dataset
    major: Dataset
    minor: Can't open object
  #002: H5VLcallback.c line 1980 in H5VL_dataset_open(): dataset open failed
    major: Virtual Object Layer
    minor: Can't open object
  #003: H5VLcallback.c line 1947 in H5VL__dataset_open(): dataset open failed
    major: Virtual Object Layer
    minor: Can't open object
  #004: H5VLnative_dataset.c line 319 in H5VL__native_dataset_open(): unable to open dataset
    major: Dataset
    minor: Can't open object
  #005: H5Dint.c line 1423 in H5D__open_name(): not found
    major: Dataset
    minor: Object not found
  #006: H5Gloc.c line 438 in H5G_loc_find(): can't find object
    major: Symbol table
    minor: Object not found
  #007: H5Gtraverse.c line 836 in H5G_traverse(): internal path traversal failed
    major: Symbol table
    minor: Object not found
  #008: H5Gtraverse.c line 613 in H5G__traverse_real(): traversal operator failed
    major: Symbol table
    minor: Callback failed
  #009: H5Gloc.c line 395 in H5G__loc_find_cb(): object 'denf' doesn't exist
    major: Symbol table
    minor: Object not found
HDF5-DIAG: Error detected in HDF5 (1.14.1-2) thread 0:
  #000: H5D.c line 402 in H5Dopen2(): unable to synchronously open dataset
    major: Dataset
    minor: Can't open object
  #001: H5D.c line 363 in H5D__open_api_common(): unable to open dataset
    major: Dataset
    minor: Can't open object
  #002: H5VLcallback.c line 1980 in H5VL_dataset_open(): dataset open failed
    major: Virtual Object Layer
    minor: Can't open object
  #003: H5VLcallback.c line 1947 in H5VL__dataset_open(): dataset open failed
    major: Virtual Object Layer
    minor: Can't open object
  #004: H5VLnative_dataset.c line 319 in H5VL__native_dataset_open(): unable to open dataset
    major: Dataset
    minor: Can't open object
  #005: H5Dint.c line 1423 in H5D__open_name(): not found
    major: Dataset
    minor: Object not found
  #006: H5Gloc.c line 438 in H5G_loc_find(): can't find object
    major: Symbol table
    minor: Object not found
  #007: H5Gtraverse.c line 836 in H5G_traverse(): internal path traversal failed
    major: Symbol table
    minor: Object not found
  #008: H5Gtraverse.c line 613 in H5G__traverse_real(): traversal operator failed
    major: Symbol table
    minor: Callback failed
  #009: H5Gloc.c line 395 in H5G__loc_find_cb(): object 'f' doesn't exist
    major: Symbol table
    minor: Object not found
 Distribution type: Non-axisymmetric Fast-ion Density Function F(energy,pitch,R,Z,Phi)
 Nenergy =  32
 Npitch  =  64
 Nr  = 128
 Nz  = 128
 Nphi  =  32
 Energy range = [ 5.63,354.38]
 Pitch  range = [-0.98, 0.98]
 R  range = [425.00,648.24]
 Z  range = [-125.00,123.05]
 Phi  range = [ 0.87, 1.21]
 Ntotal =  0.000E+00

---- FIDA/BES settings ----
At line 2813 of file fidasim.f90
Fortran runtime warning: An array temporary was created
At line 2814 of file fidasim.f90
Fortran runtime warning: An array temporary was created
 FIDA/BES System: W7X OPTEMIST FIDA SC
 Number of channels:   512

 nbi:     14:59:19 --- elapsed: 0:00:15
     # of markers:    5000000
     # of birth markers: 3 x   5000000
MC_NBI: A beam neutral has started inside the plasma.
Move the beam grid closer to the source to fix
lazersos commented 1 year ago

So v2.0.0 works OK with the files so I'm going to close this issue.