# -----------------------------------------------------------------
if ltimedep == "true":
lsflux = "ls_flux.inp"
headerlsflux1 = ["time", "wtsurf", "wqsurf", "thls", "qts", "psurf"]
headerlsflux2 = ["[s]", "[Km/s]", "[kgm/s]", "[K]", "[kg/kg]", "[Pa]"]
headerlsflux3 = ["height", "ug", "vg", "wfls", "dqdtx", "dqtdy", "dqtdtls", "dthlrad"]
# -----------------------------------------------------------------
# 3.4.1 First block of ls_flux.inp
# -----------------------------------------------------------------
nt1 = 10 # number of timesteps
dt1 = runtime / float(nt1 - 1) # timestep
lsfluxtopin = dlsin.readlsfluxinptop(case=sourcecasetitle, subcase=sourcecasesubtitle, timesteps=nt1)[
"ls_flux.inp"
] # read first block of ls_flux.inp from given case
lsflux1 = zeros((nt1, 6))
# column 0 (time)
for i in range(0, nt1):
lsflux1[i, 0] = dt1 * i
# column 3 (liquid water potential temperature at surface)
# lsflux1[i,3] = 265.+0.5*i
# column 1 (liquid water potential temperature flux at surface)
# lsflux1[:,1] = 0.2
lsflux1[:, 1] = lsfluxtopin[:, 1]
# column 2 (total humidity flux at surface)
lsflux1[:, 2] = lsfluxtopin[:, 2]
# column 3 (liquid water potential temperature at surface)
# 5.3 ls_flux.inp
#-----------------------------------------------------------------
if ltimedep == 'true':
lsflux = 'ls_flux.inp'
headerlsflux1 = ['time', 'wtsurf', 'wqsurf', 'thls', 'qts', 'psurf']
headerlsflux2 = ['[s]', '[Km/s]', '[kgm/s]', '[K]', '[kg/kg]', '[Pa]']
headerlsflux3 = ['height', 'ug', 'vg', 'wfls', 'dqdtx', 'dqtdy', 'dqtdtls', 'dthlrad']
#-----------------------------------------------------------------
# 3.4.1 First block of ls_flux.inp
#-----------------------------------------------------------------
nt1 = 10
dt1 = runtime/float(nt1-1)
lsfluxtopin = dlsin.readlsfluxinptop(case=sourcecasetitle,subcase=sourcecasesubtitle,timesteps=nt1)['ls_flux.inp']
lsflux1 = np.zeros((nt1,6))
# column 0 (time)
for i in range(0,nt1):
lsflux1[i,0] = dt1*i
# column 1 (liquid water potential temperature flux at surface)
lsflux1[:,1] = lsfluxtopin[:,1]
# column 2 (total humidity flux at surface)
lsflux1[:,2] = lsfluxtopin[:,2]
# column 3 (liquid water potential temperature at surface)
lsflux1[:,3] = lsfluxtopin[:,3]
# column 4 (total humidity at surface)
lsflux1[:,4] = lsfluxtopin[:,4]
# column 5 (psurf)
