[cig-commits] r14148 - seismo/3D/SPECFEM3D_GLOBE/branches/pluggable
tan2 at geodynamics.org
tan2 at geodynamics.org
Wed Feb 25 10:26:01 PST 2009
Author: tan2
Date: 2009-02-25 10:26:01 -0800 (Wed, 25 Feb 2009)
New Revision: 14148
Modified:
seismo/3D/SPECFEM3D_GLOBE/branches/pluggable/meshfem3D.f90
Log:
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.
Modified: seismo/3D/SPECFEM3D_GLOBE/branches/pluggable/meshfem3D.f90
===================================================================
--- seismo/3D/SPECFEM3D_GLOBE/branches/pluggable/meshfem3D.f90 2009-02-25 18:23:50 UTC (rev 14147)
+++ seismo/3D/SPECFEM3D_GLOBE/branches/pluggable/meshfem3D.f90 2009-02-25 18:26:01 UTC (rev 14148)
@@ -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
More information about the CIG-COMMITS
mailing list