[cig-commits] [commit] devel: ran the script to remove accented letters (ff4f4cb)
cig_noreply at geodynamics.org
cig_noreply at geodynamics.org
Wed Jan 21 18:35:32 PST 2015
Repository : https://github.com/geodynamics/specfem2d
On branch : devel
Link : https://github.com/geodynamics/specfem2d/compare/97ea745c241e76b7692d21be0e9f68f43d669774...d2deb1c4f3ca714ba8230050d7133f96e285962a
>---------------------------------------------------------------
commit ff4f4cb44a9dd64e788edd29ced5ef276f043a5b
Author: Dimitri Komatitsch <komatitsch at lma.cnrs-mrs.fr>
Date: Thu Jan 22 03:07:51 2015 +0100
ran the script to remove accented letters
>---------------------------------------------------------------
ff4f4cb44a9dd64e788edd29ced5ef276f043a5b
..._on_the_adjoint_tomography_workflow_dec2014.txt | 14 ++---
...current_adjoint_tomography_workflow_dec2014.txt | 2 +-
.../mesh_quality_histogram.txt | 0
.../mesh_quality_histogram.txt | 0
...d_convert_Gmsh_to_Specfem2D_Nicolas_Massacret.c | 24 ++++-----
...d_pressure_in_addition_to_displacement_2013.txt | 2 +-
doc/USER_MANUAL/manual_SPECFEM2D.tex | 0
src/cuda/write_seismograms_cuda.cu | 4 +-
src/meshfem2D/scotch_5.1.12b/LICENCE_fr.txt | 7 ---
src/meshfem2D/scotch_5.1.12b/LISEZ-MOI.txt | 8 ---
src/specfem2D/attenuation_model.f90 | 62 +++++++++++-----------
src/specfem2D/init_host_to_dev_variable.f90 | 2 +-
src/specfem2D/iterate_time.F90 | 4 +-
src/specfem2D/prepare_timerun_gpu.f90 | 22 ++++----
src/specfem2D/specfem2D_par.f90 | 2 +-
15 files changed, 69 insertions(+), 84 deletions(-)
diff --git a/Carl_Tape_his_comments_on_the_adjoint_tomography_workflow_dec2014.txt b/Carl_Tape_his_comments_on_the_adjoint_tomography_workflow_dec2014.txt
index ffadaae..1a0340e 100644
--- a/Carl_Tape_his_comments_on_the_adjoint_tomography_workflow_dec2014.txt
+++ b/Carl_Tape_his_comments_on_the_adjoint_tomography_workflow_dec2014.txt
@@ -10,19 +10,19 @@ Hi Ryan and all,
I met with Ebru, Emanuele, Frederica, Elliott, and Lion on Thursday
night at AGU. (We scheduled at the last minute and missed Min and
Hejun, unfortunately.) Below are some discussion points and questions.
-You already have plenty to discuss anyway, so you don’t have to
+You already have plenty to discuss anyway, so you don't have to
consider any of these, of course. I cced Emanuele and Frederica,
-though I don’t think they’ll be back by Monday.
+though I don't think they'll be back by Monday.
+ It seems like those in our discussion were supportive of the move
-from “legacy†utilities (perl, sac, flexwin, measure_adj) to
+from "legacy" utilities (perl, sac, flexwin, measure_adj) to
python-based utilities (obspy, pyflex, pyadj). Lion seemed quite
-willing to rewrite measure_adj into python (pyadj) – the multitaper
+willing to rewrite measure_adj into python (pyadj) - the multitaper
measurement is the key piece, and there are Slepian tapers within
python libraries, apparently.
+ Lion suggested that, in his opinion, the hdf5 file format would be
-optimal for seismograms and adjoint sources. I don’t know the details
+optimal for seismograms and adjoint sources. I don't know the details
about this. But clearly this is an important decision, since it
impacts the processing scripts, adjoint sources, the moment tensor
inversion codes, etc.
@@ -63,10 +63,10 @@ little regarding the adjoint tomography workflow since 2009. But I
have several projects that will move forward in 2015, lead by me
(southern California), my student Vipul Silwal (Alaska), my postdoc
Amir Allam (fault zones up to 5-10 Hz), and Yoshihiro Kaneko (New
-Zealand). I’d like to get these folks on the right track and am open
+Zealand). I'd like to get these folks on the right track and am open
to trying out new things. Hopefully we can help in this effort.
-I’ll talk to you soon. Thanks for including me in the call.
+I'll talk to you soon. Thanks for including me in the call.
Carl
diff --git a/Daniel_Peter_his_comments_on_the_current_adjoint_tomography_workflow_dec2014.txt b/Daniel_Peter_his_comments_on_the_current_adjoint_tomography_workflow_dec2014.txt
index 2c179f9..f04cef3 100644
--- a/Daniel_Peter_his_comments_on_the_current_adjoint_tomography_workflow_dec2014.txt
+++ b/Daniel_Peter_his_comments_on_the_current_adjoint_tomography_workflow_dec2014.txt
@@ -11,7 +11,7 @@ for 1.
i kept the routines very similar, some differences occur when reading in mesh files and e.g. where the maximum of the gradient is taken for the update step length. in future, i want to have them also merged into the same set of SPECFEM function library.
for 2.
-good point, the jacobian is used for smoothing of the kernels in the global code, but isn’t currently used for the cartesian code. that is one of the differences that has to be checked again (thus smoothing in the cartesian is based solely on point distance, not volume)
+good point, the jacobian is used for smoothing of the kernels in the global code, but isn't currently used for the cartesian code. that is one of the differences that has to be checked again (thus smoothing in the cartesian is based solely on point distance, not volume)
best wishes,
daniel
diff --git a/UTILS/Gmsh/unofficial_unsupported_convert_Gmsh_to_Specfem2D_Nicolas_Massacret.c b/UTILS/Gmsh/unofficial_unsupported_convert_Gmsh_to_Specfem2D_Nicolas_Massacret.c
index 1c551bc..6c1defe 100644
--- a/UTILS/Gmsh/unofficial_unsupported_convert_Gmsh_to_Specfem2D_Nicolas_Massacret.c
+++ b/UTILS/Gmsh/unofficial_unsupported_convert_Gmsh_to_Specfem2D_Nicolas_Massacret.c
@@ -34,7 +34,7 @@ int main(int argc, char *argv[])
while ( strcmp(&chaine[0], &nodes[0]) != 0)
{
fscanf(fichier, "%s", &chaine[0]);
- } // Lorsque l'on sort de cette boucle, on a trouvé $Nodes
+ } // Lorsque l'on sort de cette boucle, on a trouve $Nodes
/* FDEBUT DE L'ECRITURE DE NODE */
fichierNode = fopen("nodes", "r+");
@@ -44,47 +44,47 @@ int main(int argc, char *argv[])
fprintf(fichierNode, "%li \n", infos[0]); // on note le nombre de noeuds
int nbrNode = infos[0];
car = fgetc(fichier);
- while( car != 36) // Boucle permettant l'écriture du fichier contenant les coordonées des noeuds§.
+ while( car != 36) // Boucle permettant l'ecriture du fichier contenant les coordonees des noeuds?.
{
fseek(fichier, -1, SEEK_CUR);
fscanf(fichier, "%f %f %f %f", &coordNoeud[0], &coordNoeud[1], &coordNoeud[2], &coordNoeud[3]);
fprintf(fichierNode, "%f %f \n", coordNoeud[1], coordNoeud[2]);
- if(coordNoeud[2] == 0) // Calcul du nombre d'éléments suivant x.
+ if(coordNoeud[2] == 0) // Calcul du nombre d'elements suivant x.
{
x++;
}
int etat = ((coordNoeud[0]*100) / nbrNode);
-printf("\r- Ecriture des coordonées des noeuds : %d %%", etat);
+printf("\r- Ecriture des coordonees des noeuds : %d %%", etat);
car = fgetc(fichier);
car = fgetc(fichier);
}
fclose(fichierNode);
}
- else {printf("!!! Erreur : Problème à l'ouverture du fichier nodes. \n"); return 0;}
+ else {printf("!!! Erreur : Probleme a l'ouverture du fichier nodes. \n"); return 0;}
/* FIN DE L'ECRITURE DE NODE */
- printf("\n- Ecriture du fichier node terminée. \n");
+ printf("\n- Ecriture du fichier node terminee. \n");
x--;
- printf("- Nombre d'éléments suivant x = %d \n", x);
+ printf("- Nombre d'elements suivant x = %d \n", x);
/* DEBUT DE L'ECRITURE DE ELEMENT ET METERIAL */
fichierElem = fopen("element", "r+");
if (fichierElem != NULL)
{
-// printf("début écriture element et material \n");
+// printf("debut ecriture element et material \n");
fichierMat = fopen("material", "r+");
if (fichierMat != NULL)
{
fscanf(fichier, "%s", &chaine[0]); // pour sauter une ligne
fscanf(fichier, "%s", &chaine[0]);
fscanf(fichier, "%li", &infos[0]);
- fprintf(fichierElem, "%li \n", infos[0]); // on note le nombre d'éléments
+ fprintf(fichierElem, "%li \n", infos[0]); // on note le nombre d'elements
int nbrElem = infos[0];
car = fgetc(fichier);
- while( car != 36) // Pour s'arrêter au prochain $.
+ while( car != 36) // Pour s'arreter au prochain $.
{
fseek(fichier, -1, SEEK_CUR);
fscanf(fichier, "%li %li %li", &infos[0], &infos[1], &infos[2]);
@@ -96,7 +96,7 @@ int nbrElem = infos[0];
fscanf(fichier, "%li", &numNode[i]);
}
int etat = ((infos[0]*100) / nbrElem);
-printf("\r- Ecriture des éléments et materiaux : %d %%", etat);
+printf("\r- Ecriture des elements et materiaux : %d %%", etat);
fprintf(fichierMat, "%li \n", numNode[0]); //on note le numero du matos
fscanf(fichier, "%li %li %li %li", &numNode[0], &numNode[1], &numNode[2], &numNode[3]);
fprintf(fichierElem, "%li %li %li %li\n", numNode[0], numNode[1], numNode[2], numNode[3]); // on note les ID des nodes
@@ -130,7 +130,7 @@ else
{
printf("Erreur ouverture fichier \n"); return 0;
}
-printf("\n- Ecriture des fichiers element et material terminée. \n");
+printf("\n- Ecriture des fichiers element et material terminee. \n");
return 0;
}
diff --git a/UTILS/Ra_Cleave_record_pressure_in_addition_to_displacement/email_from_Ra_Cleave_about_multiple_receiver_types_in_SPECFEM2D_to_record_pressure_in_addition_to_displacement_2013.txt b/UTILS/Ra_Cleave_record_pressure_in_addition_to_displacement/email_from_Ra_Cleave_about_multiple_receiver_types_in_SPECFEM2D_to_record_pressure_in_addition_to_displacement_2013.txt
index 9da9fb5..d93590d 100644
--- a/UTILS/Ra_Cleave_record_pressure_in_addition_to_displacement/email_from_Ra_Cleave_about_multiple_receiver_types_in_SPECFEM2D_to_record_pressure_in_addition_to_displacement_2013.txt
+++ b/UTILS/Ra_Cleave_record_pressure_in_addition_to_displacement/email_from_Ra_Cleave_about_multiple_receiver_types_in_SPECFEM2D_to_record_pressure_in_addition_to_displacement_2013.txt
@@ -10,7 +10,7 @@ Here are the routines I have modified. As I mentioned the code now optionally
records pressure at each receiver. Most changes are preceded by a comment that
starts with "!RC:". The par_file needs to have the following line:
-seismo_p = .true. # Whether pressure should also be recorded
+seismo_p = .true. # Whether pressure should also be recorded
Files:
SPECFEM2D-7.0.0-RC/src/specfem2D/write_seismograms.F90
diff --git a/src/cuda/write_seismograms_cuda.cu b/src/cuda/write_seismograms_cuda.cu
index 7971fac..2368f9c 100644
--- a/src/cuda/write_seismograms_cuda.cu
+++ b/src/cuda/write_seismograms_cuda.cu
@@ -211,7 +211,7 @@ void FC_FUNC_(compute_seismograms_cuda,
switch (seismotype){
- case 1 : //Déplacement
+ case 1 : //Deplacement
if (! *any_elastic_glob) printf("\nWrong type of seismogram for a pure fluid simulation, use pressure in seismotype\n");
compute_elastic_seismogram_kernel<<<grid,threads,0,mp->compute_stream>>>( mp->nrec_local,
@@ -241,7 +241,7 @@ void FC_FUNC_(compute_seismograms_cuda,
);
break;
- case 3 : //Accélération
+ case 3 : //Acceleration
if (! *any_elastic_glob) printf("\nWrong type of seismogram for a pure fluid simulation, use pressure in seismotype\n");
compute_elastic_seismogram_kernel<<<grid,threads,0,mp->compute_stream>>>( mp->nrec_local,
mp->d_accel,
diff --git a/src/meshfem2D/scotch_5.1.12b/LICENCE_fr.txt b/src/meshfem2D/scotch_5.1.12b/LICENCE_fr.txt
deleted file mode 100644
index 89250d5..0000000
--- a/src/meshfem2D/scotch_5.1.12b/LICENCE_fr.txt
+++ /dev/null
@@ -1,7 +0,0 @@
-Cette copie de la distribution logicielle Scotch 5.1 est distribuée
-selon les termes de la licence de logiciel libre CeCILL-C.
-
-Une copie du texte de cette licence se trouve dans le sous-répertoire
-"doc", dans le fichier appelé "CeCILL-C_V1-fr.txt". Veuillez en lire
-attentivement les clauses. Vous ne pourrez utiliser cette copie de la
-distribution logicielle Scotch 5.1 que si vous les acceptez.
diff --git a/src/meshfem2D/scotch_5.1.12b/LISEZ-MOI.txt b/src/meshfem2D/scotch_5.1.12b/LISEZ-MOI.txt
deleted file mode 100644
index 1dc6dab..0000000
--- a/src/meshfem2D/scotch_5.1.12b/LISEZ-MOI.txt
+++ /dev/null
@@ -1,8 +0,0 @@
-Les clauses selon lesquelles cette copie de la distribution logicielle
-Scotch 5.1 est mise à votre disposition sont définies dans le fichier
-"LICENCE_fr.txt", situé dans le même répertoire que le présent
-fichier.
-
-Si vous les acceptez, veuillez vous référer au fichier "INSTALL.txt",
-également situé dans ce répertoire, pour consulter les instructions
-d'installation.
diff --git a/src/specfem2D/attenuation_model.f90 b/src/specfem2D/attenuation_model.f90
index e4a4c7c..ec895a6 100644
--- a/src/specfem2D/attenuation_model.f90
+++ b/src/specfem2D/attenuation_model.f90
@@ -302,16 +302,16 @@
! Lombard and Piraux, Numerical modeling of transient two-dimensional viscoelastic waves,
! Journal of Computational Physics, Volume 230, Issue 15, Pages 6099-6114 (2011)
-! Suivant les compilateurs et les options de compilation utilisées,
-! il peut y avoir des différences au 4ème chiffre significatif. C'est sans conséquences sur la précision du calcul :
-! l'erreur est de 0.015 % avec optimisation non linéaire, à comparer à 1.47 % avec Emmerich and Korn (1987).
-! Si je relance le calcul en initialisant avec le résultat précédent, ce chiffre varie à nouveau très légèrement.
+! Suivant les compilateurs et les options de compilation utilisees,
+! il peut y avoir des differences au 4eme chiffre significatif. C'est sans consequences sur la precision du calcul :
+! l'erreur est de 0.015 % avec optimisation non lineaire, a comparer a 1.47 % avec Emmerich and Korn (1987).
+! Si je relance le calcul en initialisant avec le resultat precedent, ce chiffre varie a nouveau tres legerement.
!-------------------------------------------------------------------------
! From Bruno Lombard, June 2014:
-! j'ai relu en détail :
+! j'ai relu en detail :
! [1] Carcione, Kosslof, Kosslof, "Viscoacoustic wave propagation simulation in the Earth",
! Geophysics 53-6 (1988), 769-777
@@ -322,24 +322,24 @@
! [3] Moczo, Kristek, "On the rheological models used for time-domain methods of seismic wave propagation",
! Geophysical Research Letters 32 (2005).
-! Le problème provient probablement d'une erreur récurrente dans [1,2] et datant de Liu et al 1976 :
-! l'oubli du facteur 1/N dans la fonction de relaxation d'un modèle de Zener à N éléments.
-! Il est effectivement facile de faire l'erreur. Voir l'équation (12) de [3], et le paragraphe qui suit.
+! Le probleme provient probablement d'une erreur recurrente dans [1,2] et datant de Liu et al 1976 :
+! l'oubli du facteur 1/N dans la fonction de relaxation d'un modele de Zener a N elements.
+! Il est effectivement facile de faire l'erreur. Voir l'equation (12) de [3], et le paragraphe qui suit.
-! Du coup le module de viscoélasticité est faux dans [1,2], et donc le facteur de qualité,
+! Du coup le module de viscoelasticite est faux dans [1,2], et donc le facteur de qualite,
! et donc les temps de relaxation tau_sigma...
-! Après, [2] calcule une solution analytique juste, mais avec des coefficients sans sens physique.
-! Et quand SPECFEM2D obtient un bon accord avec cette solution analytique, ça valide SPECFEM, mais pas le calcul des coefficients.
+! Apres, [2] calcule une solution analytique juste, mais avec des coefficients sans sens physique.
+! Et quand SPECFEM2D obtient un bon accord avec cette solution analytique, ca valide SPECFEM, mais pas le calcul des coefficients.
! Il y a donc une erreur dans [1,2], et [3] a raison.
-! Sa solution analytique découle d'un travail sur ses fonctions de relaxation (A4),
+! Sa solution analytique decoule d'un travail sur ses fonctions de relaxation (A4),
! qu'il injecte ensuite dans la relation de comportement (A1) et travaille en Fourier.
-! Le problème est que sa fonction de relaxation (A4) est fausse : il manque 1/N.
-! De ce fait, sa solution analytique est cohérente avec sa solution numérique.
-! Dans les deux cas, ce sont les mêmes temps de relaxation qui sont utilisés. Mais ces temps sont calculés de façon erronée.
+! Le probleme est que sa fonction de relaxation (A4) est fausse : il manque 1/N.
+! De ce fait, sa solution analytique est coherente avec sa solution numerique.
+! Dans les deux cas, ce sont les memes temps de relaxation qui sont utilises. Mais ces temps sont calcules de facon erronee.
!-------------------------------------------------------------------------
@@ -354,36 +354,36 @@
! From Emilie Blanc, April 2014:
-! le programme SolvOpt d'optimisation non-linéaire
-! avec contrainte. Ce programme prend quatre fonctions en entrée :
+! le programme SolvOpt d'optimisation non-lineaire
+! avec contrainte. Ce programme prend quatre fonctions en entree :
-! - fun() est la fonction à minimiser
+! - fun() est la fonction a minimiser
-! - grad() est le gradient de la fonction à minimiser par rapport à chaque paramètre
+! - grad() est le gradient de la fonction a minimiser par rapport a chaque parametre
-! - func() est le maximum des résidus (= 0 si toutes les contraintes sont satisfaites)
+! - func() est le maximum des residus (= 0 si toutes les contraintes sont satisfaites)
-! - gradc() est le gradient du maximum des résidus (= 0 si toutes les
+! - gradc() est le gradient du maximum des residus (= 0 si toutes les
! contraintes sont satisfaites)
-! Ce programme a été développé par Kappel et Kuntsevich. Leur article décrit l'algorithme.
+! Ce programme a ete developpe par Kappel et Kuntsevich. Leur article decrit l'algorithme.
-! J'ai utilisé ce code pour la poroélasticité haute-fréquence, et aussi en
-! viscoélasticité fractionnaire (modèle d'Andrade, avec Bruno Lombard et
-! Cédric Bellis). Nous pouvons interagir sur l'algorithme d'optimisation
-! pour votre modèle visco, et étudier l'effet des coefficients ainsi obtenus.
+! J'ai utilise ce code pour la poroelasticite haute-frequence, et aussi en
+! viscoelasticite fractionnaire (modele d'Andrade, avec Bruno Lombard et
+! Cedric Bellis). Nous pouvons interagir sur l'algorithme d'optimisation
+! pour votre modele visco, et etudier l'effet des coefficients ainsi obtenus.
!---------------------------------------------------
! From Emilie Blanc, March 2014:
-! Les entrées du programme principal sont le nombre de variables
-! diffusives, le facteur de qualité voulu Qref et la fréquence centrale f0.
+! Les entrees du programme principal sont le nombre de variables
+! diffusives, le facteur de qualite voulu Qref et la frequence centrale f0.
-! Cependant, pour l'optimisation non-linéaire, j'ai mis theta_max=100*f0
+! Cependant, pour l'optimisation non-lineaire, j'ai mis theta_max=100*f0
! et non pas theta_max=2*pi*100*f0. En effet, dans le programme, on
-! travaille sur les fréquences, et non pas sur les fréquences angulaires.
-! Cela dit, dans les deux cas j'obtiens les mêmes coefficients...
+! travaille sur les frequences, et non pas sur les frequences angulaires.
+! Cela dit, dans les deux cas j'obtiens les memes coefficients...
!---------------------------------------------------
diff --git a/src/specfem2D/init_host_to_dev_variable.f90 b/src/specfem2D/init_host_to_dev_variable.f90
index e378f36..6043818 100644
--- a/src/specfem2D/init_host_to_dev_variable.f90
+++ b/src/specfem2D/init_host_to_dev_variable.f90
@@ -456,7 +456,7 @@ num_elem_colors_acoustic(1)=0
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
-!!! Initialisation paramètres pour simulation elastique
+!!! Initialisation parametres pour simulation elastique
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
diff --git a/src/specfem2D/iterate_time.F90 b/src/specfem2D/iterate_time.F90
index 182ea48..a02657a 100644
--- a/src/specfem2D/iterate_time.F90
+++ b/src/specfem2D/iterate_time.F90
@@ -1089,8 +1089,8 @@ if (myrank == 0) write(IOUT,400)
! impose potential_gravito in order to have z displacement equal to forced
! value on the boundary
-!!!! Passe deux fois sur le même iglob
-!!!! Mais vrai pour tous les points partagés entre deux elements
+!!!! Passe deux fois sur le meme iglob
+!!!! Mais vrai pour tous les points partages entre deux elements
iglob = ibool(i,j,ispec)
displ_n = displ_x*nx + displ_z*nz
if (abs(nz) > TINYVAL) then
diff --git a/src/specfem2D/prepare_timerun_gpu.f90 b/src/specfem2D/prepare_timerun_gpu.f90
index 20ca10c..e1a7826 100644
--- a/src/specfem2D/prepare_timerun_gpu.f90
+++ b/src/specfem2D/prepare_timerun_gpu.f90
@@ -83,7 +83,7 @@
! count_nspec_acoustic : nombre local d'elements spectraux acoustiques
! prepares general fields on GPU
- !§!§ JC JC here we will need to add GPU support for the new C-PML routines
+ !! JC JC here we will need to add GPU support for the C-PML routines
call prepare_constants_device(Mesh_pointer, &
@@ -127,7 +127,7 @@
!!! Parametres fournis
-! rmass_inverse_acoustic : matrice acoustique inversée de taille nglob (nglob_acoustic = nglob s'il existe des elements acoustiques)
+! rmass_inverse_acoustic : matrice acoustique inversee de taille nglob (nglob_acoustic = nglob s'il existe des elements acoustiques)
! num_phase_ispec_acoustic : max entre nb d'element spectraux acoustiques interieur et exterieur
! phase_ispec_inner_acoustic(i,j) : i eme element spectral acoustique interieur si j=2 exterieur si j=1
! acoustic(i) : vrai si l'element spectral i est acoustique
@@ -136,14 +136,14 @@
! free_surface_ij(i,j,ispec) : i eme coordonnee du j eme point GLL de l'element spectral libre ispec
! b_reclen_potential : place en octet prise par b_nelem_acoustic_surface * GLLX
! any_elastic : vrai s'il existe des elements elastiques
-! num_fluid_solid_edges : nombre d'élements spectraux sur une frontiere elastique/acoustique
+! num_fluid_solid_edges : nombre d'elements spectraux sur une frontiere elastique/acoustique
! coupling_ac_el_ispec : tableau des elements spectraux frontiere ACOUSTIQUE
! coupling_ac_el_ij : coordonnees locales des points GLL sur la frontiere elastique/acoustique
! coupling_ac_el_normal(i,j,ispec) : i eme coordonne de la normale au point GLL j de l'element frontiere ispec
! coupling_ac_el_jacobian1Dw(i,ispec) : jacobienne ponderee du i eme point GLL de l'element frontiere ispec
-! num_colors_outer_acoustic : a initialiser plus tard quand USE_COLOR_MESH sera implementé
-! num_colors_inner_acoustic : a initialiser plus tard quand USE_COLOR_MESH sera implementé
-! num_elem_colors_acoustic : a initialiser plus tard quand USE_COLOR_MESH sera implementé
+! num_colors_outer_acoustic : a initialiser plus tard quand USE_COLOR_MESH sera implemente
+! num_colors_inner_acoustic : a initialiser plus tard quand USE_COLOR_MESH sera implemente
+! num_elem_colors_acoustic : a initialiser plus tard quand USE_COLOR_MESH sera implemente
! prepares fields on GPU for acoustic simulations
@@ -169,16 +169,16 @@
!!! Parametres fournis
-! rmass_inverse_elastic : matrice elastique inversée de taille nglob_acoustic (nglob_acoustic = nglob s'il existe des elements acoustiques)
+! rmass_inverse_elastic : matrice elastique inversee de taille nglob_acoustic (nglob_acoustic = nglob s'il existe des elements acoustiques)
! num_phase_ispec_elastic : max entre nb d'element spectraux elastiques interieur et exterieur
! phase_ispec_inner_elastic(i,j) : i eme element spectral elastique interieur si j=2 exterieur si j=1
! elastic(i) : vrai si l'element spectral i est elastique
-! num_colors_outer_elastic : a initialiser plus tard quand USE_COLOR_MESH sera implementé
-! num_colors_inner_elastic : a initialiser plus tard quand USE_COLOR_MESH sera implementé
-! num_elem_colors_elastic : a initialiser plus tard quand USE_COLOR_MESH sera implementé
+! num_colors_outer_elastic : a initialiser plus tard quand USE_COLOR_MESH sera implemente
+! num_colors_inner_elastic : a initialiser plus tard quand USE_COLOR_MESH sera implemente
+! num_elem_colors_elastic : a initialiser plus tard quand USE_COLOR_MESH sera implemente
! prepares fields on GPU for elastic simulations
- !§!§ JC JC here we will need to add GPU support for the new C-PML routines
+ !?!? JC JC here we will need to add GPU support for the new C-PML routines
if( any_elastic ) then
call prepare_fields_elastic_device(Mesh_pointer, &
rmass_inverse_elastic_one,rmass_inverse_elastic_three, &
diff --git a/src/specfem2D/specfem2D_par.f90 b/src/specfem2D/specfem2D_par.f90
index 0388d2b..21a82f7 100644
--- a/src/specfem2D/specfem2D_par.f90
+++ b/src/specfem2D/specfem2D_par.f90
@@ -713,7 +713,7 @@ module specfem_par
! parameters used in LDDRK scheme, from equation (2) of
! Berland, J., Bogey, C., & Bailly, C.
- ! Low-dissipation and low-dispersion fourth-order Runge–Kutta algorithm, Computers & Fluids, 35(10), 1459-1463.
+ ! Low-dissipation and low-dispersion fourth-order Runge-Kutta algorithm, Computers & Fluids, 35(10), 1459-1463.
Data alpha_LDDRK /0.0_CUSTOM_REAL,-0.737101392796_CUSTOM_REAL, &
-1.634740794341_CUSTOM_REAL,-0.744739003780_CUSTOM_REAL, &
-1.469897351522_CUSTOM_REAL,-2.813971388035_CUSTOM_REAL/
More information about the CIG-COMMITS
mailing list