I am trying to find the formulas for surface exchange coefficient for heat and upward moisture flux from the surface. I already know that my fluxes are calculated in the module_sf_myjsfc, but I have no idea where is AKHS calculated? And hot to obtain TLOW?PLOW?QLOW?

In the attachment is my namelist_input file.

&time_control

run_hours = 30,

start_year = 1998,

start_month = 02,

start_day = 01,

start_hour = 00,

start_minute = 00,

start_second = 00,

history_interval = 10,

frames_per_outfile = 1,

restart = .false.,

restart_interval = 60,

io_form_history = 2,

io_form_restart = 2,

io_form_input = 2,

io_form_boundary = 2,

diag_print = 0,

debug_level = 0,

nocolons = .true.,

/

&domains

time_step = 1,

use_adaptive_time_step = .false.,

step_to_output_time = .true.,

min_time_step = 1,

min_time_step_den = 2,

max_dom = 1,

s_we = 1,

e_we = 201,

s_sn = 1,

e_sn = 201,

s_vert = 1,

e_vert = 61,

dx = 100,

dy = 100,

ztop = 2000,

max_ts_locs = 20,

ts_buf_size = 200

/

&physics

mp_physics = 4,

ra_lw_physics = 4,

ra_sw_physics = 2,

radt = 10,

swint_opt = 1,

sf_sfclay_physics = 2,

sf_surface_physics = 2,

bl_pbl_physics = 0,

bldt = 0,

cu_physics = 0,

isfflx = 1,

icloud = 0,

num_soil_layers = 4,

sst_skin = 1

seaice_threshold = 0.0

fractional_seaice = 1

seaice_albedo_opt = 0

seaice_albedo_default = 0.75

seaice_snowdepth_opt = 0

seaice_snowdepth_max = 1.e10

seaice_snowdepth_min = 0.001

seaice_thickness_opt = 0

seaice_thickness_default = 1.5

/

/

&dynamics

rk_ord = 3,

diff_opt = 2,

km_opt = 2,

damp_opt = 2,

zdamp = 200.0,

dampcoef = .003,

top_lid = .false.,

rad_nudge = 0,

smdiv = 0.1,

emdiv = 0.01,

epssm = 0.1,

time_step_sound = 6,

non_hydrostatic = .true.,

moist_adv_opt = 1,

/

&bdy_control

periodic_x = .true.,

symmetric_xs = .false.,

symmetric_xe = .false.,

open_xs = .false.,

open_xe = .false.,

periodic_y = .true.,

symmetric_ys = .false.,

symmetric_ye = .false.,

open_ys = .false.,

open_ye = .false.,

/

&grib2

/

&namelist_quilt

nio_tasks_per_group = 0,

nio_groups = 1,

/

PLEASE HELP!