[CIG-SEISMO] SPECFEM3D possible bug in internal mesher

[Valentin Pavlov]

Dear SPECFEM3D Developers,

I think there is a bug in the meshfem3D program, both in the V2.0.1 
distribution, and in the SVN development version.

Please look at SVN, meshfem3D.f90, as of version 18219, line 626. It reads:

         icornerlat = int((lat - orig_y_interface_bottom) / 
spacing_*x*_interface_bottom) + 1

I think this should read

         icornerlat = int((lat - orig_y_interface_bottom) / 
spacing_*y*_interface_bottom) + 1

instead.

best regards,
- Dr. Valentin Pavlov
CTO, VPR&D
Rila Solutions EAD
1113 Sofia, Bulgaria
Acad. G. Bonchev str., bl. 27
+359 2 970 96 100
http://www.rila.com

-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://geodynamics.org/pipermail/cig-seismo/attachments/20120113/524da373/attachment.htm