[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