'fig11_JOSAA_2004' !The file names resulting from this computation will contain this chain 0.3 ! longnorm: length of the normals, computed in prof_norm.XXX. Useless! 400 !Number of sources N refering to section 5 of G. Tayeb and S. Enoch, "Combined fictitious-sources–scattering-matrix method," J. Opt. Soc. Am. A 21, 1417-1423 (2004) 13.944 !a: half width (along x) of the jacket 2.6996 !b: half height (along y) of the jacket 1. !coef4: coef4>1 put the sources farther from the profile; coef4<1 put them closer. Begin with 1 if you do not know, then experiment... 200 !nfpoly: number of Fourier coefficients used to describe the profile, called n in Eq. (10) of G. Tayeb and S. Enoch, "Combined fictitious-sources–scattering-matrix method," J. Opt. Soc. Am. A 21, 1417-1423 (2004) 1.,0.1,1 !re_eps1,inc,nbre: real part of the relative permittivity of the external medium (initial, step, number of computations) 0.,0.1,1 !im_eps1,inc,nbre: imaginary part of the relative permittivity of the external medium (initial, step, number of computations) (3.4641,0.) !cnumater2 : complex optical index of the jacket [time dependency in exp(-iwt)] 2.02,0.,1 !wlvi,incwlv,nwlv: wavelength in vacuum (initial, step, number of computations) 0 !npol: POLARISATION, 0=E// to the rods, 1=H// to the rods 2 !nchampinc: [0=plane wave] [1=wire source outside the jacket] [2=Gaussian beam] 30. !tetad : incidence angle (degres), taken from Oy (only used for plane wave or Gaussian beam) (25.,0.) !zfilsou: position (x,y) of the wire source (only used if nchampinc=1) 5. !dpaq: Gaussian beam width (only used for Gaussian beam) 3.,0. !x0paq,y0paq: waist position of the Gaussian beam 60 !nbop: number of plane waves used to build the Gaussian beam .false. !ldl: curvilinear abscissa for scaler products, see below .false. !lsvd: (true:SVD) (false:Eq normales de Gauss) 180 !nntot: Fourier coefficients of the diffracted field: from -nntot to + nntot 720 !npp: nbre of pts for g(teta). Take >= 2*(nntot+1) 3 !nntig: truncation of the diffracted field (et eventuellement interieur) for each rod (from -nntig to +nntig) 1.d-4 !relconv: recommanded value: 1.d-4. Display warnings if the bn convergence for the worse rod is less more bad than relconv 0. !centx: center of map abscissa (for field maps) 0. !centy: center of map ordinate (for field maps) 30. !lgrilx: horizontal side of the map 30. !lgrily: vertical side of the map 0.3 !quad: grid spacing (for field maps) (if=0., the map is not computed) 40 !Number of rods along x axis, on the first horizontal line of rods starting from the bottom 9 !Number of rows of rods, along the y axis -2.355589 !y coordinate of the centers of the first line of rods 0.68 !Distance between centers of rods (1.,0.) !Optical index of the rods 0.294 !Radius of the rods .false. !lacfich: intermediate file results accumulation 1 !nresu: file results (0=ser>et only)(1=all) --------------------------------------------------------------------------- Indices optiques: Pour cnumater2 et les indices lus dans nomfipos: Si Im(cnumater2)=-1. , l'indice varie en fonction de wlv, et le matériau est choisi est fonction de Re(cnumater2) selon la correspondance suivante (à compléter dans subroutine calcindice au fur et à mesure des besoins): cnumater2=(1.,-1.) => Au cnumater2=(2.,-1.) => Ag Sinon, on lit directement l'indice, qui est alors indépendant de wlv --------------------------------------------------------------------------- ATTENTION: ce programme a été conçu en pensant aux cas: nprofil = 3 ou 4, et devrait aussi fonctionner avec 1 et 2 opts = 4 Pour les autres cas, il y a probablement des ajustements à faire --------------------------------------------------------------------------- Differents profils pour la gaine: nprofil=1: Ellipse, on peut choisir la taille (a et b) et toutes options opts ci-dessous. 2: Haricot, profil défini dans subroutine haricot, opts= 4 ou 20 seulement. 3: Polygone à angles arrondis. Profil défini dans subroutine polygone, opts= 4 ou 20 seulement. Le cas actuellement programmé est un rectangle centré à l'origine. Les sommets sont calculés à partir de a et b. 4: Polygone quelconque à angles arrondis. Les sommets sont lus dans nomfipoly. opts= 4 ou 20 seulement (à vérifier). --------------------------------------------------------------------------- Differentes options de sources: opts=0: comme MVDE; e1 et e2 sont alors les coeff d'homothetie (e1<1 et e2>1) 1: selon teta; e1 et e2 sont alors des distances 2: vers centre ellipse; e1 et e2 sont alors des distances 3: sur normale; e1 et e2 sont alors des distances 4: fonction du rayon de courbure (nd non utilise; nd=2*nb). ATTENTION: visualiser la position des sources, car elles peuvent se placer du mauvais côté du profil. 10: nappes 20: pour un point du profil sur 2, on porte sur la normale coef4 * la distance entre 2 points consécutifs du profil ATTENTION: visualiser la position des sources, car elles peuvent se placer du mauvais côté du profil. --------------------------------------------------------------------------- nprosc et ldl: nprosc : le produit scalaire de la colonne (u,v) est calcule en integrant sur P: |u|**2 avec ponderation cp1 (u=champ) |v|**2 avec ponderation cp3 (v=dérivée normale) si nprosc=0 , cp3=1 si nprosc=5 , cp3 obtenu a partir de moyennes sur tous les vect. de base, et depend de la polarisation NOTA: nprosc n'est pas utilise pour les tiges IC NOTA: nprosc n'est pas utilise si la resolution se fait au moyen de SVD qui utilise son propre produit scalaire (equivalent a nprosc=0) ldl : T = on pondere les produits scalaires en tenant compte de la distance entre les pts sur P F = pas de ponderation ---------------------------------------------------------------------------