Object Info



scalar_values --> object
----------------- Properties:
scalar_values has body({[ E := ( 1 / 2
* ( 1 / h_y * ave(h_y * u ^ 2, x)
+ 1 / h_x * ave(h_x * v ^ 2, y)
)
)
, R_d := 287.04 ~ "J/kg/K"
, R_v := 461.51 ~ "J/kg/K"
, c_pd := 1004.64 ~ "J/kg/K"
, c_pv := 1869.46 ~ "J/kg/K"
, delta := c_pv / c_pd
, epsilon := R_d / R_v
, kappa := R_d / c_pd
, p := A + B * p_s
, Z := ( ( ( df(h_y * v, x)
+ f * ave(ave(h_x * h_y, x), y)
)
- df(h_x * u, y)
)
/ ave(ave(h_x * h_y * df(p, z ~ "1"), x), y)
)
, ln_p := ( df(( log(0.0 if k <= 1 \\ p)
* (0.0 if k <= 1 \\ p)
), z ~ "1")
/ df(p, z ~ "1")
)
, phi := ( phi_s
+ sum(( R_d
* T
* (q * (1 / epsilon - 1) + 1)
* ( ((1 - ( ( df(log(p), z ~ "1")
* shift(p, z, - 1 / 2)
)
/ df(p, z ~ "1")
) if
2 <= k and k <= nlev
) \\
0.693147180559945
)
+ ((( df(log(p), z ~ "1")
- ((1 - ( ( df(log(p), z ~ "1")
* shift(p, z, - 1 / 2)
)
/ df(p, z ~ "1")
) if
2 <= k and k <= nlev
) \\
0.693147180559945
)
) if
2 <= k and k <= nlev
) \\
- 0.693147180559945
)
)
), k = k + 1 .. nlev)
+ ( R_d
* T
* (q * (1 / epsilon - 1) + 1)
* ((1 - ( ( df(log(p), z ~ "1")
* shift(p, z, - 1 / 2)
)
/ df(p, z ~ "1")
) if
2 <= k and k <= nlev
) \\
0.693147180559945
)
)
)
, r := 6371.0 ~ "km"
, u_aux := u * df(p, z ~ "1")
, v_aux := v * df(p, z ~ "1")
, p_s_t := - int(( 1 / (h_x * h_y)
* ( df(h_y * u_aux, x)
+ df(h_x * v_aux, y)
)
), z ~ "1" = 1 .. nlev + 1)
, eta_p := ( int(( 1 / (h_x * h_y)
* ( df(h_y * u_aux, x)
+ df(h_x * v_aux, y)
)
), z ~ "1" = point(z) ~ "1" .. nlev + 1)
+ p_s_t * (1 - B)
)
, T_t := ( kappa * ( ( ( T
* (q * (1 / epsilon - 1) + 1)
* ( ( 1 / (h_x * h_y)
* (df(h_y * u_aux, x) + df(h_x * v_aux, y))
* ((( df(log(p), z ~ "1")
- ((1 - ( ( df(log(p), z ~
"1")
* shift(p, z, - 1 / 2)
)
/ df(p, z ~ "1")
) if
2 <= k and k <= nlev
) \\
0.693147180559945
)
) if
2 <= k and k <= nlev
) \\
- 0.693147180559945
)
)
+ ((( df(log(p), z ~ "1")
* (p_s_t + protect(int(( 1 / (h_x * h_y)
* ( df(h_y * u_aux, x)
+ df(h_x * v_aux, y)
)
), z ~
"1" = k + 1 .. nlev + 1), z))
) if
2 <= k and k <= nlev
) \\ 0.0)
)
)
+ ( 1 / (h_x * h_y)
* ( ave(( h_y
* u_aux
* ave(T * (q * (1 / epsilon - 1) + 1), x)
* df(ln_p, x)
), x)
+ ave(( h_x
* v_aux
* ave(T * (q * (1 / epsilon - 1) + 1), y)
* df(ln_p, y)
), y)
)
)
)
/ (df(p, z ~ "1") * (q * (delta - 1) + 1))
)
- ( ave(( df(T, z ~ "1")
* (eta_p if 2 <= k and k <= nlev \\ 0.0)
), z)
/ df(p, z ~ "1")
)
- ( 1 / (h_x * h_y)
* ( ( ave(h_y * u_aux * df(T, x), x)
+ ave(h_x * v_aux * df(T, y), y)
)
/ df(p, z ~ "1")
)
)
)
, q_t := ( - ( ave(( df(q, z ~ "1")
* (eta_p if 2 <= k and k <= nlev \\ 0.0)
), z)
/ df(p, z ~ "1")
)
- ( 1 / (h_x * h_y)
* ( ( ave(h_y * u_aux * df(q, x), x)
+ ave(h_x * v_aux * df(q, y), y)
)
/ df(p, z ~ "1")
)
)
)
, u_t := ( 1 / h_x * ave(ave(h_y * v_aux, x), y) * protect(Z)
- ( (R_d * T) / h_x
* df(ln_p, x)
* (q * (1 / epsilon - 1) + 1)
)
- ( ave(( df(u, z ~ "1")
* (eta_p if 2 <= k and k <= nlev \\ 0.0)
), z)
/ df(p, z ~ "1")
)
- 1 / h_x * df(E + phi, x)
)
, v_t := ( ( 1 / h_y
* ave(ave(- h_x * u_aux, x), y)
* protect(Z)
)
- ( (R_d * T) / h_y
* df(ln_p, y)
* (q * (1 / epsilon - 1) + 1)
)
- ( ave(( df(v, z ~ "1")
* (eta_p if 2 <= k and k <= nlev \\ 0.0)
), z)
/ df(p, z ~ "1")
)
- 1 / h_y * df(E + phi, y)
) ]})
----------------- Methods inherited from class object:
scalar_values can restore(_G35705)
save(_G35705, _G35706)
get_properties(_G35705, _G35706)
get_superobjects(_G35705, _G35706)
get_subobjects(_G35705, _G35706)
get_parents(_G35705, _G35706)
get_children(_G35705, _G35706)
get_instances(_G35705, _G35706)
create(_G35705, _G35706)
get_dependencies(_G35705, _G35706)