[cig-commits] [commit] devel, master: Small spelling fixes in comments. (8711be6)
cig_noreply at geodynamics.org
cig_noreply at geodynamics.org
Thu Nov 6 08:20:28 PST 2014
Repository : https://github.com/geodynamics/specfem3d_globe
On branches: devel,master
Link : https://github.com/geodynamics/specfem3d_globe/compare/bc58e579b3b0838a0968725a076f5904845437ca...be63f20cbb6f462104e949894dbe205d2398cd7f
>---------------------------------------------------------------
commit 8711be6abc8af12cbcdef6b49191dc61d53ed2a7
Author: Elliott Sales de Andrade <esalesde at physics.utoronto.ca>
Date: Sun Jun 1 03:55:44 2014 -0400
Small spelling fixes in comments.
>---------------------------------------------------------------
8711be6abc8af12cbcdef6b49191dc61d53ed2a7
src/auxiliaries/combine_AVS_DX.f90 | 2 +-
src/meshfem3D/compute_coordinates_grid.f90 | 4 ++--
src/meshfem3D/compute_volumes_and_areas.F90 | 2 +-
src/meshfem3D/model_attenuation.f90 | 6 +++---
src/meshfem3D/model_ppm.f90 | 8 ++++----
src/shared/rthetaphi_xyz.f90 | 8 ++++----
src/specfem3D/compute_kernels.F90 | 4 ++--
src/specfem3D/locate_sources.f90 | 10 +++++-----
src/specfem3D/setup_sources_receivers.f90 | 4 ++--
9 files changed, 24 insertions(+), 24 deletions(-)
diff --git a/src/auxiliaries/combine_AVS_DX.f90 b/src/auxiliaries/combine_AVS_DX.f90
index 78da285..43a93b7 100644
--- a/src/auxiliaries/combine_AVS_DX.f90
+++ b/src/auxiliaries/combine_AVS_DX.f90
@@ -701,7 +701,7 @@
! compute total area consisting of oceans
! and suppress areas that are not considered oceans if needed
-! use arbitrary threshold to suppress artefacts in ETOPO5 model
+! use arbitrary threshold to suppress artifacts in ETOPO5 model
if(xval >= -0.018) then
if(OCEANS_ONLY) xval = 0.
above_zero = above_zero + 1
diff --git a/src/meshfem3D/compute_coordinates_grid.f90 b/src/meshfem3D/compute_coordinates_grid.f90
index e8b463c..429e885 100644
--- a/src/meshfem3D/compute_coordinates_grid.f90
+++ b/src/meshfem3D/compute_coordinates_grid.f90
@@ -162,7 +162,7 @@
rgt = (r_top / R_EARTH)*gamma
rgb = (r_bottom / R_EARTH)*gamma
- ! define the mesh points on the top and the bottom in the six regions of the cubed shpere
+ ! define the mesh points on the top and the bottom in the six regions of the cubed sphere
select case (ichunk)
case(CHUNK_AB)
@@ -428,7 +428,7 @@
! keeping in mind that the code non-dimensionalizes the radius of the spherical Earth to one
rgt = R_UNIT_SPHERE*gamma
- ! define the mesh points on the top and the bottom in the six regions of the cubed shpere
+ ! define the mesh points on the top and the bottom in the six regions of the cubed sphere
select case (ichunk)
case(CHUNK_AB)
diff --git a/src/meshfem3D/compute_volumes_and_areas.F90 b/src/meshfem3D/compute_volumes_and_areas.F90
index d5eafec..0bbce07 100644
--- a/src/meshfem3D/compute_volumes_and_areas.F90
+++ b/src/meshfem3D/compute_volumes_and_areas.F90
@@ -377,7 +377,7 @@
weight = wxgll(i)*wygll(j)*wzgll(k)
- ! compute the jacobian
+ ! compute the Jacobian
xixl = xixstore(i,j,k,ispec)
xiyl = xiystore(i,j,k,ispec)
xizl = xizstore(i,j,k,ispec)
diff --git a/src/meshfem3D/model_attenuation.f90 b/src/meshfem3D/model_attenuation.f90
index 0b3abf3..a7ac26c 100644
--- a/src/meshfem3D/model_attenuation.f90
+++ b/src/meshfem3D/model_attenuation.f90
@@ -162,7 +162,7 @@
! files rather than be defined within the CODE
!
! All this subroutine does is define the Attenuation vs Radius and then Compute the Attenuation
-! Variables (tau_sigma and tau_epslion ( or tau_mu) )
+! Variables (tau_sigma and tau_epsilon ( or tau_mu) )
subroutine model_attenuation_setup(myrank,REFERENCE_1D_MODEL,RICB,RCMB, &
R670,R220,R80,AM_V,AM_S,AS_V,CRUSTAL)
@@ -749,7 +749,7 @@ subroutine attenuation_simplex_setup(nf_in,nsls_in,f_in,Q_in,tau_s_in,AS_V)
! dimension(nf)
! tau_s = Tau_sigma Stress relaxation time (see References)
! dimension(nsls)
-! tau_e = Tau_epislon Strain relaxation time (see References)
+! tau_e = Tau_epsilon Strain relaxation time (see References)
! dimension(nsls)!
! Output
! B = Real Moduli ( M2 Dahlen and Tromp pp.203 )
@@ -814,7 +814,7 @@ subroutine attenuation_simplex_setup(nf_in,nsls_in,f_in,Q_in,tau_s_in,AS_V)
!
! Uses attenuation_simplex_variables to store constant values
!
-! See atteunation_simplex_setup
+! See attenuation_simplex_setup
!
double precision function attenuation_eval(Xin,AS_V)
diff --git a/src/meshfem3D/model_ppm.f90 b/src/meshfem3D/model_ppm.f90
index 3300677..41b3571 100644
--- a/src/meshfem3D/model_ppm.f90
+++ b/src/meshfem3D/model_ppm.f90
@@ -922,7 +922,7 @@
y0 = yl(i,j,k,ispec)
z0 = zl(i,j,k,ispec)
- ! calculate weights based on gaussian smoothing
+ ! calculate weights based on Gaussian smoothing
call smoothing_weights_vec(x0,y0,z0,ispec2,sigma_h2,sigma_v2,exp_val,&
xx(:,:,:,ispec2),yy(:,:,:,ispec2),zz(:,:,:,ispec2))
@@ -941,7 +941,7 @@
tk_rhovp(i,j,k,ispec) = tk_rhovp(i,j,k,ispec) + sum(exp_val(:,:,:) * ks_rhovp(:,:,:,ispec2))
tk_rhovs(i,j,k,ispec) = tk_rhovs(i,j,k,ispec) + sum(exp_val(:,:,:) * ks_rhovs(:,:,:,ispec2))
- ! normalization, integrated values of gaussian smoothing function
+ ! normalization, integrated values of Gaussian smoothing function
bk(i,j,k,ispec) = bk(i,j,k,ispec) + sum(exp_val(:,:,:))
enddo
@@ -1086,7 +1086,7 @@
! +(zz(:,:,:,ispec2)-z0)**2 )/(2*sigma2) )*factor(:,:,:)
! from basin code smoothing:
- ! gaussian function
+ ! Gaussian function
!exp_val(:,:,:) = exp( -(xx(:,:,:,ispec2)-x0)**2/(sigma_h2) &
! -(yy(:,:,:,ispec2)-y0)**2/(sigma_h2) &
! -(zz(:,:,:,ispec2)-z0)**2/(sigma_v2) ) * factor(:,:,:)
@@ -1101,7 +1101,7 @@
call get_distance_vec(dist_h,dist_v,x0,y0,z0, &
xx_elem(ii,jj,kk),yy_elem(ii,jj,kk),zz_elem(ii,jj,kk))
- ! gaussian function
+ ! Gaussian function
exp_val(ii,jj,kk) = exp( - dist_h*dist_h/sigma_h2 &
- dist_v*dist_v/sigma_v2 ) ! * factor(ii,jj,kk)
diff --git a/src/shared/rthetaphi_xyz.f90 b/src/shared/rthetaphi_xyz.f90
index 5fc0365..16c2008 100644
--- a/src/shared/rthetaphi_xyz.f90
+++ b/src/shared/rthetaphi_xyz.f90
@@ -265,14 +265,14 @@
! - moment tensor:
! the moment tensor of the CMTSOLUTION is rotated from a spherical to a
-! cartesian representation, using the geocentric colatitude
+! Cartesian representation, using the geocentric colatitude
! - vertical component:
! the vertical component at source/receiver locations is calculated with
! respect to geocentric colatitude
! - target position in the mesh:
! to find the best source/receiver position in the mesh, the target
! location is using the geocentric colatitude and corrects the radius for
-! topography and ellipticity when calculating the x/y/z cartesian
+! topography and ellipticity when calculating the x/y/z Cartesian
! position. We then locate this x/y/z location in the mesh.
! Note that we always convert the geographic latitude to geocentric
@@ -285,11 +285,11 @@
! mesh, which should be fine. however, for the elliptical case where
! ELLIPTICITY is turn on, the mesh would be stretched out and x/y/z
! positions of the mesh nodes are then what? geocentric or geographic now?
-! and what is the cartesian position of a location given in geographic
+! and what is the Cartesian position of a location given in geographic
! coordinates? since we map topography onto geocentric colatitude and
! stretch out the ellipticity in radial direction, i think we always have
! to convert first the source/receiver positions to geocentric ones to
-! then locate the best cartesian x/y/z location in the mesh. so, i think,
+! then locate the best Cartesian x/y/z location in the mesh. so, i think,
! we are doing it the right way.
! Not sure if i made this point clear, and i would welcome some more
diff --git a/src/specfem3D/compute_kernels.F90 b/src/specfem3D/compute_kernels.F90
index c20b4dd..9fb28fc 100644
--- a/src/specfem3D/compute_kernels.F90
+++ b/src/specfem3D/compute_kernels.F90
@@ -142,7 +142,7 @@
! reason for this is that the adjoint wavefield is in general smoother
! since the adjoint sources normally are obtained for filtered traces.
! numerically, the time derivative by a finite-difference scheme should
- ! behave better for smoother wavefields, thus containing less numerical artefacts.
+ ! behave better for smoother wavefields, thus containing less numerical artifacts.
rho_kl_crust_mantle(i,j,k,ispec) = rho_kl_crust_mantle(i,j,k,ispec) &
+ deltat * (accel_crust_mantle(1,iglob) * b_displ_crust_mantle(1,iglob) &
+ accel_crust_mantle(2,iglob) * b_displ_crust_mantle(2,iglob) &
@@ -191,7 +191,7 @@
! reason for this is that the adjoint wavefield is in general smoother
! since the adjoint sources normally are obtained for filtered traces.
! numerically, the time derivative by a finite-difference scheme should
- ! behave better for smoother wavefields, thus containing less numerical artefacts.
+ ! behave better for smoother wavefields, thus containing less numerical artifacts.
rho_kl_crust_mantle(i,j,k,ispec) = rho_kl_crust_mantle(i,j,k,ispec) &
+ deltat * (accel_crust_mantle(1,iglob) * b_displ_crust_mantle(1,iglob) &
+ accel_crust_mantle(2,iglob) * b_displ_crust_mantle(2,iglob) &
diff --git a/src/specfem3D/locate_sources.f90 b/src/specfem3D/locate_sources.f90
index 24cdcf0..2143508 100644
--- a/src/specfem3D/locate_sources.f90
+++ b/src/specfem3D/locate_sources.f90
@@ -400,7 +400,7 @@
enddo
enddo
- ! calculates a gaussian mask around source point
+ ! calculates a Gaussian mask around source point
if( SAVE_SOURCE_MASK .and. SIMULATION_TYPE == 3 ) then
call calc_mask_source(mask_source,ispec,NSPEC,typical_size, &
x_target_source,y_target_source,z_target_source, &
@@ -782,7 +782,7 @@
x_target_source,y_target_source,z_target_source, &
ibool,xstore,ystore,zstore,NGLOB)
-! calculate a gaussian function mask in the crust_mantle region
+! calculate a Gaussian function mask in the crust_mantle region
! which is 0 around the source locations and 1 everywhere else
use constants
@@ -802,7 +802,7 @@
integer i,j,k,iglob
double precision dist_sq,sigma_sq
- ! standard deviation for gaussian
+ ! standard deviation for Gaussian
! (removes factor 10 added for search radius from typical_size)
sigma_sq = typical_size * typical_size / 100.0
@@ -817,7 +817,7 @@
+(y_target_source - dble(ystore(iglob)))**2 &
+(z_target_source - dble(zstore(iglob)))**2
- ! adds gaussian function value to mask
+ ! adds Gaussian function value to mask
! (mask value becomes 0 closer to source location, 1 everywhere else )
mask_source(i,j,k,ispec) = mask_source(i,j,k,ispec) &
* ( 1.0_CUSTOM_REAL - exp( - dist_sq / sigma_sq ) )
@@ -939,7 +939,7 @@
t0 = USER_T0
endif
endif
- ! convert the half duration for triangle STF to the one for gaussian STF
+ ! convert the half duration for triangle STF to the one for Gaussian STF
! note: this calculation here is only used for outputting the plot_source_time_function file
! (see setup_sources_receivers.f90)
hdur_gaussian(:) = hdur(:)/SOURCE_DECAY_MIMIC_TRIANGLE
diff --git a/src/specfem3D/setup_sources_receivers.f90 b/src/specfem3D/setup_sources_receivers.f90
index 7f53e6d..d1b79bf 100644
--- a/src/specfem3D/setup_sources_receivers.f90
+++ b/src/specfem3D/setup_sources_receivers.f90
@@ -172,7 +172,7 @@
endif
endif
- ! convert the half duration for triangle STF to the one for gaussian STF
+ ! convert the half duration for triangle STF to the one for Gaussian STF
hdur_gaussian(:) = hdur(:)/SOURCE_DECAY_MIMIC_TRIANGLE
! define t0 as the earliest start time
@@ -182,7 +182,7 @@
if( USE_FORCE_POINT_SOURCE ) then
! note: point force sources will give the dominant frequency in hdur,
! thus the main period is 1/hdur.
- ! also, these sources use a Ricker source time function instead of a gaussian.
+ ! also, these sources use a Ricker source time function instead of a Gaussian.
! for a Ricker source time function, a start time ~1.2 * main_period is a good choice
t0 = - 1.2d0 * minval(tshift_cmt(:) - 1.0d0/hdur(:))
endif
More information about the CIG-COMMITS
mailing list