|
| subroutine | most_cpu::select_bc_operators (bc_type, bc_value, q, ts, ti, kappa, utau, z0h, hi, pr) |
| | Selects different expressions for the similarity functions in MOST based on the type of bottom boundary condition for temperature.
|
| |
| subroutine | most_cpu::compute_ri_b (bc_type, g_dot_n, hi, ti, ts, magu, kappa, q, pr, ri_b) |
| | Computes the Richardson number.
|
| |
| subroutine | most_cpu::set_stability_regime (ri_b, ri_threshold) |
| | Sets the stability regime based on the Richardson number value (quite arbitrary).
|
| |
| subroutine, public | most_cpu::most_compute_cpu (u, v, w, temp, ind_r, ind_s, ind_t, ind_e, n_x, n_y, n_z, h, tau_x, tau_y, tau_z, n_nodes, lx, nelv, kappa, mu_w, rho_w, g_vec, pr, z0, z0h_in, bc_type, bc_value, tstep, ri_b_diagn, l_ob_diagn, utau_diagn, magu_diagn, ti_diagn, ts_diagn, q_diagn, h_x_idx, h_y_idx, h_z_idx) |
| | Main routine to compute the surface stresses based on MOST.
|
| |
| real(kind=rp) function | most_cpu::slaw_m_stable (z, l_ob, z0) |
| | Similarity laws and corrections for the STABLE regime: REFERENCE: Holtslag, A. A. M., & De Bruin, H. A. R. (1988). Applied Modeling of the Nighttime Surface Energy Balance over Land. Journal of Applied Meteorology, 27(6), 689–704. NOTE: This formulation is chosen for its superior behavior in very stable conditions (large z/L), avoiding the numerical decoupling found in older linear functions (e.g., Dyer).
|
| |
| real(kind=rp) function | most_cpu::slaw_h_stable (z, l_ob, z0h) |
| |
| real(kind=rp) function | most_cpu::corr_m_stable (z, l_ob) |
| |
| real(kind=rp) function | most_cpu::corr_h_stable (z, l_ob) |
| |
| real(kind=rp) function | most_cpu::slaw_m_convective (z, l_ob, z0) |
| | Similarity laws and corrections for the UNSTABLE (convective) regime: REFERENCE: Dyer, A. J. (1974), A review of flux-profile relationships, Bound.-Layer Meteorol., 7, 363-372. INTEGRATION: Paulson, C. A. (1970), The mathematical representation of wind speed and temperature profiles in the unstable atmospheric surface layer, J. Appl. Meteorol., 9, 857-861.
|
| |
| real(kind=rp) function | most_cpu::slaw_h_convective (z, l_ob, z0h) |
| |
| real(kind=rp) function | most_cpu::corr_m_convective (z, l_ob) |
| |
| real(kind=rp) function | most_cpu::corr_h_convective (z, l_ob) |
| |
| real(kind=rp) function | most_cpu::slaw_m_neutral (z, l_ob, z0) |
| | Similarity laws and corrections for the NEUTRAL regime:
|
| |
| real(kind=rp) function | most_cpu::slaw_h_neutral (z, l_ob, z0h) |
| |
| real(kind=rp) function | most_cpu::f_neumann (ri_b, z, z0, z0h, pr, l_ob, slaw_m, slaw_h) |
| | Simialrity laws (different for neumann and dirichlet bc's)
|
| |
| real(kind=rp) function | most_cpu::dfdl_neumann (l_upper, l_lower, z, z0, z0h, pr, l_ob, slaw_m, slaw_h, fd_h) |
| |
| real(kind=rp) function | most_cpu::f_dirichlet (ri_b, z, z0, z0h, pr, l_ob, slaw_m, slaw_h) |
| |
| real(kind=rp) function | most_cpu::dfdl_dirichlet (l_upper, l_lower, z, z0, z0h, pr, l_ob, slaw_m, slaw_h, fd_h) |
| |