```ln_p :: (real(- infinity .. infinity) ~ "1" field ( x(grid)                                                  , y(grid)                                                  , z(half)                                                  ) on         (i = 0 .. nlon + 1 by          j = 0 .. nlat + 1 by          k = 1 .. nlev         )        )

Scalar Form

ln_p := ( df(( log(0.0 if k <= 1 \\ p)             * (0.0 if k <= 1 \\ p)             ), z ~ "1")        / df(p, z ~ "1")        )
Dimensioned

ln_p := ( df(( log(0.0 ~ "Pa" if k <= 1 ~ "1" \\ p)             * (0.0 ~ "Pa" if k <= 1 ~ "1" \\ p)             ), z ~ "1")        / df(p, z ~ "1")        )
Discretized

ln_p := ( df_g(log(0.0 if k <= 1 \\ p) * (0.0 if k <= 1 \\ p), z)        / df_g(p, z)        )
CSE Eliminated

ln_p := ((( varidx(t_13, [i, j, 2])          / (varidx(p, [i, j, k + 1]) - varidx(p, [i, j, k]))          ) if          k <= 1         ) \\         ( (varidx(t_13, [i, j, k + 1]) - varidx(t_13, [i, j, k]))         / (varidx(p, [i, j, k + 1]) - varidx(p, [i, j, k]))         )        )```