[cig-commits] r21580 - seismo/3D/SPECFEM3D_GLOBE/trunk/src/specfem3D

[dkomati1 at geodynamics.org]

Author: dkomati1
Date: 2013-03-19 15:04:49 -0700 (Tue, 19 Mar 2013)
New Revision: 21580

Modified:
   seismo/3D/SPECFEM3D_GLOBE/trunk/src/specfem3D/compute_forces_outer_core_Dev.f90
   seismo/3D/SPECFEM3D_GLOBE/trunk/src/specfem3D/locate_regular_points.f90
   seismo/3D/SPECFEM3D_GLOBE/trunk/src/specfem3D/write_seismograms.f90
Log:
added white lines for clarity;
also removed white spaces


Modified: seismo/3D/SPECFEM3D_GLOBE/trunk/src/specfem3D/compute_forces_outer_core_Dev.f90
===================================================================
--- seismo/3D/SPECFEM3D_GLOBE/trunk/src/specfem3D/compute_forces_outer_core_Dev.f90	2013-03-19 21:48:04 UTC (rev 21579)
+++ seismo/3D/SPECFEM3D_GLOBE/trunk/src/specfem3D/compute_forces_outer_core_Dev.f90	2013-03-19 22:04:49 UTC (rev 21580)
@@ -251,6 +251,7 @@
                                 hprime_xx(i,5)*B1_m1_m2_5points(5,j)
       enddo
     enddo
+
     do k = 1,NGLLX
       do j=1,m1
         do i=1,m1
@@ -262,6 +263,7 @@
         enddo
       enddo
     enddo
+
     do j=1,m1
       do i=1,m2
         C1_mxm_m2_m1_5points(i,j) = A1_mxm_m2_m1_5points(i,1)*hprime_xxT(1,j) + &
@@ -429,6 +431,7 @@
                                 hprimewgll_xxT(i,5)*C1_m1_m2_5points(5,j)
       enddo
     enddo
+
     do k = 1,NGLLX
       do j=1,m1
         do i=1,m1
@@ -440,6 +443,7 @@
         enddo
       enddo
     enddo
+
     do j=1,m1
       do i=1,m2
         E1_mxm_m2_m1_5points(i,j) = C1_mxm_m2_m1_5points(i,1)*hprimewgll_xx(1,j) + &

Modified: seismo/3D/SPECFEM3D_GLOBE/trunk/src/specfem3D/locate_regular_points.f90
===================================================================
--- seismo/3D/SPECFEM3D_GLOBE/trunk/src/specfem3D/locate_regular_points.f90	2013-03-19 21:48:04 UTC (rev 21579)
+++ seismo/3D/SPECFEM3D_GLOBE/trunk/src/specfem3D/locate_regular_points.f90	2013-03-19 22:04:49 UTC (rev 21580)
@@ -55,7 +55,7 @@
   nlayer = 1
   do while (nlayer <= NM_KL_REG_LAYER)
     read(10,*,iostat=ios) GRID%rlayer(nlayer), GRID%ndoubling(nlayer)
-    if (ios.ne.0) exit
+    if (ios/=0) exit
     nlayer = nlayer + 1
   enddo
   close(10)

Modified: seismo/3D/SPECFEM3D_GLOBE/trunk/src/specfem3D/write_seismograms.f90
===================================================================
--- seismo/3D/SPECFEM3D_GLOBE/trunk/src/specfem3D/write_seismograms.f90	2013-03-19 21:48:04 UTC (rev 21579)
+++ seismo/3D/SPECFEM3D_GLOBE/trunk/src/specfem3D/write_seismograms.f90	2013-03-19 22:04:49 UTC (rev 21580)
@@ -585,11 +585,11 @@
   double precision DT
   character(len=2) bic
 
-  if (DT .ge. 1.0d0)  bic = 'LX'
-  if (DT .lt. 1.0d0 .and. DT .gt. 0.1d0) bic = 'MX'
-  if (DT .le. 0.1d0 .and. DT .gt. 0.0125d0) bic = 'BX'
-  if (DT .le. 0.0125d0 .and. DT .gt. 0.004d0) bic = 'HX'
-  if (DT .le. 0.004d0 .and. DT .gt. 0.001d0) bic = 'CX'
-  if (DT .le. 0.001d0) bic = 'FX'
+  if (DT >= 1.0d0)  bic = 'LX'
+  if (DT < 1.0d0 .and. DT > 0.1d0) bic = 'MX'
+  if (DT <= 0.1d0 .and. DT > 0.0125d0) bic = 'BX'
+  if (DT <= 0.0125d0 .and. DT > 0.004d0) bic = 'HX'
+  if (DT <= 0.004d0 .and. DT > 0.001d0) bic = 'CX'
+  if (DT <= 0.001d0) bic = 'FX'
 
  end subroutine band_instrument_code