[cig-commits] r19803 - in seismo/2D/SPECFEM2D/trunk: DATA EXAMPLES EXAMPLES/Abel_Brest EXAMPLES/Gmsh_example_MPI EXAMPLES/Gmsh_example_serial EXAMPLES/INDUSTRIAL_FORMAT EXAMPLES/M2_UPPA EXAMPLES/Rayleigh_wave_no_crack EXAMPLES/Tape2007 EXAMPLES/Tape2007_kernel EXAMPLES/Tromp2005 EXAMPLES/Tromp2005_kernel EXAMPLES/acoustic_poroelastic EXAMPLES/attenuation EXAMPLES/canyon EXAMPLES/fluid_solid/fluid_solid_external_mesh EXAMPLES/fluid_solid/from_2000_Geophysics_paper_flat_ocean_bottom EXAMPLES/fluid_solid/from_2000_Geophysics_paper_sinusoidal_ocean_bottom EXAMPLES/init_plane EXAMPLES/noise_layered/model_0 EXAMPLES/noise_layered/model_1 EXAMPLES/noise_layered/model_2 EXAMPLES/noise_uniform EXAMPLES/semi_infinite_homo UTILS setup
dkomati1 at geodynamics.org
dkomati1 at geodynamics.org
Sun Mar 18 18:25:42 PDT 2012
Author: dkomati1
Date: 2012-03-18 18:25:41 -0700 (Sun, 18 Mar 2012)
New Revision: 19803
Added:
seismo/2D/SPECFEM2D/trunk/EXAMPLES/change_all_2D_Par_files.csh
seismo/2D/SPECFEM2D/trunk/EXAMPLES/clean_a_Par_file.py
seismo/2D/SPECFEM2D/trunk/EXAMPLES/process_DATA_Par_files_to_update_their_format_when_new_parameters_are_added.py
Removed:
seismo/2D/SPECFEM2D/trunk/UTILS/process_DATA_Par_files_to_update_their_format_when_new_parameters_are_added.py
Modified:
seismo/2D/SPECFEM2D/trunk/DATA/Par_file
seismo/2D/SPECFEM2D/trunk/EXAMPLES/Abel_Brest/Par_file_Abel_Balanche_bathy_source_solid
seismo/2D/SPECFEM2D/trunk/EXAMPLES/Gmsh_example_MPI/Par_file_Gmsh_SqrCircles.in
seismo/2D/SPECFEM2D/trunk/EXAMPLES/Gmsh_example_serial/Par_file_Gmsh_SqrCircles.in
seismo/2D/SPECFEM2D/trunk/EXAMPLES/INDUSTRIAL_FORMAT/Par_file
seismo/2D/SPECFEM2D/trunk/EXAMPLES/M2_UPPA/Par_file_M2_UPPA
seismo/2D/SPECFEM2D/trunk/EXAMPLES/Rayleigh_wave_no_crack/Par_file_Rayleigh_2D
seismo/2D/SPECFEM2D/trunk/EXAMPLES/Tape2007/Par_file_Tape2007_132rec_checker
seismo/2D/SPECFEM2D/trunk/EXAMPLES/Tape2007/Par_file_Tape2007_onerec
seismo/2D/SPECFEM2D/trunk/EXAMPLES/Tape2007_kernel/Par_file_Tape2007_onerec
seismo/2D/SPECFEM2D/trunk/EXAMPLES/Tromp2005/Par_file_Tromp2005
seismo/2D/SPECFEM2D/trunk/EXAMPLES/Tromp2005/Par_file_Tromp2005_s100
seismo/2D/SPECFEM2D/trunk/EXAMPLES/Tromp2005_kernel/Par_file_Tromp2005
seismo/2D/SPECFEM2D/trunk/EXAMPLES/acoustic_poroelastic/Par_file_acoustic_poroelastic
seismo/2D/SPECFEM2D/trunk/EXAMPLES/attenuation/Par_file_attenuation_2D
seismo/2D/SPECFEM2D/trunk/EXAMPLES/canyon/Par_file_canyon
seismo/2D/SPECFEM2D/trunk/EXAMPLES/fluid_solid/fluid_solid_external_mesh/Par_file_fluid_solid
seismo/2D/SPECFEM2D/trunk/EXAMPLES/fluid_solid/from_2000_Geophysics_paper_flat_ocean_bottom/Par_file_fluid_solid
seismo/2D/SPECFEM2D/trunk/EXAMPLES/fluid_solid/from_2000_Geophysics_paper_sinusoidal_ocean_bottom/Par_file_fluid_solid
seismo/2D/SPECFEM2D/trunk/EXAMPLES/init_plane/Par_file_Slave
seismo/2D/SPECFEM2D/trunk/EXAMPLES/init_plane/Par_file_Slave_for
seismo/2D/SPECFEM2D/trunk/EXAMPLES/init_plane/Par_file_Slave_kernel
seismo/2D/SPECFEM2D/trunk/EXAMPLES/noise_layered/model_0/Par_file_fair
seismo/2D/SPECFEM2D/trunk/EXAMPLES/noise_layered/model_0/Par_file_good
seismo/2D/SPECFEM2D/trunk/EXAMPLES/noise_layered/model_1/Par_file_best
seismo/2D/SPECFEM2D/trunk/EXAMPLES/noise_layered/model_1/Par_file_fair
seismo/2D/SPECFEM2D/trunk/EXAMPLES/noise_layered/model_1/Par_file_good
seismo/2D/SPECFEM2D/trunk/EXAMPLES/noise_layered/model_2/Par_file_fair
seismo/2D/SPECFEM2D/trunk/EXAMPLES/noise_layered/model_2/Par_file_good
seismo/2D/SPECFEM2D/trunk/EXAMPLES/noise_uniform/Par_file_noise_1
seismo/2D/SPECFEM2D/trunk/EXAMPLES/noise_uniform/Par_file_noise_2
seismo/2D/SPECFEM2D/trunk/EXAMPLES/noise_uniform/Par_file_noise_3
seismo/2D/SPECFEM2D/trunk/EXAMPLES/semi_infinite_homo/Par_file_elastic_2D
seismo/2D/SPECFEM2D/trunk/setup/constants.h.in
Log:
added my scripts to clean all the Par_files.
set Cuthill-McKee sorting to .false. by default.
Modified: seismo/2D/SPECFEM2D/trunk/DATA/Par_file
===================================================================
--- seismo/2D/SPECFEM2D/trunk/DATA/Par_file 2012-03-18 19:59:34 UTC (rev 19802)
+++ seismo/2D/SPECFEM2D/trunk/DATA/Par_file 2012-03-19 01:25:41 UTC (rev 19803)
@@ -9,7 +9,7 @@
# parameters concerning partitioning
nproc = 1 # number of processes
partitioning_method = 3 # SCOTCH = 3, ascending order (very bad idea) = 1
-PERFORM_CUTHILL_MCKEE = .true. # perform inverse Cuthill-McKee (1969) optimization/permutation for mesh numbering
+PERFORM_CUTHILL_MCKEE = .false. # perform inverse Cuthill-McKee (1969) optimization/permutation for mesh numbering (can be very costly and not very useful)
ngnod = 9 # number of control nodes per element (4 or 9)
initialfield = .false. # use a plane wave as source or not
Modified: seismo/2D/SPECFEM2D/trunk/EXAMPLES/Abel_Brest/Par_file_Abel_Balanche_bathy_source_solid
===================================================================
--- seismo/2D/SPECFEM2D/trunk/EXAMPLES/Abel_Brest/Par_file_Abel_Balanche_bathy_source_solid 2012-03-18 19:59:34 UTC (rev 19802)
+++ seismo/2D/SPECFEM2D/trunk/EXAMPLES/Abel_Brest/Par_file_Abel_Balanche_bathy_source_solid 2012-03-19 01:25:41 UTC (rev 19803)
@@ -9,7 +9,7 @@
# parameters concerning partitioning
nproc = 1 # number of processes
partitioning_method = 3 # SCOTCH = 3, ascending order (very bad idea) = 1
-PERFORM_CUTHILL_MCKEE = .true. # perform inverse Cuthill-McKee (1969) optimization/permutation for mesh numbering
+PERFORM_CUTHILL_MCKEE = .false. # perform inverse Cuthill-McKee (1969) optimization/permutation for mesh numbering (can be very costly and not very useful)
ngnod = 9 # number of control nodes per element (4 or 9)
initialfield = .false. # use a plane wave as source or not
Modified: seismo/2D/SPECFEM2D/trunk/EXAMPLES/Gmsh_example_MPI/Par_file_Gmsh_SqrCircles.in
===================================================================
--- seismo/2D/SPECFEM2D/trunk/EXAMPLES/Gmsh_example_MPI/Par_file_Gmsh_SqrCircles.in 2012-03-18 19:59:34 UTC (rev 19802)
+++ seismo/2D/SPECFEM2D/trunk/EXAMPLES/Gmsh_example_MPI/Par_file_Gmsh_SqrCircles.in 2012-03-19 01:25:41 UTC (rev 19803)
@@ -9,7 +9,7 @@
# parameters concerning partitioning
nproc = 4 # number of processes
partitioning_method = 3 # ascending order = 1, Scotch = 3
-PERFORM_CUTHILL_MCKEE = .true. # perform inverse Cuthill-McKee (1969) optimization/permutation for mesh numbering
+PERFORM_CUTHILL_MCKEE = .false. # perform inverse Cuthill-McKee (1969) optimization/permutation for mesh numbering (can be very costly and not very useful)
ngnod = 9 # number of control nodes per element (4 or 9)
Modified: seismo/2D/SPECFEM2D/trunk/EXAMPLES/Gmsh_example_serial/Par_file_Gmsh_SqrCircles.in
===================================================================
--- seismo/2D/SPECFEM2D/trunk/EXAMPLES/Gmsh_example_serial/Par_file_Gmsh_SqrCircles.in 2012-03-18 19:59:34 UTC (rev 19802)
+++ seismo/2D/SPECFEM2D/trunk/EXAMPLES/Gmsh_example_serial/Par_file_Gmsh_SqrCircles.in 2012-03-19 01:25:41 UTC (rev 19803)
@@ -9,7 +9,7 @@
# parameters concerning partitioning
nproc = 1 # number of processes
partitioning_method = 3 # ascending order = 1, Scotch = 3
-PERFORM_CUTHILL_MCKEE = .true. # perform inverse Cuthill-McKee (1969) optimization/permutation for mesh numbering
+PERFORM_CUTHILL_MCKEE = .false. # perform inverse Cuthill-McKee (1969) optimization/permutation for mesh numbering (can be very costly and not very useful)
ngnod = 9 # number of control nodes per element (4 or 9)
Modified: seismo/2D/SPECFEM2D/trunk/EXAMPLES/INDUSTRIAL_FORMAT/Par_file
===================================================================
--- seismo/2D/SPECFEM2D/trunk/EXAMPLES/INDUSTRIAL_FORMAT/Par_file 2012-03-18 19:59:34 UTC (rev 19802)
+++ seismo/2D/SPECFEM2D/trunk/EXAMPLES/INDUSTRIAL_FORMAT/Par_file 2012-03-19 01:25:41 UTC (rev 19803)
@@ -9,7 +9,7 @@
# parameters concerning partitioning
nproc = 1 # number of processes
partitioning_method = 3 # SCOTCH = 3, ascending order (very bad idea) = 1
-PERFORM_CUTHILL_MCKEE = .true. # perform inverse Cuthill-McKee (1969) optimization/permutation for mesh numbering
+PERFORM_CUTHILL_MCKEE = .false. # perform inverse Cuthill-McKee (1969) optimization/permutation for mesh numbering (can be very costly and not very useful)
ngnod = 9 # number of control nodes per element (4 or 9)
initialfield = .false. # use a plane wave as source or not
Modified: seismo/2D/SPECFEM2D/trunk/EXAMPLES/M2_UPPA/Par_file_M2_UPPA
===================================================================
--- seismo/2D/SPECFEM2D/trunk/EXAMPLES/M2_UPPA/Par_file_M2_UPPA 2012-03-18 19:59:34 UTC (rev 19802)
+++ seismo/2D/SPECFEM2D/trunk/EXAMPLES/M2_UPPA/Par_file_M2_UPPA 2012-03-19 01:25:41 UTC (rev 19803)
@@ -9,7 +9,7 @@
# parameters concerning partitioning
nproc = 1 # number of processes
partitioning_method = 3 # SCOTCH = 3, ascending order (very bad idea) = 1
-PERFORM_CUTHILL_MCKEE = .true. # perform inverse Cuthill-McKee (1969) optimization/permutation for mesh numbering
+PERFORM_CUTHILL_MCKEE = .false. # perform inverse Cuthill-McKee (1969) optimization/permutation for mesh numbering (can be very costly and not very useful)
ngnod = 9 # number of control nodes per element (4 or 9)
initialfield = .false. # use a plane wave as source or not
Modified: seismo/2D/SPECFEM2D/trunk/EXAMPLES/Rayleigh_wave_no_crack/Par_file_Rayleigh_2D
===================================================================
--- seismo/2D/SPECFEM2D/trunk/EXAMPLES/Rayleigh_wave_no_crack/Par_file_Rayleigh_2D 2012-03-18 19:59:34 UTC (rev 19802)
+++ seismo/2D/SPECFEM2D/trunk/EXAMPLES/Rayleigh_wave_no_crack/Par_file_Rayleigh_2D 2012-03-19 01:25:41 UTC (rev 19803)
@@ -9,7 +9,7 @@
# parameters concerning partitioning
nproc = 1 # number of processes
partitioning_method = 3 # SCOTCH = 3, ascending order (very bad idea) = 1
-PERFORM_CUTHILL_MCKEE = .true. # perform inverse Cuthill-McKee (1969) optimization/permutation for mesh numbering
+PERFORM_CUTHILL_MCKEE = .false. # perform inverse Cuthill-McKee (1969) optimization/permutation for mesh numbering (can be very costly and not very useful)
ngnod = 4 # number of control nodes per element (4 or 9)
initialfield = .true. # use a plane wave as source or not
Modified: seismo/2D/SPECFEM2D/trunk/EXAMPLES/Tape2007/Par_file_Tape2007_132rec_checker
===================================================================
--- seismo/2D/SPECFEM2D/trunk/EXAMPLES/Tape2007/Par_file_Tape2007_132rec_checker 2012-03-18 19:59:34 UTC (rev 19802)
+++ seismo/2D/SPECFEM2D/trunk/EXAMPLES/Tape2007/Par_file_Tape2007_132rec_checker 2012-03-19 01:25:41 UTC (rev 19803)
@@ -9,7 +9,7 @@
# parameters concerning partitioning
nproc = 1 # number of processes
partitioning_method = 3 # SCOTCH = 3, ascending order (very bad idea) = 1
-PERFORM_CUTHILL_MCKEE = .true. # perform inverse Cuthill-McKee (1969) optimization/permutation for mesh numbering
+PERFORM_CUTHILL_MCKEE = .false. # perform inverse Cuthill-McKee (1969) optimization/permutation for mesh numbering (can be very costly and not very useful)
ngnod = 4 # number of control nodes per element (4 or 9)
initialfield = .false. # use a plane wave as source or not
Modified: seismo/2D/SPECFEM2D/trunk/EXAMPLES/Tape2007/Par_file_Tape2007_onerec
===================================================================
--- seismo/2D/SPECFEM2D/trunk/EXAMPLES/Tape2007/Par_file_Tape2007_onerec 2012-03-18 19:59:34 UTC (rev 19802)
+++ seismo/2D/SPECFEM2D/trunk/EXAMPLES/Tape2007/Par_file_Tape2007_onerec 2012-03-19 01:25:41 UTC (rev 19803)
@@ -9,7 +9,7 @@
# parameters concerning partitioning
nproc = 1 # number of processes
partitioning_method = 3 # SCOTCH = 3, ascending order (very bad idea) = 1
-PERFORM_CUTHILL_MCKEE = .true. # perform inverse Cuthill-McKee (1969) optimization/permutation for mesh numbering
+PERFORM_CUTHILL_MCKEE = .false. # perform inverse Cuthill-McKee (1969) optimization/permutation for mesh numbering (can be very costly and not very useful)
ngnod = 4 # number of control nodes per element (4 or 9)
initialfield = .false. # use a plane wave as source or not
Modified: seismo/2D/SPECFEM2D/trunk/EXAMPLES/Tape2007_kernel/Par_file_Tape2007_onerec
===================================================================
--- seismo/2D/SPECFEM2D/trunk/EXAMPLES/Tape2007_kernel/Par_file_Tape2007_onerec 2012-03-18 19:59:34 UTC (rev 19802)
+++ seismo/2D/SPECFEM2D/trunk/EXAMPLES/Tape2007_kernel/Par_file_Tape2007_onerec 2012-03-19 01:25:41 UTC (rev 19803)
@@ -9,7 +9,7 @@
# parameters concerning partitioning
nproc = 1 # number of processes
partitioning_method = 3 # SCOTCH = 3, ascending order (very bad idea) = 1
-PERFORM_CUTHILL_MCKEE = .true. # perform inverse Cuthill-McKee (1969) optimization/permutation for mesh numbering
+PERFORM_CUTHILL_MCKEE = .false. # perform inverse Cuthill-McKee (1969) optimization/permutation for mesh numbering (can be very costly and not very useful)
ngnod = 4 # number of control nodes per element (4 or 9)
initialfield = .false. # use a plane wave as source or not
Modified: seismo/2D/SPECFEM2D/trunk/EXAMPLES/Tromp2005/Par_file_Tromp2005
===================================================================
--- seismo/2D/SPECFEM2D/trunk/EXAMPLES/Tromp2005/Par_file_Tromp2005 2012-03-18 19:59:34 UTC (rev 19802)
+++ seismo/2D/SPECFEM2D/trunk/EXAMPLES/Tromp2005/Par_file_Tromp2005 2012-03-19 01:25:41 UTC (rev 19803)
@@ -9,7 +9,7 @@
# parameters concerning partitioning
nproc = 1 # number of processes
partitioning_method = 3 # SCOTCH = 3, ascending order (very bad idea) = 1
-PERFORM_CUTHILL_MCKEE = .true. # perform inverse Cuthill-McKee (1969) optimization/permutation for mesh numbering
+PERFORM_CUTHILL_MCKEE = .false. # perform inverse Cuthill-McKee (1969) optimization/permutation for mesh numbering (can be very costly and not very useful)
ngnod = 4 # number of control nodes per element (4 or 9)
initialfield = .false. # use a plane wave as source or not
Modified: seismo/2D/SPECFEM2D/trunk/EXAMPLES/Tromp2005/Par_file_Tromp2005_s100
===================================================================
--- seismo/2D/SPECFEM2D/trunk/EXAMPLES/Tromp2005/Par_file_Tromp2005_s100 2012-03-18 19:59:34 UTC (rev 19802)
+++ seismo/2D/SPECFEM2D/trunk/EXAMPLES/Tromp2005/Par_file_Tromp2005_s100 2012-03-19 01:25:41 UTC (rev 19803)
@@ -9,7 +9,7 @@
# parameters concerning partitioning
nproc = 1 # number of processes
partitioning_method = 3 # SCOTCH = 3, ascending order (very bad idea) = 1
-PERFORM_CUTHILL_MCKEE = .true. # perform inverse Cuthill-McKee (1969) optimization/permutation for mesh numbering
+PERFORM_CUTHILL_MCKEE = .false. # perform inverse Cuthill-McKee (1969) optimization/permutation for mesh numbering (can be very costly and not very useful)
ngnod = 4 # number of control nodes per element (4 or 9)
initialfield = .false. # use a plane wave as source or not
Modified: seismo/2D/SPECFEM2D/trunk/EXAMPLES/Tromp2005_kernel/Par_file_Tromp2005
===================================================================
--- seismo/2D/SPECFEM2D/trunk/EXAMPLES/Tromp2005_kernel/Par_file_Tromp2005 2012-03-18 19:59:34 UTC (rev 19802)
+++ seismo/2D/SPECFEM2D/trunk/EXAMPLES/Tromp2005_kernel/Par_file_Tromp2005 2012-03-19 01:25:41 UTC (rev 19803)
@@ -9,7 +9,7 @@
# parameters concerning partitioning
nproc = 1 # number of processes
partitioning_method = 3 # SCOTCH = 3, ascending order (very bad idea) = 1
-PERFORM_CUTHILL_MCKEE = .true. # perform inverse Cuthill-McKee (1969) optimization/permutation for mesh numbering
+PERFORM_CUTHILL_MCKEE = .false. # perform inverse Cuthill-McKee (1969) optimization/permutation for mesh numbering (can be very costly and not very useful)
ngnod = 4 # number of control nodes per element (4 or 9)
initialfield = .false. # use a plane wave as source or not
Modified: seismo/2D/SPECFEM2D/trunk/EXAMPLES/acoustic_poroelastic/Par_file_acoustic_poroelastic
===================================================================
--- seismo/2D/SPECFEM2D/trunk/EXAMPLES/acoustic_poroelastic/Par_file_acoustic_poroelastic 2012-03-18 19:59:34 UTC (rev 19802)
+++ seismo/2D/SPECFEM2D/trunk/EXAMPLES/acoustic_poroelastic/Par_file_acoustic_poroelastic 2012-03-19 01:25:41 UTC (rev 19803)
@@ -9,7 +9,7 @@
# parameters concerning partitioning
nproc = 1 # number of processes
partitioning_method = 3 # SCOTCH = 3, ascending order (very bad idea) = 1
-PERFORM_CUTHILL_MCKEE = .true. # perform inverse Cuthill-McKee (1969) optimization/permutation for mesh numbering
+PERFORM_CUTHILL_MCKEE = .false. # perform inverse Cuthill-McKee (1969) optimization/permutation for mesh numbering (can be very costly and not very useful)
ngnod = 9 # number of control nodes per element (4 or 9)
initialfield = .false. # use a plane wave as source or not
Modified: seismo/2D/SPECFEM2D/trunk/EXAMPLES/attenuation/Par_file_attenuation_2D
===================================================================
--- seismo/2D/SPECFEM2D/trunk/EXAMPLES/attenuation/Par_file_attenuation_2D 2012-03-18 19:59:34 UTC (rev 19802)
+++ seismo/2D/SPECFEM2D/trunk/EXAMPLES/attenuation/Par_file_attenuation_2D 2012-03-19 01:25:41 UTC (rev 19803)
@@ -9,7 +9,7 @@
# parameters concerning partitioning
nproc = 1 # number of processes
partitioning_method = 3 # SCOTCH = 3, ascending order (very bad idea) = 1
-PERFORM_CUTHILL_MCKEE = .true. # perform inverse Cuthill-McKee (1969) optimization/permutation for mesh numbering
+PERFORM_CUTHILL_MCKEE = .false. # perform inverse Cuthill-McKee (1969) optimization/permutation for mesh numbering (can be very costly and not very useful)
ngnod = 4 # number of control nodes per element (4 or 9)
initialfield = .false. # use a plane wave as source or not
Modified: seismo/2D/SPECFEM2D/trunk/EXAMPLES/canyon/Par_file_canyon
===================================================================
--- seismo/2D/SPECFEM2D/trunk/EXAMPLES/canyon/Par_file_canyon 2012-03-18 19:59:34 UTC (rev 19802)
+++ seismo/2D/SPECFEM2D/trunk/EXAMPLES/canyon/Par_file_canyon 2012-03-19 01:25:41 UTC (rev 19803)
@@ -9,7 +9,7 @@
# parameters concerning partitioning
nproc = 1 # number of processes
partitioning_method = 3 # SCOTCH = 3, ascending order (very bad idea) = 1
-PERFORM_CUTHILL_MCKEE = .true. # perform inverse Cuthill-McKee (1969) optimization/permutation for mesh numbering
+PERFORM_CUTHILL_MCKEE = .false. # perform inverse Cuthill-McKee (1969) optimization/permutation for mesh numbering (can be very costly and not very useful)
ngnod = 4 # number of control nodes per element (4 or 9)
initialfield = .true. # use a plane wave as source or not
Added: seismo/2D/SPECFEM2D/trunk/EXAMPLES/change_all_2D_Par_files.csh
===================================================================
--- seismo/2D/SPECFEM2D/trunk/EXAMPLES/change_all_2D_Par_files.csh (rev 0)
+++ seismo/2D/SPECFEM2D/trunk/EXAMPLES/change_all_2D_Par_files.csh 2012-03-19 01:25:41 UTC (rev 19803)
@@ -0,0 +1,21 @@
+#!/bin/csh
+
+# I get the list of files using
+# alias list_all_2D_Par_files_to_change " ls DATA/Par_file EXAMPLES/*/Par_fi* EXAMPLES/*/*/Par_fi* "
+#
+# and then I manually remove EXAMPLES/unused_older_examples_DATA_to_sort , which is by definition in an obsolete/incompatible format
+#
+
+foreach file ( DATA/Par_file EXAMPLES/Abel_Brest/Par_file_Abel_Balanche_bathy_source_solid EXAMPLES/acoustic_poroelastic/Par_file_acoustic_poroelastic EXAMPLES/attenuation/Par_file_attenuation_2D EXAMPLES/canyon/Par_file_canyon EXAMPLES/fluid_solid/fluid_solid_external_mesh/Par_file_fluid_solid EXAMPLES/fluid_solid/from_2000_Geophysics_paper_flat_ocean_bottom/Par_file_fluid_solid EXAMPLES/fluid_solid/from_2000_Geophysics_paper_sinusoidal_ocean_bottom/Par_file_fluid_solid EXAMPLES/Gmsh_example_MPI/Par_file_Gmsh_SqrCircles.in EXAMPLES/Gmsh_example_serial/Par_file_Gmsh_SqrCircles.in EXAMPLES/INDUSTRIAL_FORMAT/Par_file EXAMPLES/init_plane/Par_file_Slave EXAMPLES/init_plane/Par_file_Slave_for EXAMPLES/init_plane/Par_file_Slave_kernel EXAMPLES/M2_UPPA/Par_file_M2_UPPA EXAMPLES/noise_layered/model_0/Par_file_fair EXAMPLES/noise_layered/model_0/Par_file_good EXAMPLES/noise_layered/model_1/Par_file_best EXAMPLES/noise_layered/model_1/Par_file_fair EXAMPLES/noise_layered/model_1/Par_file_good EXAMPLES/noise_layered/model_2/Par_file_fair EXAMPLES/noise_layered/model_2/Par_file_good EXAMPLES/noise_uniform/Par_file_noise_1 EXAMPLES/noise_uniform/Par_file_noise_2 EXAMPLES/noise_uniform/Par_file_noise_3 EXAMPLES/Rayleigh_wave_no_crack/Par_file_Rayleigh_2D EXAMPLES/Rayleigh_wave_with_crack/Par_file_Rayleigh_2D EXAMPLES/semi_infinite_homo/Par_file_elastic_2D EXAMPLES/Tape2007_kernel/Par_file_Tape2007_onerec EXAMPLES/Tape2007/Par_file_Tape2007_132rec_checker EXAMPLES/Tape2007/Par_file_Tape2007_onerec EXAMPLES/Tromp2005_kernel/Par_file_Tromp2005 EXAMPLES/Tromp2005/Par_file_Tromp2005 EXAMPLES/Tromp2005/Par_file_Tromp2005_s100 )
+
+set newfile = "../"$file
+
+# DK DK leave the white spaces, in order to have the right alignment with other variables
+sed -e '1,$s/PERFORM_CUTHILL_MCKEE = .true. # perform inverse Cuthill-McKee (1969) optimization\/permutation for mesh numbering/PERFORM_CUTHILL_MCKEE = .false. # perform inverse Cuthill-McKee (1969) optimization\/permutation for mesh numbering (can be very costly and not very useful)/g' < $newfile > __________temp_27_zzzyyyyxxxx__________
+
+clean_a_Par_file.py __________temp_27_zzzyyyyxxxx__________ $file
+
+mv __________temp_27_zzzyyyyxxxx__________ $newfile
+
+end
+
Property changes on: seismo/2D/SPECFEM2D/trunk/EXAMPLES/change_all_2D_Par_files.csh
___________________________________________________________________
Name: svn:executable
+ *
Added: seismo/2D/SPECFEM2D/trunk/EXAMPLES/clean_a_Par_file.py
===================================================================
--- seismo/2D/SPECFEM2D/trunk/EXAMPLES/clean_a_Par_file.py (rev 0)
+++ seismo/2D/SPECFEM2D/trunk/EXAMPLES/clean_a_Par_file.py 2012-03-19 01:25:41 UTC (rev 19803)
@@ -0,0 +1,35 @@
+#!/usr/bin/perl
+
+#
+# Script to clean a given file
+#
+# Author : Dimitri Komatitsch, EPS - Harvard University, May 1998
+#
+
+ $name = $ARGV[0];
+ $nametoprint = $ARGV[1];
+
+# change tabs to white spaces
+ system("expand -6 < $name > _____testzzzXXXyyy_____");
+ $f90name = $name;
+ print STDOUT "Cleaning file $nametoprint ...\n";
+
+ open(FILEF77,"<_____testzzzXXXyyy_____");
+ open(FILEF90,">$f90name");
+
+# read the input file
+ while($line = <FILEF77>) {
+
+# suppress trailing white spaces and carriage return
+ $line =~ s/\s*$//;
+
+# write the output line
+ print FILEF90 "$line\n";
+
+ }
+
+ close(FILEF77);
+ close(FILEF90);
+
+ system("rm -f _____testzzzXXXyyy_____");
+
Property changes on: seismo/2D/SPECFEM2D/trunk/EXAMPLES/clean_a_Par_file.py
___________________________________________________________________
Name: svn:executable
+ *
Modified: seismo/2D/SPECFEM2D/trunk/EXAMPLES/fluid_solid/fluid_solid_external_mesh/Par_file_fluid_solid
===================================================================
--- seismo/2D/SPECFEM2D/trunk/EXAMPLES/fluid_solid/fluid_solid_external_mesh/Par_file_fluid_solid 2012-03-18 19:59:34 UTC (rev 19802)
+++ seismo/2D/SPECFEM2D/trunk/EXAMPLES/fluid_solid/fluid_solid_external_mesh/Par_file_fluid_solid 2012-03-19 01:25:41 UTC (rev 19803)
@@ -9,7 +9,7 @@
# parameters concerning partitioning
nproc = 1 # number of processes
partitioning_method = 3 # SCOTCH = 3, ascending order (very bad idea) = 1
-PERFORM_CUTHILL_MCKEE = .true. # perform inverse Cuthill-McKee (1969) optimization/permutation for mesh numbering
+PERFORM_CUTHILL_MCKEE = .false. # perform inverse Cuthill-McKee (1969) optimization/permutation for mesh numbering (can be very costly and not very useful)
ngnod = 4 # number of control nodes per element (4 or 9)
initialfield = .false. # use a plane wave as source or not
Modified: seismo/2D/SPECFEM2D/trunk/EXAMPLES/fluid_solid/from_2000_Geophysics_paper_flat_ocean_bottom/Par_file_fluid_solid
===================================================================
--- seismo/2D/SPECFEM2D/trunk/EXAMPLES/fluid_solid/from_2000_Geophysics_paper_flat_ocean_bottom/Par_file_fluid_solid 2012-03-18 19:59:34 UTC (rev 19802)
+++ seismo/2D/SPECFEM2D/trunk/EXAMPLES/fluid_solid/from_2000_Geophysics_paper_flat_ocean_bottom/Par_file_fluid_solid 2012-03-19 01:25:41 UTC (rev 19803)
@@ -9,7 +9,7 @@
# parameters concerning partitioning
nproc = 1 # number of processes
partitioning_method = 3 # SCOTCH = 3, ascending order (very bad idea) = 1
-PERFORM_CUTHILL_MCKEE = .true. # perform inverse Cuthill-McKee (1969) optimization/permutation for mesh numbering
+PERFORM_CUTHILL_MCKEE = .false. # perform inverse Cuthill-McKee (1969) optimization/permutation for mesh numbering (can be very costly and not very useful)
ngnod = 4 # number of control nodes per element (4 or 9)
initialfield = .false. # use a plane wave as source or not
Modified: seismo/2D/SPECFEM2D/trunk/EXAMPLES/fluid_solid/from_2000_Geophysics_paper_sinusoidal_ocean_bottom/Par_file_fluid_solid
===================================================================
--- seismo/2D/SPECFEM2D/trunk/EXAMPLES/fluid_solid/from_2000_Geophysics_paper_sinusoidal_ocean_bottom/Par_file_fluid_solid 2012-03-18 19:59:34 UTC (rev 19802)
+++ seismo/2D/SPECFEM2D/trunk/EXAMPLES/fluid_solid/from_2000_Geophysics_paper_sinusoidal_ocean_bottom/Par_file_fluid_solid 2012-03-19 01:25:41 UTC (rev 19803)
@@ -9,7 +9,7 @@
# parameters concerning partitioning
nproc = 1 # number of processes
partitioning_method = 3 # SCOTCH = 3, ascending order (very bad idea) = 1
-PERFORM_CUTHILL_MCKEE = .true. # perform inverse Cuthill-McKee (1969) optimization/permutation for mesh numbering
+PERFORM_CUTHILL_MCKEE = .false. # perform inverse Cuthill-McKee (1969) optimization/permutation for mesh numbering (can be very costly and not very useful)
ngnod = 9 # number of control nodes per element (4 or 9)
initialfield = .false. # use a plane wave as source or not
Modified: seismo/2D/SPECFEM2D/trunk/EXAMPLES/init_plane/Par_file_Slave
===================================================================
--- seismo/2D/SPECFEM2D/trunk/EXAMPLES/init_plane/Par_file_Slave 2012-03-18 19:59:34 UTC (rev 19802)
+++ seismo/2D/SPECFEM2D/trunk/EXAMPLES/init_plane/Par_file_Slave 2012-03-19 01:25:41 UTC (rev 19803)
@@ -9,7 +9,7 @@
# parameters concerning partitioning
nproc = 1 # number of processes
partitioning_method = 3 # SCOTCH = 3, ascending order (very bad idea) = 1
-PERFORM_CUTHILL_MCKEE = .true. # perform inverse Cuthill-McKee (1969) optimization/permutation for mesh numbering
+PERFORM_CUTHILL_MCKEE = .false. # perform inverse Cuthill-McKee (1969) optimization/permutation for mesh numbering (can be very costly and not very useful)
ngnod = 4 # number of control nodes per element (4 or 9)
initialfield = .true. # use a plane wave as source or not
Modified: seismo/2D/SPECFEM2D/trunk/EXAMPLES/init_plane/Par_file_Slave_for
===================================================================
--- seismo/2D/SPECFEM2D/trunk/EXAMPLES/init_plane/Par_file_Slave_for 2012-03-18 19:59:34 UTC (rev 19802)
+++ seismo/2D/SPECFEM2D/trunk/EXAMPLES/init_plane/Par_file_Slave_for 2012-03-19 01:25:41 UTC (rev 19803)
@@ -9,7 +9,7 @@
# parameters concerning partitioning
nproc = 1 # number of processes
partitioning_method = 3 # SCOTCH = 3, ascending order (very bad idea) = 1
-PERFORM_CUTHILL_MCKEE = .true. # perform inverse Cuthill-McKee (1969) optimization/permutation for mesh numbering
+PERFORM_CUTHILL_MCKEE = .false. # perform inverse Cuthill-McKee (1969) optimization/permutation for mesh numbering (can be very costly and not very useful)
ngnod = 4 # number of control nodes per element (4 or 9)
initialfield = .true. # use a plane wave as source or not
Modified: seismo/2D/SPECFEM2D/trunk/EXAMPLES/init_plane/Par_file_Slave_kernel
===================================================================
--- seismo/2D/SPECFEM2D/trunk/EXAMPLES/init_plane/Par_file_Slave_kernel 2012-03-18 19:59:34 UTC (rev 19802)
+++ seismo/2D/SPECFEM2D/trunk/EXAMPLES/init_plane/Par_file_Slave_kernel 2012-03-19 01:25:41 UTC (rev 19803)
@@ -9,7 +9,7 @@
# parameters concerning partitioning
nproc = 1 # number of processes
partitioning_method = 3 # SCOTCH = 3, ascending order (very bad idea) = 1
-PERFORM_CUTHILL_MCKEE = .true. # perform inverse Cuthill-McKee (1969) optimization/permutation for mesh numbering
+PERFORM_CUTHILL_MCKEE = .false. # perform inverse Cuthill-McKee (1969) optimization/permutation for mesh numbering (can be very costly and not very useful)
ngnod = 4 # number of control nodes per element (4 or 9)
initialfield = .false. # use a plane wave as source or not
Modified: seismo/2D/SPECFEM2D/trunk/EXAMPLES/noise_layered/model_0/Par_file_fair
===================================================================
--- seismo/2D/SPECFEM2D/trunk/EXAMPLES/noise_layered/model_0/Par_file_fair 2012-03-18 19:59:34 UTC (rev 19802)
+++ seismo/2D/SPECFEM2D/trunk/EXAMPLES/noise_layered/model_0/Par_file_fair 2012-03-19 01:25:41 UTC (rev 19803)
@@ -9,7 +9,7 @@
# parameters concerning partitioning
nproc = 1 # number of processes
partitioning_method = 3 # SCOTCH = 3, ascending order (very bad idea) = 1
-PERFORM_CUTHILL_MCKEE = .true. # perform inverse Cuthill-McKee (1969) optimization/permutation for mesh numbering
+PERFORM_CUTHILL_MCKEE = .false. # perform inverse Cuthill-McKee (1969) optimization/permutation for mesh numbering (can be very costly and not very useful)
ngnod = 9 # number of control nodes per element (4 or 9)
initialfield = .false. # use a plane wave as source or not
Modified: seismo/2D/SPECFEM2D/trunk/EXAMPLES/noise_layered/model_0/Par_file_good
===================================================================
--- seismo/2D/SPECFEM2D/trunk/EXAMPLES/noise_layered/model_0/Par_file_good 2012-03-18 19:59:34 UTC (rev 19802)
+++ seismo/2D/SPECFEM2D/trunk/EXAMPLES/noise_layered/model_0/Par_file_good 2012-03-19 01:25:41 UTC (rev 19803)
@@ -9,7 +9,7 @@
# parameters concerning partitioning
nproc = 1 # number of processes
partitioning_method = 3 # SCOTCH = 3, ascending order (very bad idea) = 1
-PERFORM_CUTHILL_MCKEE = .true. # perform inverse Cuthill-McKee (1969) optimization/permutation for mesh numbering
+PERFORM_CUTHILL_MCKEE = .false. # perform inverse Cuthill-McKee (1969) optimization/permutation for mesh numbering (can be very costly and not very useful)
ngnod = 9 # number of control nodes per element (4 or 9)
initialfield = .false. # use a plane wave as source or not
Modified: seismo/2D/SPECFEM2D/trunk/EXAMPLES/noise_layered/model_1/Par_file_best
===================================================================
--- seismo/2D/SPECFEM2D/trunk/EXAMPLES/noise_layered/model_1/Par_file_best 2012-03-18 19:59:34 UTC (rev 19802)
+++ seismo/2D/SPECFEM2D/trunk/EXAMPLES/noise_layered/model_1/Par_file_best 2012-03-19 01:25:41 UTC (rev 19803)
@@ -9,7 +9,7 @@
# parameters concerning partitioning
nproc = 1 # number of processes
partitioning_method = 3 # SCOTCH = 3, ascending order (very bad idea) = 1
-PERFORM_CUTHILL_MCKEE = .true. # perform inverse Cuthill-McKee (1969) optimization/permutation for mesh numbering
+PERFORM_CUTHILL_MCKEE = .false. # perform inverse Cuthill-McKee (1969) optimization/permutation for mesh numbering (can be very costly and not very useful)
ngnod = 9 # number of control nodes per element (4 or 9)
initialfield = .false. # use a plane wave as source or not
Modified: seismo/2D/SPECFEM2D/trunk/EXAMPLES/noise_layered/model_1/Par_file_fair
===================================================================
--- seismo/2D/SPECFEM2D/trunk/EXAMPLES/noise_layered/model_1/Par_file_fair 2012-03-18 19:59:34 UTC (rev 19802)
+++ seismo/2D/SPECFEM2D/trunk/EXAMPLES/noise_layered/model_1/Par_file_fair 2012-03-19 01:25:41 UTC (rev 19803)
@@ -9,7 +9,7 @@
# parameters concerning partitioning
nproc = 1 # number of processes
partitioning_method = 3 # SCOTCH = 3, ascending order (very bad idea) = 1
-PERFORM_CUTHILL_MCKEE = .true. # perform inverse Cuthill-McKee (1969) optimization/permutation for mesh numbering
+PERFORM_CUTHILL_MCKEE = .false. # perform inverse Cuthill-McKee (1969) optimization/permutation for mesh numbering (can be very costly and not very useful)
ngnod = 9 # number of control nodes per element (4 or 9)
initialfield = .false. # use a plane wave as source or not
Modified: seismo/2D/SPECFEM2D/trunk/EXAMPLES/noise_layered/model_1/Par_file_good
===================================================================
--- seismo/2D/SPECFEM2D/trunk/EXAMPLES/noise_layered/model_1/Par_file_good 2012-03-18 19:59:34 UTC (rev 19802)
+++ seismo/2D/SPECFEM2D/trunk/EXAMPLES/noise_layered/model_1/Par_file_good 2012-03-19 01:25:41 UTC (rev 19803)
@@ -9,7 +9,7 @@
# parameters concerning partitioning
nproc = 1 # number of processes
partitioning_method = 3 # SCOTCH = 3, ascending order (very bad idea) = 1
-PERFORM_CUTHILL_MCKEE = .true. # perform inverse Cuthill-McKee (1969) optimization/permutation for mesh numbering
+PERFORM_CUTHILL_MCKEE = .false. # perform inverse Cuthill-McKee (1969) optimization/permutation for mesh numbering (can be very costly and not very useful)
ngnod = 9 # number of control nodes per element (4 or 9)
initialfield = .false. # use a plane wave as source or not
Modified: seismo/2D/SPECFEM2D/trunk/EXAMPLES/noise_layered/model_2/Par_file_fair
===================================================================
--- seismo/2D/SPECFEM2D/trunk/EXAMPLES/noise_layered/model_2/Par_file_fair 2012-03-18 19:59:34 UTC (rev 19802)
+++ seismo/2D/SPECFEM2D/trunk/EXAMPLES/noise_layered/model_2/Par_file_fair 2012-03-19 01:25:41 UTC (rev 19803)
@@ -9,7 +9,7 @@
# parameters concerning partitioning
nproc = 1 # number of processes
partitioning_method = 3 # SCOTCH = 3, ascending order (very bad idea) = 1
-PERFORM_CUTHILL_MCKEE = .true. # perform inverse Cuthill-McKee (1969) optimization/permutation for mesh numbering
+PERFORM_CUTHILL_MCKEE = .false. # perform inverse Cuthill-McKee (1969) optimization/permutation for mesh numbering (can be very costly and not very useful)
ngnod = 9 # number of control nodes per element (4 or 9)
initialfield = .false. # use a plane wave as source or not
Modified: seismo/2D/SPECFEM2D/trunk/EXAMPLES/noise_layered/model_2/Par_file_good
===================================================================
--- seismo/2D/SPECFEM2D/trunk/EXAMPLES/noise_layered/model_2/Par_file_good 2012-03-18 19:59:34 UTC (rev 19802)
+++ seismo/2D/SPECFEM2D/trunk/EXAMPLES/noise_layered/model_2/Par_file_good 2012-03-19 01:25:41 UTC (rev 19803)
@@ -9,7 +9,7 @@
# parameters concerning partitioning
nproc = 1 # number of processes
partitioning_method = 3 # SCOTCH = 3, ascending order (very bad idea) = 1
-PERFORM_CUTHILL_MCKEE = .true. # perform inverse Cuthill-McKee (1969) optimization/permutation for mesh numbering
+PERFORM_CUTHILL_MCKEE = .false. # perform inverse Cuthill-McKee (1969) optimization/permutation for mesh numbering (can be very costly and not very useful)
ngnod = 9 # number of control nodes per element (4 or 9)
initialfield = .false. # use a plane wave as source or not
Modified: seismo/2D/SPECFEM2D/trunk/EXAMPLES/noise_uniform/Par_file_noise_1
===================================================================
--- seismo/2D/SPECFEM2D/trunk/EXAMPLES/noise_uniform/Par_file_noise_1 2012-03-18 19:59:34 UTC (rev 19802)
+++ seismo/2D/SPECFEM2D/trunk/EXAMPLES/noise_uniform/Par_file_noise_1 2012-03-19 01:25:41 UTC (rev 19803)
@@ -9,7 +9,7 @@
# parameters concerning partitioning
nproc = 1 # number of processes
partitioning_method = 3 # SCOTCH = 3, ascending order (very bad idea) = 1
-PERFORM_CUTHILL_MCKEE = .true. # perform inverse Cuthill-McKee (1969) optimization/permutation for mesh numbering
+PERFORM_CUTHILL_MCKEE = .false. # perform inverse Cuthill-McKee (1969) optimization/permutation for mesh numbering (can be very costly and not very useful)
ngnod = 9 # number of control nodes per element (4 or 9)
initialfield = .false. # use a plane wave as source or not
Modified: seismo/2D/SPECFEM2D/trunk/EXAMPLES/noise_uniform/Par_file_noise_2
===================================================================
--- seismo/2D/SPECFEM2D/trunk/EXAMPLES/noise_uniform/Par_file_noise_2 2012-03-18 19:59:34 UTC (rev 19802)
+++ seismo/2D/SPECFEM2D/trunk/EXAMPLES/noise_uniform/Par_file_noise_2 2012-03-19 01:25:41 UTC (rev 19803)
@@ -9,7 +9,7 @@
# parameters concerning partitioning
nproc = 1 # number of processes
partitioning_method = 3 # SCOTCH = 3, ascending order (very bad idea) = 1
-PERFORM_CUTHILL_MCKEE = .true. # perform inverse Cuthill-McKee (1969) optimization/permutation for mesh numbering
+PERFORM_CUTHILL_MCKEE = .false. # perform inverse Cuthill-McKee (1969) optimization/permutation for mesh numbering (can be very costly and not very useful)
ngnod = 9 # number of control nodes per element (4 or 9)
initialfield = .false. # use a plane wave as source or not
Modified: seismo/2D/SPECFEM2D/trunk/EXAMPLES/noise_uniform/Par_file_noise_3
===================================================================
--- seismo/2D/SPECFEM2D/trunk/EXAMPLES/noise_uniform/Par_file_noise_3 2012-03-18 19:59:34 UTC (rev 19802)
+++ seismo/2D/SPECFEM2D/trunk/EXAMPLES/noise_uniform/Par_file_noise_3 2012-03-19 01:25:41 UTC (rev 19803)
@@ -9,7 +9,7 @@
# parameters concerning partitioning
nproc = 1 # number of processes
partitioning_method = 3 # SCOTCH = 3, ascending order (very bad idea) = 1
-PERFORM_CUTHILL_MCKEE = .true. # perform inverse Cuthill-McKee (1969) optimization/permutation for mesh numbering
+PERFORM_CUTHILL_MCKEE = .false. # perform inverse Cuthill-McKee (1969) optimization/permutation for mesh numbering (can be very costly and not very useful)
ngnod = 9 # number of control nodes per element (4 or 9)
initialfield = .false. # use a plane wave as source or not
Copied: seismo/2D/SPECFEM2D/trunk/EXAMPLES/process_DATA_Par_files_to_update_their_format_when_new_parameters_are_added.py (from rev 19802, seismo/2D/SPECFEM2D/trunk/UTILS/process_DATA_Par_files_to_update_their_format_when_new_parameters_are_added.py)
===================================================================
--- seismo/2D/SPECFEM2D/trunk/EXAMPLES/process_DATA_Par_files_to_update_their_format_when_new_parameters_are_added.py (rev 0)
+++ seismo/2D/SPECFEM2D/trunk/EXAMPLES/process_DATA_Par_files_to_update_their_format_when_new_parameters_are_added.py 2012-03-19 01:25:41 UTC (rev 19803)
@@ -0,0 +1,191 @@
+# -*- coding: utf-8 -*-
+"""
+Created on Tue Nov 15 09:00:00 2011
+Updated on Fri Jan 13 2012
+
+Processing of Par_file to update them to new format
+
+Usage : "python PathTo/EXAMPLES/ProcessParFileParametersToNewRelease.py"
+This will process all Par_file starting from current directory
+
+ at author: Cristini Paul, Laboratoire de Mecanique et d'Acoustique, CNRS, Marseille, France
+"""
+import os, string, sys
+from os import listdir, walk
+from shutil import copy, move
+from os.path import exists
+#------------------------------------------------------------------------------
+def LoadLig(Fichier):
+ f = open(Fichier,'r')
+ ligs= f.readlines()
+ f.close()
+ return ligs
+#------------------------------------------------------------------------------
+def mylister(currdir):
+ for file in os.listdir(currdir):
+ path=os.path.join(currdir, file)
+ if not os.path.isdir(path):
+ #print path
+ Fichiers.append(path)
+ else:
+ mylister(path)
+#------------------------------------------------------------------------------
+def ProcessParfile_r19201(fic):
+ # Open the file and get all lines from Par_file
+ ligs= LoadLig(fic)
+ # Test pour voir si le traitement a déjà été fait
+ for lig in ligs:
+ if 'ADD_PERIODIC_CONDITIONS' in lig:
+ print '----> '+fic+' already processed to r19201'
+ return
+ # New additions to the Par_file
+ a1='PERFORM_CUTHILL_MCKEE = .false. # perform inverse Cuthill-McKee (1969) optimization/permutation for mesh numbering (can be very costly and not very useful)\n'
+ a2='USER_T0 = 0.0d0 # use this t0 as earliest starting time rather than the automatically calculated one\n'
+ a3='SU_FORMAT = .false. # output seismograms in Seismic Unix format (adjoint traces will be read in the same format)\n'
+ a4='factor_subsample_image = 1 # factor to subsample color images output by the code (useful for very large models)\n'+ \
+ 'POWER_DISPLAY_COLOR = 0.30d0 # non linear display to enhance small amplitudes in color images\n'+ \
+ 'DRAW_WATER_CONSTANT_BLUE_IN_JPG = .true. # display acoustic layers as constant blue in JPEG images, because they likely correspond to water\n'
+ a5='US_LETTER = .false. # US letter paper or European A4\n'+ \
+ 'USE_SNAPSHOT_NUMBER_IN_FILENAME = .false. # use snapshot number in the file name of JPEG color snapshots instead of the time step\n'
+ a6='\n# for horizontal periodic conditions: detect common points between left and right edges\n'+ \
+ 'ADD_PERIODIC_CONDITIONS = .false.\n\n'+ \
+ '# horizontal periodicity distance for periodic conditions\n'+ \
+ 'PERIODIC_horiz_dist = 0.3597d0\n\n'+ \
+ '# grid point detection tolerance for periodic conditions\n'+ \
+ 'PERIODIC_DETECT_TOL = 3.3334d-6\n'
+ #--------------------------------------------------------------------------
+ # Ajout des parametres supplementaires
+ #
+ for ilg, lig in enumerate(ligs):
+ if lig.startswith('partitioning'):
+ ligs.insert(ilg+1,a1)
+
+ if lig.startswith('deltat'):
+ ligs.insert(ilg+1,a2)
+
+ if lig.startswith('rec_normal'):
+ ligs.insert(ilg+1,a3)
+
+ if lig.startswith('subsamp'):
+ ligs[ilg]=string.replace(ligs[ilg],'subsamp ','subsamp_postscript',1)
+ ligs.insert(ilg+1,a4)
+
+ if lig.startswith('sizemax'):
+ ligs.insert(ilg+1,a5)
+
+ if lig.startswith('absorbing_conditions'):
+ ligs.insert(ilg+1,a6)
+ #
+ move(fic,fic+'.before_update_to_r19201')
+ fm = open(fic,'w')
+ fm.writelines(ligs)
+ fm.close()
+ #
+ print 'xxxxx------> '+fic+' processed to r19201'
+ return
+#------------------------------------------------------------------------------
+def ProcessParfile_r19340(fic):
+ # Open the file and get all lines from Par_file
+ ligs= LoadLig(fic)
+ # Teste si le traitement a déjà été fait
+ for lig in ligs:
+ if 'nreceiversets' in lig:
+ print '----> '+fic+' already processed to r19340'
+ return
+ #
+ # Ajout des parametres supplementaires
+ #
+ for ilg, lig in enumerate(ligs):
+ if lig.startswith('nreceiverlines'):
+ ligs[ilg]=ligs[ilg].replace('lines','sets ')
+ #
+ move(fic,fic+'.before_update_to_r19340')
+ #
+ fm = open(fic,'w')
+ fm.writelines(ligs)
+ fm.close()
+ #
+ print 'xxxxx------> '+fic+' processed to r19340'
+ return
+#------------------------------------------------------------------------------
+def ProcessParfile_r19346(fic):
+ # Open the file and get all lines from Par_file
+ ligs= LoadLig(fic)
+ # Teste si le traitement a déjà été fait
+ for lig in ligs:
+ if 'ATTENUATION_PORO_FLUID_PART' in lig:
+ print '----> '+fic+' already processed to r19346'
+ return
+ #--------------------------------------------------------------------------
+ # Ajout des parametres supplementaires
+ #
+ for ilg, lig in enumerate(ligs):
+ if lig.startswith('TURN_ATTENUATION_ON'):
+ ligs[ilg]=ligs[ilg].replace('TURN_ATTENUATION_ON ', \
+ 'ATTENUATION_VISCOELASTIC_SOLID')
+ if lig.startswith('TURN_VISCATTENUATION_ON'):
+ ligs[ilg]=ligs[ilg].replace('TURN_VISCATTENUATION_ON ', \
+ 'ATTENUATION_PORO_FLUID_PART')
+ #
+ move(fic,fic+'.before_update_to_r19346')
+ #
+ fm = open(fic,'w')
+ fm.writelines(ligs)
+ fm.close()
+ #
+ print 'xxxxx------> '+fic+' processed to r19346'
+ return
+#------------------------------------------------------------------------------
+def ProcessParfile_r19xxx(fic):
+ # Open the file and get all lines from Par_file
+ ligs= LoadLig(fic)
+ # Teste si le traitement a déjà été fait
+ for lig in ligs:
+ if 'time_stepping_scheme' in lig:
+ print '----> '+fic+' already processed to r19xxx'
+ return
+ #
+ a1='time_stepping_scheme = 1 # 1 = Newmark (2nd order), \
+ 2 = LDDRK4-6 (4th-order 6-stage low storage Runge-Kutta), \
+ 3 = classical 4th-order 4-stage Runge-Kutta\n'
+
+ #--------------------------------------------------------------------------
+ # Ajout des parametres supplementaires
+ #
+ for ilg, lig in enumerate(ligs):
+ if lig.startswith('USER_T0'):
+ ligs.insert(ilg+1,a1)
+ #
+ move(fic,fic+'.before_update_to_r19xxx')
+ #
+ fm = open(fic,'w')
+ fm.writelines(ligs)
+ fm.close()
+ #
+ print 'xxxxx------> '+fic+' processed to r19xxx'
+ return
+#------------------------------------------------------------------------------
+if __name__=='__main__':
+ ## Liste de tous les fichiers à partir du répertoire courant
+ Fichiers=[]
+ mylister('.')
+ #
+ print '~'*80
+ Ct_Par_file=0
+ for fic in Fichiers:
+ repert, ficname = os.path.split(fic)
+ if not( ('.svn' in repert) or ('unused' in repert) or \
+ '.before_update_to_' in ficname):
+ if ficname.startswith('Par_file'):
+ print 'Analysis of file : '+fic
+ if not (ficname.endswith('~')):
+ Ct_Par_file+=1
+ ProcessParfile_r19201(fic)
+ ProcessParfile_r19340(fic)
+ ProcessParfile_r19346(fic)
+ ProcessParfile_r19xxx(fic)
+ print '~'*80
+ #
+ print 'Number of Par_file analysed : ', Ct_Par_file
+ print 'END OF Par_file PROCESSING'
+
Modified: seismo/2D/SPECFEM2D/trunk/EXAMPLES/semi_infinite_homo/Par_file_elastic_2D
===================================================================
--- seismo/2D/SPECFEM2D/trunk/EXAMPLES/semi_infinite_homo/Par_file_elastic_2D 2012-03-18 19:59:34 UTC (rev 19802)
+++ seismo/2D/SPECFEM2D/trunk/EXAMPLES/semi_infinite_homo/Par_file_elastic_2D 2012-03-19 01:25:41 UTC (rev 19803)
@@ -9,7 +9,7 @@
# parameters concerning partitioning
nproc = 1 # number of processes
partitioning_method = 3 # SCOTCH = 3, ascending order (very bad idea) = 1
-PERFORM_CUTHILL_MCKEE = .true. # perform inverse Cuthill-McKee (1969) optimization/permutation for mesh numbering
+PERFORM_CUTHILL_MCKEE = .false. # perform inverse Cuthill-McKee (1969) optimization/permutation for mesh numbering (can be very costly and not very useful)
ngnod = 4 # number of control nodes per element (4 or 9)
initialfield = .false. # use a plane wave as source or not
Deleted: seismo/2D/SPECFEM2D/trunk/UTILS/process_DATA_Par_files_to_update_their_format_when_new_parameters_are_added.py
===================================================================
--- seismo/2D/SPECFEM2D/trunk/UTILS/process_DATA_Par_files_to_update_their_format_when_new_parameters_are_added.py 2012-03-18 19:59:34 UTC (rev 19802)
+++ seismo/2D/SPECFEM2D/trunk/UTILS/process_DATA_Par_files_to_update_their_format_when_new_parameters_are_added.py 2012-03-19 01:25:41 UTC (rev 19803)
@@ -1,191 +0,0 @@
-# -*- coding: utf-8 -*-
-"""
-Created on Tue Nov 15 09:00:00 2011
-Updated on Fri Jan 13 2012
-
-Processing of Par_file to update them to new format
-
-Usage : "python PathTo/SPECFEM2D/UTILS/ProcessParFileParametersToNewRelease.py"
-This will process all Par_file starting from current directory
-
- at author: Cristini Paul, Laboratoire de Mecanique et d'Acoustique, CNRS, Marseille, France
-"""
-import os, string, sys
-from os import listdir, walk
-from shutil import copy, move
-from os.path import exists
-#------------------------------------------------------------------------------
-def LoadLig(Fichier):
- f = open(Fichier,'r')
- ligs= f.readlines()
- f.close()
- return ligs
-#------------------------------------------------------------------------------
-def mylister(currdir):
- for file in os.listdir(currdir):
- path=os.path.join(currdir, file)
- if not os.path.isdir(path):
- #print path
- Fichiers.append(path)
- else:
- mylister(path)
-#------------------------------------------------------------------------------
-def ProcessParfile_r19201(fic):
- # Open the file and get all lines from Par_file
- ligs= LoadLig(fic)
- # Test pour voir si le traitement a déjà été fait
- for lig in ligs:
- if 'ADD_PERIODIC_CONDITIONS' in lig:
- print '----> '+fic+' already processed to r19201'
- return
- # New additions to the Par_file
- a1='PERFORM_CUTHILL_MCKEE = .true. # perform inverse Cuthill-McKee (1969) optimization/permutation for mesh numbering\n'
- a2='USER_T0 = 0.0d0 # use this t0 as earliest starting time rather than the automatically calculated one\n'
- a3='SU_FORMAT = .false. # output seismograms in Seismic Unix format (adjoint traces will be read in the same format)\n'
- a4='factor_subsample_image = 1 # factor to subsample color images output by the code (useful for very large models)\n'+ \
- 'POWER_DISPLAY_COLOR = 0.30d0 # non linear display to enhance small amplitudes in color images\n'+ \
- 'DRAW_WATER_CONSTANT_BLUE_IN_JPG = .true. # display acoustic layers as constant blue in JPEG images, because they likely correspond to water\n'
- a5='US_LETTER = .false. # US letter paper or European A4\n'+ \
- 'USE_SNAPSHOT_NUMBER_IN_FILENAME = .false. # use snapshot number in the file name of JPEG color snapshots instead of the time step\n'
- a6='\n# for horizontal periodic conditions: detect common points between left and right edges\n'+ \
- 'ADD_PERIODIC_CONDITIONS = .false.\n\n'+ \
- '# horizontal periodicity distance for periodic conditions\n'+ \
- 'PERIODIC_horiz_dist = 0.3597d0\n\n'+ \
- '# grid point detection tolerance for periodic conditions\n'+ \
- 'PERIODIC_DETECT_TOL = 3.3334d-6\n'
- #--------------------------------------------------------------------------
- # Ajout des parametres supplementaires
- #
- for ilg, lig in enumerate(ligs):
- if lig.startswith('partitioning'):
- ligs.insert(ilg+1,a1)
-
- if lig.startswith('deltat'):
- ligs.insert(ilg+1,a2)
-
- if lig.startswith('rec_normal'):
- ligs.insert(ilg+1,a3)
-
- if lig.startswith('subsamp'):
- ligs[ilg]=string.replace(ligs[ilg],'subsamp ','subsamp_postscript',1)
- ligs.insert(ilg+1,a4)
-
- if lig.startswith('sizemax'):
- ligs.insert(ilg+1,a5)
-
- if lig.startswith('absorbing_conditions'):
- ligs.insert(ilg+1,a6)
- #
- move(fic,fic+'.before_update_to_r19201')
- fm = open(fic,'w')
- fm.writelines(ligs)
- fm.close()
- #
- print 'xxxxx------> '+fic+' processed to r19201'
- return
-#------------------------------------------------------------------------------
-def ProcessParfile_r19340(fic):
- # Open the file and get all lines from Par_file
- ligs= LoadLig(fic)
- # Teste si le traitement a déjà été fait
- for lig in ligs:
- if 'nreceiversets' in lig:
- print '----> '+fic+' already processed to r19340'
- return
- #
- # Ajout des parametres supplementaires
- #
- for ilg, lig in enumerate(ligs):
- if lig.startswith('nreceiverlines'):
- ligs[ilg]=ligs[ilg].replace('lines','sets ')
- #
- move(fic,fic+'.before_update_to_r19340')
- #
- fm = open(fic,'w')
- fm.writelines(ligs)
- fm.close()
- #
- print 'xxxxx------> '+fic+' processed to r19340'
- return
-#------------------------------------------------------------------------------
-def ProcessParfile_r19346(fic):
- # Open the file and get all lines from Par_file
- ligs= LoadLig(fic)
- # Teste si le traitement a déjà été fait
- for lig in ligs:
- if 'ATTENUATION_PORO_FLUID_PART' in lig:
- print '----> '+fic+' already processed to r19346'
- return
- #--------------------------------------------------------------------------
- # Ajout des parametres supplementaires
- #
- for ilg, lig in enumerate(ligs):
- if lig.startswith('TURN_ATTENUATION_ON'):
- ligs[ilg]=ligs[ilg].replace('TURN_ATTENUATION_ON ', \
- 'ATTENUATION_VISCOELASTIC_SOLID')
- if lig.startswith('TURN_VISCATTENUATION_ON'):
- ligs[ilg]=ligs[ilg].replace('TURN_VISCATTENUATION_ON ', \
- 'ATTENUATION_PORO_FLUID_PART')
- #
- move(fic,fic+'.before_update_to_r19346')
- #
- fm = open(fic,'w')
- fm.writelines(ligs)
- fm.close()
- #
- print 'xxxxx------> '+fic+' processed to r19346'
- return
-#------------------------------------------------------------------------------
-def ProcessParfile_r19xxx(fic):
- # Open the file and get all lines from Par_file
- ligs= LoadLig(fic)
- # Teste si le traitement a déjà été fait
- for lig in ligs:
- if 'time_stepping_scheme' in lig:
- print '----> '+fic+' already processed to r19xxx'
- return
- #
- a1='time_stepping_scheme = 1 # 1 = Newmark (2nd order), \
- 2 = LDDRK4-6 (4th-order 6-stage low storage Runge-Kutta), \
- 3 = classical 4th-order 4-stage Runge-Kutta\n'
-
- #--------------------------------------------------------------------------
- # Ajout des parametres supplementaires
- #
- for ilg, lig in enumerate(ligs):
- if lig.startswith('USER_T0'):
- ligs.insert(ilg+1,a1)
- #
- move(fic,fic+'.before_update_to_r19xxx')
- #
- fm = open(fic,'w')
- fm.writelines(ligs)
- fm.close()
- #
- print 'xxxxx------> '+fic+' processed to r19xxx'
- return
-#------------------------------------------------------------------------------
-if __name__=='__main__':
- ## Liste de tous les fichiers à partir du répertoire courant
- Fichiers=[]
- mylister('.')
- #
- print '~'*80
- Ct_Par_file=0
- for fic in Fichiers:
- repert, ficname = os.path.split(fic)
- if not( ('.svn' in repert) or ('unused' in repert) or \
- '.before_update_to_' in ficname):
- if ficname.startswith('Par_file'):
- print 'Analysis of file : '+fic
- if not (ficname.endswith('~')):
- Ct_Par_file+=1
- ProcessParfile_r19201(fic)
- ProcessParfile_r19340(fic)
- ProcessParfile_r19346(fic)
- ProcessParfile_r19xxx(fic)
- print '~'*80
- #
- print 'Number of Par_file analysed : ', Ct_Par_file
- print 'END OF Par_file PROCESSING'
-
\ No newline at end of file
Modified: seismo/2D/SPECFEM2D/trunk/setup/constants.h.in
===================================================================
--- seismo/2D/SPECFEM2D/trunk/setup/constants.h.in 2012-03-18 19:59:34 UTC (rev 19802)
+++ seismo/2D/SPECFEM2D/trunk/setup/constants.h.in 2012-03-19 01:25:41 UTC (rev 19803)
@@ -46,7 +46,7 @@
! implement Cuthill-McKee or replace with identity permutation
logical, parameter :: ACTUALLY_IMPLEMENT_PERM_OUT = .false.
logical, parameter :: ACTUALLY_IMPLEMENT_PERM_INN = .false.
- logical, parameter :: ACTUALLY_IMPLEMENT_PERM_WHOLE = .false. ! .true.
+ logical, parameter :: ACTUALLY_IMPLEMENT_PERM_WHOLE = .true.
! create file DATA/model_velocity.dat_output or not (can be huge and slow, thus off by default)
logical, parameter :: OUTPUT_MODEL_VELOCITY_FILE = .false.
More information about the CIG-COMMITS
mailing list