[cig-commits] r22959 - seismo/3D/SPECFEM3D_GLOBE/trunk/src/meshfem3D

lefebvre at geodynamics.org lefebvre at geodynamics.org
Thu Oct 10 09:10:18 PDT 2013


Author: lefebvre
Date: 2013-10-10 09:10:17 -0700 (Thu, 10 Oct 2013)
New Revision: 22959

Modified:
   seismo/3D/SPECFEM3D_GLOBE/trunk/src/meshfem3D/create_regions_mesh.F90
Log:
ipass allocation for global simulation fixed in create_regions_mesh.

Modified: seismo/3D/SPECFEM3D_GLOBE/trunk/src/meshfem3D/create_regions_mesh.F90
===================================================================
--- seismo/3D/SPECFEM3D_GLOBE/trunk/src/meshfem3D/create_regions_mesh.F90	2013-10-10 14:58:33 UTC (rev 22958)
+++ seismo/3D/SPECFEM3D_GLOBE/trunk/src/meshfem3D/create_regions_mesh.F90	2013-10-10 16:10:17 UTC (rev 22959)
@@ -631,22 +631,24 @@
   normal_bottom(:,:,:,:) = 0.0; normal_top(:,:,:,:) = 0.0
 
   ! Stacey
-  if( ipass == 1 .and. NCHUNKS /= 6 ) then
-    allocate(nimin(2,NSPEC2DMAX_YMIN_YMAX), &
-             nimax(2,NSPEC2DMAX_YMIN_YMAX), &
-             njmin(2,NSPEC2DMAX_XMIN_XMAX), &
-             njmax(2,NSPEC2DMAX_XMIN_XMAX), &
-             nkmin_xi(2,NSPEC2DMAX_XMIN_XMAX), &
-             nkmin_eta(2,NSPEC2DMAX_YMIN_YMAX),stat=ier)
-    if(ier /= 0) stop 'error in allocate 14'
-  else
-    allocate(nimin(1,1), &
-             nimax(1,1), &
-             njmin(1,1), &
-             njmax(1,1), &
-             nkmin_xi(1,1), &
-             nkmin_eta(1,1),stat=ier)
-    if(ier /= 0) stop 'error in allocate 14'
+  if (ipass == 1) then
+    if (NCHUNKS /= 6 ) then
+      allocate(nimin(2,NSPEC2DMAX_YMIN_YMAX), &
+               nimax(2,NSPEC2DMAX_YMIN_YMAX), &
+               njmin(2,NSPEC2DMAX_XMIN_XMAX), &
+               njmax(2,NSPEC2DMAX_XMIN_XMAX), &
+               nkmin_xi(2,NSPEC2DMAX_XMIN_XMAX), &
+               nkmin_eta(2,NSPEC2DMAX_YMIN_YMAX),stat=ier)
+      if(ier /= 0) stop 'error in allocate 14'
+    else
+      allocate(nimin(1,1), &
+               nimax(1,1), &
+               njmin(1,1), &
+               njmax(1,1), &
+               nkmin_xi(1,1), &
+               nkmin_eta(1,1),stat=ier)
+      if(ier /= 0) stop 'error in allocate 14'
+    endif
   endif
 
   nimin(:,:) = 0; nimax(:,:) = 0



More information about the CIG-COMMITS mailing list