Object Info



indexed_values --> object
------------------ Properties:
indexed_values has body({[ (varidx(E, [i, j, k]) :=
( 1 / 2
* ( ( 1 / 2
* ( ( ( varidx(v, [i, j, k]) ^ 2
* varidx(hxv, [i, j])
)
+ ( varidx(v, [i, j - 1, k]) ^ 2
* varidx(hxv, [i, j - 1])
)
)
/ varidx(hxt, [i, j])
)
)
+ ( 1 / 2
* ( ( ( varidx(u, [i, j, k]) ^ 2
* varidx(hyu, [i, j])
)
+ ( varidx(u, [i - 1, j, k]) ^ 2
* varidx(hyu, [i - 1, j])
)
)
/ varidx(hyt, [i, j])
)
)
)
)
)
, R_d := 287.04
, R_v := 461.51
, c_pd := 1004.64
, c_pv := 1869.46
, delta := 1.86082576843446
, epsilon := 0.621958354098503
, kappa := 0.285714285714286
, (varidx(p, [i, j, k]) :=
varidx(A, [k]) + varidx(B, [k]) * varidx(p_s, [i, j])
)
, (varidx(ln_p, [i, j, k]) :=
( ( ( log(varidx(p, [i, j, k + 1]))
* varidx(p, [i, j, k + 1])
)
- ( log(0.0 if k <= 1 \\ varidx(p, [i, j, k]))
* (0.0 if k <= 1 \\ varidx(p, [i, j, k]))
)
)
/ (varidx(p, [i, j, k + 1]) - varidx(p, [i, j, k]))
)
)
, (varidx(phi, [i, j, k]) :=
( sum(( 287.04
* varidx(T, [i, j, k])
* (( 0.607824693422519
* varidx(q, [i, j, k])
) + 1)
* ( (1 - ( ( varidx(p, [i, j, k])
* ( log(varidx(p, [i, j, k + 1]))
- log(varidx(p, [i, j, k]))
)
)
/ ( varidx(p, [i, j, k + 1])
- varidx(p, [i, j, k])
)
) \\
0.693147180559945
)
+ (( log(varidx(p, [i, j, k + 1]))
- log(varidx(p, [i, j, k]))
- (1 - ( ( varidx(p, [i, j, k])
* ( log(varidx(p, [i, j, k + 1]))
- log(varidx(p, [i, j, k]))
)
)
/ ( varidx(p, [i, j, k + 1])
- varidx(p, [i, j, k])
)
) \\
0.693147180559945
)
) \\
- 0.693147180559945
)
)
), k = k + 1 .. nlev)
+ varidx(phi_s, [i, j])
+ ( 287.04
* varidx(T, [i, j, k])
* (0.607824693422519 * varidx(q, [i, j, k]) + 1)
* ((1 - ( ( varidx(p, [i, j, k])
* ( log(varidx(p, [i, j, k + 1]))
- log(varidx(p, [i, j, k]))
)
)
/ ( varidx(p, [i, j, k + 1])
- varidx(p, [i, j, k])
)
) if
2 <= k
) \\
0.693147180559945
)
)
)
)
, r := 6371000.0
, (varidx(Z, [i, j, k]) :=
2 * ( ( ( ( 1.56961230576048E-7
* rdth
* ( varidx(hxu, [i, j]) * varidx(u, [i, j, k])
- ( varidx(hxu, [i, j + 1])
* varidx(u, [i, j + 1, k])
)
)
)
+ ( 1 / 2
* varidx(f, [i, j])
* ( ( 1 / 2
* ( varidx(hxt, [i, j]) * varidx(hyt, [i, j])
+ ( varidx(hxt, [i + 1, j])
* varidx(hyt, [i + 1, j])
)
)
)
+ ( 1 / 2
* ( ( varidx(hxt, [i, j + 1])
* varidx(hyt, [i, j + 1])
)
+ ( varidx(hxt, [i + 1, j + 1])
* varidx(hyt, [i + 1, j + 1])
)
)
)
)
)
)
- ( 1.56961230576048E-7
* rdlam
* ( varidx(hyv, [i, j]) * varidx(v, [i, j, k])
- ( varidx(hyv, [i + 1, j])
* varidx(v, [i + 1, j, k])
)
)
)
)
/ ( ( 1 / 2
* ( ( varidx(hxt, [i, j])
* varidx(hyt, [i, j])
* ( varidx(p, [i, j, k + 1])
- varidx(p, [i, j, k])
)
)
+ ( varidx(hxt, [i + 1, j])
* varidx(hyt, [i + 1, j])
* ( varidx(p, [i + 1, j, k + 1])
- varidx(p, [i + 1, j, k])
)
)
)
)
+ ( 1 / 2
* ( ( varidx(hxt, [i, j + 1])
* varidx(hyt, [i, j + 1])
* ( varidx(p, [i, j + 1, k + 1])
- varidx(p, [i, j + 1, k])
)
)
+ ( varidx(hxt, [i + 1, j + 1])
* varidx(hyt, [i + 1, j + 1])
* ( varidx(p, [i + 1, j + 1, k + 1])
- varidx(p, [i + 1, j + 1, k])
)
)
)
)
)
)
)
, (varidx(u_aux, [i, j, k]) :=
( varidx(u, [i, j, k])
* ( ( 1 / 2
* ( varidx(p, [i, j, k + 1])
+ varidx(p, [i + 1, j, k + 1])
)
)
- ( 1 / 2
* (varidx(p, [i, j, k]) + varidx(p, [i + 1, j, k]))
)
)
)
)
, (varidx(v_aux, [i, j, k]) :=
( varidx(v, [i, j, k])
* ( ( 1 / 2
* ( varidx(p, [i, j, k + 1])
+ varidx(p, [i, j + 1, k + 1])
)
)
- ( 1 / 2
* (varidx(p, [i, j, k]) + varidx(p, [i, j + 1, k]))
)
)
)
)
, (varidx(p_s_t, [i, j]) :=
( ( - ( 1.56961230576048E-7
* rdlam
* ( ( sum(varidx(u_aux, [ i
, j
, k ]), k = 1 .. nlev)
* varidx(hyu, [i, j])
)
- ( sum(varidx(u_aux, [i - 1, j, k]), k = 1 .. nlev)
* varidx(hyu, [i - 1, j])
)
)
)
- ( 1.56961230576048E-7
* rdth
* ( ( sum(varidx(v_aux, [ i
, j
, k ]), k = 1 .. nlev)
* varidx(hxv, [i, j])
)
- ( sum(varidx(v_aux, [i, j - 1, k]), k = 1 .. nlev)
* varidx(hxv, [i, j - 1])
)
)
)
)
/ (varidx(hxt, [i, j]) * varidx(hyt, [i, j]))
)
)
, (varidx(eta_p, [i, j, k]) :=
( ( ( ( 1.56961230576048E-7
* rdlam
* ( ( sum(varidx(u_aux, [ i
, j
, k ]), k = k .. nlev)
* varidx(hyu, [i, j])
)
- ( sum(varidx(u_aux, [i - 1, j, k]), k = k .. nlev)
* varidx(hyu, [i - 1, j])
)
)
)
+ ( 1.56961230576048E-7
* rdth
* ( ( sum(varidx(v_aux, [ i
, j
, k ]), k = k .. nlev)
* varidx(hxv, [i, j])
)
- ( sum(varidx(v_aux, [i, j - 1, k]), k = k .. nlev)
* varidx(hxv, [i, j - 1])
)
)
)
)
/ (varidx(hxt, [i, j]) * varidx(hyt, [i, j]))
)
+ varidx(p_s_t, [i, j]) * (1 - varidx(B, [k]))
)
)
, (varidx(T_t, [i, j, k]) :=
( ( ( ( 0.285714285714286
* varidx(T, [i, j, k])
* ( 1 / (( 0.860825768434464
* varidx(q, [i, j, k])
) + 1)
+ ( 0.607824693422519
* ( varidx(q, [i, j, k])
/ (( 0.860825768434464
* varidx(q, [i, j, k])
) + 1)
)
)
)
* ( ( ( ( ( 1.56961230576048E-7
* rdlam
* ( ( varidx(hyu, [i, j])
* varidx(u_aux, [i, j, k])
)
- ( varidx(hyu, [i - 1, j])
* varidx(u_aux, [i - 1, j, k])
)
)
)
+ ( 1.56961230576048E-7
* rdth
* ( ( varidx(hxv, [i, j])
* varidx(v_aux, [i, j, k])
)
- ( varidx(hxv, [i, j - 1])
* varidx(v_aux, [i, j - 1, k])
)
)
)
)
* ((( log(varidx(p, [i, j, k + 1]))
- log(varidx(p, [i, j, k]))
- (1 - ( ( varidx(p, [i, j, k])
* ( log(varidx(p, [i, j, k + 1]))
- log(varidx(p, [ i
, j
, k ]))
)
)
/ ( varidx(p, [i, j, k + 1])
- varidx(p, [i, j, k])
)
) \\
0.693147180559945
)
) if
2 <= k
) \\
- 0.693147180559945
)
)
/ (varidx(hxt, [i, j]) * varidx(hyt, [i, j]))
)
+ (( ( log(varidx(p, [i, j, k + 1]))
- log(varidx(p, [i, j, k]))
)
* ( varidx(p_s_t, [i, j])
+ ( ( ( 1.56961230576048E-7
* rdlam
* ( sum(varidx(u_aux, [ i
, j
, k ]), k = k + 1 .. nlev) * varidx(hyu, [i, j])
- ( sum(varidx(u_aux, [i - 1, j, k]), k = k + 1 .. nlev)
* varidx(hyu, [i - 1, j])
)
)
)
+ ( 1.56961230576048E-7
* rdth
* ( sum(varidx(v_aux, [ i
, j
, k ]), k = k + 1 .. nlev) * varidx(hxv, [i, j])
- ( sum(varidx(v_aux, [i, j - 1, k]), k = k + 1 .. nlev)
* varidx(hxv, [i, j - 1])
)
)
)
)
/ ( varidx(hxt, [i, j])
* varidx(hyt, [i, j])
)
)
)
) if 2 <= k \\ 0.0)
)
)
+ ( ( ( ( 1 / 2
* ( ( 1.56961230576048E-7
* rdlam
* varidx(hyu, [i, j])
* varidx(u_aux, [i, j, k])
* ( varidx(T, [i, j, k])
- varidx(T, [i + 1, j, k])
)
)
- ( 1.56961230576048E-7
* rdlam
* varidx(hyu, [i - 1, j])
* varidx(u_aux, [i - 1, j, k])
* ( varidx(T, [i, j, k])
- varidx(T, [i - 1, j, k])
)
)
)
)
+ ( 1 / 2
* ( ( 1.56961230576048E-7
* rdth
* varidx(hxv, [i, j])
* varidx(v_aux, [i, j, k])
* ( varidx(T, [i, j, k])
- varidx(T, [i, j + 1, k])
)
)
- ( 1.56961230576048E-7
* rdth
* varidx(hxv, [i, j - 1])
* varidx(v_aux, [i, j - 1, k])
* ( varidx(T, [i, j, k])
- varidx(T, [i, j - 1, k])
)
)
)
)
)
- ( 0.285714285714286
* ( ( ( 1.60782469342252
* ( ( 1 / 2
* ( ( 7.84806152880239E-8
* rdlam
* varidx(hyu, [i, j])
* varidx(u_aux, [i, j, k])
* ( varidx(ln_p, [i, j, k])
- varidx(ln_p, [i + 1, j, k])
)
* ( ( varidx(T, [i, j, k])
* varidx(q, [i, j, k])
)
+ ( varidx(T, [i + 1, j, k])
* varidx(q, [i + 1, j, k])
)
)
)
- ( 7.84806152880239E-8
* rdlam
* varidx(hyu, [i - 1, j])
* varidx(u_aux, [i - 1, j, k])
* ( varidx(ln_p, [i, j, k])
- varidx(ln_p, [i - 1, j, k])
)
* ( ( varidx(T, [i, j, k])
* varidx(q, [i, j, k])
)
+ ( varidx(T, [i - 1, j, k])
* varidx(q, [i - 1, j, k])
)
)
)
)
)
+ ( 1 / 2
* ( ( 7.84806152880239E-8
* rdth
* varidx(hxv, [i, j])
* varidx(v_aux, [i, j, k])
* ( varidx(ln_p, [i, j, k])
- varidx(ln_p, [i, j + 1, k])
)
* ( ( varidx(T, [i, j, k])
* varidx(q, [i, j, k])
)
+ ( varidx(T, [i, j + 1, k])
* varidx(q, [i, j + 1, k])
)
)
)
- ( 7.84806152880239E-8
* rdth
* varidx(hxv, [i, j - 1])
* varidx(v_aux, [i, j - 1, k])
* ( varidx(ln_p, [i, j, k])
- varidx(ln_p, [i, j - 1, k])
)
* ( ( varidx(T, [i, j, k])
* varidx(q, [i, j, k])
)
+ ( varidx(T, [i, j - 1, k])
* varidx(q, [i, j - 1, k])
)
)
)
)
)
)
)
+ ( 1 / 2
* ( ( 7.84806152880239E-8
* rdlam
* varidx(hyu, [i, j])
* varidx(u_aux, [i, j, k])
* ( varidx(ln_p, [i, j, k])
- varidx(ln_p, [i + 1, j, k])
)
* ( ( varidx(T, [i, j, k])
* (1 - varidx(q, [i, j, k]))
)
+ ( varidx(T, [i + 1, j, k])
* (1 - varidx(q, [i + 1, j, k]))
)
)
)
- ( 7.84806152880239E-8
* rdlam
* varidx(hyu, [i - 1, j])
* varidx(u_aux, [i - 1, j, k])
* ( varidx(ln_p, [i, j, k])
- varidx(ln_p, [i - 1, j, k])
)
* ( ( varidx(T, [i, j, k])
* (1 - varidx(q, [i, j, k]))
)
+ ( varidx(T, [i - 1, j, k])
* (1 - varidx(q, [i - 1, j, k]))
)
)
)
)
)
+ ( 1 / 2
* ( ( 7.84806152880239E-8
* rdth
* varidx(hxv, [i, j])
* varidx(v_aux, [i, j, k])
* ( varidx(ln_p, [i, j, k])
- varidx(ln_p, [i, j + 1, k])
)
* ( ( varidx(T, [i, j, k])
* (1 - varidx(q, [i, j, k]))
)
+ ( varidx(T, [i, j + 1, k])
* (1 - varidx(q, [i, j + 1, k]))
)
)
)
- ( 7.84806152880239E-8
* rdth
* varidx(hxv, [i, j - 1])
* varidx(v_aux, [i, j - 1, k])
* ( varidx(ln_p, [i, j, k])
- varidx(ln_p, [i, j - 1, k])
)
* ( ( varidx(T, [i, j, k])
* (1 - varidx(q, [i, j, k]))
)
+ ( varidx(T, [i, j - 1, k])
* (1 - varidx(q, [i, j - 1, k]))
)
)
)
)
)
)
/ (( 0.860825768434464
* varidx(q, [i, j, k])
) + 1)
)
)
)
/ (varidx(hxt, [i, j]) * varidx(hyt, [i, j]))
)
)
- ( 1 / 2
* ( ( ( varidx(T, [i, j, k + 1])
- varidx(T, [i, j, k])
)
* ((varidx(eta_p, [i, j, k + 1]) if
k + 1 <= nlev
) \\ 0.0)
)
+ ( ( varidx(T, [i, j, k])
- varidx(T, [i, j, k - 1])
)
* (varidx(eta_p, [ i
, j
, k ]) if 2 <= k \\ 0.0)
)
)
)
)
/ (varidx(p, [i, j, k + 1]) - varidx(p, [i, j, k]))
)
)
, (varidx(q_t, [i, j, k]) :=
( ( ( ( ( 1 / 2
* ( ( 1.56961230576048E-7
* rdlam
* varidx(hyu, [i, j])
* varidx(u_aux, [i, j, k])
* ( varidx(q, [i, j, k])
- varidx(q, [i + 1, j, k])
)
)
- ( 1.56961230576048E-7
* rdlam
* varidx(hyu, [i - 1, j])
* varidx(u_aux, [i - 1, j, k])
* ( varidx(q, [i, j, k])
- varidx(q, [i - 1, j, k])
)
)
)
)
+ ( 1 / 2
* ( ( 1.56961230576048E-7
* rdth
* varidx(hxv, [i, j])
* varidx(v_aux, [i, j, k])
* ( varidx(q, [i, j, k])
- varidx(q, [i, j + 1, k])
)
)
- ( 1.56961230576048E-7
* rdth
* varidx(hxv, [i, j - 1])
* varidx(v_aux, [i, j - 1, k])
* ( varidx(q, [i, j, k])
- varidx(q, [i, j - 1, k])
)
)
)
)
)
/ (varidx(hxt, [i, j]) * varidx(hyt, [i, j]))
)
- ( 1 / 2
* ( ( ( varidx(q, [i, j, k + 1])
- varidx(q, [i, j, k])
)
* ((varidx(eta_p, [i, j, k + 1]) if
k + 1 <= nlev
) \\ 0.0)
)
+ ( ( varidx(q, [i, j, k])
- varidx(q, [i, j, k - 1])
)
* (varidx(eta_p, [ i
, j
, k ]) if 2 <= k \\ 0.0)
)
)
)
)
/ (varidx(p, [i, j, k + 1]) - varidx(p, [i, j, k]))
)
)
, (varidx(u_t, [i, j, k]) :=
( 1 / 2 * ( ( ( ( varidx(u, [i, j, k + 1])
- varidx(u, [i, j, k])
)
* ((( 1 / 2
* ( varidx(eta_p, [i, j, k + 1])
+ varidx(eta_p, [ i + 1
, j
, k + 1 ])
)
) if
k + 1 <= nlev
) \\ 0.0)
)
+ ( ( varidx(u, [i, j, k])
- varidx(u, [i, j, k - 1])
)
* ((( 1 / 2
* ( varidx(eta_p, [i, j, k])
+ varidx(eta_p, [i + 1, j, k])
)
) if
2 <= k
) \\ 0.0)
)
)
/ ( ( 1 / 2
* ( varidx(p, [i, j, k])
+ varidx(p, [i + 1, j, k])
)
)
- ( 1 / 2
* ( varidx(p, [i, j, k + 1])
+ varidx(p, [i + 1, j, k + 1])
)
)
)
)
+ ( ( ( 1.56961230576048E-7
* rdlam
* ( (varidx(E, [i, j, k]) + varidx(phi, [i, j, k]))
- varidx(E, [i + 1, j, k])
- varidx(phi, [i + 1, j, k])
)
)
+ ( 2.25270758122744E-5
* rdlam
* (( 0.30391234671126
* ( varidx(q, [i, j, k])
+ varidx(q, [i + 1, j, k])
)
) + 1)
* ( varidx(ln_p, [i, j, k])
- varidx(ln_p, [i + 1, j, k])
)
* ( varidx(T, [i, j, k])
+ varidx(T, [i + 1, j, k])
)
)
+ ( 1 / 4
* ( varidx(Z, [i, j, k])
+ varidx(Z, [i, j - 1, k])
)
* ( ( 1 / 2
* ( ( varidx(hyv, [i, j])
* varidx(v_aux, [i, j, k])
)
+ ( varidx(hyv, [i + 1, j])
* varidx(v_aux, [i + 1, j, k])
)
)
)
+ ( 1 / 2
* ( ( varidx(hyv, [i, j - 1])
* varidx(v_aux, [i, j - 1, k])
)
+ ( varidx(hyv, [i + 1, j - 1])
* varidx(v_aux, [i + 1, j - 1, k])
)
)
)
)
)
)
/ varidx(hxu, [i, j])
)
)
)
, (varidx(v_t, [i, j, k]) :=
( 1 / 2 * ( ( ( ( varidx(v, [i, j, k + 1])
- varidx(v, [i, j, k])
)
* ((( 1 / 2
* ( varidx(eta_p, [i, j, k + 1])
+ varidx(eta_p, [ i
, j + 1
, k + 1 ])
)
) if
k + 1 <= nlev
) \\ 0.0)
)
+ ( ( varidx(v, [i, j, k])
- varidx(v, [i, j, k - 1])
)
* ((( 1 / 2
* ( varidx(eta_p, [i, j, k])
+ varidx(eta_p, [i, j + 1, k])
)
) if
2 <= k
) \\ 0.0)
)
)
/ ( ( 1 / 2
* ( varidx(p, [i, j, k])
+ varidx(p, [i, j + 1, k])
)
)
- ( 1 / 2
* ( varidx(p, [i, j, k + 1])
+ varidx(p, [i, j + 1, k + 1])
)
)
)
)
+ ( ( ( ( 1.56961230576048E-7
* rdth
* ( ( varidx(E, [i, j, k])
+ varidx(phi, [i, j, k])
)
- varidx(E, [i, j + 1, k])
- varidx(phi, [i, j + 1, k])
)
)
+ ( 2.25270758122744E-5
* rdth
* (( 0.30391234671126
* ( varidx(q, [i, j, k])
+ varidx(q, [i, j + 1, k])
)
) + 1)
* ( varidx(ln_p, [i, j, k])
- varidx(ln_p, [i, j + 1, k])
)
* ( varidx(T, [i, j, k])
+ varidx(T, [i, j + 1, k])
)
)
)
- ( 1 / 4
* ( varidx(Z, [i, j, k])
+ varidx(Z, [i - 1, j, k])
)
* ( ( 1 / 2
* ( ( varidx(hxu, [i, j])
* varidx(u_aux, [i, j, k])
)
+ ( varidx(hxu, [i - 1, j])
* varidx(u_aux, [i - 1, j, k])
)
)
)
+ ( 1 / 2
* ( ( varidx(hxu, [i, j + 1])
* varidx(u_aux, [i, j + 1, k])
)
+ ( varidx(hxu, [i - 1, j + 1])
* varidx(u_aux, [i - 1, j + 1, k])
)
)
)
)
)
)
/ varidx(hyv, [i, j])
)
)
) ]})
------------------ Methods inherited from class object:
indexed_values can restore(_G26729)
save(_G26729, _G26730)
get_properties(_G26729, _G26730)
get_superobjects(_G26729, _G26730)
get_subobjects(_G26729, _G26730)
get_parents(_G26729, _G26730)
get_children(_G26729, _G26730)
get_instances(_G26729, _G26730)
create(_G26729, _G26730)
get_dependencies(_G26729, _G26730)