[cig-commits] [commit] pluggable: Moved the computation of rotation matrix forward. The rotation matrix is needed when calling compute_coord_main_mesh(). (94d54c9)

[cig_noreply at geodynamics.org]

Repository : ssh://geoshell/specfem3d_globe

On branch  : pluggable
Link       : https://github.com/geodynamics/specfem3d_globe/compare/64e1b38f0c5ebb4056cce0b15d41c0b9f94ab6e5...099a4d330d5b173b21e51ad441f9f429e5d37842

>---------------------------------------------------------------

commit 94d54c9be6b1f036890f1457b8d0b89112208d80
Author: Eh Tan <tan2 at earth.sinica.edu.tw>
Date:   Wed Feb 25 18:26:01 2009 +0000

    Moved the computation of rotation matrix forward. The rotation matrix is needed when calling compute_coord_main_mesh().
    
    Also created 3 common blocks for passing parameters to citcoms earth model. The common blocks are temporary and will be removed later.


>---------------------------------------------------------------

94d54c9be6b1f036890f1457b8d0b89112208d80
 meshfem3D.f90 | 16 +++++++++++-----
 1 file changed, 11 insertions(+), 5 deletions(-)

diff --git a/meshfem3D.f90 b/meshfem3D.f90
index c2d8f86..e5128a9 100644
--- a/meshfem3D.f90
+++ b/meshfem3D.f90
@@ -380,6 +380,11 @@
   logical :: CUT_SUPERBRICK_XI,CUT_SUPERBRICK_ETA
   integer, dimension(MAX_NUM_REGIONS) :: NGLOB1D_RADIAL_TEMP
 
+  ! parameters for passing information to citcoms mantle model
+  common /for_citcoms1/ NPROC_XI, NPROC_ETA, NCHUNKS
+  common /for_citcoms2/ ANGULAR_WIDTH_XI_RAD, ANGULAR_WIDTH_ETA_RAD
+  common /for_citcoms3/ rotation_matrix
+
 ! ************** PROGRAM STARTS HERE **************
 
 ! initialize the MPI communicator and start the NPROCTOT MPI processes.
@@ -855,6 +860,12 @@
   write(IMAIN,*)
 
   endif
+
+! compute rotation matrix from Euler angles
+  ANGULAR_WIDTH_XI_RAD = ANGULAR_WIDTH_XI_IN_DEGREES * PI / 180.d0
+  ANGULAR_WIDTH_ETA_RAD = ANGULAR_WIDTH_ETA_IN_DEGREES * PI / 180.d0
+  if(NCHUNKS /= 6) call euler_angles(rotation_matrix,CENTER_LONGITUDE_IN_DEGREES,CENTER_LATITUDE_IN_DEGREES,GAMMA_ROTATION_AZIMUTH)
+
   if(ELLIPTICITY) call make_ellipticity(nspl,rspl,espl,espl2,ONE_CRUST)
 
   call bcast_model(myrank, trim(LOCAL_PATH))
@@ -911,11 +922,6 @@
     write(IMAIN,*) 'Central cube is at a radius of ',R_CENTRAL_CUBE/1000.d0,' km'
   endif
 
-! compute rotation matrix from Euler angles
-  ANGULAR_WIDTH_XI_RAD = ANGULAR_WIDTH_XI_IN_DEGREES * PI / 180.d0
-  ANGULAR_WIDTH_ETA_RAD = ANGULAR_WIDTH_ETA_IN_DEGREES * PI / 180.d0
-  if(NCHUNKS /= 6) call euler_angles(rotation_matrix,CENTER_LONGITUDE_IN_DEGREES,CENTER_LATITUDE_IN_DEGREES,GAMMA_ROTATION_AZIMUTH)
-
 ! volume of the slice
   volume_total = ZERO