[cig-commits] r8463 - in seismo/2D/SPECFEM2D/trunk: . DATA
walter at geodynamics.org
walter at geodynamics.org
Fri Dec 7 15:48:08 PST 2007
Author: walter
Date: 2007-12-07 15:48:07 -0800 (Fri, 07 Dec 2007)
New Revision: 8463
Modified:
seismo/2D/SPECFEM2D/trunk/DATA/Par_file
seismo/2D/SPECFEM2D/trunk/DATA/Par_file_Paul_acoustic
seismo/2D/SPECFEM2D/trunk/DATA/Par_file_Paul_elastic
seismo/2D/SPECFEM2D/trunk/DATA/Par_file_cours_M2_UPPA
seismo/2D/SPECFEM2D/trunk/meshfem2D.f90
seismo/2D/SPECFEM2D/trunk/specfem2D.f90
Log:
added output_postscript_image output_PNM_image flags to Par_file in 2D code to suppress
postscript snapshots and/or PNM images to make code run faster when doing numerous runs
Modified: seismo/2D/SPECFEM2D/trunk/DATA/Par_file
===================================================================
--- seismo/2D/SPECFEM2D/trunk/DATA/Par_file 2005-02-19 14:43:28 UTC (rev 8462)
+++ seismo/2D/SPECFEM2D/trunk/DATA/Par_file 2007-12-07 23:48:07 UTC (rev 8463)
@@ -62,6 +62,8 @@
# display parameters
#
itaff = 100 ! display frequency in time steps
+output_postscript_image = .true. ! output Postscript image of the results
+output_PNM_image = .true. ! output PNM image of the results
vecttype = 1 ! display 1=displ 2=veloc 3=accel
cutvect = 1. ! amplitude min en % pour vector plots
meshvect = .true. ! display mesh on vector plots or not
Modified: seismo/2D/SPECFEM2D/trunk/DATA/Par_file_Paul_acoustic
===================================================================
--- seismo/2D/SPECFEM2D/trunk/DATA/Par_file_Paul_acoustic 2005-02-19 14:43:28 UTC (rev 8462)
+++ seismo/2D/SPECFEM2D/trunk/DATA/Par_file_Paul_acoustic 2007-12-07 23:48:07 UTC (rev 8463)
@@ -62,6 +62,8 @@
# display parameters
#
itaff = 100 ! display frequency in time steps
+output_postscript_image = .true. ! output Postscript image of the results
+output_PNM_image = .true. ! output PNM image of the results
vecttype = 2 ! display 1=displ 2=veloc 3=accel
cutvect = 1. ! amplitude min en % pour vector plots
meshvect = .true. ! display mesh on vector plots or not
Modified: seismo/2D/SPECFEM2D/trunk/DATA/Par_file_Paul_elastic
===================================================================
--- seismo/2D/SPECFEM2D/trunk/DATA/Par_file_Paul_elastic 2005-02-19 14:43:28 UTC (rev 8462)
+++ seismo/2D/SPECFEM2D/trunk/DATA/Par_file_Paul_elastic 2007-12-07 23:48:07 UTC (rev 8463)
@@ -62,6 +62,8 @@
# display parameters
#
itaff = 100 ! display frequency in time steps
+output_postscript_image = .true. ! output Postscript image of the results
+output_PNM_image = .true. ! output PNM image of the results
vecttype = 1 ! display 1=displ 2=veloc 3=accel
cutvect = 1. ! amplitude min en % pour vector plots
meshvect = .true. ! display mesh on vector plots or not
Modified: seismo/2D/SPECFEM2D/trunk/DATA/Par_file_cours_M2_UPPA
===================================================================
--- seismo/2D/SPECFEM2D/trunk/DATA/Par_file_cours_M2_UPPA 2005-02-19 14:43:28 UTC (rev 8462)
+++ seismo/2D/SPECFEM2D/trunk/DATA/Par_file_cours_M2_UPPA 2007-12-07 23:48:07 UTC (rev 8463)
@@ -62,6 +62,8 @@
# display parameters
#
itaff = 100 ! display frequency in time steps
+output_postscript_image = .true. ! output Postscript image of the results
+output_PNM_image = .true. ! output PNM image of the results
vecttype = 1 ! display 1=displ 2=veloc 3=accel
cutvect = 1. ! amplitude min en % pour vector plots
meshvect = .true. ! display mesh on vector plots or not
Modified: seismo/2D/SPECFEM2D/trunk/meshfem2D.f90
===================================================================
--- seismo/2D/SPECFEM2D/trunk/meshfem2D.f90 2005-02-19 14:43:28 UTC (rev 8462)
+++ seismo/2D/SPECFEM2D/trunk/meshfem2D.f90 2007-12-07 23:48:07 UTC (rev 8463)
@@ -61,7 +61,7 @@
integer ngnod,nt,nx,nz,nxread,nzread
integer icodematread
- logical codehaut,codebas,codegauche,codedroite
+ logical codehaut,codebas,codegauche,codedroite,output_postscript_image,output_PNM_image
double precision tang1,tangN,vpzone,vszone,poisson_ratio
double precision cutvect,xspacerec,zspacerec
@@ -241,6 +241,8 @@
! read display parameters
call read_value_integer(IIN_PAR,IGNORE_JUNK,itaff)
+ call read_value_logical(IIN_PAR,IGNORE_JUNK,output_postscript_image)
+ call read_value_logical(IIN_PAR,IGNORE_JUNK,output_PNM_image)
call read_value_integer(IIN_PAR,IGNORE_JUNK,vecttype)
call read_value_double_precision(IIN_PAR,IGNORE_JUNK,cutvect)
call read_value_logical(IIN_PAR,IGNORE_JUNK,meshvect)
@@ -605,8 +607,8 @@
write(15,*) 'gnuplot interpol'
write(15,*) gnuplot,interpol
- write(15,*) 'itaff colors numbers'
- write(15,*) itaff,' 1 0'
+ write(15,*) 'itaff output_postscript_image output_PNM_image colors numbers'
+ write(15,*) itaff,output_postscript_image,output_PNM_image,' 1 0'
write(15,*) 'meshvect modelvect boundvect cutvect subsamp nx_sem_PNM'
write(15,*) meshvect,modelvect,boundvect,cutvect,subsamp,nxread
Modified: seismo/2D/SPECFEM2D/trunk/specfem2D.f90
===================================================================
--- seismo/2D/SPECFEM2D/trunk/specfem2D.f90 2005-02-19 14:43:28 UTC (rev 8462)
+++ seismo/2D/SPECFEM2D/trunk/specfem2D.f90 2007-12-07 23:48:07 UTC (rev 8463)
@@ -133,7 +133,7 @@
integer numat,ngnod,nspec,iptsdisp,nelemabs,nelemsurface
logical interpol,meshvect,modelvect,boundvect,read_external_model,initialfield,abshaut, &
- outputgrid,gnuplot,ELASTIC,TURN_ANISOTROPY_ON,TURN_ATTENUATION_ON
+ outputgrid,gnuplot,ELASTIC,TURN_ANISOTROPY_ON,TURN_ATTENUATION_ON,output_postscript_image,output_PNM_image
double precision cutvect,anglerec,xirec,gammarec
@@ -230,7 +230,7 @@
read(IIN,*) gnuplot,interpol
read(IIN,40) datlin
- read(IIN,*) itaff,colors,numbers
+ read(IIN,*) itaff,output_postscript_image,output_PNM_image,colors,numbers
read(IIN,40) datlin
read(IIN,*) meshvect,modelvect,boundvect,cutvect,subsamp,nx_sem_PNM
@@ -1628,6 +1628,8 @@
!
!---- affichage postscript
!
+ if(output_postscript_image) then
+
write(IOUT,*) 'Dump PostScript'
! for elastic medium
@@ -1691,9 +1693,13 @@
endif
write(IOUT,*) 'Fin dump PostScript'
+ endif
+
!
!---- affichage image PNM
!
+ if(output_PNM_image) then
+
write(IOUT,*) 'Creation image PNM de taille ',NX_IMAGE_PNM,' x ',NZ_IMAGE_PNM
donnees_image_PNM_2D(:,:) = 0.d0
@@ -1722,6 +1728,8 @@
write(IOUT,*) 'Fin creation image PNM'
+ endif
+
!---- save temporary seismograms
call write_seismograms(sisux,sisuz,station_name,network_name,NSTEP,nrec,deltat,sismostype,st_xval,it,t0)
More information about the cig-commits
mailing list