Thanks to visit codestin.com
Credit goes to github.com

Skip to content

Commit 6280047

Browse files
author
Daniel Peter
committed
updates headers and version numbers to version 5.1.1
1 parent 771b54b commit 6280047

File tree

184 files changed

+548
-471
lines changed

Some content is hidden

Large Commits have some content hidden by default. Use the searchbox below for content that may be hidden.

184 files changed

+548
-471
lines changed

Makefile.in

Lines changed: 6 additions & 6 deletions
Original file line numberDiff line numberDiff line change
@@ -7,7 +7,7 @@
77
# Princeton University, USA
88
# and University of Pau / CNRS / INRIA, France
99
# (c) Princeton University / California Institute of Technology and University of Pau / CNRS / INRIA
10-
# February 2011
10+
# April 2011
1111
#
1212
# This program is free software; you can redistribute it and/or modify
1313
# it under the terms of the GNU General Public License as published by
@@ -77,13 +77,13 @@ obj:
7777
mkdir -p obj
7878

7979

80-
xmeshfem3D: required
80+
xmeshfem3D: required
8181
(cd src/meshfem3D; make)
8282

83-
xspecfem3D: required
83+
xspecfem3D: required
8484
(cd src/specfem3D; make)
8585

86-
xcreate_header_file: required
86+
xcreate_header_file: required
8787
(cd src/create_header_file; make)
8888

8989
xcheck_buffers_1D: required
@@ -101,10 +101,10 @@ xcheck_buffers_faces_chunks: required
101101
xcombine_AVS_DX: required
102102
(cd src/auxiliaries ; make xcombine_AVS_DX)
103103

104-
xcombine_vol_data: required
104+
xcombine_vol_data: required
105105
(cd src/auxiliaries ; make xcombine_vol_data)
106106

107-
xcombine_surf_data: required
107+
xcombine_surf_data: required
108108
(cd src/auxiliaries ; make xcombine_surf_data)
109109

110110
xconvolve_source_timefunction: required

UTILS/Cluster/estimate_best_values_runs/estimate_best_values_runs.f90

Lines changed: 2 additions & 2 deletions
Original file line numberDiff line numberDiff line change
@@ -1,13 +1,13 @@
11
!=====================================================================
22
!
3-
! S p e c f e m 3 D G l o b e V e r s i o n 5 . 0
3+
! S p e c f e m 3 D G l o b e V e r s i o n 5 . 1
44
! --------------------------------------------------
55
!
66
! Main authors: Dimitri Komatitsch and Jeroen Tromp
77
! Princeton University, USA
88
! and University of Pau / CNRS / INRIA, France
99
! (c) Princeton University / California Institute of Technology and University of Pau / CNRS / INRIA
10-
! March 2010
10+
! April 2011
1111
!
1212
! This program is free software; you can redistribute it and/or modify
1313
! it under the terms of the GNU General Public License as published by

UTILS/Cluster/test_cache_size.f90

Lines changed: 1 addition & 1 deletion
Original file line numberDiff line numberDiff line change
@@ -1,6 +1,6 @@
11
!=====================================================================
22
!
3-
! S p e c f e m 3 D G l o b e V e r s i o n 5 . 0
3+
! S p e c f e m 3 D G l o b e V e r s i o n 5 . 1
44
! --------------------------------------------------
55
!
66
! Main authors: Dimitri Komatitsch and Jeroen Tromp

UTILS/Profiles/write_profile.f90

Lines changed: 10 additions & 12 deletions
Original file line numberDiff line numberDiff line change
@@ -1,13 +1,13 @@
11
!=====================================================================
22
!
3-
! S p e c f e m 3 D G l o b e V e r s i o n 5 . 0
3+
! S p e c f e m 3 D G l o b e V e r s i o n 5 . 1
44
! --------------------------------------------------
55
!
66
! Main authors: Dimitri Komatitsch and Jeroen Tromp
77
! Princeton University, USA
88
! and University of Pau / CNRS / INRIA, France
99
! (c) Princeton University / California Institute of Technology and University of Pau / CNRS / INRIA
10-
! March 2010
10+
! April 2011
1111
!
1212
! This program is free software; you can redistribute it and/or modify
1313
! it under the terms of the GNU General Public License as published by
@@ -24,8 +24,6 @@
2424
! 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA.
2525
!
2626
!=====================================================================
27-
!
28-
! United States Government Sponsorship Acknowledged.
2927

3028
!=====================================================================
3129
! write_profile.f90
@@ -35,13 +33,13 @@
3533
! intervals, with discontinuities honored by the mesh and those of the MOHO and
3634
! crust/ocean indicated by repeating points with different values.
3735

38-
! The code shortcuts through meshfem3D -> create_regions_mesh ->
39-
! create_regular_elements -> compute_element_properties -> compute_element_properties ->
40-
! get_model, cutting and pasting the relevant parts.
36+
! The code shortcuts through meshfem3D -> create_regions_mesh ->
37+
! create_regular_elements -> compute_element_properties -> compute_element_properties ->
38+
! get_model, cutting and pasting the relevant parts.
4139

4240
! This code IS NOT MAINTAINED and most likely will not work without modification as
4341
! it calls functions that may have changed since it was last updated.
44-
!
42+
!
4543
! I have tried to indicate where different parts come from to facilitate updates.
4644
! - vala hjorleifsdottir ([email protected])
4745
!=====================================================================
@@ -123,7 +121,7 @@ program xwrite_profile
123121
logical, dimension(MAX_NUMBER_OF_MESH_LAYERS) :: this_region_has_a_doubling
124122
double precision, dimension(MAX_NUMBER_OF_MESH_LAYERS) :: rmins,rmaxs
125123

126-
124+
127125
integer, dimension(NB_SQUARE_CORNERS,NB_CUT_CASE) :: DIFF_NSPEC1D_RADIAL
128126
integer, dimension(NB_SQUARE_EDGES_ONEDIR,NB_CUT_CASE) :: DIFF_NSPEC2D_XI,DIFF_NSPEC2D_ETA
129127
logical :: CUT_SUPERBRICK_XI,CUT_SUPERBRICK_ETA
@@ -160,7 +158,7 @@ program xwrite_profile
160158
logical:: elem_in_crust,elem_in_mantle
161159
!---
162160

163-
!--- from
161+
!--- from
164162
! local parameters
165163
double precision xmesh,ymesh,zmesh
166164
!---
@@ -284,7 +282,7 @@ program xwrite_profile
284282
if(ilayer == 1) then
285283
rmin = 0.0d0
286284
rmax = rmins(NUMBER_OF_MESH_LAYERS-1)
287-
idoubling = IFLAG_INNER_CORE_NORMAL
285+
idoubling = IFLAG_INNER_CORE_NORMAL
288286
else
289287
rmin = rmins(NUMBER_OF_MESH_LAYERS-ilayer+1)
290288
rmax = rmaxs(NUMBER_OF_MESH_LAYERS-ilayer+1)
@@ -368,7 +366,7 @@ program xwrite_profile
368366
if(r >= rmax*0.999999d0) r_prem = rmax*0.999999d0
369367

370368
! convert from rthetaphi to xyz to use in function calls.
371-
369+
372370
call rthetaphi_2_xyz(xmesh,ymesh,zmesh,r_prem,theta,phi)
373371

374372
!!!!-------------------------------------

UTILS/Visualization/Paraview/create_slice_VTK.f90

Lines changed: 40 additions & 39 deletions
Original file line numberDiff line numberDiff line change
@@ -1,13 +1,13 @@
11
!=====================================================================
22
!
3-
! S p e c f e m 3 D G l o b e V e r s i o n 5 . 0
3+
! S p e c f e m 3 D G l o b e V e r s i o n 5 . 1
44
! --------------------------------------------------
55
!
66
! Main authors: Dimitri Komatitsch and Jeroen Tromp
77
! Princeton University, USA
88
! and University of Pau / CNRS / INRIA, France
99
! (c) Princeton University / California Institute of Technology and University of Pau / CNRS / INRIA
10-
! March 2010
10+
! April 2011
1111
!
1212
! This program is free software; you can redistribute it and/or modify
1313
! it under the terms of the GNU General Public License as published by
@@ -25,6 +25,7 @@
2525
!
2626
!=====================================================================
2727

28+
2829
program create_slice_VTK
2930

3031
! this programs creates a vtk file that specifies the (velocity) model values on each element point,
@@ -39,7 +40,7 @@ program create_slice_VTK
3940
! cd to your SPECFEM3D_GLOBE root directory:
4041
! > cd SPECFEM3D_GLOBE/
4142
! create symbolic link:
42-
! > ln -s UTILS/Visualization/Paraview/create_slice_VTK.f90
43+
! > ln -s UTILS/Visualization/Paraview/create_slice_VTK.f90
4344
! compile with:
4445
! > f90 -o xcreate_slice_VTK create_slice_VTK.f90
4546
! run :
@@ -69,7 +70,7 @@ program create_slice_VTK
6970
real x, y, z, dat
7071
integer numpoin, iglob, n1, n2, n3, n4, n5, n6, n7, n8
7172
integer iglob1, iglob2, iglob3, iglob4, iglob5, iglob6, iglob7, iglob8
72-
73+
7374

7475
! starts here--------------------------------------------------------------------------------------------------
7576
do i = 1, 6
@@ -114,7 +115,7 @@ program create_slice_VTK
114115
print*,'no file: ',trim(arg(1))
115116
stop 'Error opening slices file'
116117
endif
117-
118+
118119
do while (1 == 1)
119120
read(20,'(a)',iostat=ios) sline
120121
if (ios /= 0) exit
@@ -155,11 +156,11 @@ program create_slice_VTK
155156
print*,'error ',ios
156157
print*,'file:',trim(dimension_file)
157158
stop 'Error opening file'
158-
endif
159+
endif
159160
read(27,*) nspec(it)
160161
read(27,*) nglob(it)
161162
close(27)
162-
163+
163164

164165
enddo
165166

@@ -173,8 +174,8 @@ program create_slice_VTK
173174
do it = 1, num_node
174175

175176
iproc = node_list(it)
176-
177-
177+
178+
178179
print *, ' '
179180
print *, 'Reading slice ', iproc
180181
write(prname_topo,'(a,i6.6,a,i1,a)') trim(in_topo_dir)//'/proc',iproc,'_reg',ir,'_'
@@ -192,11 +193,11 @@ program create_slice_VTK
192193
data(:,:,:,:) = -1.e10
193194
read(27) data(:,:,:,1:nspec(it))
194195
close(27)
195-
196+
196197
print *,trim(data_file)
197198
print *,' min/max value: ',minval(data(:,:,:,1:nspec(it))),maxval(data(:,:,:,1:nspec(it)))
198199
print *
199-
200+
200201
! topology file
201202
topo_file = trim(prname_topo) // 'solver_data_2' // '.bin'
202203
open(unit = 28,file = trim(topo_file),status='old',action='read', iostat = ios, form='unformatted')
@@ -214,7 +215,7 @@ program create_slice_VTK
214215
read(28) zstore(1:nglob(it))
215216
read(28) ibool(:,:,:,1:nspec(it))
216217
close(28)
217-
218+
218219

219220
write(mesh_file,'(a,i1,a)') trim(outdir)//'/' // 'reg_',ir,'_'//trim(filename)
220221
print *, trim(mesh_file)
@@ -250,12 +251,12 @@ subroutine write_VTK_data_gll_cr(nspec,nglob, &
250251
include "constants.h"
251252

252253
integer :: nspec,nglob
253-
254-
! global coordinates
254+
255+
! global coordinates
255256
integer, dimension(NGLLX,NGLLY,NGLLZ,nspec) :: ibool
256257
real(kind=CUSTOM_REAL), dimension(nglob) :: xstore_dummy,ystore_dummy,zstore_dummy
257258

258-
! gll data values array
259+
! gll data values array
259260
real(kind=CUSTOM_REAL), dimension(NGLLX,NGLLY,NGLLZ,nspec) :: gll_data
260261

261262
! file name
@@ -266,81 +267,81 @@ subroutine write_VTK_data_gll_cr(nspec,nglob, &
266267
! write source and receiver VTK files for Paraview
267268
write(IMAIN,*) ' vtk file: '
268269
write(IMAIN,*) ' ',prname_file(1:len_trim(prname_file))//'.vtk'
269-
270+
270271
open(IOVTK,file=prname_file(1:len_trim(prname_file))//'.vtk',status='unknown')
271272
write(IOVTK,'(a)') '# vtk DataFile Version 3.1'
272273
write(IOVTK,'(a)') 'material model VTK file'
273274
write(IOVTK,'(a)') 'ASCII'
274275
write(IOVTK,'(a)') 'DATASET UNSTRUCTURED_GRID'
275-
276+
276277
! writes out all points for each element, not just global ones
277278
write(IOVTK, '(a,i12,a)') 'POINTS ', nspec*8, ' float'
278279
do ispec=1,nspec
279-
i = ibool(1,1,1,ispec)
280+
i = ibool(1,1,1,ispec)
280281
write(IOVTK,'(3e18.6)') xstore_dummy(i),ystore_dummy(i),zstore_dummy(i)
281-
282+
282283
i = ibool(NGLLX,1,1,ispec)
283284
write(IOVTK,'(3e18.6)') xstore_dummy(i),ystore_dummy(i),zstore_dummy(i)
284-
285+
285286
i = ibool(NGLLX,NGLLY,1,ispec)
286287
write(IOVTK,'(3e18.6)') xstore_dummy(i),ystore_dummy(i),zstore_dummy(i)
287-
288+
288289
i = ibool(1,NGLLY,1,ispec)
289-
write(IOVTK,'(3e18.6)') xstore_dummy(i),ystore_dummy(i),zstore_dummy(i)
290-
290+
write(IOVTK,'(3e18.6)') xstore_dummy(i),ystore_dummy(i),zstore_dummy(i)
291+
291292
i = ibool(1,1,NGLLZ,ispec)
292293
write(IOVTK,'(3e18.6)') xstore_dummy(i),ystore_dummy(i),zstore_dummy(i)
293-
294+
294295
i = ibool(NGLLX,1,NGLLZ,ispec)
295296
write(IOVTK,'(3e18.6)') xstore_dummy(i),ystore_dummy(i),zstore_dummy(i)
296-
297+
297298
i = ibool(NGLLX,NGLLY,NGLLZ,ispec)
298299
write(IOVTK,'(3e18.6)') xstore_dummy(i),ystore_dummy(i),zstore_dummy(i)
299-
300+
300301
i = ibool(1,NGLLY,NGLLZ,ispec)
301-
write(IOVTK,'(3e18.6)') xstore_dummy(i),ystore_dummy(i),zstore_dummy(i)
302+
write(IOVTK,'(3e18.6)') xstore_dummy(i),ystore_dummy(i),zstore_dummy(i)
302303
enddo
303304
write(IOVTK,*) ""
304305

305306
! note: indices for vtk start at 0
306-
write(IOVTK,'(a,i12,i12)') "CELLS ",nspec,nspec*9
307+
write(IOVTK,'(a,i12,i12)') "CELLS ",nspec,nspec*9
307308
do ispec=1,nspec
308309
write(IOVTK,'(9i12)') 8,(ispec-1)*8,(ispec-1)*8+1,(ispec-1)*8+2,(ispec-1)*8+3,&
309310
(ispec-1)*8+4,(ispec-1)*8+5,(ispec-1)*8+6,(ispec-1)*8+7
310311
enddo
311312
write(IOVTK,*) ""
312-
313+
313314
! type: hexahedrons
314315
write(IOVTK,'(a,i12)') "CELL_TYPES ",nspec
315316
write(IOVTK,*) (12,ispec=1,nspec)
316317
write(IOVTK,*) ""
317-
318+
318319
! writes out gll-data (velocity) for each element point
319320
write(IOVTK,'(a,i12)') "POINT_DATA ",nspec*8
320321
write(IOVTK,'(a)') "SCALARS gll_data float"
321322
write(IOVTK,'(a)') "LOOKUP_TABLE default"
322323
do ispec = 1,nspec
323-
i = ibool(1,1,1,ispec)
324+
i = ibool(1,1,1,ispec)
324325
write(IOVTK,'(3e18.6)') gll_data(1,1,1,ispec)
325-
326+
326327
i = ibool(NGLLX,1,1,ispec)
327328
write(IOVTK,'(3e18.6)') gll_data(NGLLX,1,1,ispec)
328-
329+
329330
i = ibool(NGLLX,NGLLY,1,ispec)
330331
write(IOVTK,'(3e18.6)') gll_data(NGLLX,NGLLY,1,ispec)
331-
332+
332333
i = ibool(1,NGLLY,1,ispec)
333-
write(IOVTK,'(3e18.6)') gll_data(1,NGLLY,1,ispec)
334-
334+
write(IOVTK,'(3e18.6)') gll_data(1,NGLLY,1,ispec)
335+
335336
i = ibool(1,1,NGLLZ,ispec)
336337
write(IOVTK,'(3e18.6)') gll_data(1,1,NGLLZ,ispec)
337-
338+
338339
i = ibool(NGLLX,1,NGLLZ,ispec)
339340
write(IOVTK,'(3e18.6)') gll_data(NGLLX,1,NGLLZ,ispec)
340-
341+
341342
i = ibool(NGLLX,NGLLY,NGLLZ,ispec)
342343
write(IOVTK,'(3e18.6)') gll_data(NGLLX,NGLLY,NGLLZ,ispec)
343-
344+
344345
i = ibool(1,NGLLY,NGLLZ,ispec)-1
345346
write(IOVTK,'(3e18.6)') gll_data(1,NGLLY,NGLLZ,ispec)
346347
enddo

UTILS/Visualization/opendx_AVS/constants_topo.h

Lines changed: 1 addition & 1 deletion
Original file line numberDiff line numberDiff line change
@@ -1,6 +1,6 @@
11
!=====================================================================
22
!
3-
! S p e c f e m 3 D G l o b e V e r s i o n 5 . 0
3+
! S p e c f e m 3 D G l o b e V e r s i o n 5 . 1
44
! --------------------------------------------------
55
!
66
! Main authors: Dimitri Komatitsch and Jeroen Tromp

UTILS/Visualization/opendx_AVS/convert_DX_to_AVS_slices.f90

Lines changed: 1 addition & 1 deletion
Original file line numberDiff line numberDiff line change
@@ -1,6 +1,6 @@
11
!=====================================================================
22
!
3-
! S p e c f e m 3 D G l o b e V e r s i o n 5 . 0
3+
! S p e c f e m 3 D G l o b e V e r s i o n 5 . 1
44
! --------------------------------------------------
55
!
66
! Main authors: Dimitri Komatitsch and Jeroen Tromp

UTILS/Visualization/opendx_AVS/convert_DX_to_AVS_surface.f90

Lines changed: 1 addition & 1 deletion
Original file line numberDiff line numberDiff line change
@@ -1,6 +1,6 @@
11
!=====================================================================
22
!
3-
! S p e c f e m 3 D G l o b e V e r s i o n 5 . 0
3+
! S p e c f e m 3 D G l o b e V e r s i o n 5 . 1
44
! --------------------------------------------------
55
!
66
! Main authors: Dimitri Komatitsch and Jeroen Tromp

UTILS/Visualization/opendx_AVS/convert_avs2dx_quads.f90

Lines changed: 1 addition & 1 deletion
Original file line numberDiff line numberDiff line change
@@ -1,6 +1,6 @@
11
!=====================================================================
22
!
3-
! S p e c f e m 3 D G l o b e V e r s i o n 5 . 0
3+
! S p e c f e m 3 D G l o b e V e r s i o n 5 . 1
44
! --------------------------------------------------
55
!
66
! Main authors: Dimitri Komatitsch and Jeroen Tromp

0 commit comments

Comments
 (0)