https://github.com/geodynamics/citcoms
Revision 08ab04ba711da9da93c8ac383be32c2b0df1a213 authored by Eh Tan on 19 June 2008, 01:04:06 UTC, committed by Eh Tan on 19 June 2008, 01:04:06 UTC
Use E->surf_det, instead of the surface jacobian returned by get_global_1d_shape_fn(), to expand spherical harmonics. 

The surface jacobian returned by get_global_1d_shape_fn() is defined at the top radius of the current mesh. When nprocz>1, the top radius is not the same as the outer radius of the sphere for processors. As a result, the spherical harmonics expansion gave incorrect result when nprocz>1 in the previous revisions.

1 parent d3025a8
History
Tip revision: 08ab04ba711da9da93c8ac383be32c2b0df1a213 authored by Eh Tan on 19 June 2008, 01:04:06 UTC
Compute E->surf_det (top surface area, ie. surface jacobian, of an element) on a unit sphere.
Tip revision: 08ab04b
File Mode Size
checkpoint
Makefile.am -rw-r--r-- 1.0 KB
array2d.cc -rw-r--r-- 3.9 KB
coupled.cfg -rw-r--r-- 1.8 KB
exchange.py -rwxr-xr-x 4.9 KB
gnomonic.c -rw-r--r-- 2.2 KB
multicoupled.cfg -rw-r--r-- 2.3 KB
signon.py -rwxr-xr-x 2.4 KB
test1.sh -rwxr-xr-x 2.3 KB
test2.sh -rwxr-xr-x 4.2 KB
test5.sh -rwxr-xr-x 2.6 KB

back to top