[cig-commits] r16285 - in seismo/3D/SPECFEM3D_GLOBE/trunk: . DATA UTILS UTILS/Visualization/Paraview

danielpeter at geodynamics.org danielpeter at geodynamics.org
Sat Feb 20 12:40:53 PST 2010


Author: danielpeter
Date: 2010-02-20 12:40:52 -0800 (Sat, 20 Feb 2010)
New Revision: 16285

Added:
   seismo/3D/SPECFEM3D_GLOBE/trunk/DATA/CMTSOLUTION.default
   seismo/3D/SPECFEM3D_GLOBE/trunk/DATA/Par_file.default
   seismo/3D/SPECFEM3D_GLOBE/trunk/DATA/STATIONS.default
   seismo/3D/SPECFEM3D_GLOBE/trunk/INSTALL
Removed:
   seismo/3D/SPECFEM3D_GLOBE/trunk/DATA/CMTSOLUTION
   seismo/3D/SPECFEM3D_GLOBE/trunk/DATA/Par_file
   seismo/3D/SPECFEM3D_GLOBE/trunk/DATA/STATIONS
   seismo/3D/SPECFEM3D_GLOBE/trunk/UTILS/Paraview/
   seismo/3D/SPECFEM3D_GLOBE/trunk/UTILS/Visualization/Paraview/global_slice_number.daniel.pl
Modified:
   seismo/3D/SPECFEM3D_GLOBE/trunk/README_SPECFEM3D_GLOBE
Log:
further changes to move towards the new version 'Tiger':

the files DATA/Par_file, DATA/CMTSOLUTION, DATA/STATIONS are
renamed:
DATA/Par_file       -> DATA/Par_file.default
DATA/CMTSOLUTION    -> DATA/CMTSOLUTION.default
DATA/STATIONS       -> DATA/STATIONS.default

to avoid having the user changes in his own files being overwritten
when updating with SVN to the new version.

the new default simulation has slightly changed: 
it will run on 150 CPUs, using a higher NEX resolution and a
more realistic 3D reference model run.

a step-by-step tutorial describing how to run this default simulation, 
is added in: 
README_SPECFEM3D_GLOBE

similar with the installation process, described now in a new file:
INSTALL

these procedures repeat in short what is already described in the user manual.




Deleted: seismo/3D/SPECFEM3D_GLOBE/trunk/DATA/CMTSOLUTION
===================================================================
--- seismo/3D/SPECFEM3D_GLOBE/trunk/DATA/CMTSOLUTION	2010-02-20 19:58:48 UTC (rev 16284)
+++ seismo/3D/SPECFEM3D_GLOBE/trunk/DATA/CMTSOLUTION	2010-02-20 20:40:52 UTC (rev 16285)
@@ -1,13 +0,0 @@
-PDE 1994  6  9  0 33 16.40 -13.8300  -67.5600 637.0 6.9 6.8 NORTHERNBOLIVIA
-event name:     060994A
-time shift:      0.0000
-half duration:  35.0000
-latitude:      -13.8200
-longitude:     -67.2500
-depth:         647.1000
-Mrr:      -7.600000e+27
-Mtt:       7.700000e+27
-Mpp:      -2.000000e+26
-Mrt:      -2.500000e+28
-Mrp:       4.000000e+26
-Mtp:      -2.500000e+27

Added: seismo/3D/SPECFEM3D_GLOBE/trunk/DATA/CMTSOLUTION.default
===================================================================
--- seismo/3D/SPECFEM3D_GLOBE/trunk/DATA/CMTSOLUTION.default	                        (rev 0)
+++ seismo/3D/SPECFEM3D_GLOBE/trunk/DATA/CMTSOLUTION.default	2010-02-20 20:40:52 UTC (rev 16285)
@@ -0,0 +1,13 @@
+PDE 1994  6  9  0 33 16.40 -13.8300  -67.5600 637.0 6.9 6.8 NORTHERNBOLIVIA
+event name:     060994A
+time shift:      0.0000
+half duration:  35.0000
+latitude:      -13.8200
+longitude:     -67.2500
+depth:         647.1000
+Mrr:      -7.600000e+27
+Mtt:       7.700000e+27
+Mpp:      -2.000000e+26
+Mrt:      -2.500000e+28
+Mrp:       4.000000e+26
+Mtp:      -2.500000e+27

Deleted: seismo/3D/SPECFEM3D_GLOBE/trunk/DATA/Par_file
===================================================================
--- seismo/3D/SPECFEM3D_GLOBE/trunk/DATA/Par_file	2010-02-20 19:58:48 UTC (rev 16284)
+++ seismo/3D/SPECFEM3D_GLOBE/trunk/DATA/Par_file	2010-02-20 20:40:52 UTC (rev 16285)
@@ -1,114 +0,0 @@
-
-# forward or adjoint simulation
-SIMULATION_TYPE                 = 1
-SAVE_FORWARD                    = .false.  # save last frame of forward simulation or not
-
-# number of chunks (1,2,3 or 6)
-NCHUNKS                         = 6
-
-# angular width of the first chunk (not used if full sphere with six chunks)
-ANGULAR_WIDTH_XI_IN_DEGREES   = 90.d0      # angular size of a chunk
-ANGULAR_WIDTH_ETA_IN_DEGREES  = 90.d0
-CENTER_LATITUDE_IN_DEGREES    = 40.d0
-CENTER_LONGITUDE_IN_DEGREES   = 10.d0
-GAMMA_ROTATION_AZIMUTH        = 20.d0
-
-# number of elements at the surface along the two sides of the first chunk
-# (must be multiple of 16 and 8 * multiple of NPROC below)
-NEX_XI                          = 96
-NEX_ETA                         = 96
-
-# number of MPI processors along the two sides of the first chunk
-NPROC_XI                        = 6
-NPROC_ETA                       = 6
-
-# 1D models with real structure:
-# 1D_isotropic_prem, 1D_transversely_isotropic_prem, 1D_iasp91, 1D_1066a, 1D_ak135, 1D_ref, 1D_ref_iso, 1D_jp3d,1D_sea99
-#
-# 1D models with only one fictitious averaged crustal layer:
-# 1D_isotropic_prem_onecrust, 1D_transversely_isotropic_prem_onecrust, 1D_iasp91_onecrust, 1D_1066a_onecrust, 1D_ak135_onecrust
-#
-# fully 3D models:
-# transversely_isotropic_prem_plus_3D_crust_2.0, 3D_anisotropic, 3D_attenuation,
-# s20rts, s362ani, s362iso, s362wmani, s362ani_prem, s362ani_3DQ, s362iso_3DQ,
-# s29ea, s29ea,sea99_jp3d1994,sea99,jp3d1994,heterogen
-MODEL                           = 1D_isotropic_prem
-
-# parameters describing the Earth model
-OCEANS                          = .false.
-ELLIPTICITY                     = .false.
-TOPOGRAPHY                      = .false.
-GRAVITY                         = .false.
-ROTATION                        = .false.
-ATTENUATION                     = .false.
-
-# absorbing boundary conditions for a regional simulation
-ABSORBING_CONDITIONS            = .false.
-
-# record length in minutes
-RECORD_LENGTH_IN_MINUTES        = 75.0d0
-
-# save AVS or OpenDX movies
-#MOVIE_COARSE saves movie only at corners of elements (SURFACE OR VOLUME) 
-#MOVIE_COARSE does not work with create_movie_AVS_DX
-MOVIE_SURFACE                   = .false.
-MOVIE_VOLUME                    = .false.
-MOVIE_COARSE                    = .false.
-NTSTEP_BETWEEN_FRAMES           = 100
-HDUR_MOVIE                      = 0.d0
-
-# save movie in volume.  Will save element if center of element is in prescribed volume
-# top/bottom: depth in KM, use MOVIE_TOP = -100 to make sure the surface is stored.
-# west/east: longitude, degrees East [-180/180] top/bottom: latitute, degrees North [-90/90]
-# start/stop: frames will be stored at MOVIE_START + i*NSTEP_BETWEEN_FRAMES, where i=(0,1,2..) and iNSTEP_BETWEEN_FRAMES <= MOVIE_STOP
-# movie_volume_type: 1=strain, 2=time integral of strain, 3=\mu*time integral of strain
-# type 4 saves the trace and deviatoric stress in the whole volume, 5=displacement, 6=velocity
-MOVIE_VOLUME_TYPE               = 2
-MOVIE_TOP_KM                    = -100.0
-MOVIE_BOTTOM_KM                 = 1000.0
-MOVIE_WEST_DEG                  = -90.0
-MOVIE_EAST_DEG                  = 90.0
-MOVIE_NORTH_DEG                 = 90.0
-MOVIE_SOUTH_DEG                 = -90.0
-MOVIE_START                     = 0
-MOVIE_STOP                      = 40000
-
-# save mesh files to check the mesh
-SAVE_MESH_FILES                 = .false.
-
-# restart files (number of runs can be 1, 2 or 3, choose 1 for no restart files)
-NUMBER_OF_RUNS                  = 1
-NUMBER_OF_THIS_RUN              = 1
-
-# path to store the local database files on each node
-LOCAL_PATH                      = /scratch/
-
-# interval at which we output time step info and max of norm of displacement
-NTSTEP_BETWEEN_OUTPUT_INFO      = 1000
-
-# interval in time steps for temporary writing of seismograms
-NTSTEP_BETWEEN_OUTPUT_SEISMOS   = 5000000
-NTSTEP_BETWEEN_READ_ADJSRC      = 1000
-
-# output format for the seismograms (one can use either or all of the three formats)
-OUTPUT_SEISMOS_ASCII_TEXT       = .true.
-OUTPUT_SEISMOS_SAC_ALPHANUM     = .false.
-OUTPUT_SEISMOS_SAC_BINARY       = .false.
-
-# rotate seismograms to Radial-Transverse-Z or use default North-East-Z reference frame
-ROTATE_SEISMOGRAMS_RT           = .false.
-
-# decide if master process writes all the seismograms or if all processes do it in parallel
-WRITE_SEISMOGRAMS_BY_MASTER     = .true.
-
-# save all seismograms in one large combined file instead of one file per seismogram
-# to avoid overloading shared non-local file systems such as GPFS for instance
-SAVE_ALL_SEISMOS_IN_ONE_FILE    = .false.
-USE_BINARY_FOR_LARGE_FILE       = .false.
-
-# flag to impose receivers at the surface or allow them to be buried
-RECEIVERS_CAN_BE_BURIED         = .true.
-
-# print source time function
-PRINT_SOURCE_TIME_FUNCTION      = .false.
-

Added: seismo/3D/SPECFEM3D_GLOBE/trunk/DATA/Par_file.default
===================================================================
--- seismo/3D/SPECFEM3D_GLOBE/trunk/DATA/Par_file.default	                        (rev 0)
+++ seismo/3D/SPECFEM3D_GLOBE/trunk/DATA/Par_file.default	2010-02-20 20:40:52 UTC (rev 16285)
@@ -0,0 +1,114 @@
+
+# forward or adjoint simulation
+SIMULATION_TYPE                 = 1
+SAVE_FORWARD                    = .false.  # save last frame of forward simulation or not
+
+# number of chunks (1,2,3 or 6)
+NCHUNKS                         = 6
+
+# angular width of the first chunk (not used if full sphere with six chunks)
+ANGULAR_WIDTH_XI_IN_DEGREES   = 90.d0      # angular size of a chunk
+ANGULAR_WIDTH_ETA_IN_DEGREES  = 90.d0
+CENTER_LATITUDE_IN_DEGREES    = 40.d0
+CENTER_LONGITUDE_IN_DEGREES   = 10.d0
+GAMMA_ROTATION_AZIMUTH        = 20.d0
+
+# number of elements at the surface along the two sides of the first chunk
+# (must be multiple of 16 and 8 * multiple of NPROC below)
+NEX_XI                          = 240
+NEX_ETA                         = 240
+
+# number of MPI processors along the two sides of the first chunk
+NPROC_XI                        = 5
+NPROC_ETA                       = 5
+
+# 1D models with real structure:
+# 1D_isotropic_prem, 1D_transversely_isotropic_prem, 1D_iasp91, 1D_1066a, 1D_ak135, 1D_ref, 1D_ref_iso, 1D_jp3d,1D_sea99
+#
+# 1D models with only one fictitious averaged crustal layer:
+# 1D_isotropic_prem_onecrust, 1D_transversely_isotropic_prem_onecrust, 1D_iasp91_onecrust, 1D_1066a_onecrust, 1D_ak135_onecrust
+#
+# fully 3D models:
+# transversely_isotropic_prem_plus_3D_crust_2.0, 3D_anisotropic, 3D_attenuation,
+# s20rts, s362ani, s362iso, s362wmani, s362ani_prem, s362ani_3DQ, s362iso_3DQ,
+# s29ea, s29ea,sea99_jp3d1994,sea99,jp3d1994,heterogen
+MODEL                           = s362ani
+
+# parameters describing the Earth model
+OCEANS                          = .true.
+ELLIPTICITY                     = .true.
+TOPOGRAPHY                      = .true.
+GRAVITY                         = .true.
+ROTATION                        = .true.
+ATTENUATION                     = .true.
+
+# absorbing boundary conditions for a regional simulation
+ABSORBING_CONDITIONS            = .false.
+
+# record length in minutes
+RECORD_LENGTH_IN_MINUTES        = 15.0d0
+
+# save AVS or OpenDX movies
+#MOVIE_COARSE saves movie only at corners of elements (SURFACE OR VOLUME) 
+#MOVIE_COARSE does not work with create_movie_AVS_DX
+MOVIE_SURFACE                   = .false.
+MOVIE_VOLUME                    = .false.
+MOVIE_COARSE                    = .false.
+NTSTEP_BETWEEN_FRAMES           = 100
+HDUR_MOVIE                      = 0.d0
+
+# save movie in volume.  Will save element if center of element is in prescribed volume
+# top/bottom: depth in KM, use MOVIE_TOP = -100 to make sure the surface is stored.
+# west/east: longitude, degrees East [-180/180] top/bottom: latitute, degrees North [-90/90]
+# start/stop: frames will be stored at MOVIE_START + i*NSTEP_BETWEEN_FRAMES, where i=(0,1,2..) and iNSTEP_BETWEEN_FRAMES <= MOVIE_STOP
+# movie_volume_type: 1=strain, 2=time integral of strain, 3=\mu*time integral of strain
+# type 4 saves the trace and deviatoric stress in the whole volume, 5=displacement, 6=velocity
+MOVIE_VOLUME_TYPE               = 2
+MOVIE_TOP_KM                    = -100.0
+MOVIE_BOTTOM_KM                 = 1000.0
+MOVIE_WEST_DEG                  = -90.0
+MOVIE_EAST_DEG                  = 90.0
+MOVIE_NORTH_DEG                 = 90.0
+MOVIE_SOUTH_DEG                 = -90.0
+MOVIE_START                     = 0
+MOVIE_STOP                      = 40000
+
+# save mesh files to check the mesh
+SAVE_MESH_FILES                 = .false.
+
+# restart files (number of runs can be 1, 2 or 3, choose 1 for no restart files)
+NUMBER_OF_RUNS                  = 1
+NUMBER_OF_THIS_RUN              = 1
+
+# path to store the local database files on each node
+LOCAL_PATH                      = ./DATABASES_MPI
+
+# interval at which we output time step info and max of norm of displacement
+NTSTEP_BETWEEN_OUTPUT_INFO      = 1000
+
+# interval in time steps for temporary writing of seismograms
+NTSTEP_BETWEEN_OUTPUT_SEISMOS   = 5000000
+NTSTEP_BETWEEN_READ_ADJSRC      = 1000
+
+# output format for the seismograms (one can use either or all of the three formats)
+OUTPUT_SEISMOS_ASCII_TEXT       = .true.
+OUTPUT_SEISMOS_SAC_ALPHANUM     = .false.
+OUTPUT_SEISMOS_SAC_BINARY       = .false.
+
+# rotate seismograms to Radial-Transverse-Z or use default North-East-Z reference frame
+ROTATE_SEISMOGRAMS_RT           = .false.
+
+# decide if master process writes all the seismograms or if all processes do it in parallel
+WRITE_SEISMOGRAMS_BY_MASTER     = .true.
+
+# save all seismograms in one large combined file instead of one file per seismogram
+# to avoid overloading shared non-local file systems such as GPFS for instance
+SAVE_ALL_SEISMOS_IN_ONE_FILE    = .false.
+USE_BINARY_FOR_LARGE_FILE       = .false.
+
+# flag to impose receivers at the surface or allow them to be buried
+RECEIVERS_CAN_BE_BURIED         = .true.
+
+# print source time function
+PRINT_SOURCE_TIME_FUNCTION      = .false.
+

Deleted: seismo/3D/SPECFEM3D_GLOBE/trunk/DATA/STATIONS
===================================================================
--- seismo/3D/SPECFEM3D_GLOBE/trunk/DATA/STATIONS	2010-02-20 19:58:48 UTC (rev 16284)
+++ seismo/3D/SPECFEM3D_GLOBE/trunk/DATA/STATIONS	2010-02-20 20:40:52 UTC (rev 16285)
@@ -1,129 +0,0 @@
-AAK        II       42.6390     74.4940    1645.0    30.0
-ABKT       II       37.9304     58.1189     678.0     7.0
-ABPO       II      -19.0180     47.2290    1528.0     5.3
-ALE        II       82.5033    -62.3500      60.0     0.0
-ARU        II       56.4302     58.5625     250.0     0.0
-ASCN       II       -7.9327    -14.3601     173.0   100.0
-BFO        II       48.3319      8.3311     589.0     0.0
-BORG       II       64.7474    -21.3268     110.0    95.0
-BRVK       II       53.0581     70.2828     330.0    15.0
-CMLA       II       37.7637    -25.5243     429.0     7.0
-COCO       II      -12.1901     96.8349       1.0    70.0
-DGAR       II       -7.4121     72.4525       1.0     2.0
-EFI        II      -51.6753    -58.0637     110.0    80.0
-ERM        II       42.0150    143.1572      40.0     0.0
-ESK        II       55.3167     -3.2050     242.0     0.0
-FFC        II       54.7250   -101.9783     338.0     0.0
-GAR        II       39.0000     70.3167    1300.0     0.0
-HOPE       II      -54.2836    -36.4879      20.0     0.0
-JTS        II       10.2908    -84.9525     340.0     0.0
-KAPI       II       -5.0142    119.7517     300.0   100.0
-KDAK       II       57.7828   -152.5835     152.0     5.5
-KIV        II       43.9562     42.6888    1210.0     0.0
-KURK       II       50.7154     78.6202     184.0    25.0
-KWAJ       II        8.8019    167.6130       0.0     0.0
-LVZ        II       67.8979     34.6514     630.0   200.0
-MBAR       II       -0.6019     30.7382    1390.0   100.0
-MSEY       II       -4.6737     55.4792     475.0    91.0
-MSVF       II      -17.7448    178.0528     801.1   100.0
-NIL        II       33.6506     73.2686     629.0    68.0
-NNA        II      -11.9875    -76.8422     575.0    40.0
-NRIL       II       69.5049     88.4414      92.0   506.0
-NVS        II       54.8404     83.2346     150.0     0.0
-OBN        II       55.1146     36.5674     160.0    30.0
-PALK       II        7.2728     80.7022     460.0    90.0
-RAYN       II       23.5225     45.5032     631.0     2.0
-RPN        II      -27.1267   -109.3344     110.0     0.0
-SACV       II       14.9702    -23.6085     387.0    97.0
-SHEL       II      -15.9588     -5.7457     537.0    60.0
-SUR        II      -32.3797     20.8117    1770.0     0.0
-TAU        II      -42.9099    147.3204     132.0     0.0
-TLY        II       51.6807    103.6438     579.0    20.0
-UOSS       II       24.9453     56.2042     284.4     0.0
-WRAB       II      -19.9336    134.3600     366.0   100.0
-XPFO       II       33.6107   -116.4555    1280.0     0.0
-AAE        IU        9.0292     38.7656    2442.0     0.0
-ADK        IU       51.8823   -176.6842     130.0     0.0
-AFI        IU      -13.9093   -171.7773     705.0     1.0
-ANMO       IU       34.9459   -106.4572    1720.0   100.0
-ANTO       IU       39.8680     32.7934     895.0   195.0
-BBSR       IU       32.3713    -64.6963      -1.3    31.4
-BILL       IU       68.0653    166.4531     320.0     0.0
-BOCO       IU        4.5869    -74.0432    3137.0    23.0
-CASY       IU      -66.2792    110.5354       5.0     5.0
-CCM        IU       38.0557    -91.2446     171.0    51.0
-CHTO       IU       18.8141     98.9443     420.0     0.0
-COLA       IU       64.8736   -147.8616     200.0     0.0
-COR        IU       44.5855   -123.3046     110.0     0.0
-CTAO       IU      -20.0882    146.2545     320.0    37.0
-DAV        IU        7.0697    125.5791     149.0     1.0
-DWPF       IU       28.1103    -81.4327    -132.0   162.0
-FUNA       IU       -8.5259    179.1966      19.0     1.0
-FURI       IU        8.8952     38.6798    2565.0     5.0
-GNI        IU       40.1480     44.7410    1509.0   100.0
-GRFO       IU       49.6909     11.2203     384.0   116.0
-GUMO       IU       13.5893    144.8684      61.0   109.0
-HKT        IU       29.9618    -95.8384    -413.0   450.0
-HNR        IU       -9.4387    159.9475       0.0   100.0
-HRV        IU       42.5064    -71.5583     200.0     0.0
-INCN       IU       37.4776    126.6239      79.0     1.0
-JOHN       IU       16.7329   -169.5292     -37.0    39.0
-KBS        IU       78.9154     11.9385      90.0     3.0
-KEV        IU       69.7565     27.0035      85.0    15.0
-KIEV       IU       50.7012     29.2242     140.0    40.0
-KIP        IU       21.4233   -158.0150      37.0    33.0
-KMBO       IU       -1.1271     37.2525    1930.0    20.0
-KNTN       IU       -2.7744   -171.7186      18.0     2.0
-KONO       IU       59.6491      9.5982     216.0   340.0
-KOWA       IU       14.4967     -4.0140     316.0     5.0
-LCO        IU      -29.0110    -70.7004    2300.0     0.0
-LSZ        IU      -15.2779     28.1882    1200.0     0.0
-LVC        IU      -22.6127    -68.9111    2930.0    30.0
-MA2        IU       59.5756    150.7700     337.0     2.0
-MAJO       IU       36.5457    138.2041     405.0     0.0
-MAKZ       IU       46.8080     81.9770     590.0    10.0
-MBWA       IU      -21.1590    119.7313     181.0     9.0
-MIDW       IU       28.2155   -177.3697      17.8     1.0
-MSKU       IU       -1.6557     13.6116     287.0    25.0
-NAI        IU       -1.2739     36.8037    1692.0     0.0
-NWAO       IU      -32.9277    117.2390     370.9     9.1
-OTAV       IU        0.2376    -78.4508    3495.0    15.0
-PAB        IU       39.5446     -4.3499     950.0     0.0
-PAYG       IU       -0.6742    -90.2861     170.0   100.0
-PET        IU       53.0233    158.6499     105.0     5.0
-PMG        IU       -9.4047    147.1597      90.0     0.0
-PMSA       IU      -64.7744    -64.0489      40.0     0.0
-POHA       IU       19.7573   -155.5326    1910.0    80.0
-PTCN       IU      -25.0713   -130.0953     218.0     2.0
-PTGA       IU       -0.7308    -59.9666     141.0    96.0
-QSPA       IU      -89.9289    144.4382    2847.0     3.0
-RAIO       IU       46.0403   -122.8851      11.0     0.0
-RAO        IU      -29.2450   -177.9290      59.5     0.5
-RAR        IU      -21.2125   -159.7733     -72.0   100.0
-RCBR       IU       -5.8274    -35.9014     291.0   109.0
-RSSD       IU       44.1212   -104.0359    2022.7    67.3
-SAML       IU       -8.9489    -63.1831     120.0     0.0
-SBA        IU      -77.8492    166.7572      48.0     2.0
-SDV        IU        8.8839    -70.6340    1588.0    32.0
-SFJD       IU       66.9961    -50.6208     329.0     1.0
-SJG        IU       18.1091    -66.1500     420.0     0.0
-SLBS       IU       23.6858   -109.9443     825.0     0.0
-SNZO       IU      -41.3087    174.7043      20.0   100.0
-SSPA       IU       40.6358    -77.8876     170.0   100.0
-TARA       IU        1.3549    172.9229      19.0     1.0
-TATO       IU       24.9735    121.4971      77.1    82.9
-TBT        IU       28.6794    -17.9145     180.0    40.0
-TEIG       IU       20.2263    -88.2763      15.0    25.0
-TIXI       IU       71.6341    128.8667      40.0     0.0
-TOL        IU       39.8814     -4.0485     480.0     0.0
-TRIS       IU      -37.0681    -12.3152      58.0     2.0
-TRQA       IU      -38.0568    -61.9787     439.0   101.0
-TSUM       IU      -19.2022     17.5838    1260.0     0.0
-TUC        IU       32.3098   -110.7847     909.0     1.0
-ULN        IU       47.8651    107.0532    1610.0     0.0
-WAKE       IU       19.2834    166.6520      19.0     1.0
-WCI        IU       38.2289    -86.2939      78.0   132.0
-WVT        IU       36.1297    -87.8300     170.0     0.0
-XMAS       IU        2.0448   -157.4457      19.0     1.0
-YAK        IU       62.0310    129.6805     110.0    14.0
-YSS        IU       46.9587    142.7604     148.0     2.0

Added: seismo/3D/SPECFEM3D_GLOBE/trunk/DATA/STATIONS.default
===================================================================
--- seismo/3D/SPECFEM3D_GLOBE/trunk/DATA/STATIONS.default	                        (rev 0)
+++ seismo/3D/SPECFEM3D_GLOBE/trunk/DATA/STATIONS.default	2010-02-20 20:40:52 UTC (rev 16285)
@@ -0,0 +1,129 @@
+AAK        II       42.6390     74.4940    1645.0    30.0
+ABKT       II       37.9304     58.1189     678.0     7.0
+ABPO       II      -19.0180     47.2290    1528.0     5.3
+ALE        II       82.5033    -62.3500      60.0     0.0
+ARU        II       56.4302     58.5625     250.0     0.0
+ASCN       II       -7.9327    -14.3601     173.0   100.0
+BFO        II       48.3319      8.3311     589.0     0.0
+BORG       II       64.7474    -21.3268     110.0    95.0
+BRVK       II       53.0581     70.2828     330.0    15.0
+CMLA       II       37.7637    -25.5243     429.0     7.0
+COCO       II      -12.1901     96.8349       1.0    70.0
+DGAR       II       -7.4121     72.4525       1.0     2.0
+EFI        II      -51.6753    -58.0637     110.0    80.0
+ERM        II       42.0150    143.1572      40.0     0.0
+ESK        II       55.3167     -3.2050     242.0     0.0
+FFC        II       54.7250   -101.9783     338.0     0.0
+GAR        II       39.0000     70.3167    1300.0     0.0
+HOPE       II      -54.2836    -36.4879      20.0     0.0
+JTS        II       10.2908    -84.9525     340.0     0.0
+KAPI       II       -5.0142    119.7517     300.0   100.0
+KDAK       II       57.7828   -152.5835     152.0     5.5
+KIV        II       43.9562     42.6888    1210.0     0.0
+KURK       II       50.7154     78.6202     184.0    25.0
+KWAJ       II        8.8019    167.6130       0.0     0.0
+LVZ        II       67.8979     34.6514     630.0   200.0
+MBAR       II       -0.6019     30.7382    1390.0   100.0
+MSEY       II       -4.6737     55.4792     475.0    91.0
+MSVF       II      -17.7448    178.0528     801.1   100.0
+NIL        II       33.6506     73.2686     629.0    68.0
+NNA        II      -11.9875    -76.8422     575.0    40.0
+NRIL       II       69.5049     88.4414      92.0   506.0
+NVS        II       54.8404     83.2346     150.0     0.0
+OBN        II       55.1146     36.5674     160.0    30.0
+PALK       II        7.2728     80.7022     460.0    90.0
+RAYN       II       23.5225     45.5032     631.0     2.0
+RPN        II      -27.1267   -109.3344     110.0     0.0
+SACV       II       14.9702    -23.6085     387.0    97.0
+SHEL       II      -15.9588     -5.7457     537.0    60.0
+SUR        II      -32.3797     20.8117    1770.0     0.0
+TAU        II      -42.9099    147.3204     132.0     0.0
+TLY        II       51.6807    103.6438     579.0    20.0
+UOSS       II       24.9453     56.2042     284.4     0.0
+WRAB       II      -19.9336    134.3600     366.0   100.0
+XPFO       II       33.6107   -116.4555    1280.0     0.0
+AAE        IU        9.0292     38.7656    2442.0     0.0
+ADK        IU       51.8823   -176.6842     130.0     0.0
+AFI        IU      -13.9093   -171.7773     705.0     1.0
+ANMO       IU       34.9459   -106.4572    1720.0   100.0
+ANTO       IU       39.8680     32.7934     895.0   195.0
+BBSR       IU       32.3713    -64.6963      -1.3    31.4
+BILL       IU       68.0653    166.4531     320.0     0.0
+BOCO       IU        4.5869    -74.0432    3137.0    23.0
+CASY       IU      -66.2792    110.5354       5.0     5.0
+CCM        IU       38.0557    -91.2446     171.0    51.0
+CHTO       IU       18.8141     98.9443     420.0     0.0
+COLA       IU       64.8736   -147.8616     200.0     0.0
+COR        IU       44.5855   -123.3046     110.0     0.0
+CTAO       IU      -20.0882    146.2545     320.0    37.0
+DAV        IU        7.0697    125.5791     149.0     1.0
+DWPF       IU       28.1103    -81.4327    -132.0   162.0
+FUNA       IU       -8.5259    179.1966      19.0     1.0
+FURI       IU        8.8952     38.6798    2565.0     5.0
+GNI        IU       40.1480     44.7410    1509.0   100.0
+GRFO       IU       49.6909     11.2203     384.0   116.0
+GUMO       IU       13.5893    144.8684      61.0   109.0
+HKT        IU       29.9618    -95.8384    -413.0   450.0
+HNR        IU       -9.4387    159.9475       0.0   100.0
+HRV        IU       42.5064    -71.5583     200.0     0.0
+INCN       IU       37.4776    126.6239      79.0     1.0
+JOHN       IU       16.7329   -169.5292     -37.0    39.0
+KBS        IU       78.9154     11.9385      90.0     3.0
+KEV        IU       69.7565     27.0035      85.0    15.0
+KIEV       IU       50.7012     29.2242     140.0    40.0
+KIP        IU       21.4233   -158.0150      37.0    33.0
+KMBO       IU       -1.1271     37.2525    1930.0    20.0
+KNTN       IU       -2.7744   -171.7186      18.0     2.0
+KONO       IU       59.6491      9.5982     216.0   340.0
+KOWA       IU       14.4967     -4.0140     316.0     5.0
+LCO        IU      -29.0110    -70.7004    2300.0     0.0
+LSZ        IU      -15.2779     28.1882    1200.0     0.0
+LVC        IU      -22.6127    -68.9111    2930.0    30.0
+MA2        IU       59.5756    150.7700     337.0     2.0
+MAJO       IU       36.5457    138.2041     405.0     0.0
+MAKZ       IU       46.8080     81.9770     590.0    10.0
+MBWA       IU      -21.1590    119.7313     181.0     9.0
+MIDW       IU       28.2155   -177.3697      17.8     1.0
+MSKU       IU       -1.6557     13.6116     287.0    25.0
+NAI        IU       -1.2739     36.8037    1692.0     0.0
+NWAO       IU      -32.9277    117.2390     370.9     9.1
+OTAV       IU        0.2376    -78.4508    3495.0    15.0
+PAB        IU       39.5446     -4.3499     950.0     0.0
+PAYG       IU       -0.6742    -90.2861     170.0   100.0
+PET        IU       53.0233    158.6499     105.0     5.0
+PMG        IU       -9.4047    147.1597      90.0     0.0
+PMSA       IU      -64.7744    -64.0489      40.0     0.0
+POHA       IU       19.7573   -155.5326    1910.0    80.0
+PTCN       IU      -25.0713   -130.0953     218.0     2.0
+PTGA       IU       -0.7308    -59.9666     141.0    96.0
+QSPA       IU      -89.9289    144.4382    2847.0     3.0
+RAIO       IU       46.0403   -122.8851      11.0     0.0
+RAO        IU      -29.2450   -177.9290      59.5     0.5
+RAR        IU      -21.2125   -159.7733     -72.0   100.0
+RCBR       IU       -5.8274    -35.9014     291.0   109.0
+RSSD       IU       44.1212   -104.0359    2022.7    67.3
+SAML       IU       -8.9489    -63.1831     120.0     0.0
+SBA        IU      -77.8492    166.7572      48.0     2.0
+SDV        IU        8.8839    -70.6340    1588.0    32.0
+SFJD       IU       66.9961    -50.6208     329.0     1.0
+SJG        IU       18.1091    -66.1500     420.0     0.0
+SLBS       IU       23.6858   -109.9443     825.0     0.0
+SNZO       IU      -41.3087    174.7043      20.0   100.0
+SSPA       IU       40.6358    -77.8876     170.0   100.0
+TARA       IU        1.3549    172.9229      19.0     1.0
+TATO       IU       24.9735    121.4971      77.1    82.9
+TBT        IU       28.6794    -17.9145     180.0    40.0
+TEIG       IU       20.2263    -88.2763      15.0    25.0
+TIXI       IU       71.6341    128.8667      40.0     0.0
+TOL        IU       39.8814     -4.0485     480.0     0.0
+TRIS       IU      -37.0681    -12.3152      58.0     2.0
+TRQA       IU      -38.0568    -61.9787     439.0   101.0
+TSUM       IU      -19.2022     17.5838    1260.0     0.0
+TUC        IU       32.3098   -110.7847     909.0     1.0
+ULN        IU       47.8651    107.0532    1610.0     0.0
+WAKE       IU       19.2834    166.6520      19.0     1.0
+WCI        IU       38.2289    -86.2939      78.0   132.0
+WVT        IU       36.1297    -87.8300     170.0     0.0
+XMAS       IU        2.0448   -157.4457      19.0     1.0
+YAK        IU       62.0310    129.6805     110.0    14.0
+YSS        IU       46.9587    142.7604     148.0     2.0

Added: seismo/3D/SPECFEM3D_GLOBE/trunk/INSTALL
===================================================================
--- seismo/3D/SPECFEM3D_GLOBE/trunk/INSTALL	                        (rev 0)
+++ seismo/3D/SPECFEM3D_GLOBE/trunk/INSTALL	2010-02-20 20:40:52 UTC (rev 16285)
@@ -0,0 +1,188 @@
+
+!=======================================================================!
+!                                                                       !
+!   specfem3D is a 3-D spectral-element solver for the Earth.           !
+!   It uses a mesh generated by program meshfem3D                       !
+!                                                                       !
+!=======================================================================!
+
+Instructions on how to install and use SPECFEM3D_GLOBE are
+available in the manual located in directory USERS_MANUAL.
+
+Main developers: Dimitri Komatitsch, Jeroen Tromp, Qinya Liu and David Michea.
+
+
+
+Basic installation
+------------------
+
+1. untar the SPECFEM3D_GLOBE package: 
+
+   > tar -zxvf SPECFEM3D_GLOBE*.tar.gz 
+
+   in the following, we will assume you set the root directory of the 
+   package name to:
+   
+   SPECFEM3D_GLOBE/
+   
+2. configure the software for your system:
+
+    > cd SPECFEM3D_GLOBE/
+    > ./configure
+
+    by default, it uses gfortran as a Fortran compiler, mpif90 for MPI compilation
+    and gcc as a C compiler. 
+    
+    in order to use a different compiler, use for example:
+    
+        > ./configure FC=ifort
+  
+    in this case, it would make use of
+    
+    - ifort, Intel Fortran Compiler 
+      (see http://software.intel.com/en-us/intel-compilers/ )
+      with which we see very good performance results.   
+        
+    - a default MPI installation which provides mpif90 as fortran wrapper command.     
+      An excellent package is provided by Open MPI (see http://www.open-mpi.org/).
+      
+      for the example above, you would make sure that the mpif90 command 
+      would also use ifort as a Fortran compiler.
+    
+    - a C compiler provided from the GNU compiler collection (see http://gcc.gnu.org/).
+    
+    in case `configure` was successful, it will create the files:
+
+    ./Makefile
+    ./constants.h
+    ./precision.h
+    
+    please make sure, your installation is working and that 
+    the created files 'constant.h' and 'Makefile' satisfy
+    you needs. 
+
+    more information is given in the manual provided in USERS_MANUAL.
+
+3. compile the package:
+
+   > make
+   
+   
+well done!
+
+
+
+
+Configure options
+-----------------
+
+'configure' accepts the following options:
+
+FC                          Fortran90 compiler command name, 
+                            default is 'gfortran'
+ 
+MPIFC                       MPI Fortran90 command name, 
+                            default is 'mpif90'
+
+CC                          C compiler command name,
+                            default is 'gcc'
+
+FLAGS_CHECK                 Compiler flags for non-critical subroutines, 
+                            default is '-O3'
+                  
+FLAGS_NO_CHECK              Compiler flags for creating fast, production-run code for 
+                            critical subroutines,
+                            default is '-O3'
+
+LOCAL_PATH_IS_ALSO_GLOBAL   Directory LOCAL_PATH (specified in DATA/Par_file) on the
+                            cluster can be accessed by all compute nodes ('true') or not ('false'),
+                            default is 'true'
+
+
+and flags:
+
+
+--enable-double-precision   The package can run either in single or in double precision.
+                            default is single precision 
+                            
+--help                      Directs configure to print a usage screen 
+                            
+
+run with 
+  ./configure <option>=<my_option> --<flag>
+  
+
+Compiler specific flags are set in file 'flags.guess' in order to provide
+suitable default options for different compilers. 
+
+you might however want to modify the created 'Makefile' and specify your best 
+system compiler flags in order to ensure optimal performance of the code. 
+
+
+
+
+
+Troobleshooting
+---------------
+
+FAQ
+
+1. configuration fails:
+
+   Examine the log file 'config.log'. It contains detailed informations.
+   in many cases, the path's to these specific compiler commands F90, 
+   CC and MPIF90 won't be correct if `configure` fails.
+   
+   please make sure that you have a working installation of a Fortran compiler, 
+   a C compiler and an MPI implementation. you should be able to compile this
+   little program code:
+   
+      program main
+        include 'mpif.h'
+        integer, parameter :: CUSTOM_MPI_TYPE = MPI_REAL
+        integer ier
+        call MPI_INIT(ier)
+        call MPI_BARRIER(MPI_COMM_WORLD,ier)
+        call MPI_FINALIZE(ier)
+      end
+
+    
+2. compilation fails stating :
+    ... 
+    In file ./compute_element_properties.f90:44
+
+      use meshfem3D_models_par
+                         1
+    Fatal Error: File 'meshfem3d_models_par.mod' opened at (1) is not a GFORTRAN module file
+    ...
+    
+  make sure, you're pointing to the right 'mpif90' wrapper command. 
+  
+  normally, this message will appear when you're mixing two different fortran
+  compilers. that is, using e.g. gfortran to compile non-MPI files 
+  and mpif90, wrapper provided for e.g. ifort, to compile MPI-files. 
+  the module will be created by the wrapper, thus ifort, while the gfortran compiler 
+  is trying to read that module for the compilation of 'compute_element_properties.f90'.
+  
+  fix: e.g. specify > ./configure FC=gfortran MPIF90=/usr/local/openmpi-gfortran/bin/mpif90
+  
+  
+3. compilation fails with:
+    ...
+    obj/specfem3D.o: In function `MAIN__':
+    specfem3D.f90:(.text+0xb66): relocation truncated to fit: R_X86_64_32S against `.bss'
+    ...
+    
+  you're probably using some resolution settings in 'DATA/Par_file' which are 
+  too big for a single processor on your system. the solver tries to statically allocate arrays,
+  which can not be handled by a 32-bit address anymore in this case.
+  
+  fix: check the static memory needed by the solver, which is outputted
+       when you run: > make xcreate_header_file
+                     > ./xcreate_header_file
+       
+       the size of static arrays per slice has to fit on to a single processor.
+  
+       you can either try to use e.g. -mcmodel=medium as additional compiler flag,
+       or better, change your resolution settings for NPROC_XI,NPROC_ETA and NEX_XI,NEX_ETA.
+       

Modified: seismo/3D/SPECFEM3D_GLOBE/trunk/README_SPECFEM3D_GLOBE
===================================================================
--- seismo/3D/SPECFEM3D_GLOBE/trunk/README_SPECFEM3D_GLOBE	2010-02-20 19:58:48 UTC (rev 16284)
+++ seismo/3D/SPECFEM3D_GLOBE/trunk/README_SPECFEM3D_GLOBE	2010-02-20 20:40:52 UTC (rev 16285)
@@ -11,3 +11,108 @@
 
 Main developers: Dimitri Komatitsch, Jeroen Tromp, Qinya Liu and David Michea.
 
+
+Installation
+------------
+
+In order to compile the package on your system, please consult the
+manual located in directory USERS_MANUAL and the file './INSTALL' within this
+package directory.
+
+
+Default simulation
+------------------
+
+A short step-by-step procedure to run a default simulation is described
+here. 
+
+1. To set up a default simulation, you need to have the files:
+
+   DATA/Par_file 
+   DATA/CMTSOLUTION   
+   DATA/STATIONS
+
+   we will use the default simulation files '*.default' provided in the package:
+   > cd SPECFEM3D_GLOBE/
+   > cp DATA/Par_file.default DATA/Par_file
+   > cp DATA/CMTSOLUTION.default DATA/CMTSOLUTION
+   > cp DATA/STATIONS.default DATA/STATIONS
+   
+   The default simulation is an earthquake that occurred in 
+   Bolivia 1994 with a moment magnitude of ~6.8.
+   
+   The simulation will run over the entire globe and require a total of 150 CPUs.
+   The seismogram resolution will be accurate down to ~18 seconds.
+   
+
+2. compile the mesher:
+
+   > make xmeshfem3D
+
+
+3. run the mesher on your cluster:
+
+   > qsub go_mesher_pbs.bash
+   
+   - Different cluster scripts are provided as examples in the directory
+      UTILS/Cluster. Locate your script and copy it the package
+      root directory: 
+      e.g. > cp UTILS/Cluster/pbs/go_mesher_pbs.bash ./
+      
+      The example above uses a version of the Portable Batch System (PBS).
+      Example scripts are also provided for Load Sharing Facility (LSF)
+      and Sun Grid Engine (SGE). make sure, you modify this script to
+      fit to your specific cluster setup.
+  
+   - please, before you submit your job, make sure to have created the directory specified as 
+      'LOCAL_PATH' in your DATA/Par_file. by default the DATA/Par_file specifies:
+          ..
+          # path to store the local database files on each node
+          LOCAL_PATH                      = ./DATABASES_MPI
+          ..
+          
+      make sure this points to the right location for your cluster and that
+      you created the directory before submitting your job:
+      e.g. > mkdir ./DATABASES_MPI  
+      
+      (at this point, you should also know if your directory for LOCAL_PATH
+       is accessible by all compute nodes or not. set the flag
+       LOCAL_PATH_IS_ALSO_GLOBAL in constants.h accordingly )
+       
+   after a successful run, the mesher will have created a perfectly, load-balanced 
+   hexahedral spectral-element mesh stored in this directory.
+
+   note: for the default example, the mesher will take about 30 minutes
+         to create the mesh.
+         the file disk space needed for this global mesh is about ~51 GB.
+         
+     
+4. compile the solver:
+
+  > make xspecfem3D
+  
+  re-compilation is needed since the solver tries to use static memory allocation
+  and compiler optimizations that will work out much better once the details of
+  your simulation mesh are known. 
+
+
+5. run the solver on your cluster:
+
+  > qsub go_solver_pbs.bash
+  
+  - Again, different cluster scripts are provided in the directory 
+    UTILS/Cluster. modify the example scripts to fit to your specific 
+    cluster setup.
+        
+  after a successful simulation, all seismograms are stored into the directory:
+  OUTPUT_FILES/
+    
+  note: the default simulation will take about 1 h 40 m.
+        the file disk space for all seismograms in OUTPUT_FILES/ is about ~77 MB.
+
+
+
+well done!
+
+
+   
\ No newline at end of file

Deleted: seismo/3D/SPECFEM3D_GLOBE/trunk/UTILS/Visualization/Paraview/global_slice_number.daniel.pl
===================================================================
--- seismo/3D/SPECFEM3D_GLOBE/trunk/UTILS/Visualization/Paraview/global_slice_number.daniel.pl	2010-02-20 19:58:48 UTC (rev 16284)
+++ seismo/3D/SPECFEM3D_GLOBE/trunk/UTILS/Visualization/Paraview/global_slice_number.daniel.pl	2010-02-20 20:40:52 UTC (rev 16285)
@@ -1,179 +0,0 @@
-#!/usr/bin/perl -w
-
-#  This program figures out the global slice number for given simulation parameters
-#  uses xglobal_slice_number and xnormal_plane, compile first from global_slice_util/
-# Qinya Liu, May 2007, Caltech
-
-# modify the following line for the correct location of perl libs (UTILS/lib)
-use lib '/opt/seismo-util/lib/perl';
-use CMT_TOOLS;
-
-if (@ARGV != 3) {die("Usage: global_slice_number.pl CMTSOLUTION STATIONS_ADJOINT Par_file\n");}
-$cmt = $ARGV[0];
-$sta = $ARGV[1];
-$par_file = $ARGV[2];
-
-# obtain event location
-($elat,$elon) = get_cmt_location($cmt);
-
-# obtain station location
- at lat = ();
- at lon = ();
-open(IN,$sta);
- at stations = <IN>;
-close(IN);
-$nstations = @stations;
-
-for($i=0;$i<$nstations;$i++){
-  $line = $stations[$i];
-  (undef,undef,$slat,$slon,undef,undef) = split(" ",$line);
-  push(@lat,$slat);
-  push(@lon,$slon);
-
-  print "event = ($elat,$elon);  station = ($lat[$i],$lon[$i]) \n";
-
-}
-#($slat,$slon) = split(" ",`awk 'NR == 1 {print \$3, \$4}' $sta`);
-#print "event = ($elat,$elon);  station = ($slat,$slon) \n";
-
-# obtain simulation parameters
-if (not -f $par_file) {die("Check if $par_file exists or not\n");}
-($nchunk) = split(" ",`grep NCHUNKS $par_file|awk '{print \$3}'`);
-($nproc) = split(" ",`grep NPROC_XI $par_file|awk '{print \$3}'`);
-print "NCHUNKS = $nchunk, nproc = $nproc\n";
-if ($nchunk != 6) {
-  ($xi_width) = split(" ",`grep ANGULAR_WIDTH_XI_IN_DEGREES $par_file|awk '{print \$3}'`);
-  ($eta_width) = split(" ",`grep ANGULAR_WIDTH_ETA_IN_DEGREES $par_file|awk '{print \$3}'`);
-  ($clat) = split(" ",`grep CENTER_LATITUDE_IN_DEGREES $par_file|awk '{print \$3}'`);
-  ($clon) = split(" ",`grep CENTER_LONGITUDE_IN_DEGREES $par_file|awk '{print \$3}'`);
-  ($grot) = split(" ",`grep GAMMA_ROTATION_AZIMUTH $par_file|awk '{print \$3}'`);
-  print "xi_width = $xi_width, eta_width = $eta_width; clat = $clat, clon = $clon; grot = $grot\n";
-}
-
-# all
-open(SLICE_ALL,">slices_all");
-for ($i = 0; $i < $nproc*$nproc*$nchunk; $i++ ) {
-    print SLICE_ALL "$i\n";
-}
-close(SLICE_ALL);  
-
-# minor 
-print "compute slices along minor arc ...\n";
-
-for($istation=0;$istation<$nstations;$istation++){
-  $slat = $lat[$istation];
-  $slon = $lon[$istation];
-  
-  if ($nchunk == 6) {
-    $command = "./xglobal_slice_number $elon $elat $slon $slat $nproc 0";
-  }else {
-    $command = "./xglobal_slice_number $elon $elat $slon $slat $nproc 0 $nchunk $xi_width $eta_width $clon $clat $grot";}
-  
-  if( $istation == 0 ){
-    system(" $command > slices_minor ");
-  }else {
-    system(" $command >> slices_minor ");  
-  }
-}
-
-#sorts out doubles
-open(IN,"slices_minor");
- at slices = <IN>;
-close(IN);
-$nslices = @slices;
- at uni = ();
-for($i=0;$i<$nslices;$i++){
-  $sl = $slices[$i];
-  $found = 0;
-  for($j=0;$j<@uni;$j++){
-    if( $sl == $uni[$j] ){ 
-      $found = 1; 
-      last;
-    }
-  }  
-  if( $found == 0 ){
-    push(@uni,$sl);
-  }
-}
-open(OUT,">slices_minor");
-for($i=0;$i<@uni;$i++){
-  print OUT "$uni[$i]";
-}
-close(OUT);
-
-# major
-print "compute slices along major arc ...\n";
-
-for($istation=0;$istation<$nstations;$istation++){
-  $slat = $lat[$istation];
-  $slon = $lon[$istation];
-  
-  if ($nchunk == 6) {
-    $command = "./xglobal_slice_number $elon $elat $slon $slat $nproc 1";
-  }else {
-    $command = "./xglobal_slice_number $elon $elat $slon $slat $nproc 1 $nchunk $xi_width $eta_width $clon $clat $grot";}
-
-  if( $istation == 0 ){
-    system(" $command > slices_major ");
-  }else {
-    system(" $command >> slices_major ");  
-  }
-}
-
-#sorts out doubles
-open(IN,"slices_major");
- at slices = <IN>;
-close(IN);
-$nslices = @slices;
- at uni = ();
-for($i=0;$i<$nslices;$i++){
-  $sl = $slices[$i];
-  $found = 0;
-  for($j=0;$j<@uni;$j++){
-    if( $sl == $uni[$j] ){ 
-      $found = 1; 
-      last;
-    }
-  }  
-  if( $found == 0 ){
-    push(@uni,$sl);
-  }
-}
-open(OUT,">slices_major");
-for($i=0;$i<@uni;$i++){
-  print OUT "$uni[$i]";
-}
-close(OUT);
-
-
-
-
-
-if ($nchunk == 6) {
-  print "output slices_ic file\n";
-  $nslice = $nproc * $nproc;
-  `cat slices_minor slices_major | sort -g | awk '\$1 < $nslice {print \$1}' > slice_ab_old`;
-
-  open(SLICE,"slice_ab_old");
-  @tmp=<SLICE>;
-  close(SLICE);
-
-  open(SLICE_IC,">slices_ic");
-  for ($i = 0; $i < $nslice; $i ++ ) {
-    for ($j = 0; $j < @tmp; $j ++ ){
-      if ($i == $tmp[$j]) {last;}
-      if ($j == @tmp - 1) {print SLICE_IC "$i\n";}
-    }
-  }
-  close(SLICE_IC);
-}
-
-# figure out the normal to the source-receiver plane
-print "calculate normal to source and receiver plane\n";
-$command = "./xnormal_plane $elat $elon $slat $slon";
-$result = `$command`;
-print "$result\n";
-
-# clean up temp files
-system(" rm -f slice_ab_old gcarc_station.txt xsection_translate.txt");
-



More information about the CIG-COMMITS mailing list