[cig-commits] r18831 - in seismo/3D/SPECFEM3D/tags/v1.4.4_last_BASIN: . DATA/m16_model

[carltape at geodynamics.org]

Author: carltape
Date: 2011-08-16 17:22:05 -0700 (Tue, 16 Aug 2011)
New Revision: 18831

Modified:
   seismo/3D/SPECFEM3D/tags/v1.4.4_last_BASIN/DATA/m16_model/Par_file_m16
   seismo/3D/SPECFEM3D/tags/v1.4.4_last_BASIN/DATA/m16_model/README_m16_model
   seismo/3D/SPECFEM3D/tags/v1.4.4_last_BASIN/combine_vol_data.f90
Log:
updated notes for generating m16 socal model


Modified: seismo/3D/SPECFEM3D/tags/v1.4.4_last_BASIN/DATA/m16_model/Par_file_m16
===================================================================
--- seismo/3D/SPECFEM3D/tags/v1.4.4_last_BASIN/DATA/m16_model/Par_file_m16	2011-08-15 23:38:27 UTC (rev 18830)
+++ seismo/3D/SPECFEM3D/tags/v1.4.4_last_BASIN/DATA/m16_model/Par_file_m16	2011-08-17 00:22:05 UTC (rev 18831)
@@ -47,7 +47,7 @@
 HDUR_MOVIE                      = 0.0
 
 # save AVS or OpenDX mesh files to check the mesh
-SAVE_MESH_FILES                 = .false.
+SAVE_MESH_FILES                 = .true.
 
 # path to store the local database file on each node
 LOCAL_PATH                      =  /scratch/ctape.326257

Modified: seismo/3D/SPECFEM3D/tags/v1.4.4_last_BASIN/DATA/m16_model/README_m16_model
===================================================================
--- seismo/3D/SPECFEM3D/tags/v1.4.4_last_BASIN/DATA/m16_model/README_m16_model	2011-08-15 23:38:27 UTC (rev 18830)
+++ seismo/3D/SPECFEM3D/tags/v1.4.4_last_BASIN/DATA/m16_model/README_m16_model	2011-08-17 00:22:05 UTC (rev 18831)
@@ -15,6 +15,7 @@
 
    The Par_file has the option set to check the mesh:
    SAVE_MESH_FILES                 = .true.
+   (For many simulations, you only need to run the mesher once anyway.)
 
 2. Check that 366 (= 168 x 2) .bin files are present in /bin_files/
 

Modified: seismo/3D/SPECFEM3D/tags/v1.4.4_last_BASIN/combine_vol_data.f90
===================================================================
--- seismo/3D/SPECFEM3D/tags/v1.4.4_last_BASIN/combine_vol_data.f90	2011-08-15 23:38:27 UTC (rev 18830)
+++ seismo/3D/SPECFEM3D/tags/v1.4.4_last_BASIN/combine_vol_data.f90	2011-08-17 00:22:05 UTC (rev 18831)
@@ -67,11 +67,13 @@
       print *, ' files have been collected in input_dir, output mesh file goes to output_dir '
       print *, ' give 0 for low resolution and 1 for high resolution'
       stop ' Reenter command line options'
+    else
+      print *, i, trim(arg(i))
     endif
   enddo
 
 ! get slice list
-  if (trim(arg(6)) == '') then
+!  if (trim(arg(6)) == '') then
     num_node = 0
     open(unit = 20, file = trim(arg(1)), status = 'unknown',iostat = ios)
     if (ios /= 0) then
@@ -91,18 +93,18 @@
     indir= arg(3)
     outdir = arg(4)
     read(arg(5),*) ires
-  else
-    read(arg(1),*) proc1
-    read(arg(2),*) proc2
-    do iproc = proc1, proc2
-      node_list(iproc - proc1 + 1) = iproc
-    enddo
-    num_node = proc2 - proc1 + 1
-    filename = arg(3)
-    indir = arg(4)
-    outdir = arg(5)
-    read(arg(6),*) ires
-  endif
+!  else
+!    read(arg(1),*) proc1
+!    read(arg(2),*) proc2
+!    do iproc = proc1, proc2
+!      node_list(iproc - proc1 + 1) = iproc
+!    enddo
+!    num_node = proc2 - proc1 + 1
+!    filename = arg(3)
+!    indir = arg(4)
+!    outdir = arg(5)
+!    read(arg(6),*) ires
+!  endif
 
   if (ires == 0) then
     HIGH_RESOLUTION_MESH = .false.
@@ -314,7 +316,11 @@
       enddo !ispec
     endif
 
-    if (numpoin /= npoint) stop 'Error: number of points are not consistent'
+    if (numpoin /= npoint) then
+       print *, 'numpoin = ',numpoin
+       print *, 'npoint = ',npoint
+       stop 'Error: number of points are not consistent'
+    endif
     np = np + npoint
 
   enddo  ! all slices for points