[Dart-dev] [3640] DART/trunk/assim_tools/assim_tools_mod.f90: Corrected trapezoidal option in boxcar2 that was not in use.

[nancy at ucar.edu]

An HTML attachment was scrubbed...
URL: http://mailman.ucar.edu/pipermail/dart-dev/attachments/20081103/347c1c64/attachment.html
-------------- next part --------------
Modified: DART/trunk/assim_tools/assim_tools_mod.f90
===================================================================
--- DART/trunk/assim_tools/assim_tools_mod.f90	2008-10-30 21:47:34 UTC (rev 3639)
+++ DART/trunk/assim_tools/assim_tools_mod.f90	2008-11-03 16:55:23 UTC (rev 3640)
@@ -1706,7 +1706,7 @@
    ! Height of prior in this bin is mass/width; Only needed for trapezoidal
    ! If two ensemble members are the same, set height to -1 as flag
    if(x(i) == x(i - 1)) then 
-      height = -1.0_r8
+      height(i) = -1.0_r8
    else
       height(i) = 1.0_r8 / ((ens_size + 1.0_r8) * (x(i) - x(i-1)))
    endif
@@ -1813,14 +1813,9 @@
    endif
 end do
 
-! Can now compute sorted increments
-do i = 1, ens_size
-   sort_inc(i) = new_ens(i) - x(i)
-end do
-
 ! Now, need to convert to increments for unsorted
 do i = 1, ens_size
-   obs_inc(e_ind(i)) = sort_inc(i)
+   obs_inc(e_ind(i)) = new_ens(i) - x(i)
 end do
 
 end subroutine obs_increment_boxcar2