grid_view_part: grid_view_part

Description Usage Arguments Examples

Description

grid_view_part

Usage

1
grid_view_part(grid, part = 0.01, alpha = 1, color_p = PT.Color[["O"]], color_n = PT.Color[["N"]], ...)

Arguments

grid
part
alpha
color_p
color_n
...

Examples

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
##---- Should be DIRECTLY executable !! ----
##-- ==>  Define data, use random,
##--	or do  help(data=index)  for the standard data sets.

## The function is currently defined as
function (grid, part = 0.01, alpha = 1, color_p = PT.Color[["O"]], 
    color_n = PT.Color[["N"]], ...) 
{
    npoints <- grid$ngridx * grid$ngridy * grid$ngridz
    
    ivalue <- 1
    for (igridx in 1:grid$ngridx) {
        for (igridy in 1:grid$ngridy) {
            for (igridz in 1:grid$ngridz) {
                values[ivalue] <- abs(grid$val[igridx, igridy, 
                  igridz])
                ivalue <- ivalue + 1
            }
        }
    }
    svalues <- sort(values, decreasing = TRUE)
    level <- svalues[ceiling(npoints * part)]
    cat(sprintf("level=%g\n", level))
    grid_view_level(grid, level, alpha, color_p, color_n, ...)
  }

conmolfields documentation built on May 2, 2019, 4:18 p.m.