running build
running config_cc
INFO: unifing config_cc, config, build_clib, build_ext, build commands --compiler options
running config_fc
INFO: unifing config_fc, config, build_clib, build_ext, build commands --fcompiler options
running build_src
INFO: build_src
INFO: building extension "mstm" sources
INFO: f2py options: []
INFO: f2py:> /tmp/tmp8y23f3q4/src.linux-x86_64-3.10/mstmmodule.c
creating /tmp/tmp8y23f3q4/src.linux-x86_64-3.10
Reading f2cmap from '.f2py_f2cmap' ...
        Mapping "complex(kind= kind = 8 )" to "complex_double"
        Mapping "complex(kind=4)" to "complex_float"
        Warning: redefinition of {'complex':{'8':'complex_float'->'complex_double'}}
        Mapping "complex(kind=8)" to "complex_double"
        Mapping "complex(kind=dp)" to "complex_double"
        Mapping "complex(kind=sp)" to "complex_float"
Successfully applied user defined f2cmap changes
Reading fortran codes...
        Reading file './code/mstm-intrinsics.f90' (format:free)
rmbadname1: Replacing "char" with "char_bn".
rmbadname1: Replacing "char" with "char_bn".
rmbadname1: Replacing "char" with "char_bn".
        Reading file './code/mpidefs-serial.f90' (format:free)
rmbadname1: Replacing "type" with "type_bn".
rmbadname1: Replacing "size" with "size_bn".
rmbadname1: Replacing "rank" with "rank_bn".
        Reading file './code/mstm-v4.0.f90' (format:free)
rmbadname1: Replacing "const" with "const_bn".
analyzeline: intent(callback) f is already in argument list
analyzeline: intent(callback) qsub is already in argument list
rmbadname1: Replacing "const" with "const_bn".
rmbadname1: Replacing "const" with "const_bn".
linked_sphere_list <re.Match object; span=(0, 18), match='linked_sphere_list'>
host_list <re.Match object; span=(0, 9), match='host_list'>
translation_data <re.Match object; span=(0, 16), match='translation_data'>
surface_ref_data <re.Match object; span=(0, 16), match='surface_ref_data'>
pl_translation_data <re.Match object; span=(0, 19), match='pl_translation_data'>
rmbadname1: Replacing "rank" with "rank_bn".
rmbadname1: Replacing "rank" with "rank_bn".
rmbadname1: Replacing "rank" with "rank_bn".
rmbadname1: Replacing "rank" with "rank_bn".
node_data <re.Match object; span=(0, 9), match='node_data'>
linked_ilist <re.Match object; span=(0, 12), match='linked_ilist'>
rmbadname1: Replacing "index" with "index_bn".
coefficient_list <re.Match object; span=(0, 16), match='coefficient_list'>
rmbadname1: Replacing "rank" with "rank_bn".
rmbadname1: Replacing "rank" with "rank_bn".
rmbadname1: Replacing "rank" with "rank_bn".
rmbadname1: Replacing "rank" with "rank_bn".
rmbadname1: Replacing "rank" with "rank_bn".
rmbadname1: Replacing "rank" with "rank_bn".
rmbadname1: Replacing "rank" with "rank_bn".
rmbadname1: Replacing "const" with "const_bn".
rmbadname1: Replacing "rank" with "rank_bn".
rmbadname1: Replacing "rank" with "rank_bn".
rmbadname1: Replacing "rank" with "rank_bn".
rmbadname1: Replacing "rank" with "rank_bn".
grid_info <re.Match object; span=(0, 9), match='grid_info'>
cell_info <re.Match object; span=(0, 9), match='cell_info'>
linked_cell_list <re.Match object; span=(0, 16), match='linked_cell_list'>
vector_storage <re.Match object; span=(0, 14), match='vector_storage'>
linked_sphere_data <re.Match object; span=(0, 18), match='linked_sphere_data'>
rmbadname1: Replacing "rank" with "rank_bn".
rmbadname1: Replacing "rank" with "rank_bn".
rmbadname1: Replacing "rank" with "rank_bn".
rmbadname1: Replacing "rank" with "rank_bn".
coll_list <re.Match object; span=(0, 9), match='coll_list'>
l_list <re.Match object; span=(0, 6), match='l_list'>
rmbadname1: Replacing "index" with "index_bn".
c_list <re.Match object; span=(0, 6), match='c_list'>
rmbadname1: Replacing "rank" with "rank_bn".
rmbadname1: Replacing "rank" with "rank_bn".
rmbadname1: Replacing "rank" with "rank_bn".
rmbadname1: Replacing "rank" with "rank_bn".
rmbadname1: Replacing "rank" with "rank_bn".
rmbadname1: Replacing "rank" with "rank_bn".
rmbadname1: Replacing "rank" with "rank_bn".
rmbadname1: Replacing "rank" with "rank_bn".
rmbadname1: Replacing "rank" with "rank_bn".
rmbadname1: Replacing "rank" with "rank_bn".
rmbadname1: Replacing "rank" with "rank_bn".
rmbadname1: Replacing "rank" with "rank_bn".
rmbadname1: Replacing "rank" with "rank_bn".
Post-processing...
        Block: mstm
                        Block: intrinsics
                                Block: mytime
In: :mstm:./code/mstm-intrinsics.f90:intrinsics:mytime
get_parameters: got "eval() arg 1 must be a string, bytes or code object""character*1 loop_var_type(5)" is considered as "character loop_var_type(5,1)"; "intent(c)" is forced.
analyzevars: character array "character*256 loop_var_label(5)" is considered as "character loop_var_label(5,256)"; "intent(c)" is forced.
                                Block: variable_list_operation
                                Block: inputdata
analyzevars: character array "character*256 inputfiledata(*)" is considered as "character inputfiledata(*,256)"; "intent(c)" is forced.
                                Block: main_calling_program
                                Block: configuration_average_calling_program
                                Block: incidence_average_calling_program
                                Block: subtract_1_from_0
                                Block: surface_absorptance_calculation
                                Block: sample_incident_direction
                                Block: effective_extinction_coefficient_ratio
                                Block: read_sphere_data_input_file
                                Block: generate_random_configuration
                                Block: scattering_matrix_calculation
                                Block: output_header
                                Block: print_run_variables
                                Block: print_calculation_results
analyzevars: character array "character*2 smlabel(16)" is considered as "character smlabel(16,2)"; "intent(c)" is forced.
                                Block: scaled_scat_mat
                                Block: scat_mat_to_phase_mat
                                Block: checkpositions
                                Block: set_string_to_int_variable
                                Block: set_string_to_real_variable
                                Block: set_string_to_real_array_variable
                                Block: set_string_to_cmplx_variable
                                Block: set_string_to_logical_variable
                                Block: set_string_to_logical_array_variable
                        Block: main
analyzevars: character array "character*256 inputfiledata(:)" is considered as "character inputfiledata(:,256)"; "intent(c)""cb_f_in_qng__user__routines"
append_needs: unknown need 'int'
append_needs: unknown need 'double'
      def f(centr,fcentr,[n]): return
    Building module "mstm"...
                Constructing COMMON block support for "qkernelcommon"...
                  nodr,integrationmodel,x,y,z,w,k0y,kz,ri
                Constructing F90 module support for "intrinsics"...
                Creating wrapper for Fortran function "mytime"("mytime")...
            Constructing wrapper function "intrinsics.mytime"...
              mytime = mytime()
                Creating wrapper for Fortran function "mstm_nargs"("mstm_nargs")...
            Constructing wrapper function "intrinsics.mstm_nargs"...
              mstm_nargs = mstm_nargs()
            Constructing wrapper function "intrinsics.mstm_getarg"...
              mstm_getarg(char_bn)
            Constructing wrapper function "intrinsics.mstm_fseek"...
              mstm_fseek(unit,position,whence,ierr)
                Constructing F90 module support for "mpidefs"...
                  Variables: mpi_comm_world mstm_mpi_comm_world mstm_mpi_sum mstm_mpi_max mstm_mpi_min mpi_comm_null mstm_global_rank
                Creating wrapper for Fortran function "mstm_mpi_wtime"("mstm_mpi_wtime")...
            Constructing wrapper function "mpidefs.mstm_mpi_wtime"...
              mstm_mpi_wtime = mstm_mpi_wtime()
            Constructing wrapper function "mpidefs.mstm_mpi"...
getarrdims:warning: assumed shape array, using 0 instead of '*'
getarrdims:warning: assumed shape array, using 0 instead of '*'
getarrdims:warning: assumed shape array, using 0 instead of '*'
getarrdims:warning: assumed shape array, using 0 instead of '*'
getarrdims:warning: assumed shape array, using 0 instead of '*'
getarrdims:warning: assumed shape array, using 0 instead of '*'
getarrdims:warning: assumed shape array, using 0 instead of '*'
getarrdims:warning: assumed shape array, using 0 instead of '*'
getarrdims:warning: assumed shape array, using 0 instead of '*'
getarrdims:warning: assumed shape array, using 0 instead of '*'
getarrdims:warning: assumed shape array, using 0 instead of '*'
              mstm_mpi(mpi_command,[mpi_recv_buf_i,mpi_recv_buf_r,mpi_recv_buf_c,mpi_recv_buf_dp,mpi_recv_buf_dc,mpi_send_buf_i,mpi_send_buf_r,mpi_send_buf_c,mpi_send_buf_dp,mpi_send_buf_dc,mpi_number,mpi_comm,mpi_group,mpi_rank,mpi_size,mpi_new_comm,mpi_new_group,mpi_new_group_list,mpi_operation,mpi_color,mpi_key,mpi_tag,mpi_flag])
                Constructing F90 module support for "numconstants"...
                  Variables: light_up print_intermediate_results global_rank monen pi bcof fnr vwh_coef vcc_const fnm1_const fn_const fnp1_const tran_coef
getarrdims:warning: assumed shape array, using 0 instead of ':'
getarrdims:warning: assumed shape array, using 0 instead of ':'
getarrdims:warning: assumed shape array, using 0 instead of ':'
getarrdims:warning: assumed shape array, using 0 instead of ':'
getarrdims:warning: assumed shape array, using 0 instead of ':'
getarrdims:warning: assumed shape array, using 0 instead of ':'
getarrdims:warning: assumed shape array, using 0 instead of ':'
getarrdims:warning: assumed shape array, using 0 instead of ':'
getarrdims:warning: assumed shape array, using 0 instead of ':'
getarrdims:warning: assumed shape array, using 0 instead of ':'
getarrdims:warning: assumed shape array, using 0 instead of ':'
getarrdims:warning: assumed shape array, using 0 instead of ':'
getarrdims:warning: assumed shape array, using 0 instead of ':'
getarrdims:warning: assumed shape array, using 0 instead of ':'
getarrdims:warning: assumed shape array, using 0 instead of ':'
getarrdims:warning: assumed shape array, using 0 instead of ':'
getarrdims:warning: assumed shape array, using 0 instead of ':'
getarrdims:warning: assumed shape array, using 0 instead of ':'
getarrdims:warning: assumed shape array, using 0 instead of ':'
getarrdims:warning: assumed shape array, using 0 instead of ':'
            Constructing wrapper function "numconstants.init"...
              init(notd)
                Constructing F90 module support for "specialfuncs"...
            Constructing wrapper function "specialfuncs.timewrite"...
              timewrite(iunit,char1,time,[line_break])
            Constructing wrapper function "specialfuncs.ricbessel"...
              ricbessel(ds,eps,nmax,psi,[n])
            Constructing wrapper function "specialfuncs.richankel"...
              richankel(ds,xi,[n])
            Constructing wrapper function "specialfuncs.cricbessel"...
              cricbessel(ds,psi,[n])
            Constructing wrapper function "specialfuncs.crichankel"...
              crichankel(ds,xi,[n])
            Constructing wrapper function "specialfuncs.cspherebessel"...
              cspherebessel(z,csj,csy,[n])
            Constructing wrapper function "specialfuncs.ch12n"...
              ch12n(z,nm,chf1,[n])
            Constructing wrapper function "specialfuncs.ciknb"...
              ciknb(z,nm,cbi,cdi,cbk,cdk,[n])
            Constructing wrapper function "specialfuncs.bessel_integer_complex"...
              bessel_integer_complex(z,nmax,b,[n])
            Constructing wrapper function "specialfuncs.cjynb"...
              cjynb(z,nm,cbj,cdj,cby,cdy,[n])
                Creating wrapper for Fortran function "msta1"("msta1")...
            Constructing wrapper function "specialfuncs.msta1"...
              msta1 = msta1(x,mp)
                Creating wrapper for Fortran function "msta2"("msta2")...
            Constructing wrapper function "specialfuncs.msta2"...
              msta2 = msta2(x,n,mp)
                Creating wrapper for Fortran function "envj"("envj")...
            Constructing wrapper function "specialfuncs.envj"...
              envj = envj(n,x)
            Constructing wrapper function "specialfuncs.vcfunc"...
              vcfunc(m,n,k,l,vcn)
            Constructing wrapper function "specialfuncs.vcfuncuprec"...
              vcfuncuprec(m,n,k,l,wmax,vcn)
            Constructing wrapper function "specialfuncs.normalizedlegendre"...
              normalizedlegendre(cbe,dc,[mmax,nmax])
            Constructing wrapper function "specialfuncs.rotcoef"...
              rotcoef(cbe,nmax,dc,[kmax])
            Constructing wrapper function "specialfuncs.crotcoef"...
              crotcoef(cbe,nmax,dc,[kmax,sin_beta])
            Constructing wrapper function "specialfuncs.complexpivec"...
              complexpivec(cb,nodr,pivec,[icon,lr_model,azimuth_angle,index_model])
            Constructing wrapper function "specialfuncs.taufunc"...
              taufunc(cb,tau,[nmax])
            Constructing wrapper function "specialfuncs.rotvec"...
              rotvec(alpha,beta,gamma,mmax,amn,idir,[nmax])
            Constructing wrapper function "specialfuncs.genplanewavecoef"...
              genplanewavecoef(alpha,cb,pmnp0,[nodr,lr_tran])
            Constructing wrapper function "specialfuncs.gaussianbeamcoef"...
              gaussianbeamcoef(alpha,cbeta,cbeam,pmnp0,[nodr,lr_tran])
            Constructing wrapper function "specialfuncs.axialtrancoefrecurrence"...
              axialtrancoefrecurrence(itype,r,ri,nmax,lmax,ac,[ndim])
            Constructing wrapper function "specialfuncs.axialtrancoefinit"...
              axialtrancoefinit(nmax)
            Constructing wrapper function "specialfuncs.gentrancoefconstants"...
              gentrancoefconstants(nodrmax)
            Constructing wrapper function "specialfuncs.gentranmatrix"...
              gentranmatrix(nodr_s,nodr_t,translation_vector,ac_matrix,[refractive_index,vswf_type,mode_s,mode_t,index_model])
            Constructing wrapper function "specialfuncs.tranordertest"...
              tranordertest(r,ri,lmax,eps,nmax)
                Creating wrapper for Fortran function "atcadd"("atcadd")...
            Constructing wrapper function "specialfuncs.atcadd"...
              atcadd = atcadd(m,n,ntot)
                Creating wrapper for Fortran function "atcdim"("atcdim")...
            Constructing wrapper function "specialfuncs.atcdim"...
              atcdim = atcdim(ntot,ltot)
                Creating wrapper for Fortran function "moffset"("moffset")...
            Constructing wrapper function "specialfuncs.moffset"...
              moffset = moffset(m,ntot,ltot)
            Constructing wrapper function "specialfuncs.cartosphere"...
              cartosphere(xp,r,ct,ep)
            Constructing wrapper function "specialfuncs.eulerrotation"...
getarrdims:warning: assumed shape array, using 0 instead of '*'
getarrdims:warning: assumed shape array, using 0 instead of '*'
              eulerrotation(xp,eulerangf,dir,xprot,[num])
            Constructing wrapper function "specialfuncs.ephicoef"...
              ephicoef(ep,epm,[nodr])
            Constructing wrapper function "specialfuncs.planewavetruncationorder"...
              planewavetruncationorder(r,rimedium,eps,nodr)
            Constructing wrapper function "specialfuncs.vwhcalc"...
              vwhcalc(rpos,ri,nodr,itype,vwh,[index_model,lr_to_mode])
            Constructing wrapper function "specialfuncs.scalar_wave_function"...
              scalar_wave_function(nodr,itype,x,y,z,ri,swf)
            Constructing wrapper function "specialfuncs.reciprocal_scalar_wave_function"...
              reciprocal_scalar_wave_function(nodr,kx,ky,x,y,z,ri,swf)
            Constructing wrapper function "specialfuncs.twobytwoinverse"...
              twobytwoinverse(mat,imat)
            Constructing wrapper function "specialfuncs.mtransfer"...
              mtransfer(cin,cout,[nin,nout])
                Creating wrapper for Fortran function "amnaddress"("amnaddress")...
            Constructing wrapper function "specialfuncs.amnaddress"...
              amnaddress = amnaddress(m,n,l,model)
                Creating wrapper for Fortran function "amnpaddress"("amnpaddress")...
            Constructing wrapper function "specialfuncs.amnpaddress"...
              amnpaddress = amnpaddress(m,n,p,l,model)
            Constructing wrapper function "specialfuncs.lr_mode_transformation"...
              lr_mode_transformation(nodr,alr,amode,[lr_to_mode])
            Constructing wrapper function "specialfuncs.degree_transformation"...
              degree_transformation(nodr,ain,aout)
            Constructing wrapper function "specialfuncs.lu_decomposition"...
              lu_decomposition(a,indx,d,ierr,[n])
            Constructing wrapper function "specialfuncs.lu_backsubstitution"...
              lu_backsubstitution(a,indx,b,[n])
            Constructing wrapper function "specialfuncs.groupfilename"...
              groupfilename(firststring,number,laststring,newstring)
            Constructing wrapper function "specialfuncs.qng"...
              qng(a,b,epsabs,epsrel,resultf,abserr,neval,ier,[n])
            Constructing wrapper function "specialfuncs.gkintegrate"...
routsign2map: Confused: function gkintegrate has externals ['qsub'] but no "use" statement.
sign2map: Confused: external qsub is not in lcb_map[].
append_needs: unknown need 'qsub'
append_needs: unknown need 'qsub'
              qint = gkintegrate(ntot,t0,t1,subdiv,errorcodes,inteps,mindiv,maxnumdiv)
            Constructing wrapper function "specialfuncs.realsort"...
              realsort(limits,eps,nlimits,[nlimits0])
                Constructing F90 module support for "surface_subroutines"...
                  Variables: source_sum include_direct_source pole_integration max_number_plane_boundaries max_singular_points source_order target_order max_azimuth_mode number_layers source_layer target_layer source_layer2 number_limits max_gf_iterations number_gf_iterations source_order2 error_codes energy_kernel_region number_singular_points singular_point_polarization number_plane_boundaries maximum_integration_subdivisions source_z target_z azimuth_angle incident_field_boundary radial_distance max_s gf_error_epsilon s_sc1 s_sc2 max_gf max_bf max_pi max_picon top_boundary bot_boundary source_z2 integration_error singular_points incident_field_scale pole_integration_radius singular_gf_value incident_lateral_vector g_cut g_sing_mag layer_thickness integration_limit_epsilon integration_error_epsilon real_axis_integration_limit gf_switch_factor s_scale_constant minimum_integration_spacing plane_boundary_position real_axis_limits source_ri target_ri pole_integration_s layer_ref_index source_coefficient source_coefficient_1 source_coefficient_2 layer_ref_index(0) layer_ref_index(1:max_number_plane_boundaries)
getarrdims:warning: assumed shape array, using 0 instead of ':'
getarrdims:warning: assumed shape array, using 0 instead of ':'
getarrdims:warning: assumed shape array, using 0 instead of ':'
getarrdims:warning: assumed shape array, using 0 instead of ':'
getarrdims:warning: assumed shape array, using 0 instead of ':'
getarrdims:warning: assumed shape array, using 0 instead of ':'
getarrdims:warning: assumed shape array, using 0 instead of ':'
getarrdims:warning: assumed shape array, using 0 instead of ':'
getarrdims:warning: assumed shape array, using 0 instead of ':'
getarrdims:warning: assumed shape array, using 0 instead of ':'
getctype: No C-type found in "{'=': '(1.d0', 'attrspec': ['optional']}", assuming void.
Traceback (most recent call last):
  File "/home/user/.local/bin/f2py3.10", line 8, in <module>
    sys.exit(main())
  File "/home/user/.local/lib/python3.10/site-packages/numpy/f2py/f2py2e.py", line 702, in main
    run_compile()
  File "/home/user/.local/lib/python3.10/site-packages/numpy/f2py/f2py2e.py", line 669, in run_compile
    setup(ext_modules=[ext])
  File "/home/user/.local/lib/python3.10/site-packages/numpy/distutils/core.py", line 169, in setup
    return old_setup(**new_attr)
  File "/usr/lib64/python3.10/distutils/core.py", line 148, in setup
    dist.run_commands()
  File "/usr/lib64/python3.10/distutils/dist.py", line 966, in run_commands
    self.run_command(cmd)
  File "/usr/lib64/python3.10/distutils/dist.py", line 985, in run_command
    cmd_obj.run()
  File "/home/user/.local/lib/python3.10/site-packages/numpy/distutils/command/build.py", line 62, in run
    old_build.run(self)
  File "/usr/lib64/python3.10/distutils/command/build.py", line 135, in run
    self.run_command(cmd_name)
  File "/usr/lib64/python3.10/distutils/cmd.py", line 313, in run_command
    self.distribution.run_command(command)
  File "/usr/lib64/python3.10/distutils/dist.py", line 985, in run_command
    cmd_obj.run()
  File "/home/user/.local/lib/python3.10/site-packages/numpy/distutils/command/build_src.py", line 144, in run
    self.build_sources()
  File "/home/user/.local/lib/python3.10/site-packages/numpy/distutils/command/build_src.py", line 161, in build_sources
    self.build_extension_sources(ext)
  File "/home/user/.local/lib/python3.10/site-packages/numpy/distutils/command/build_src.py", line 321, in build_extension_sources
    sources = self.f2py_sources(sources, ext)
  File "/home/user/.local/lib/python3.10/site-packages/numpy/distutils/command/build_src.py", line 562, in f2py_sources
    numpy.f2py.run_main(f2py_options + ['--lower',
  File "/home/user/.local/lib/python3.10/site-packages/numpy/f2py/f2py2e.py", line 475, in run_main
    ret = buildmodules(postlist)
  File "/home/user/.local/lib/python3.10/site-packages/numpy/f2py/f2py2e.py", line 397, in buildmodules
    dict_append(ret[name], rules.buildmodule(module, um))
  File "/home/user/.local/lib/python3.10/site-packages/numpy/f2py/rules.py", line 1254, in buildmodule
    mr, wrap = f90mod_rules.buildhooks(m)
  File "/home/user/.local/lib/python3.10/site-packages/numpy/f2py/f90mod_rules.py", line 143, in buildhooks
    at = capi_maps.c2capi_map[ct]
KeyError: 'void'