[cig-commits] [commit] devel, master: switch to direct access i/o for writing noise wavefields (7e7b86a)

[cig_noreply at geodynamics.org]

Repository : https://github.com/geodynamics/specfem2d

On branches: devel,master
Link       : https://github.com/geodynamics/specfem2d/compare/fc67e6fd7ad890705b2b72b4b3c509accb22249e...e9ca46c40131588d89d7b0883250bc6584ce6b4c

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

commit 7e7b86a1b4ab3ecdf2242d2677cbe7e090df24c6
Author: Ryan Modrak <rmodrak at princeton.edu>
Date:   Thu Oct 6 06:56:17 2011 +0000

    switch to direct access i/o for writing noise wavefields


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

7e7b86a1b4ab3ecdf2242d2677cbe7e090df24c6
 noise_uniform/adj_cc.f90 | 15 ++++++++++++---
 1 file changed, 12 insertions(+), 3 deletions(-)

diff --git a/noise_uniform/adj_cc.f90 b/noise_uniform/adj_cc.f90
index c00bd9f..a30a293 100644
--- a/noise_uniform/adj_cc.f90
+++ b/noise_uniform/adj_cc.f90
@@ -9,12 +9,12 @@ integer, parameter :: differentiate = 1
 logical, parameter :: use_filtering = .true.
 
 !choose exactly one of the following window options
-logical, parameter :: use_negative_branch = .true.
-logical, parameter :: use_positive_branch = .false.
+logical, parameter :: use_negative_branch = .false.
+logical, parameter :: use_positive_branch = .true.
 logical, parameter :: use_custom_window = .false.
 
 !choose whether to time reverse, carried out subsequent to all other processing
-logical, parameter :: time_reverse = .true.
+logical, parameter :: time_reverse = .false.
 
 
 
@@ -135,17 +135,25 @@ if (use_custom_window) then
   it_end   = ceiling((t_end - t(1))/dt)
   if (it_begin < 1) it_begin = 1
   if (it_end > nt) it_end = nt
+
 elseif (use_positive_branch) then
   write(*,*) 'Choosing positive branch'
   it_begin = nthalf+1
   it_end   = nt
+  if (it_begin < nthalf) it_begin = nthalf
+  if (it_end > nt) it_end = nt
+
 elseif (use_negative_branch) then
   write(*,*) 'Choosing negative branch'
   it_begin = 1
   it_end   = nthalf
+  if (it_begin < 1) it_begin = 1
+  if (it_end > nthalf) it_end = nthalf
+
 else
   write(*,*) 'Must select one of the following: positive_branch, &
               negative_branch, custom_window.'
+
 endif
 
 write(*,'(a,2f10.3)') ' Time range: ', t(1), t(nt)
@@ -153,6 +161,7 @@ write(*,'(a,2f10.3)') ' Window:     ', t(it_begin), t(it_end)
 write(*,'(a,f10.3,f10.3)') ' Filtering:  ', 1./freq_high, 1./freq_low
 
 !! Tukey taper
+alpha = w_tukey
 k=0
 do it = it_begin,it_end
   k=k+1