[cig-commits] r17090 - seismo/3D/SPECFEM3D/trunk/meshfem3D

[pieyre at geodynamics.org]

Author: pieyre
Date: 2010-08-16 00:49:05 -0700 (Mon, 16 Aug 2010)
New Revision: 17090

Modified:
   seismo/3D/SPECFEM3D/trunk/meshfem3D/read_parameter_file.f90
Log:
minor modification in the mesher concerning the doublings positions


Modified: seismo/3D/SPECFEM3D/trunk/meshfem3D/read_parameter_file.f90
===================================================================
--- seismo/3D/SPECFEM3D/trunk/meshfem3D/read_parameter_file.f90	2010-08-16 05:04:25 UTC (rev 17089)
+++ seismo/3D/SPECFEM3D/trunk/meshfem3D/read_parameter_file.f90	2010-08-16 07:49:05 UTC (rev 17090)
@@ -79,7 +79,7 @@
 ! second dimension : #rho  #vp  #vs  #Q_flag  #anisotropy_flag #domain_id
   double precision, dimension(:,:), pointer :: material_properties
 
-  integer i,ireg,imat
+  integer i,ireg,imat,idoubl
 
 ! open parameter file
   open(unit=IIN,file='DATA/Par_file',status='old',action='read')
@@ -136,6 +136,14 @@
   call read_value_integer(IIN,IGNORE_JUNK,ner_doublings(2), 'mesher.NZ_DOUGLING_2')
   if(err_occurred() /= 0) return
 
+  if(ner_doublings(1) < ner_doublings(2) .and. NDOUBLINGS == 2) then
+    idoubl = ner_doublings(1)
+    ner_doublings(1) = ner_doublings(2)
+    ner_doublings(2) = idoubl
+  end if
+
+
+
   call read_value_logical(IIN,IGNORE_JUNK,CREATE_ABAQUS_FILES, 'mesher.CREATE_ABAQUS_FILES')
   if(err_occurred() /= 0) return
   call read_value_logical(IIN,IGNORE_JUNK,CREATE_DX_FILES, 'mesher.CREATE_DX_FILES')