[cig-commits] r8455 - in seismo/2D/SPECFEM2D/trunk: . DATA

[walter at geodynamics.org]

Author: walter
Date: 2007-12-07 15:47:26 -0800 (Fri, 07 Dec 2007)
New Revision: 8455

Modified:
   seismo/2D/SPECFEM2D/trunk/DATA/Par_file_Paul_acoustic
   seismo/2D/SPECFEM2D/trunk/Makefile
   seismo/2D/SPECFEM2D/trunk/meshfem2D.f90
   seismo/2D/SPECFEM2D/trunk/write_seismograms.f90
Log:
fixed ngnod = 9 bug in 2D code.
also fixed "Ux_file.bin" binary file too large problem in 2D code
(now first erase the file, then recreate it)


Modified: seismo/2D/SPECFEM2D/trunk/DATA/Par_file_Paul_acoustic
===================================================================
--- seismo/2D/SPECFEM2D/trunk/DATA/Par_file_Paul_acoustic	2005-01-25 22:04:11 UTC (rev 8454)
+++ seismo/2D/SPECFEM2D/trunk/DATA/Par_file_Paul_acoustic	2007-12-07 23:47:26 UTC (rev 8455)
@@ -50,9 +50,9 @@
 #
 # receiver line parameters
 #
-enreg_surf                      = .false.         ! enregistrement volume ou surface
+enreg_surf                      = .false.        ! enregistrement volume ou surface
 sismostype                      = 2              ! record 1=displ 2=veloc 3=accel
-nrec                            = 11            ! number of receivers
+nrec                            = 101            ! number of receivers
 xdeb                            = 300.           ! first receiver x in meters
 zdeb                            = 1700.          ! first receiver z in meters
 xfin                            = 3700.          ! last receiver x in meters

Modified: seismo/2D/SPECFEM2D/trunk/Makefile
===================================================================
--- seismo/2D/SPECFEM2D/trunk/Makefile	2005-01-25 22:04:11 UTC (rev 8454)
+++ seismo/2D/SPECFEM2D/trunk/Makefile	2007-12-07 23:47:26 UTC (rev 8455)
@@ -46,7 +46,7 @@
 all: default
 
 clean:
-	/bin/rm -f xmeshfem2D xmeshfem2D.trace xspecfem2D xspecfem2D.trace $O/*.o *.o $O/*.il *.mod core *.gnu *.ps Ux*.bin Uz*.bin image*.pnm xconvolve_source_timefunction *receiver_line_* plotgnu source.txt *.sem*
+	/bin/rm -r -f xmeshfem2D xmeshfem2D.trace xspecfem2D xspecfem2D.trace $O/*.o *.o $O/*.il *.mod core *.gnu *.ps Ux*.bin Uz*.bin image*.pnm xconvolve_source_timefunction *receiver_line_* plotgnu source.txt *.sem* OUTPUT_FILES/*
 
 meshfem2D: $(OBJS_MESHFEM2D)
 	$(LINK) $(FLAGS_CHECK) -o xmeshfem2D $(OBJS_MESHFEM2D)

Modified: seismo/2D/SPECFEM2D/trunk/meshfem2D.f90
===================================================================
--- seismo/2D/SPECFEM2D/trunk/meshfem2D.f90	2005-01-25 22:04:11 UTC (rev 8454)
+++ seismo/2D/SPECFEM2D/trunk/meshfem2D.f90	2007-12-07 23:47:26 UTC (rev 8455)
@@ -173,6 +173,7 @@
   if(ngnod == 9) then
     nx = nx * 2
     nz = nz * 2
+    nz_layer = nz_layer * 2
   endif
 
 ! read absorbing boundaries parameters

Modified: seismo/2D/SPECFEM2D/trunk/write_seismograms.f90
===================================================================
--- seismo/2D/SPECFEM2D/trunk/write_seismograms.f90	2005-01-25 22:04:11 UTC (rev 8454)
+++ seismo/2D/SPECFEM2D/trunk/write_seismograms.f90	2007-12-07 23:47:26 UTC (rev 8455)
@@ -117,6 +117,11 @@
     enddo
   enddo
 
+! delete the old file
+  open(unit=11,file='OUTPUT_FILES/Ux_file.bin',status='unknown')
+  close(11,status='delete')
+
+! write the new file
   open(unit=11,file='OUTPUT_FILES/Ux_file.bin',status='unknown',access='direct',recl=4*NSTEP*nrec)
   write(11,rec=1) buffer_SEP_binary
   close(11)
@@ -130,6 +135,11 @@
     enddo
   enddo
 
+! delete the old file
+  open(unit=11,file='OUTPUT_FILES/Uz_file.bin',status='unknown')
+  close(11,status='delete')
+
+! write the new file
   open(unit=11,file='OUTPUT_FILES/Uz_file.bin',status='unknown',access='direct',recl=4*NSTEP*nrec)
   write(11,rec=1) buffer_SEP_binary
   close(11)