LAYTYP, AMFORM, LAYHT = 1, 1, 0.0
H_base, P_base, INTERTYP = None, None, 1
# layer properties
Layer = Layer(RADIUS=RADIUS,
LAYTYP=LAYTYP,
NLAY=NLAY,
LAYINT=LAYINT,
NINT=101,
AMFORM=AMFORM,
INTERTYP=INTERTYP,
H_base=H_base,
P_base=P_base)
BASEH, BASEP, BASET, HEIGHT, PRESS, TEMP, TOTAM, AMOUNT, PP, LAYSF, DELH\
= Layer.integrate(H=H,P=P,T=T, LAYANG=LAYANG, ID=ID,VMR=VMR)
print('{:<8} {:<8} {:<8} {:<8} {:<8} {:<8} {:<8} '.format(
'BASEH', 'DELH', "BASEP", 'BASET', 'TOTAM', 'PRESS', 'TEMP'),
end="")
for j in range(NVMR):
print('{:<8} {:<8}'.format('AMOUNT', 'PP'), end='')
for i in range(len(BASEH)):
print(
'\n{} {:<8.2f} {:<8.2f} {:<8.3E} {:<8.0f} {:<8.3e} {:<8.3e} {:<8.0f} '.
format(i + 1, BASEH[i] / 1e3, DELH[i] / 1e3, BASEP[i] / 101325,
BASET[i], TOTAM[i] / 1e4, PRESS[i] / 101325, TEMP[i]),
end="")
for j in range(NVMR):
print('{:<10.3E} {:<10.3E}'.format(AMOUNT[i, j] / 1e4,
PP[i, j] / 101325),
