Posted to tcl by mig at Fri Jan 09 02:55:01 GMT 2015view raw

  1. for {set i5 0} {$i5 < $ncells} {incr i5} {
  2. set calck [expr {$i5/$grdsize}]
  3. set grdindex [expr {$i5-$calck*$grdsize}]
  4. set volpt [[contGeom GetOutput] GetPoint $i5];# HERE!
  5. set volz [lindex $volpt 2]
  6. if {$calck < $volzdim-1} {
  7. set i5b [expr {$i5 + $grdsize}]
  8. } else {
  9. set i5b [expr {$i5 - $grdsize}]
  10. }
  11. set volptb [[contGeom GetOutput] GetPoint $i5b];# HERE!
  12. set volzb [lindex $volptb 2]
  13. set zgriddim [expr {abs($volz - $volzb)}]
  14.  
  15. set cellvol [expr {$dx*$dy*$zgriddim}]
  16.  
  17. set ptscal [[[[contGeom GetOutput] GetPointData] GetScalars] GetValue $i5];# HERE!
  18.  
  19. if {$volz>=$elev1 && $volz<=$elev2} {
  20. if {$threshtype==1 && $ptscal>= $volcontvalU} {
  21. sumcells
  22. } elseif {$threshtype==0 && $ptscal>=$volcontvalL && $ptscal<=$volcontvalU} {
  23. sumcells
  24. } elseif {$threshtype==-1 && $ptscal<=$volcontvalL} {
  25. sumcells
  26. }
  27. } else {
  28. eval [[[[contGeom GetOutput] GetPointData] GetScalars] SetValue $i5 $repor];# HERE!!!
  29. }
  30. }
  31.