[cig-commits] r13384 - seismo/2D/SPECFEM2D/branches/BIOT

[cmorency at geodynamics.org]

Author: cmorency
Date: 2008-11-24 12:12:41 -0800 (Mon, 24 Nov 2008)
New Revision: 13384

Modified:
   seismo/2D/SPECFEM2D/branches/BIOT/meshfem2D.F90
   seismo/2D/SPECFEM2D/branches/BIOT/specfem2D.F90
Log:
bug fixed with nxread/nzread



Modified: seismo/2D/SPECFEM2D/branches/BIOT/meshfem2D.F90
===================================================================
--- seismo/2D/SPECFEM2D/branches/BIOT/meshfem2D.F90	2008-11-24 17:07:54 UTC (rev 13383)
+++ seismo/2D/SPECFEM2D/branches/BIOT/meshfem2D.F90	2008-11-24 20:12:41 UTC (rev 13384)
@@ -1257,9 +1257,6 @@
      write(15,*) 'nelemabs nelem_acoustic_surface num_fluid_solid_edges num_fluid_poro_edges num_solid_poro_edges'
      write(15,*) nelemabs_loc,nelem_acoustic_surface_loc,nedges_coupled_loc,nedges_acporo_coupled_loc,nedges_elporo_coupled_loc
 
-     write(15,*) 'nxread, nzread'
-     write(15,*) nxread,nzread
-
      write(15,*) 'Material sets Isotropic (Anisotropic: to be defined)'
      do i=1,nb_materials
     write(15,*) i,icodemat(i),rho_s(i),rho_f(i),phi(i),tortuosity(i),permxx(i),permxz(i),permzz(i),kappa_s(i),&

Modified: seismo/2D/SPECFEM2D/branches/BIOT/specfem2D.F90
===================================================================
--- seismo/2D/SPECFEM2D/branches/BIOT/specfem2D.F90	2008-11-24 17:07:54 UTC (rev 13383)
+++ seismo/2D/SPECFEM2D/branches/BIOT/specfem2D.F90	2008-11-24 20:12:41 UTC (rev 13384)
@@ -256,7 +256,7 @@
   double precision :: cutsnaps,sizemax_arrows,anglerec,xirec,gammarec
 
 ! for absorbing and acoustic free surface conditions
-  integer :: ispec_acoustic_surface,inum,numabsread,nxread,nzread
+  integer :: ispec_acoustic_surface,inum,numabsread
   logical :: codeabsread(4)
   real(kind=CUSTOM_REAL) :: nx,nz,weight,xxi,zgamma
 
@@ -679,8 +679,6 @@
   read(IIN,*) numat,ngnod,nspec,pointsdisp,plot_lowerleft_corner_only
   read(IIN,"(a80)") datlin
   read(IIN,*) nelemabs,nelem_acoustic_surface,num_fluid_solid_edges,num_fluid_poro_edges,num_solid_poro_edges
-  read(IIN,"(a80)") datlin
-  read(IIN,*) nxread,nzread
 !
 !---- allocate arrays
 !
@@ -723,25 +721,25 @@
   allocate(numabs(nelemabs))
   allocate(codeabs(4,nelemabs))
 
-  allocate(ibegin_bottom(nxread))
-  allocate(iend_bottom(nxread))
-  allocate(ibegin_top(nxread))
-  allocate(iend_top(nxread))
+  allocate(ibegin_bottom(nelemabs))
+  allocate(iend_bottom(nelemabs))
+  allocate(ibegin_top(nelemabs))
+  allocate(iend_top(nelemabs))
 
-  allocate(jbegin_left(nzread))
-  allocate(jend_left(nzread))
-  allocate(jbegin_right(nzread))
-  allocate(jend_right(nzread))
+  allocate(jbegin_left(nelemabs))
+  allocate(jend_left(nelemabs))
+  allocate(jbegin_right(nelemabs))
+  allocate(jend_right(nelemabs))
 
-  allocate(ibegin_bottom_poro(nxread))
-  allocate(iend_bottom_poro(nxread))
-  allocate(ibegin_top_poro(nxread))
-  allocate(iend_top_poro(nxread))
+  allocate(ibegin_bottom_poro(nelemabs))
+  allocate(iend_bottom_poro(nelemabs))
+  allocate(ibegin_top_poro(nelemabs))
+  allocate(iend_top_poro(nelemabs))
 
-  allocate(jbegin_left_poro(nzread))
-  allocate(jend_left_poro(nzread))
-  allocate(jbegin_right_poro(nzread))
-  allocate(jend_right_poro(nzread))
+  allocate(jbegin_left_poro(nelemabs))
+  allocate(jend_left_poro(nelemabs))
+  allocate(jbegin_right_poro(nelemabs))
+  allocate(jend_right_poro(nelemabs))
 
 !
 !---- print element group main parameters
@@ -880,10 +878,10 @@
     nspec_xmax = ZERO
     nspec_zmin = ZERO
     nspec_zmax = ZERO
-    allocate(ib_xmin(nzread))
-    allocate(ib_xmax(nzread))
-    allocate(ib_zmin(nxread))
-    allocate(ib_zmax(nxread))
+    allocate(ib_xmin(nelemabs))
+    allocate(ib_xmax(nelemabs))
+    allocate(ib_zmin(nelemabs))
+    allocate(ib_zmax(nelemabs))
     do inum = 1,nelemabs
        if (codeabs(IBOTTOM,inum)) then
          nspec_zmin = nspec_zmin + 1