[cig-commits] r8609 - seismo/2D/SPECFEM2D/trunk

[walter at geodynamics.org]

Author: walter
Date: 2007-12-07 16:00:52 -0800 (Fri, 07 Dec 2007)
New Revision: 8609

Modified:
   seismo/2D/SPECFEM2D/trunk/Makefile
   seismo/2D/SPECFEM2D/trunk/create_color_image.f90
   seismo/2D/SPECFEM2D/trunk/datim.f90
   seismo/2D/SPECFEM2D/trunk/locate_receivers.F90
Log:
added IBM xlf compiler.

Modified: seismo/2D/SPECFEM2D/trunk/Makefile
===================================================================
--- seismo/2D/SPECFEM2D/trunk/Makefile	2007-12-04 02:59:46 UTC (rev 8608)
+++ seismo/2D/SPECFEM2D/trunk/Makefile	2007-12-08 00:00:52 UTC (rev 8609)
@@ -65,6 +65,12 @@
 FLAGS_NOCHECK = -std=gnu -fimplicit-none -frange-check -O2 -Wunused-labels -Waliasing -Wampersand -Wsurprising -Wline-truncation -Wunderflow
 FLAGS_CHECK = $(FLAGS_NOCHECK) -fbounds-check
 
+# IBM
+#F90 = xlf_r
+#CC = xlc -q64
+#FLAGS_NOCHECK = -qextname=attenuation_compute_param -O3 -qsave -qstrict -q64 -qtune=ppc970 -qarch=ppc64v -qcache=auto -qfree=f90 -Q -qsuffix=f=f90 -qhalt=w -qflttrap=en:ov:zero:inv -qfullpath -qsigtrap
+#FLAGS_CHECK = $(FLAGS_NOCHECK) -qddim
+
 LINK = $(F90)
 
 #LIB = /opt/metis-4.0/gcc64/lib/libmetis.a /opt/scotch-4.0/gcc64/lib/libscotch.a  /opt/scotch-4.0/gcc64/lib/libscotcherr.a

Modified: seismo/2D/SPECFEM2D/trunk/create_color_image.f90
===================================================================
--- seismo/2D/SPECFEM2D/trunk/create_color_image.f90	2007-12-04 02:59:46 UTC (rev 8608)
+++ seismo/2D/SPECFEM2D/trunk/create_color_image.f90	2007-12-08 00:00:52 UTC (rev 8609)
@@ -151,9 +151,18 @@
 
 ! compute maximum amplitude
   amplitude_max = maxval(abs(color_image_2D_data))
-  vpmin = minval(image_color_vp_display)
-  vpmax = maxval(image_color_vp_display)
-
+  vpmin = HUGEVAL
+  vpmax = TINYVAL
+  do iy=1,NY
+    do ix=1,NX
+      if ( iglob_image_color_2D(ix,iy) > -1 ) then
+        vpmin = min(vpmin,image_color_vp_display(ix,iy))
+        vpmax = max(vpmax,image_color_vp_display(ix,iy))
+      endif
+      
+    enddo
+  enddo
+  
 ! in the PNM format, the image starts in the upper-left corner
   do iy=NY,1,-1
     do ix=1,NX

Modified: seismo/2D/SPECFEM2D/trunk/datim.f90
===================================================================
--- seismo/2D/SPECFEM2D/trunk/datim.f90	2007-12-04 02:59:46 UTC (rev 8608)
+++ seismo/2D/SPECFEM2D/trunk/datim.f90	2007-12-08 00:00:52 UTC (rev 8609)
@@ -54,8 +54,8 @@
   character(len=16) dateprint
   character(len=8) timeprint
 
-  datein = ''
-  timein = ''
+  datein = ' '
+  timein = ' '
 
   call date_and_time(datein,timein)
 

Modified: seismo/2D/SPECFEM2D/trunk/locate_receivers.F90
===================================================================
--- seismo/2D/SPECFEM2D/trunk/locate_receivers.F90	2007-12-04 02:59:46 UTC (rev 8608)
+++ seismo/2D/SPECFEM2D/trunk/locate_receivers.F90	2007-12-08 00:00:52 UTC (rev 8609)
@@ -238,11 +238,9 @@
 call MPI_BCAST(which_proc_receiver(1),nrec,MPI_INTEGER,0,MPI_COMM_WORLD,ierror)
 
 #else
-print *, 'POY', nproc
-call flush(6)
+
 gather_final_distance(:,1) = final_distance(:)
-print *, 'POY'
-call flush(6)
+
 gather_xi_receiver(:,1) = xi_receiver(:)
 gather_gamma_receiver(:,1) = gamma_receiver(:)
 gather_ispec_selected_rec(:,1) = ispec_selected_rec(:)