[cig-commits] r15633 - seismo/3D/SPECFEM3D_GLOBE/trunk
dkomati1 at geodynamics.org
dkomati1 at geodynamics.org
Sat Aug 29 09:14:02 PDT 2009
Author: dkomati1
Date: 2009-08-29 09:14:01 -0700 (Sat, 29 Aug 2009)
New Revision: 15633
Modified:
seismo/3D/SPECFEM3D_GLOBE/trunk/create_movie_GMT_global.f90
Log:
fixed a line length problem in a comment that generated a compiler warning
Modified: seismo/3D/SPECFEM3D_GLOBE/trunk/create_movie_GMT_global.f90
===================================================================
--- seismo/3D/SPECFEM3D_GLOBE/trunk/create_movie_GMT_global.f90 2009-08-29 16:08:47 UTC (rev 15632)
+++ seismo/3D/SPECFEM3D_GLOBE/trunk/create_movie_GMT_global.f90 2009-08-29 16:14:01 UTC (rev 15633)
@@ -171,7 +171,7 @@
if(MOVIE_COARSE) then
!daniel: note that
! nex_per_proc_xi*nex_per_proc_eta = nex_xi*nex_eta/nproc = nspec2d_top(iregion_crust_mantle) used in specfem3D.f90
- ! and ilocnum = nmovie_points = 2 * 2 * NEX_XI * NEX_ETA / NPROC
+ ! and ilocnum = nmovie_points = 2 * 2 * NEX_XI * NEX_ETA / NPROC
ilocnum = 2 * 2 * NEX_PER_PROC_XI*NEX_PER_PROC_ETA
NIT =NGLLX-1
else
@@ -375,11 +375,11 @@
ieoff = (ielm+(i-1)+(j-1)*(NGLLX-1))+1
endif
- !daniel: for movie_coarse e.g. x(i,j) is defined at x(1,1), x(1,NGLLY), x(NGLLX,1) and x(NGLLX,NGLLY)
- ! be aware that for the cubed sphere, the mapping changes for different chunks,
- ! i.e. e.g. x(1,1) and x(5,5) flip left and right sides of the elements in geographical coordinates
- if( MOVIE_COARSE) then
- if( NCHUNKS == 6 ) then
+! Daniel: for movie_coarse e.g. x(i,j) is defined at x(1,1), x(1,NGLLY), x(NGLLX,1) and x(NGLLX,NGLLY)
+! be aware that for the cubed sphere, the mapping changes for different chunks,
+! i.e. e.g. x(1,1) and x(5,5) flip left and right sides of the elements in geographical coordinates
+ if(MOVIE_COARSE) then
+ if(NCHUNKS == 6) then
! chunks mapped such that element corners increase in long/lat
select case (iproc/NPROC+1)
case(CHUNK_AB)
@@ -415,11 +415,11 @@
case default
stop 'incorrect chunk number'
end select
- else
+ else
xp(ieoff) = dble(x(1,1))
yp(ieoff) = dble(y(1,1))
zp(ieoff) = dble(z(1,1))
- field_display(ieoff) = dble(displn(1,1))
+ field_display(ieoff) = dble(displn(1,1))
endif ! NCHUNKS
else
xp(ieoff) = dble(x(i,j))
@@ -427,7 +427,7 @@
zp(ieoff) = dble(z(i,j))
field_display(ieoff) = dble(displn(i,j))
endif ! MOVIE_COARSE
-
+
enddo !i
enddo !j
More information about the CIG-COMMITS
mailing list