from numpy import log10 
from water2 import MorelMaritorena

wav = [412., 443., 490., 510., 565., 670., 865., 1020.]    # wavelength (nm)
chl = 0.8     # chlorophyll (mg/m3)
bbs = 0.      # 
sz  = 30.0     # solar zenith (deg)
vz  = 20.0     # view zenith (deg)
ra  = 90.0     # relative azimuth (deg)
wind = 5.0    # wind speed (m/s)

wmod = MorelMaritorena('../auxdata/common/')

rho_w = wmod.calc(wav, log10(chl), bbs, sza=sz, vza=vz, raa=ra, ws=wind)

print('Water Reflectance(wave): ')
print(rho_w)

'''
$ python run_morel.py
Water Reflectance(wave): 
[0.009709887206554413, 0.009241287596523762, 0.010147521272301674, 0.008623463101685047, 
 0.006389400456100702, 0.0009943331824615598, 0.00010384673805674538, 0.0]
'''