&namparam ! method :: 0 : bartlett : méthode linéaire (valeur par défaut) ! 1 ! mvdr : méthode adaptative - inversion sur K (voir option nivBruit) ! ln_csdm :: 1 : build CSDM and compute as conj(replica(:))*CSDM(:,:)*replica(:) ! 0 : z*conj(z) with z = sum (exp (j*(phi_data-phi_model(:)))) ! nivBruit :: ajout de bruit sur la diagonale de K (si nivBruit non nul, on ajoute norm(K)*10^(-nivBruit/10)) (0 par défaut) ! normalisation :: 0 : noNorm : pas de normalisation (valeur par défaut) ! 1 : normSp : normalisation fréquentielle (blanchiement: "prewhitening") ! 2 : normTps : normalisation en temps (norm. par l'énergie) ! azimuth :: 1 : beam(c,theta) see namtheta section ! 0 : beam(c,x,y) see namxy section method = 0 ln_csdm = 0 nivBruit = 0 normalisation = 1 azimuth = 0 Fs = 500 ! input data sampling rate (Hz), must be an integer cn_ifile = "<<>>" ! input hdf5 file jbeg = 0 ! The starting point (in nb of samples, not in sec) (Caution: index starts at 0) jlen = 500 ! Window length (in nb of samples) joverlap = 250 ! overlap (in nb of samples), joverlap=jlen means no overlapping nit = 5 ! 43200 ! 5 ! 172799 ! 43200 ! 18000 ! 86399 / &namcoh ! ln_coh :: 0 : separate bartlett contribution for each frequency ! 1 : compute bartlett coherently (sum all station and frequencies both is a large sum) only with bartlett method ln_csdm=0 ln_coh = 1 / &namio nn_io = 5 ! 1 ! 7200 ! 1 ! 7200 ! read input every nn_io iterations / &velmodel ln_vm = 0 cn_vmfile = "GridTrav3.h5" Ngc = 3154000 ! 4000000 ! number of grid nodes for velocity model / &namtheta ! Range of azimuth (degree) [min:step:max] ! case azimuth=1 tmin = 0 tmax = 360 tstep = 10 / &namxy ! Range of x and y (m) [min:step:max] ! case azimuth=0 ! if nn_optim=2, use [xyz]min/max for bounds and [xyz]step for starting point ! if nn_optim=4 (RESOLVE Argentière ZO_2018 experiment), ! [xyz]min/max are not used ! and [xyz]step are used as [xyz]start coordinates of starting point xmin = -3000.0 xmax = 3000.0 xstep = <<>> ! used as xstart ymin = -3000.0 ymax = 3000.0 ystep = <<>> ! used as ystart / &namHsource ! elevation of the source (m above sea level) [min:step:max] ! case azimuth=0 ! if nn_optim=2, use [xyz]min/max for bounds and [xyz]step for starting point ! if nn_optim=4 (RESOLVE Argentière ZO_2018 experiment), ! [xyz]min/max are not used ! and [xyz]step are used as [xyz]start coordinates of starting point zmin = 1613.0 zmax = 2391.0 zstep = <<>> ! used as zstart / &namspeed ! Range of speed (m/s) [min:step:max] ! if nn_optim=4 (RESOLVE Argentière ZO_2018 experiment), ! [c]min/max are not used ! and cstep is used as cstart cmin = 1300.0 cmax = 3600.0 cstep = 1800.0 ! used as cstart / &namindice_freq ! Range of frequencies (Hz) [min:step:max] fmin = 15.0 fmax = 19.0 fstep = 0.1 ! if you want non regular fstep, use negative integer (example fstep = -10), thus fmin and fmax are ignored and ! user should provide freq.txt ASCII 10-line file with 10 frequencies / &namcoord nn_latlon = 0 ! [0] : input station coords are given in x/y (same as starting point search grid) ! [1] : input station coords are given in latitude/longitude, and ! the code defines its own centered grid (depending on active stations) / &naminput nn_removezero = 0 ! [0] : keep input data even if zero everywhere (during the whole jlen) ! [1] : remove station if zero everywhere during jlen, thus number of station may change at each iteration / &namrandom ln_randomdata = 0 ! [0] : ! [1] : replace phidata (from input data) by random [0,2pi] (Note that we still read input data anyway) / &namdiags nn_diags = 0 ! [0] : ! [1] : save some useful diagnostics about the optimisation path (should be used only if nit=1) / &namoptim nn_optim = 4 ! flag for optimisation parameters ! [1] : xyc (Chloe) : use namxy and namspeed sections to iterate over start for optimisation ! [2] : xyz (Gosia) : fixed cel = cstep (cmin and cmax not used) ! use namxy and namHsource to init starting point and for boundings ! [3] : xyc (Argentiere) : use namxy and namspeed to init starting point and for boundings ! fixed z = zstep ! [4] : xyzc (Argentiere) : use namxy and namspeed and namHsource to init starting point and for boundings nn_nlopt = 0 ! nlopt [1] or asa047 [0] library for optimisation dn_reqmin = 1.0D-02 ! the terminating limit for the variance of function values dn_istep_x = 10.0D+00 ! determines the size and shape of the initial simplex. ! The relative magnitudes of its elements should reflect the units of the variables. dn_istep_y = 10.0D+00 dn_istep_z = 10.0D+00 ! dn_istep_c = 10.0D+00 ! nn_konvge = 100 ! the convergence check is carried out every KONVGE iterations. nn_kcount = 3000 ! the maximum number of function evaluations. / &nambound ! namelist section to define boundary conditions for the optimisation ! z penalty only nn_zpen = 1 ! [1] : Apply penalty on Bartlett value if function converges above parametric surface (only for nn_optim=4) dn_c0 = 2367.0 dn_cx = 0.08378 dn_cxx = -4.837E-05 dn_cy = -0.0425 dn_cyy = -0.0001038 dn_cxy = 1.698E-05 ! exemple argentiere ! real(kind=8), parameter :: p00 = 2367 , p10 = 0.08378 , p01 = -0.0425 , p20 = -4.837E-05 , p11 = 1.698E-05 , p02 = -0.0001038 ! zz = p00 + p10*xx + p01*yy + p20*xx*xx + p11*xx*yy + p02*yy*yy ! exemple sggs ! zz = 2028.0 and all other values = 0 / &namsa ! only for beamforming_simulated_annealing.f90 code NSA = 15000 ! 30000 i_restart = 1250 Step_param = 50 N_pdf = 7000 T0 = 6.0 ! 15.0 T0min = 0.6 T_param_beg = 1.0 T_param_end = 1.0 N = 4 ! numbr of params for grid search : xyc or xyzc /