[cig-commits] [commit] : renamed create_movie_AVS_DX.f90 to create_movie_shakemap_AVS_DX_GMT.f90. Updated the acknowledgements section of the manual. (ba5d5e9)
cig_noreply at geodynamics.org
cig_noreply at geodynamics.org
Thu Nov 14 20:16:02 PST 2013
Repository : ssh://geoshell/specfem3d
On branch :
Link : https://github.com/geodynamics/specfem2d/compare/1e201257d91c794056b990a43329e05d04f77454...0000000000000000000000000000000000000000
>---------------------------------------------------------------
commit ba5d5e984009e1f6ddf900bfc484900d22be3a10
Author: Dimitri Komatitsch <komatitsch at lma.cnrs-mrs.fr>
Date: Wed Apr 29 22:50:18 2009 +0000
renamed create_movie_AVS_DX.f90 to create_movie_shakemap_AVS_DX_GMT.f90.
Updated the acknowledgements section of the manual.
>---------------------------------------------------------------
ba5d5e984009e1f6ddf900bfc484900d22be3a10
create_movie_shakemap_AVS_DX_GMT.f90 | 68 +++++++++---------------------------
1 file changed, 16 insertions(+), 52 deletions(-)
diff --git a/create_movie_shakemap_AVS_DX_GMT.f90 b/create_movie_shakemap_AVS_DX_GMT.f90
index 6e155e7..074e9b0 100644
--- a/create_movie_shakemap_AVS_DX_GMT.f90
+++ b/create_movie_shakemap_AVS_DX_GMT.f90
@@ -24,11 +24,12 @@
!=====================================================================
!
-!--- create a movie of vertical component of surface displacement or velocity
-!--- in AVS or OpenDX format
+!--- create a movie of the vertical component of surface displacement or velocity
+!--- or a ShakeMap(R) (i.e. map of the maximum absolute value of the two horizontal components
+!--- of the velocity vector) in AVS, OpenDX or GMT format
!
- program create_movie_AVS_DX
+ program create_movie_shakemap
implicit none
@@ -46,7 +47,7 @@
integer iformat,nframes,iframe,inumber,inorm,iscaling_shake
integer ibool_number,ibool_number1,ibool_number2,ibool_number3,ibool_number4
- logical USE_OPENDX,USE_AVS,USE_GMT,UNIQUE_FILE,plot_shaking_map
+ logical USE_OPENDX,USE_AVS,USE_GMT,plot_shaking_map
real(kind=CUSTOM_REAL), dimension(:,:), allocatable :: x,y,z,display
real(kind=CUSTOM_REAL) xcoord,ycoord,zcoord
@@ -164,34 +165,28 @@
allocate(store_val_uz(ilocnum,0:NPROC-1))
print *,'1 = create files in OpenDX format'
- print *,'2 = create files in AVS UCD format with individual files'
- print *,'3 = create files in AVS UCD format with one time-dependent file'
- print *,'4 = create files in GMT xyz Ascii long/lat/Uz format'
+ print *,'2 = create files in AVS UCD format'
+ print *,'3 = create files in GMT xyz Ascii long/lat/Uz format'
print *,'any other value = exit'
print *
+
print *,'enter value:'
read(5,*) iformat
- if(iformat<1 .or. iformat>4) stop 'exiting...'
+
+ if(iformat<1 .or. iformat>3) stop 'exiting...'
+
if(iformat == 1) then
USE_OPENDX = .true.
USE_AVS = .false.
USE_GMT = .false.
- UNIQUE_FILE = .false.
else if(iformat == 2) then
USE_OPENDX = .false.
USE_AVS = .true.
USE_GMT = .false.
- UNIQUE_FILE = .false.
- else if(iformat == 3) then
- USE_OPENDX = .false.
- USE_AVS = .true.
- USE_GMT = .false.
- UNIQUE_FILE = .true.
else
USE_OPENDX = .false.
USE_AVS = .false.
USE_GMT = .true.
- UNIQUE_FILE = .false.
endif
if(.not. USE_GMT) then
@@ -577,7 +572,6 @@
open(unit=11,file=trim(OUTPUT_FILES)//outputname,status='unknown')
write(11,*) 'object 1 class array type float rank 1 shape 3 items ',nglob,' data follows'
else if(USE_AVS) then
- if(UNIQUE_FILE) stop 'cannot use unique file AVS option for shaking map'
write(outputname,"('/AVS_shaking_map.inp')")
open(unit=11,file=trim(OUTPUT_FILES)//outputname,status='unknown')
write(11,*) nglob,' ',nspectot_AVS_max,' 1 0 0'
@@ -595,17 +589,9 @@
open(unit=11,file=trim(OUTPUT_FILES)//outputname,status='unknown')
write(11,*) 'object 1 class array type float rank 1 shape 3 items ',nglob,' data follows'
else if(USE_AVS) then
- if(UNIQUE_FILE .and. iframe == 1) then
- open(unit=11,file=trim(OUTPUT_FILES)//'/AVS_movie_all.inp',status='unknown')
- write(11,*) nframes
- write(11,*) 'data'
- write(11,"('step',i1,' image',i1)") 1,1
- write(11,*) nglob,' ',nspectot_AVS_max
- else if(.not. UNIQUE_FILE) then
- write(outputname,"('/AVS_movie_',i6.6,'.inp')") ivalue
- open(unit=11,file=trim(OUTPUT_FILES)//outputname,status='unknown')
- write(11,*) nglob,' ',nspectot_AVS_max,' 1 0 0'
- endif
+ write(outputname,"('/AVS_movie_',i6.6,'.inp')") ivalue
+ open(unit=11,file=trim(OUTPUT_FILES)//outputname,status='unknown')
+ write(11,*) nglob,' ',nspectot_AVS_max,' 1 0 0'
else if(USE_GMT) then
write(outputname,"('/gmt_movie_',i6.6,'.xyz')") ivalue
open(unit=11,file=trim(OUTPUT_FILES)//outputname,status='unknown')
@@ -635,9 +621,6 @@
else
-! if unique file, output geometry only once
- if(.not. UNIQUE_FILE .or. iframe == 1) then
-
! output list of points
mask_point = .false.
ipoin = 0
@@ -681,28 +664,11 @@
endif
enddo
- endif
-
if(USE_OPENDX) then
write(11,*) 'attribute "element type" string "quads"'
write(11,*) 'attribute "ref" string "positions"'
write(11,*) 'object 3 class array type float rank 0 items ',nglob,' data follows'
else
- if(UNIQUE_FILE) then
-! step number for AVS multistep file
- if(iframe > 1) then
- if(iframe < 10) then
- write(11,"('step',i1,' image',i1)") iframe,iframe
- else if(iframe < 100) then
- write(11,"('step',i2,' image',i2)") iframe,iframe
- else if(iframe < 1000) then
- write(11,"('step',i3,' image',i3)") iframe,iframe
- else
- write(11,"('step',i4,' image',i4)") iframe,iframe
- endif
- endif
- write(11,*) '1 0'
- endif
! dummy text for labels
write(11,*) '1 1'
write(11,*) 'a, b'
@@ -749,14 +715,12 @@
! end of test for GMT format
endif
- if(.not. UNIQUE_FILE) close(11)
+ close(11)
! end of loop and test on all the time steps for all the movie images
endif
enddo
- if(UNIQUE_FILE) close(11)
-
print *
print *,'done creating movie or shaking map'
print *
@@ -789,7 +753,7 @@
deallocate(display)
endif
- end program create_movie_AVS_DX
+ end program create_movie_shakemap
!
!=====================================================================
More information about the CIG-COMMITS
mailing list