Object Info


u_t is defined in /home/robert/projects/Ctadel/dyn.s:315
u_t --> declared_var(script("dyn")) --> declared_object(script("dyn")) --> object
derived_var --> dependent_var --> variable --> object
object
------- Properties:
u_t has cse_value(((( ( 1 / 2
* varidx(t_192, [i, j, nlev]) / varidx(t_47, [i, j, nlev])
)
+ ( ( rdlam * ( ( 1.56961230576048E-7
* ( (varidx(E, [i, j, k]) + varidx(phi, [i, j, nlev]))
- varidx(E, [i + 1, j, k])
- varidx(phi, [i + 1, j, nlev])
)
)
+ ( 2.25270758122744E-5
* (( 0.30391234671126
* ( varidx(q, [i, j, nlev])
+ varidx(q, [i + 1, j, nlev])
)
) + 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 / 16
* (varidx(Z, [i, j, k]) + varidx(Z, [i, j - 1, k]))
* ( ( varidx(hyv, [i, j])
* varidx(t_49, [i, j, nlev])
* varidx(v, [i, j, nlev])
)
+ ( varidx(hyv, [i, j - 1])
* varidx(t_49, [i, j - 1, nlev])
* varidx(v, [i, j - 1, nlev])
)
+ ( varidx(hyv, [i + 1, j])
* varidx(t_49, [i + 1, j, nlev])
* varidx(v, [i + 1, j, nlev])
)
+ ( varidx(hyv, [i + 1, j - 1])
* varidx(t_49, [i + 1, j - 1, nlev])
* varidx(v, [i + 1, j - 1, nlev])
)
)
)
)
/ varidx(hxu, [i, j])
)
) if
nlev <= k
) \\
( 1 / 2 * varidx(t_192, [i, j, 2]) / varidx(t_47, [i, j, 1])
+ ( ( rdlam * ( ( 1.56961230576048E-7
* ( (varidx(E, [i, j, k]) + varidx(phi, [i, j, 1]))
- varidx(E, [i + 1, j, k])
- varidx(phi, [i + 1, j, 1])
)
)
+ ( 2.25270758122744E-5
* (( 0.30391234671126
* ( varidx(q, [i, j, 1])
+ varidx(q, [i + 1, j, 1])
)
) + 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 / 16
* (varidx(Z, [i, j, k]) + varidx(Z, [i, j - 1, k]))
* ( ( varidx(hyv, [i, j])
* varidx(t_49, [i, j, 1])
* varidx(v, [i, j, 1])
)
+ ( varidx(hyv, [i, j - 1])
* varidx(t_49, [i, j - 1, 1])
* varidx(v, [i, j - 1, 1])
)
+ ( varidx(hyv, [i + 1, j])
* varidx(t_49, [i + 1, j, 1])
* varidx(v, [i + 1, j, 1])
)
+ ( varidx(hyv, [i + 1, j - 1])
* varidx(t_49, [i + 1, j - 1, 1])
* varidx(v, [i + 1, j - 1, 1])
)
)
)
)
/ varidx(hxu, [i, j])
)
) if k <= 1 \\
( ( 1 / 2
* ( (varidx(t_192, [i, j, k]) + varidx(t_192, [i, j, k + 1]))
/ varidx(t_47, [i, j, k])
)
)
+ ( ( rdlam * ( ( 1.56961230576048E-7
* ( (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
* (( 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 / 16
* (varidx(Z, [i, j, k]) + varidx(Z, [i, j - 1, k]))
* ( ( varidx(hyv, [i, j])
* varidx(t_49, [i, j, k])
* varidx(v, [i, j, k])
)
+ ( varidx(hyv, [i, j - 1])
* varidx(t_49, [i, j - 1, k])
* varidx(v, [i, j - 1, k])
)
+ ( varidx(hyv, [i + 1, j])
* varidx(t_49, [i + 1, j, k])
* varidx(v, [i + 1, j, k])
)
+ ( varidx(hyv, [i + 1, j - 1])
* varidx(t_49, [i + 1, j - 1, k])
* varidx(v, [i + 1, j - 1, k])
)
)
)
)
/ varidx(hxu, [i, j])
)
)
))
point([i, j, k])
indexed_value(( 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])
)
))
derived_stencil([0 .. 0, 0 .. 0, 0 .. 0])
derived_domain([1 .. nlon, 1 .. nlat, 1 .. nlev])
discrete_value(( ( 1 / hxu
* ( ( ( R_d
* ave_g(T, x)
* df_g(ln_p, x)
* (ave_g(q, x) * (1 - 1 / epsilon) - 1)
)
+ ( ave_h(ave_g(hyv * v_aux, x), y)
* protect(ave_h(Z, y))
)
)
- df_g(E + phi, x)
)
)
- ( ave_g(( df_h(u, z)
* ((ave_g(eta_p, x) if
2 <= k and k <= nlev
) \\ 0.0)
), z)
/ df_g(ave_g(p, x), z)
)
))
grid(grid([half, grid, half, _G23950]))
dimensioned_value(( 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") + 1 ~ "1")
)
- ( ave(( df(u, z ~ "1")
* (eta_p if 2 ~ "1" <= k and k <= nlev \\
0.0 ~ " kg s^-3 / m"
)
), z)
/ df(p, z ~ "1")
)
- 1 / h_x * df(E + phi, x)
))
unit(unit(0, [0, 0, 1, -2, 0, 0, 0, 0]))
scalar_value(( 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)
))
fortran_name(PDUDT)
value(( 1 / h_x * ave(ave(h_y * v_aux, x), y) * protect(Z)
- ( R_d / h_x
* df(ln_p, x)
* (T * (q * (1 / epsilon - 1) + 1) \\ T)
)
- ( ave(( df(u, z ~ "1")
* (eta_p if on (k = 2 .. nlev) \\ real(0))
), z)
/ df(p, z ~ "1")
)
- 1 / h_x * df(E + phi, x)
))
type(real)
declared_domain([1 .. nlon, 1 .. nlat, 1 .. nlev])
dependency([x, y, z])
------- Properties inherited from class dependent_var:
u_t has instance_name(v)
------- Methods inherited from class object:
u_t can restore(_G24164)
save(_G24164, _G24165)
get_properties(_G24164, _G24165)
get_superobjects(_G24164, _G24165)
get_subobjects(_G24164, _G24165)
get_parents(_G24164, _G24165)
get_children(_G24164, _G24165)
get_instances(_G24164, _G24165)
create(_G24164, _G24165)
get_dependencies(_G24164, _G24165)
------- Methods inherited from class variable:
u_t can get_store_cost(_G24175, 10)
get_load_cost(_G24175, 1)